The chemical defensome of five model teleost fish

Organ-specific biotransformation in salmonids: Insight into intrinsic enzyme activity and biotransformation of three micropollutants

Aquatic ecosystems continue to be threatened by chemical pollution. To what extent organisms are able to cope with chemical exposure depends on their ability to display mechanisms of defense across different organs. Among these mechanisms, biotransformation processes represent key physiological responses that facilitate detoxification and reduce the bioaccumulation potential of chemicals. Biotransformation does not only depend on the ability of different organs to display biotransformation enzymes but also on the affinity of chemicals towards these enzymes. In the present study, we explored the ability of different organs and of two freshwater fish to support biotransformation processes through the determination of in vitro phase I and II biotransformation enzyme activity, and their role in supporting intrinsic clearance and the formation of biotransformation products. Three environmentally relevant pollutants were evaluated: the polycyclic aromatic hydrocarbon (PAH) pyrene (as recommended by the OECD 319b test guideline), the fungicide azoxystrobin, and the pharmaceutical propranolol. Comparative studies using S9 sub-cellular fractions derived from the liver, intestine, gills, and brain of brown trout (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) revealed significant phase I and II enzyme activity in all organs. However, organ- and species-specific differences were found. In brown trout, significant extrahepatic biotransformation was observed for pyrene but not for azoxystrobin and propranolol. In rainbow trout, the brain appeared to biotransform azoxystrobin. In this same species, propranolol appeared to be biotransformed by the intestine and gills. Biotransformation products could be detected only from hepatic biotransformation, and their profiles and formation rates displayed species-specific patterns and occurred at different magnitudes. Altogether, our findings further contribute to the current understanding of organ-specific biotransformation capacity, beyond the expression and activity of enzymes, and its dependence on specific enzyme-chemical interactions to support mechanisms of defense against exposure.

Integrative Inducer Intervention and Transcriptomic Analyses Reveal the Metabolism of Paralytic Shellfish Toxins in Azumapecten farreri

Paralytic shellfish toxins (PSTs) are widely distributed neurotoxins, and the PST metabolic detoxification mechanism in bivalves has received increasing attention. To reveal the effect of phase I (cytochrome P450)-II (GST)-III (ABC transport) metabolic systems on the PST metabolism in Azumapecten farreri, this study amplified stress on the target systems using rifampicin, dl-α-tocopherol, and colchicine; measured PST levels; and conducted transcriptomic analyses. The highest toxin content reached 1623.48 μg STX eq/kg in the hepatopancreas and only 8.8% of that in the gills. Inducer intervention significantly decreased hepatopancreatic PST accumulation. The proportional reductions in the rifampicin-, dl-α-tocopherol-, and colchicine-induced groups were 55.3%, 50.4%, and 36.1%, respectively. Transcriptome analysis showed that 11 modules were significantly correlated with PST metabolism (six positive/five negative), with phase I CYP450 and phase II glutathione metabolism significantly enriched in negatively correlated pathways. Twenty-three phase I-II-III core genes were further validated using qRT-PCR and correlated with PST metabolism, revealing that CYP46A1, CYP4F6, GSTM1, and ABCF2 were significantly correlated, while CYP4F11 and ABCB1 were indirectly correlated. In conclusion, phase I-II-III detoxification enzyme systems jointly participate in the metabolic detoxification of PSTs in A. farreri. This study provides key data support to profoundly elucidate the PST metabolic detoxification mechanism in bivalves.

Biomarker selection depends on gene function and organ: the case of the cytochrome P450 family genes in freshwater fish exposed to chronic pollution

