Variation in variance means more than mean variations: what does variability tell us about population health status?

Effects of captive-breeding conditions on metabolic and performance traits in an endangered, endemic cyprinidontiform fish

Captive breeding and stocking are commonly employed strategies for enhancing fisheries and conserving endangered fish species. However, hatchery-raised fish often exhibit reduced performance in the wild, displaying alterations in physiological, morphological, and behavioral traits. We tested for differences in swimming capacity and metabolic traits between wild and hatchery-reared individuals of the Spanish toothcarp (Aphanius iberus) from 2 different populations. Furthermore, we experimentally tested if these changes translated into fitness differences after their stocking into the wild. There were significant differences in swimming capacity and metabolic traits between wild and hatchery-reared individuals and also between the 2 populations. Captive-bred individuals displayed consistently lower metabolic rates than wild individuals from the same population (30-76% lower). Critical swimming speed rather differed between the 2 populations. Sex-specific differences were observed in maximum and standard metabolic rates, with wild individuals and females generally exhibiting higher values but with some exceptions. During a 3-month experiment, survival rates did not significantly differ between wild and captive-bred fish. Captive-bred individuals started smaller but exhibited rapid growth during the experiment. Initially, larger captive-bred fish had lower body conditions than their wild counterparts, but these differences progressively diminished. In summary, captive-bred individuals of this fish species showed lower metabolic rates, although the differences with wild individuals slightly depended on sex and size.

First quantitative biomonitoring study of two ports (marina, commerce) in French littoral area: Evaluation of metals released into the marine environment and resulting from galvanic anodes

Coastal areas are important interface environments between marine and terrestrial ecosystems and are also areas of high economic activity related, among others, to trade, fishing and tourism. The port areas of La Rochelle (France) are part of these areas with strong anthropic activities and are therefore subject to potential contamination, particularly with trace metals. Among the various sources of discharges, galvanic (or sacrificial) cathodic protection can be questioned. Indeed, few studies have been carried out on the potential impact of the degradation of sacrificial anodes in ports on the health of bivalves. A study was therefore carried out on the black scallop Mimachlamys varia in the commercial port and the marina of La Rochelle using the caging method. Several biomarkers of the health status of individuals were observed in the laboratory (SOD, GST, MDA, AChE, and LAC). Among the different results obtained over two exposure periods, it appears that the potential effects of sacrificial anodes, whether new or present on site for several years, are masked by meteorological conditions, as well as by harbor activities.

Environmental pollution impacts: Are p values over-valued?

An examination revealed the dominance of the published literature of environmental science by p values. Meanwhile, the use of effect size has been neglected in publications reporting primary data, yet the size of effect is often more informative than p values inference in assessing the effects of pollution on living organisms, comparing susceptibility/resistance among organisms, and ranking pollutants according to their potency, among others. Statistical significance does not necessarily mean biological, practical, or scientific significance, and its use based on (often misinterpreted) p values reflects the average response or effect at average conditions based on an assumed linear model fit to the entire sample. However, pollution impacts and organismal responses are rarely characterized by linear and symmetric features, and dichotomous 'statistical significance' based on p values is inadequate to fully describe data and findings. Considering 'the fallacy of the average', variance, and differential response of different population percentiles in new studies would provide otherwise wasted biologically, practically, or scientifically significant information. Since p values often inform as to whether some findings warrant further examination, journals should consider mandating the reporting of effect sizes and confidence intervals, together with p values (should they be used), to provide more integrated information regarding pollution impacts. Moreover, replacing 'statistical significance' with language of evidence, especially in key components of publications, such as abstracts and conclusions, could help preventing potential misleading of the public and decision and policy makers.

First study on the determination of baseline biomarkers in Mimachlamys varia for an intra-port environmental biomonitoring in French Atlantic coastline (La Rochelle)

The characterization of organic and inorganic environmental pollution in coastal ecosystems, such as port areas, is complex and difficult to carry out due to the effect of environmental variables, as well as anthropic activities. For this study, the objective was to define a statistical method, taking into account the confounding factors influence, to define reference values for biomarkers in the black scallop (Mimachlamys varia). Thus, for five biomarkers (SOD (Superoxyde Dismutase), GST (Glutathion-S Transferase), MDA (Malondialdehyde), AChE (Acetylcholinesterase) and LAC (Laccase)), reference data could be described for individuals placed on sites more or less strongly impacted by specific environmental contaminations in port areas for more than two years, which had never been done before. All these results enabled us to calibrate and validate our approach in terms of active biomonitoring for the evaluation of a good ecological status of the environment of a port located on the French Atlantic coast.

