Study area and sample collection
Specimens of Mycteroperca rosacea were obtained in the port of Santa Rosalía (27° 20.353 'N; 112° 15.797' W), Baja California Sur, in the Gulf of California, Mexico (Fig. 1). Since 1885, this port has been extracting copper through the French mining company “Compagnie du Boleo, S.A.”. Currently, the company operates under the name “Minera y Metalúrgica del Boleo S.A. de C.V.” (Huerta-Díaz et al. 2014). About three million tons of accumulated slag are accumulated in the mine site, but unquantified amounts are already deposited on the adjacent marine environment (Shumilin et al. 2013; Huerta-Díaz et al. 2014). In addition to the mining industry, Santa Rosalía is height port for fishing and tourism vessels, and cargo and passenger ferries, and numerous small boats (Huerta-Díaz et al. 2014). The marine environment is, therefore, potentially polluted for heavy metals and marine organisms that inhabit the area, including the leopard grouper, being impacted.
Organisms were captured monthly (March 2014 to May 2015) by spearfishing while free diving, trying to obtain a good representation of different sizes. Specimens were stored on ice until processing at the Fish Ecology Laboratory at CICIMAR-IPN, where the total length (TL, cm) and weight (g) were recorded; sex was identified by direct observation of the gonads and later corroborated with histological analyses (Nikolsky 1963; Pérez-Olivas et al. 2018). Four maturity stages were considered according to Pérez-Olivas et al. (2018).
A total of 345 M. rosacea individuals was caught, with sizes ranging from 21 to 74 cm total length; 93 specimens were male and 185 were female. A total of 15 bisexual immature organisms, as well as, 52 individuals for which sex could not be histologically identified. These specimens were excluded in the comparative analyses by sex and maturity stages. According to fish size, the specimens were binned in small (n = 161, 47%), medium (n = 148, 43%), and large (n = 36, 10%). Besides, the specimens were categorized by caught season, in cold (n = 244, 71%) and warm season (n = 101, 21%), by year: 2014 (n = 180, 72%) and 2015 (n = 68, 27%).
Fish were dissected and the skin was removed to extract 5.0 g of muscle tissue from the anterior dorsal part of each individual. Each sample was tagged, stored in plastic bags, and frozen at -20ºC.
Heavy metal analysis
Muscle tissues samples were lyophilize at a 0.120 mBar pressure and -40ºC for 72 h (Labconco, FreeZone 2.5). The water content (%) was calculated by weigh differences between fresh frozen and dried sample. The dried samples were ground using an agate mortar, homogenized, and packed into trace metal cleaned plastic vials.
Then, a lot of 248 composite muscle samples were sent to the Stable Isotope Laboratory at the ICMYL-UNAM at Mazatlán for the elemental analysis. Samples were processed and analyzed in HEPA filtered air (Class 1000), trace metal clean laboratory using high-purity reagents (trace metal grade) and water (18 MW cm-1; Milli-Q academic). Aliquots of fish muscle samples, blanks and certified reference material were digested in Teflon vials (Savillex) with 10 mL concentrated nitric acid (HNO3). Containers were then placed on mod-block unit (120°C) for 4 h. After digestion, samples were transferred to polyethylene vials and diluted with Milli-Q water to known volume (about 25 mL).
Elements were analyzed in a Thermo Scientific Element XR magnetic sector high-resolution inductively coupled plasma source mass spectrometry (HR-ICP-MS) (Soto-Jiménez et al. 2008). Accuracy of the analyses was determined with concurrent processing and analysis of certified reference material (CRM) Dog-Fish muscle (DORM-3) of Institute for National Measurement Standards of the National Research Council Canada. Recovery values of CRM for Cd was 92%, Cu 94%, Pb 88%, and Zn 92%. The detection limits were <6 ng g-1 dry weight for Cd, <10 ng g-1 for Pb, and <20 ng g-1 for Cu and Zn. Concentrations of metals in fish were calculated in μg/g dry weight but reported as averages ± standard deviations in μg/g wet weight.
Toxicological evaluation
The provisional tolerable weekly intake (PTWI) of Cu, Cd, Pb, and Zn (PTWI, in μg week-1 kg-1 of body weight) to establish the amount of leopard grouper fillet that the human population could consume without a health risk, was calculated by PTWI = [FAW] [MT] [PAW].
Where FAW is the average weight of fish fillet consumed by an adult (12 kg per person per year; De la O-Burrola et al. 2015), MT is the metal concentration (μg /g, w.w.) in muscle tissue of leopard grouper, and PAW is the average weight of the general population, considering women (pregnant and lactating) and children (4 to 6 years old).
Average worldwide weights established by the World Health Organization (e.g., 70 kg for men, 60 kg for women, and 16 kg for 4-6 year old children) and the average weights for the Santa Rosalía population (75 kg for men, 69 kg for women, and 18 kg for 4-6 year old children) (INEGI 2015) were taken into account. These average weights were also used for further evaluations.
The daily mineral intake of essential elements (Cu, Zn) was also calculated based on the consumption recommended by the Institute of Medicine of the National Academies of the United States (www.iom.edu), using the following equation: DMI = C * 100 / RDI
Where DMI = daily mineral intake of essential elements, C = Cu or Zn concentration in 100 g of fresh fish weight, RDI = recommended daily intake, estimated as safe and adequate for each group in the population. CuRDI: men, women, and children = 900 μg /day, pregnant and lactating women =1000 μg /day; ZnRDI: men, women, pregnant, and lactating women = 40 μg /day, children = 34μg /day. These values were estimated for a healthy population.
Additionally, the potential health risk due to harmful effects from long-term leopard grouper consumption (e.g., months or years) was calculated as the target hazard quotient, where high THQ values (>1) would represent a health risk. The THQ was calculated as follows:
Where C is the concentration of Cu, Pb, Cd, or Zn (μg /g); EF is the exposure frequency (days/year); ED is the exposure duration (years/time); FIR is the fish intake rate (g/day); RfD is the daily fish intake rate (g/day); BW is the average body weight of the human population (kg); and AT is the average exposure time (days). All parameters except C, FIR, and BW were obtained from real data for the Santa Rosalía, BCS, Mexico population. The remaining data used in the formula were obtained from the FAO (2017), US EPA (2015), and Yi et al. (2011) and are shown in Table 1.
Table 1 Parameters and values used for the THQ formula for each heavy metal
Factor
|
Definition
|
Unit
|
Value
|
Reference
|
EF
|
Exposure frequency
|
Days/year
|
365
|
Yi et al. (2011)
|
ED
|
Duration
|
Years
|
78
|
Yi et al. (2011)
|
FIR
|
Fish intake rate
|
g/day
|
32.8
|
|
RfD
|
Reference dose
|
μg/kg/day
|
1.0E-01 (Cd), 2.0E+01 (Pb), 4.0E+01 (Cu), 3.0E+02 (Zn)
|
US EPA (2015)
|
BW
|
Average weight
|
Kg
|
70 ± 75
|
WHO; INEGI (2015)
|
AT
|
Average exposure time
|
Days
|
365*78 = 28.470
|
Yi et al. (2011)
|
Because heavy metal (Cu, Pb, Cd, and Zn) interactions can cause multiple effects (Gu et al., 2017), we considered the added effect of these elements in the THQ as was proposed by Chien et al. (2002) and Gu et al. (2017), as shown in the following equation:
Total THQ (TTHQ) = ∑THQ Cu,Pb,Cd,Zn
Data analysis
Data were grouped by sex (males and females), maturity stages (stage 1, stage 2, stage 3 and stage 4) and size (small: <36 cm, medium: >36 cm and ≤51cm, large: >51 cm). Seasonality was defined according to temperature records for the study period obtained from MODIS-AQUA satellite images with 1.1 km resolution. Anomalies were recorded based on the annual average of 23ºC (Moreno-Sánchez et al. 2019). Months with positive anomalies were assigned to the warm season and negative anomalies were assigned to the cold season. Simple linear regression (LR) analysis was carried to assess the association between the concentration of each heavy metal and the size of leopard groupers. Data base was split into categories (sex, maturity stages and season) and a LR of the heavy metal concentration as a function of size was fitted for each heavy metal and each category. With this analysis, the hypothesis that the b coefficient of the LR model was zero (e.g. no association between heavy metal concentration and size) was tested. All analyses were performed using packages found in R (R Core Team 2017).