Pollution and its effects have been of major concern in recent decades. Many strategies and markers have been developed to assess their effects on biota. Cytochrome P450 (CYP) genes have received significant attention in this context because of their relationship with detoxification and activation of exogenous compounds. While their expression has been identified as a pollution exposure biomarker, in most cases, it has been tested only after acute exposures and for CYP genes associated with exogenous compounds. To elucidate CYP gene expression patterns under chronic pollution exposure, we have used the silverside Basilichthys microlepidotus as a model, which inhabits the Maipo River Basin, a freshwater system with different pollution levels. We performed next-generation RNA sequencing of liver and gill tissues from polluted and non-polluted populations. We found most CYP genes were not dysregulated by pollution, and the seven genes that were present and differentially expressed in liver and gill were mainly downregulated. Three CYP genes associated with exogenous compounds showed differential expression in the gill, while four CYP genes associated with endogenous compounds showed differential expression in the liver. The findings presented here highlight the importance of CYP genes, his family, tissues and his interaction in the context of pollution biomarkers use.

Ivermectin Toxicokinetics in Rainbow Trout (Oncorhynchus mykiss) following P-glycoprotein Induction

Alterations in ivermectin (IVM, 22,23-dihydro avermectin B1a+22,23-dihydro avermectin B1b) toxicokinetics following P-glycoprotein (P-gp) induction by clotrimazole (CTZ) were examined in rainbow trout (Oncorhynchus mykiss) to assess the potential importance of P-gp activity levels in xenobiotic distribution and kinetics in fish. Control and fish pretreated with CTZ (30 µmol/kg) were administered 175 µg/kg 3H-IVM into the caudal vasculature. At various time points (0.25, 0.5, 1, 3, 24, 48, 96, and 168 h) following injection, tissues (blood, liver, kidney, gill, intestines, brain [5 regions], eye, gonad and fat) were removed analyzed for IVM-derived radioactivity. IVM concentration declined in blood, liver, kidney and gill, and concentrations in other tissues remained constant over the sampling period. The highest measured concentrations were found in kidney, followed by liver, with the lowest values found in brain, eye and gonad. The highest % of the administered dose was found in the liver and kidney in the immediate hours post-administration, and in the intestines and fat at 24 h post-administration. P-gp induction by CTZ did not alter IVM distribution or any calculated toxicokinetic parameter (AUC, mean residence time, T1/2, clearance rate, volume of distribution), suggesting that P-gp induction may be limited or that P-gp plays a lesser role in xenobiotic kinetics in fish compared to mammals.

Multi-tissue proteogenomic analysis for mechanistic toxicology studies in non-model species

New approach methodologies (NAM), including omics and in vitro approaches, are contributing to the implementation of 3R (reduction, refinement and replacement) strategies in regulatory science and risk assessment. In this study, we present an integrative transcriptomics and proteomics analysis workflow for the validation and revision of complex fish genomes and demonstrate how proteogenomics expression matrices can be used to support multi-level omics data integration in non-model species in vivo and in vitro. Using Atlantic salmon as an example, we constructed proteogenomic databases from publicly available transcriptomic data and in-house generated RNA-Seq and LC-MS/MS data. Our analysis identified ∼80,000 peptides, providing direct evidence of translation for over 40,000 RefSeq structures. The data also highlighted 183 co-located peptide groups that supported a single transcript each, and in each case, either corrected a previous annotation, supported Ensembl annotations not present in RefSeq, or identified novel previously unannotated genes. Proteogenomics data-derived expression matrices revealed distinct profiles for the different tissue types analyzed. Focusing on proteins involved in defense against xenobiotics, we detected distinct expression patterns across different salmon tissues and observed homology in the expression of chemical defense proteins between in vivo and in vitro liver systems. Our study demonstrates the potential of proteogenomic analyses in extending our understanding of complex fish genomes and provides an advanced bioinformatic toolkit to support the further development of NAMs and their application in regulatory science and (eco)toxicological studies of non-model species.

Pesticides bioassays using neotropical aquatic species: Trends during the last twenty years and future challenges in Argentina