Newly Hatched Stage I American Lobster (Homarus americanus) Survival Following Exposure to Physically and Chemically Dispersed Crude Oil

Standard model species are commonly used in toxicity tests due to their biological and technical advantages but studying native species increases the specificity and relevance of results generated for the potential risk assessment to an ecosystem. Accounting for intraspecies variability and other factors, such as chemical and physical characterization of test medium, is necessary to develop a reproducible bioassay for toxicity testing with native species. In this study, larval stage I American lobster (Homarus americanus), a commercially important and native species of Atlantic Canada, was used as the test species. Toxicity tests were first conducted by exposing lobster larvae to a reference toxicant of copper sulphate (CuSO₄) and then to physically and chemically (using Corexit 9500A) dispersed oil (WAF and CEWAF, respectively). The effect on larval survival was estimated by calculating the 24-h median effect concentration (24-h EC50), and there was no difference between WAF or CEWAF exposure when the results are reported on a total petroleum hydrocarbon (TPH) basis. The 24-h EC50s ranged from 2.54 to 9.73 mg TPH/L when all trials (n = 19) are considered together. The HC5 (hazardous concentration for 5 per cent of the population) value was 2.52 mg TPH/L and similar to the EC50 value when all trials were pooled. To evaluate the reproducibility of the lobster toxicity tests, inter-trial variability was determined, and the resultant coefficients of variation (%CV) were compared to those reported for two standard test species, mysid shrimp (Americamysis bahia) and inland silverside (Menidia beryillina). This comparison showed that the %CV for the lobster toxicity tests were lower than those for the standard species tests indicating that the described larval lobster toxicity bioassay produces reliable and repeatable results.

Sex-specific thermal tolerance limits in the ditch shrimp Palaemon varians: Eco-evolutionary implications under a warming ocean

As global temperatures continue to rise due to global change, marine heatwaves are also becoming more frequent and intense, impacting marine biodiversity patterns worldwide. Organisms inhabiting shallow water environments, such as the commercially relevant ditch shrimp Palaemon varians, are expected to be the most affected by rising temperatures. Thus, addressing species' thermal ecology and climate extinction-risk is crucial to foster climate-smart conservation strategies for shallow water ecosystems. Here, we estimated sex-specific upper thermal tolerance limits for P. varians via the Critical Thermal Maximum method (CTmax), using loss of equilibrium as endpoint. We further calculated thermal safety margins for males and females and tested for correlations between upper thermal limits and shrimps' body size. To determine sex-biased variation in P. varians' traits (CTmax, weight and length), we compared trait variation between females and males through the coefficient of variation ratio (lnCVR). Females displayed an average CTmax value 1.8% lower than males (CTmax_females = 37.0 °C vs CTmax_males = 37.7 °C). This finding may be related to the larger body size exhibited by females (156% heavier and 39% larger than males), as both length and weight had a significant effect on CTmax. The high energetic investment of females in offspring may also contribute to the differences recorded in thermal tolerance. Overall, organisms with a smaller body-size displayed a greater tolerance to elevated temperature, thus suggesting that smaller individuals may be positively selected in warmer environments. This selection may result in a reduction of size-at-maturity and shifts in sex ratio, given the sexual dimorphism in body size of shrimps. The thermal safety margin of P. varians was narrow (∼2.2 °C for males and ∼1.5 °C for females), revealing the vulnerability of this species to ocean warming and heatwaves.

Application in a biomonitoring context of three-spined stickleback immunomarker reference ranges

The relevance of a biomarker for biomonitoring programs was influenced both by the knowledge on biomarker natural inter-individual and site variabilities and by the sensitivity of the biomarker towards environmental perturbations. To minimize data misinterpretation, robustness reference values for biomarkers were important in biomonitoring programs. Specific three-spined stickleback, Gasterosteus aculeatus, immune reference ranges for field studies had been determined based on laboratory data and one reference station (Contentieuse river at Houdancourt). In this study, data obtained in one uncontaminated and three contaminated sites were compared to these reference ranges as a validation step before considering them for larger scale biomonitoring programs. When the field reference range were compared to data from the uncontaminated station (Béronelle), only few deviations were shown. In this way, data coming from uncontaminated station (Béronelle) was integrated in the field reference ranges to improve the evaluation of site variability. The new field reference ranges provided better discrimination of sites and spanned a larger range of fish lengths than the initial reference ranges. Furthermore, the results suggest lysosomal presence during several months and phagocytosis capacity in autumn may be the most relevant immunomarkers towards identifying contaminated sites. In the future, combining this reference value approach with active biomonitoring could facilitate the obtention of data in multiple stream conditions.