The presence of pesticides in aquatic ecosystems is one of the most relevant stressors which biota usually face. Laboratory tests using model organisms for pesticides toxicity assessment are employed worldwide. The use of these species has been encouraged in the scientific community due to their advantageous features and their acceptation by regulatory and standardization organizations. However, non-model species as well as those belonging particular ecosystems could contribute in the laboratory-field toxicity extrapolation. In this context, this work aims on exploring the state of the ecotoxicological studies of pesticides in neotropical aquatic species, focusing on bioassays performed in Argentina over the last 20 years as a case of study. Furthermore, we analyzed the possible advantages and disadvantages of these studies, possible differential sensitivities among native and model species, and future challenges to be faced. The analysis of more than 150 publications allowed identify the chemical identity of tested compounds, organisms used for the bioassays, characteristics of the experimental designs, and the toxicity endpoints. Particularly, the studied cases showed that the tested chemicals are related to those most used in the agricultural activity in Argentina, the predilection for particular species in some taxonomic groups (e.g. amphibians), and the wide election of biochemical biomarkers in the studies. Regarding the sensitivity comparison between native and non-native species, the amount of data available indicates that there is not a clear difference beyond some particular cases. However, deeper understanding of toxic effects of pesticides on non-model species could help in a more comprehensive ecological risk assessment in different ecosystems.

No observed developmental effects in early life stages of capelin (Mallotus villosus) exposed to a water-soluble fraction of crude oil during embryonic development

The rise in offshore oil and gas operations, maritime shipping, and tourism in northern latitudes enhances the risk of oil spills to sub-Arctic and Arctic coastal environments. Therefore, there is a need to understand the potential adverse effects of petroleum on key species in these areas. Here, we investigated the effects of oil exposure on the early life stages of capelin (Mallotus villosus), an ecologically and commercially important Barents Sea forage fish species that spawns along the coast of Northern Norway. Capelin embryos were exposed to five different concentrations (corresponding to 0.5-19 µg/L total PAHs) of water-soluble fraction (WSF) of crude oil from 6 days post fertilization (dpf) until hatch (25 dpf), and development of larvae in clean seawater was monitored until 52 dpf. None of the investigated endpoints (embryo development, larval length, heart rate, arrhythmia, and larval mortality) showed any effects. Our results suggest that the early life stages of capelin may be more robust to crude oil exposure than similar life stages of other fish species.

Thermal imprinting during embryogenesis modifies skin repair in juvenile European sea bass (Dicentrarchus labrax)

Fish skin is a multifunctional tissue that develops during embryogenesis, a developmental stage highly susceptible to epigenetic marks. In this study, the impact of egg incubation temperature on the regeneration of a cutaneous wound caused by scale removal in juvenile European sea bass was evaluated. Sea bass eggs were incubated at 11, 13.5 and 16 °C until hatching and then were reared at a common temperature until 9 months when the skin was damaged and sampled at 0, 1 and 3 days after scale removal and compared to the intact skin from the other flank. Skin damage elicited an immediate significant (p < 0.001) up-regulation of pcna in fish from eggs incubated at higher temperatures. In fish from eggs incubated at 11 °C there was a significant (p < 0.001) up-regulation of krt2 compared to fish from higher thermal backgrounds 1 day after skin damage. Damaged epidermis was regenerated after 3 days in all fish irrespective of the thermal background, but in fish from eggs incubated at 11 °C the epidermis was significantly (p < 0.01) thinner compared to other groups, had less goblet cells and less melanomacrophages. The thickness of the dermis increased during regeneration of wounded skin irrespective of the thermal background and by 3 days was significantly (p < 0.01) thicker than the dermis from the intact flank. The expression of genes for ECM remodelling (mmp9, colXα, col1α1, sparc, and angptl2b) and innate immunity (lyg1, lalba, sod1, csf-1r and pparγ) changed during regeneration but were not affected by egg thermal regime. Overall, the results indicate that thermal imprinting of eggs modifies the damage-repair response in juvenile sea bass skin.

Enniatin B and beauvericin affect intestinal cell function and hematological processes in Atlantic salmon (Salmo salar) after acute exposure