Multiple thresholds and trajectories of microbial biodiversity predicted across browning gradients by neural networks and decision tree learning

Ecological association studies often assume monotonicity such as between biodiversity and environmental properties although there is growing evidence that nonmonotonic relations dominate in nature. Here, we apply machine-learning algorithms to reveal the nonmonotonic association between microbial diversity and an anthropogenic-induced large-scale change, the browning of freshwaters, along a longitudinal gradient covering 70 boreal lakes in Scandinavia. Measures of bacterial richness and evenness (alpha-diversity) showed nonmonotonic trends in relation to environmental gradients, peaking at intermediate levels of browning. Depending on the statistical methods, variables indicative for browning could explain 5% of the variance in bacterial community composition (beta-diversity) when applying standard methods assuming monotonic relations and up to 45% with machine-learning methods taking non-monotonicity into account. This non-monotonicity observed at the community level was explained by the complex interchangeable nature of individual taxa responses as shown by a high degree of nonmonotonic responses of individual bacterial sequence variants to browning. Furthermore, the nonmonotonic models provide the position of thresholds and predict alternative bacterial diversity trajectories in boreal freshwater as a result of ongoing climate and land-use changes, which in turn will affect entire ecosystem metabolism and likely greenhouse gas production.

The inter-individual variance can provide additional information for the ecotoxicologists beside the mean

The hypothesis that the inter-individual parameter variability is an unexploited area of ecotoxicology was proposed several decades ago. Although some illustrative examples were presented to support this hypothesis in the last decades, it has never been tested on an extensive, coherent database. In this study, variance changes of 105 dose-response curves were analysed. All data originated from the same experiment, where the effects of the insecticide Trebon EC were investigated in a dose-response manner on 15 traits of the collembolan Folsomia candida in four subsequent generations and two types of insecticide treatments. A consistent relationship between inter-individual variance and insecticide application was found in 2 (first clutch size and growth-reproduction trade-off) out of the 15 of the parameters. Contrary to the mean, the variance of the first clutch size showed consistent differences compared to the control. Furthermore, the variance of the growth-reproduction trade-off was consistently different from the control except in one case (F3 generation of the transgenerational treatment). Higher first clutch size variances were found in F1 and a lower one in the F2 and F3 generations than in that of the control. This overall pattern of the variance changes of the first clutch size and the trade-off seems to be a quick response to the insecticide application. In the short term, we have found that variance increased with insecticide treatment (P and F1 generation), because phenotypic variance generally increases due to environmental stress. Disruptive selection could be another mechanism between the more detoxification less reproduction strategy and the more reproduction less detoxification strategy. However, in the later generations (F2-F3) the variance decreases compared to the control, which could be because on short term selection stronger on the viability parameters and in long-term selection on reproduction becomes stronger. According to our results, analysis of the variance changes of some parameters may give information about the effects of the pesticide even when the mean does not predict any impact. Testing variance changes are important in ecotoxicology because variance change can signalise toxicant impact even when the mean does not change in certain cases.

Molecular response of a sub-antarctic population of the blue mussel (Mytilus edulis platensis) to a moderate thermal stress

The Kerguelen Islands (49°26'S, 69°50'E) represent a unique environment due to their geographical isolation, which protects them from anthropogenic pollution. The ability of the endemic mussel, part of the Mytilus complex, to cope with moderate heat stress was explored using omic tools. Transcripts involved in six major metabolic functions were selected and the qRT-PCR data indicated mainly changes in aerobic and anaerobic energy metabolism and stress response. Proteomic comparisons revealed a typical stress response pattern with cytoskeleton modifications and elements suggesting increased energy metabolism. Results also suggest conservation of protein homeostasis by the long-lasting presence of HSP while a general decrease in transcription is observed. The overall findings are consistent with an adaptive response to moderate stresses in mussels in good physiological condition, i.e. living in a low-impact site, and with the literature concerning this model species. Therefore, local blue mussels could be advantageously integrated into biomonitoring strategies, especially in the context of Global Change.

Microplastics as a potential risk for aquatic environment organisms – a review