Unintentional use of mold-infested plant-based feed ingredients are sources of mycotoxins in fish feeds. The presence of the emerging mycotoxins ENNB and BEA in Norwegian commercial fish feeds and plant-based feed ingredients has raised concerns regarding the health effects on farmed Atlantic salmon (Salmon salar). Atlantic salmon pre-smolts were exposed to non-lethal doses of BEA and ENNB (ctrl, 50 and 500 μg/kg feed for 12 h), after which total RNA sequencing of the intestine and liver was carried out to evaluate gut health and identify possible hepatological changes after acute dietary exposure. ENNB and BEA did not trigger acute toxicity, however ENNB caused the onset of pathways linked to acute intestinal inflammation and BEA exposures caused the onset of hepatic hematological disruption. The prevalence and concentration of ENNB found in today's commercial feed could affect the fish health if consumed over a longer time-period.

Sodium arsenite but not aluminum chloride stimulates ABC transporter activity in renal proximal tubules of killifish (Fundulus heteroclitus)

ABC export proteins including Multidrug resistance-related protein 2 (Mrp2) serve as detoxification mechanism in renal proximal tubules due to active transport of xenobiotics and metabolic waste products into primary urine. The environmental pollutants aluminum and arsenic interfere with a multitude of regulatory mechanisms in the body and here their impact on ABC transporter function was studied. NaAsO₂ but not AlCl₃ rapidly stimulated Mrp2-mediated Texas Red (TR) transport in isolated renal proximal tubules from killifish, a well-established laboratory model for the determination of efflux transporter activity by utilizing fluorescent substrates for the ABC transporters of interest and confocal microscopy followed by image analysis. This observed stimulation remained unaffected by the translation inhibitor cycloheximide (CHX), but it was abrogated by antagonists and inhibitors of the endothelin receptor type B (ET_B)/nitric oxide synthase (NOS)/protein kinase C (PKC) signaling pathway. NaAsO₂-triggered effects were abolished as a consequence of PKCα inhibition through Gö6976 and PKCα inhibitor peptide C2-4. Phosphatidylinositol 3-kinase (PI3K) inhibitor LY 294,002 as well as the mammalian target of rapamycin (mTOR) inhibitor rapamycin suppressed NaAsO₂-triggered stimulation of luminal TR transport. In addition, the stimulatory effect of NaAsO₂ was abolished by GSK650394, an inhibitor of serum- and glucocorticoid-inducible kinase 1 (SGK1), which is an important downstream target. Environmentally relevant concentrations of NaAsO₂ further stimulated transport function of P-glycoprotein (P-gp), Multidrug resistance-related protein 4 (Mrp4) and Breast cancer resistance protein (Bcrp) while AlCl₃ was ineffective. To our knowledge, this is the first report engaging in the impact of NaAsO₂ on efflux transporter signaling and it may contribute to the understanding of defense mechanisms versus this worrying pollutant.

Variation in the Transcriptome Response and Detoxification Gene Diversity Drives Pesticide Tolerance in Fishes

Pesticides are critical for invasive species management but often have negative effects on nontarget native biota. Tolerance to pesticides should have an evolutionary basis, but this is poorly understood. Invasive sea lamprey (Petromyzon marinus) populations in North America have been controlled with a pesticide lethal to them at lower concentrations than native fishes. We addressed how interspecific variation in gene expression and detoxification gene diversity confer differential pesticide sensitivity in two fish species. We exposed sea lamprey and bluegill (Lepomis macrochirus), a tolerant native species, to 3-trifluoromethyl-4-nitrophenol (TFM), a pesticide commonly used in sea lamprey control. We then used whole-transcriptome sequencing of gill and liver to characterize the cellular response in both species. Comparatively, bluegill exhibited a larger number of detoxification genes expressed and a larger number of responsive transcripts overall, which likely contributes to greater tolerance to TFM. Understanding the genetic and physiological basis for pesticide tolerance is crucial for managing invasive species.

Medaka embryos as a model for metabolism of anabolic steroids

In anti-doping science, the knowledge of drug metabolism is a prerequisite to identify analytical targets for the detection of misused prohibited substances. As the most obvious way to study xenobiotic metabolism, the administration to human volunteers, faces ethical concerns, there is a need for model systems. In the present study, we investigated whether Oryzias latipes (medaka) embryos might be an alternative, non-animal test model to study human-like metabolism. In the present study, we exposed medaka embryos at the morula stage to the anabolic steroid metandienone (10 µM or 50 µM) for a period of 2 or 8 days. According to the fish embryo toxicity test (OECD test), we assessed the developmental status of the embryos. We further investigated metandienone metabolites by high-performance liquid chromatography- and gas chromatography-mass spectrometry. Medaka embryos produced three mono-hydroxylated and one reduced metabolite known from human biotransformation. Developmental malformations were observed for the exposition to 50 µM metandienone, while a significant elevation of the heart beat was also present in those individuals exposed to the lower dose for 8 days. The present study demonstrates that the medaka embryo represents a promising model to study human-like metabolism. Moreover, the judgement of developmental parameters of the fish embryos enables for the simultaneous assessment of toxicity.

A genome assembly of the Atlantic chub mackerel (Scomber colias): a valuable teleost fishing resource

The Atlantic chub mackerel, Scomber colias (Gmelin, 1789), is a medium-sized pelagic fish with substantial importance in the fisheries of the Atlantic Ocean and the Mediterranean Sea. Over the past decade, this species has gained special relevance, being one of the main targets of pelagic fisheries in the NE Atlantic. Here, we sequenced and annotated the first high-quality draft genome assembly of S. colias, produced with PacBio HiFi long reads and Illumina paired-end short reads. The estimated genome size is 814 Mbp, distributed into 2,028 scaffolds and 2,093 contigs with an N50 length of 4.19 and 3.34 Mbp, respectively. We annotated 27,675 protein-coding genes and the BUSCO analyses indicated high completeness, with 97.3% of the single-copy orthologs in the Actinopterygii library profile. The present genome assembly represents a valuable resource to address the biology and management of this relevant fishery. Finally, this genome assembly ranks fourth in high-quality genome assemblies within the order Scombriformes and first in the genus Scomber.

In vitro-in vivo biotransformation and phase I metabolite profiling of benzo[a]pyrene in Gulf killifish (Fundulus grandis) populations with different exposure histories

Chronic exposure to pollution may lead populations to display evolutionary adaptations associated with cellular and physiological mechanisms of defense against xenobiotics. This could result in differences in the way individuals of the same species, but inhabiting different areas, cope with chemical exposure. In the present study, we explore two Gulf killifish (Fundulus grandis) populations with different exposure histories for potential differences in the biotransformation of benzo[a]pyrene (BaP), and conduct a comparative evaluation of in vitro and in vivo approaches to describe the applicability of new approach methodologies (NAMs) for biotransformation assessments. Pollution-adapted and non-adapted F. grandis were subjected to intraperitoneal (IP) injections of BaP in time-course exposures, prior to measurements of CYP biotransformation activity, BaP liver concentrations, and the identification and quantification of phase I metabolites. Additionally, substrate depletion bioassays using liver S9 fractions were employed for measurements of intrinsic hepatic clearance and to evaluate the production of metabolites in vitro. Pollution-adapted F. grandis presented significantly lower CYP1A activity and intrinsic clearance rates that were 3 to 4 times lower than non-adapted fish. The metabolite profiling of BaP showed the presence of 1‑hydroxy-benzo[a]pyrene in both the in vitro and in vivo approaches but with no significant population differences. Contrarily, 9‑hydroxy-benzo[a]pyrene and benzo[a]pyrene-4,5-dihydrodiol, only identified through the in vivo approach, presented higher concentrations in the bile of pollution-adapted fish relative to non-adapted individuals. These observations further the understanding of the evolutionary adaptation of F. grandis inhabiting heavily polluted environments in the Houston Ship Channel, TX, USA, and highlight the need to consider the evolutionary history of populations of interest during the implementation of NAMs.