In the modern world, plastic has become a major commodity on global scale and is presented in all life aspects. The production of plastics is increasing dramatically throughout the world and is considered to be a serious threat for the aquatic environment. Scientists started to raise questions as to where all the disappeared plastic debris had gone, as the plastic material is a persistent synthetic polymer. The disappearance of plastic suggests that microplastics are generated by fragmentation of larger plastic debris. They also enter the aquatic environment from the cosmetics, clothing, and other industrial manufacturing. Microplastic particles spread in the aquatic, terrestrial, and atmospheric environments due to their small size and ubiquity, showing their high bioavailability. The presence of microplastics has been demonstrated in the digestive tracts of most aquatic organisms at various trophic levels. Evidence of negative effects of microplastics has been found in relation to fertility, mortality, oxidative stress, immune system or metabolic processes.

Calcium homeostasis and stable fatty acid composition underpin heatwave tolerance of the keystone polychaete Hediste diversicolor

Extreme weather events, such as heatwaves, are becoming increasingly frequent, long-lasting and severe as global climate change continues, shaping marine biodiversity patterns worldwide. Increased risk of overheating and mortality across major taxa have been recurrently observed, jeopardizing the sustainability of ecosystem services. Molecular responses of species, which scale up to physiological and population responses, are determinant processes that modulate species sensitivity or tolerance to extreme weather events. Here, by integrating proteomic, fatty acid profiling and physiological approaches, we show that the tolerance of the intertidal ragworm Hediste diversicolor, a keystone species in estuarine ecosystems and an emergent blue bio-resource, to long-lasting heatwaves (24  vs 30 °C for 30 days) is shaped by calcium homeostasis, immune function and stability of fatty acid profiles. These features potentially enabled H. diversicolor to increase its thermal tolerance limit by 0.81 °C under the heatwave scenario and maintain survival. No growth trade-offs were detected, as wet weight remained stable across conditions. Biological variation of physiological parameters was lower when compared to molecular measures. Proteins showed an overall elevated coefficient of variation, although decreasing molecular variance under the heatwave scenario was observed for both proteins and fatty acids. This finding is consistent with the phenomenon of physiological canalization in extreme environments and contradicts the theory that novel conditions increase trait variation. Our results show that keystone highly valued marine polychaetes are tolerant to heatwaves, confirming the potential of H. diversicolor as a blue bio-resource and opening new avenues for sustainable marine aquaculture development.

Invertebrate responses to microplastic ingestion: Reviewing the role of the antioxidant system

Microplastic ingestion in invertebrates reduces somatic and reproductive growth. This could be caused by energy reserves being detracted from growth processes and redistributed to maintenance processes that preserve life. A potential sink for this diverted energy is the antioxidant system, which minimises oxidative damage and reinstates redox homeostasis following disturbances caused by exposure to pollution. Several microplastic studies have used genetic and molecular redox biomarkers to assess how microplastic ingestion affects the functioning of the antioxidant system. This systematic review synthesises the current understanding of redox biomarker responses in invertebrates that have ingested microplastics. We found that biomarker response information exists for only seven invertebrate taxa, and early life stages have received little scientific attention. The microplastics used by most studies were polystyrene (45% of studies), spherical (51% of studies), and were < 10 μm in diameter (31% of studies). We found multiple examples of microplastic ingestion posing an oxidative challenge to invertebrates, which required upregulation of antioxidant system components. However, the lack of systematic experiments prevented us from clearly identifying which characteristic of microplastics caused these responses. We identify several areas for consideration when investigating biomarker responses to microplastic ingestion and offer research priorities for future studies.

Toxicity of manganese to various life stages of selected marine cnidarian species

Manganese (Mn) pollution in marine waters is increasing and sensitivities to this metal vary widely among marine species. The aims of this study were to characterise Mn chemistry in seawater, and evaluate the toxic effects of Mn on various life stages of two scleractinian corals - the branching sp. Acropora spathulata and massive sp. Platygyra daedalea, and the anemone Exaiptasia pallida. Analytical and theoretical characterisation experiments showed that 97-100% of Mn (II) additions ≤ 200 mg/L in seawater were soluble over 72 h and largely assumed labile complexes. Concentrations estimated to reduce coral fertilisation success by 50% (5.5-h EC50) were 237 mg/L for A. spathulata and 164 mg/L for P. daedalea. A relatively low 72-h LC50 of 7 mg/L was calculated for A. spathulata larvae. In a pilot test using fragments of adult A. spathulata, intact coral tissue rapidly sloughed away from the underlying skeleton at very low concentrations with a 48-h EC50 of just 0.7 mg/L. For E. pallida, survival, tentacle retraction and reproduction were unaffected by prolonged high exposures (12-d NOEC 54 mg/L). This study provides important data supporting the derivation of separate water quality guidelines for Mn in systems with and without coral - a decision recently considered by Australian and New Zealand authorities. It demonstrates the high sensitivity of coral larvae and adult colonies to Mn and the potential risks associated with relying on other early life stage tests and/or E. pallida as ecotoxicological representatives of critically important scleractinian corals.

Beyond the mean: A comparison of trace- and macroelement correlation profiles of two lacustrine populations of the crayfish Procambarus clarkii

In invertebrate biomonitors of chemical pollution, emphasis has been generally given to mean accumulation patterns and how they reflect varying environmental levels of contamination. Intra-population variability, and how it relates with individual phenotypic traits, has received less attention. Here, a set of analytes including trace elements (B, Ba, Cd, Cr, Cu, Fe, Li, Mn, Ni, Pb, Sr, V, and Zn), macroelements (C, Ca, K, Mg, N, Na), and carbon and nitrogen stable isotopes (δ¹³C and δ¹⁵N) was measured in two populations of the crayfish Procambarus clarkii from Lake Trasimeno and Lake Bolsena (Central Italy). The influence of location, sex, body size, and condition factor was assessed; in addition, the analyte correlation profiles of the two populations were compared to verify their congruence. In general, significant inter-lake differences were observed in the concentration of both trace- and macroelements in crayfish tissues, generally mirroring the local chemistry of water and of benthic non-living matrices (sediment and plant detritus). Crayfish CN isotopic signatures excluded the occurrence of inter-lake variations in their omnivorous trophic habits. Correlation profiles varied considerably between the two populations in the nature and strength of bivariate relationships. However, Mantel tests and procrustean analyses indicated a general, significant congruence; C, N, and, to a lesser extent K, Li, Ni, Pb, and δ¹³C showed the highest procrustean residuals, suggesting that their associations with other analytes may be partially influenced by inter-population differences in growing phases. Our study indicates that the local geochemistry of the lacustrine environment influences the elemental fingerprint of Procambarus clarkii; the considerable inter-individual variability in the concentration of analytes, however, does not significantly reflect on their association, thus corroborating its effectiveness as an indicator species.

The characterization of iron (III) in seawater and related toxicity to early life stages of scleractinian corals

Currently toxicity data for iron (Fe) in seawater are limited; furthermore, these data are of poor quality as a result of the importance of Fe solubility in test solutions being overlooked. The present study characterized the solubility and lability of Fe(III) in seawater and then examined the effects of Fe(III) on the fertilization success and larval survival of the tropical marine scleractinian corals Acropora spathulata and Platygyra daedalea. We present the first assessment of the effects of Fe on the early life stages of scleractinian corals. Concentrations of both soluble and labile forms of Fe were very low, with dissolved Fe concentrations ≤0.195 mg/L in bioassay test solutions and chemical determinations revealing labile Fe concentrations ≤1.21 mg/L. For fertilization experiments, the median effect concentration (EC50) value for total Fe was 25 mg/L for the most sensitive species, P. daedalea, whereas the EC50 values for A. spathulata ranged between 40 and 66 mg/L. The median lethal concentration value for P. daedalea larval survival was 47 mg/L Fe after 72-h exposure. We provide Fe toxicity data for tropical marine keystone species that could be used to help generate more reliable guideline values for Fe in marine waters. Environ Toxicol Chem 2018;37:1104-1114. © 2017 SETAC.

A novel methodology for the determination of biomarker baseline levels in the marine polychaete Hediste diversicolor

Identifying environmental damage due to anthropogenic activities is a focal point for scientists and policy makers like those involved in the European Water Framework Directive (WFD). Many of these approaches focus on ecological endpoints for assessing environmental perturbations, which lead to policies emphasizing mitigation rather than prevention. Biomarkers provide early-warning indicators of stress but it is necessary to distinguish their natural variations from those induced by chemical stress. The global aim of this study was to establish a baseline assessment criterion (BAC) using historical data in a reference site to define toxicity thresholds. We have developed a multiple polynomial regression model (MPR) accounting the influence of salinity, temperature and size of individual on energetic reserves (glycogen and lipids) in the marine polychaete Hediste diversicolor. The model identified a complex, orthogonal relationship between confounding factors and glycogen and a linear relationship between lipids and size of individuals.

Detecting the exposure to Cd and PCBs by means of a non-invasive transcriptomic approach in laboratory and wild contaminated European eels (Anguilla anguilla)

Detecting and separating specific effects of contaminants in a multi-stress field context remain a major challenge in ecotoxicology. In this context, the aim of this study was to assess the usefulness of a non-invasive transcriptomic method, by means of a complementary DNA (cDNA) microarray comprising 1000 candidate genes, on caudal fin clips. Fin gene transcription patterns of European eels (Anguilla anguilla) exposed in the laboratory to cadmium (Cd) or a polychloro-biphenyl (PCBs) mixture but also of wild eels from three sampling sites with differing contamination levels were compared to test whether fin clips may be used to detect and discriminate the exposure to these contaminants. Also, transcriptomic profiles from the liver and caudal fin of eels experimentally exposed to Cd were compared to assess the detection sensitivity of the fin transcriptomic response. A similar number of genes were differentially transcribed in the fin and liver in response to Cd exposure, highlighting the detection sensitivity of fin clips. Moreover, distinct fin transcription profiles were observed in response to Cd or PCB exposure. Finally, the transcription profiles of eels from the most contaminated site clustered with those from laboratory-exposed fish. This study thus highlights the applicability and usefulness of performing gene transcription assays on non-invasive tissue sampling in order to detect the in situ exposure to Cd and PCBs in fish.

Ecological relevance of biomarkers in monitoring studies of macro-invertebrates and fish in Mediterranean rivers

Mediterranean rivers are probably one of the most singular and endangered ecosystems worldwide due to the presence of many endemic species and a long history of anthropogenic impacts. Besides a conservation value per se, biodiversity is related to the services that ecosystems provide to society and the ability of these to cope with stressors, including climate change. Using macro-invertebrates and fish as sentinel organisms, this overview presents a synthesis of the state of the art in the application of biomarkers (stress and enzymatic responses, endocrine disruptors, trophic tracers, energy and bile metabolites, genotoxic indicators, histopathological and behavioural alterations, and genetic and cutting edge omic markers) to determine the causes and effects of anthropogenic stressors on the biodiversity of European Mediterranean rivers. We also discuss how a careful selection of sentinel species according to their ecological traits and the food-web structure of Mediterranean rivers could increase the ecological relevance of biomarker responses. Further, we provide suggestions to better harmonise ecological realism with experimental design in biomarker studies, including statistical analyses, which may also deliver a more comprehensible message to managers and policy makers. By keeping on the safe side the health status of populations of multiple-species in a community, we advocate to increase the resilience of fluvial ecosystems to face present and forecasted stressors. In conclusion, this review provides evidence that multi-biomarker approaches detect early signs of impairment in populations, and supports their incorporation in the standardised procedures of the Water Frame Work Directive to better appraise the status of European water bodies.

Insight into the transgenerational effect of benzo[a]pyrene on bone formation in a teleost fish (Oryzias latipes)

Recent cross-generational studies in teleost fish have raised the awareness that high levels of benzo[a]pyrene (BaP) could affect skeletal integrity in the directly exposed F0 and their F1-F2. However, no further details were provided about the causes for abnormalities on the molecular and cellular level and the persistence of such sub-organismal impairments at the transgenerational scale (beyond F2). Adult Oryzias latipes were exposed to 1μg/L BaP for 21days. The F1-F3 were examined for skeletal deformities, histopathological alterations of vertebral bodies and differential expression of key genes of bone metabolism. Significant increase of dorsal-ventral vertebral compression was evident in ancestrally exposed larvae. Histopathological analysis revealed abnormal loss of notochord sheath, a lack of notochord epithelial integrity, reduced bone tissue and decreased osteoblast abundance. A significant downregulation of ATF4 and/or osterix and a high biological variability of COL10, coupled with a significant deregulation of SOX9a/b in the F1-F3 suggest that ancestral BaP exposure most likely perturbed chordoblasts, chondroblast and osteoblast differentiation, resulting in defective notochord sheath repair and rendering the vertebral column more vulnerable to compression. The present findings provide novel molecular and cellular insights into BaP-induced transgenerational bone impairment in the unexposed F3. From the ecological risk assessment perspective, BaP needs to be regarded as a transgenerational skeletal toxicant, which exerts a far-reaching impact on fish survival and fitness. Given that basic mechanisms of cartilage/bone formation are conserved between medaka and mammals, the results may also shed light on the potential transgenerational effect of BaP on the genesis of skeletal diseases in humans.