Ancestors' Gift: Parental Early Exposure to the Environmentally Realistic Pesticide Mixture Drives Offspring Phenotype in a Larger Extent Than Direct Exposure in the Pacific Oyster, Crassostrea gigas.

Marine organisms are threatened by the presence of pesticides in coastal waters. Among them, the Pacific oyster is one of the most studied invertebrates in marine ecotoxicology where numerous studies highlighted the multiscale impacts of pesticides. In the past few years, a growing body of literature has reported the epigenetic outcomes of xenobiotics. Because DNA methylation is an epigenetic mark implicated in organism development and is meiotically heritable, it raises the question of the multigenerational implications of xenobiotic-induced epigenetic alterations. Therefore, we performed a multigenerational exposure to an environmentally relevant mixture of 18 pesticides (nominal sum concentration: 2.85 µg·L-1) during embryo-larval stages (0-48 hpf) of a second generation (F1) for which parents where already exposed or not in F0. Gene expression, DNA methylation, and physiological end points were assessed throughout the life cycle of individuals. Overall, the multigenerational effect has a greater influence on the phenotype than the exposure itself. Thus, multigenerational phenotypic effects were observed: individuals descending from exposed parents exhibited lower epinephrine-induced metamorphosis and field survival rates. At the molecular level, RNA-seq and Methyl-seq data analyses performed in gastrula embryos and metamorphosis-competent pediveliger (MCP) larvae revealed a clear F0 treatment-dependent discrimination. Some genes implicated into shell secretion and immunity exhibited F1:F0 treatment interaction patterns (e.g., Calm and Myd88). Those results suggest that low chronic environmental pesticide contamination can alter organisms beyond the individual scale level and have long-term adaptive implications.

Oyster reproduction is affected by exposure to polystyrene microplastics.

Plastics are persistent synthetic polymers that accumulate as waste in the marine environment. Microplastic (MP) particles are derived from the breakdown of larger debris or can enter the environment as microscopic fragments. Because filter-feeder organisms ingest MP while feeding, they are likely to be impacted by MP pollution. To assess the impact of polystyrene microspheres (micro-PS) on the physiology of the Pacific oyster, adult oysters were experimentally exposed to virgin micro-PS (2 and 6 µm in diameter; 0.023 mg·L(-1)) for 2 mo during a reproductive cycle. Effects were investigated on ecophysiological parameters; cellular, transcriptomic, and proteomic responses; fecundity; and offspring development. Oysters preferentially ingested the 6-µm micro-PS over the 2-µm-diameter particles. Consumption of microalgae and absorption efficiency were significantly higher in exposed oysters, suggesting compensatory and physical effects on both digestive parameters. After 2 mo, exposed oysters had significant decreases in oocyte number (-38%), diameter (-5%), and sperm velocity (-23%). The D-larval yield and larval development of offspring derived from exposed parents decreased by 41% and 18%, respectively, compared with control offspring. Dynamic energy budget modeling, supported by transcriptomic profiles, suggested a significant shift of energy allocation from reproduction to structural growth, and elevated maintenance costs in exposed oysters, which is thought to be caused by interference with energy uptake. Molecular signatures of endocrine disruption were also revealed, but no endocrine disruptors were found in the biological samples. This study provides evidence that micro-PS cause feeding modifications and reproductive disruption in oysters, with significant impacts on offspring.

Additive transcriptomic variation associated with reproductive traits suggest local adaptation in a recently settled population of the Pacific oyster, Crassostrea gigas.

BACKGROUND: Originating from Northeast Asia, the Pacific oyster Crassostrea gigas has been introduced into a large number of countries for aquaculture purpose. Following introduction, the Pacific oyster has turned into an invasive species in an increasing number of coastal areas, notably recently in Northern Europe. METHODS: To explore potential adaptation of reproductive traits in populations with different histories, we set up a common garden experiment based on the comparison of progenies from two populations of Pacific oyster sampled in France and Denmark and their hybrids. Sex ratio, condition index and microarray gene expression in gonads, were analyzed in each progeny (n = 60). RESULTS: A female-biased sex-ratio and a higher condition index were observed in the Danish progeny, possibly reflecting an evolutionary reproductive strategy to increase the potential success of natural recruitment in recently settled population. Using multifarious statistical approaches and accounting for sex differences we identified several transcripts differentially expressed between the Danish and French progenies, for which additive genetic basis is suspected (showing intermediate expression levels in hybrids, and therefore additivity). Candidate transcripts included mRNA coding for sperm quality and insulin metabolism, known to be implicated in coordinated control and success of reproduction. CONCLUSIONS: Observed differences suggest that adaptation of invasive populations might have occurred during expansion acting on reproductive traits, and in particular on a female-biased sex-ratio, gamete quality and fertility.

The Pacific oyster reproduction is affected by early-life exposure to environmental pesticide mixture: A multigenerational study.

Pesticides threat marine organisms worldwide. Among them, the Pacific oyster is a bivalve mollusc model in marine ecotoxicology. A large body of literature already stated on the multiple-scale effects pesticides can trigger in the Pacific oyster, throughout its life cycle and in a delayed manner. In particular, reproductive toxicity is of major concern because of its influence on population dynamics. However, past studies mostly investigated pesticide reprotoxicity as a direct effect of exposure during gametogenesis or directly on gametes and little is known about the influence of an early embryo exposure on the breed capacity. Therefore, we studied delayed and multigenerational consequences through gametogenesis features (i.e. sex ratio, glycogen content, gene expression) and reproductive success in two consecutive oyster generations (F0 and F1) exposed to an environmentally-relevant pesticide mixture (sum nominal concentration: 2.85 µg.L-1) during embryo-larval development (0-48 h post fertilization, hpf). In the first generation, glycogen content increased in exposed individuals and the expression of some gametogenesis target genes was modified. The reproductive success measured 48 hpf was higher in exposed individuals. A multigenerational influence was observed in the second generation, with feminisation, acceleration of gametogenesis processes and the sex-specific modification of glycogen metabolism in individuals from exposed parents. This study is the first to highlight the delayed effects on reproduction induced by an early exposure to pesticides, and its multigenerational implications in the Pacific oyster. It suggests that environmental pesticide contamination can have impacts on the recruitment and the dynamics of natural oyster populations exposed during their embryo-larval phase.

Rapid mitochondrial adjustments in response to short-term hypoxia and re-oxygenation in the Pacific oyster, Crassostrea gigas.

As oxygen concentrations in marine coastal habitats can fluctuate rapidly and drastically, sessile marine organisms such as the oyster Crassostrea gigas can experience marked and rapid oxygen variations. In this study, we investigated the responses of oyster gill mitochondria to short-term hypoxia (3 and 12 h, at 1.7 mg O2 l(-1)) and subsequent re-oxygenation. Mitochondrial respiratory rates (states 3 and 4 stimulated by glutamate) and phosphorylation efficiency [respiratory control ratio (RCR) and the relationship between ADP and oxygen consumption (ADP/O)] were measured. Cytochrome c oxidase (CCO) activity and cytochrome concentrations (a, b, c1 and c) were measured to investigate the rearrangements of respiratory chain subunits. The potential implication of an alternative oxidase (AOX) was investigated using an inhibitor of the respiratory chain (antimycin A) and through gene expression analysis in gills and digestive gland. Results indicate a downregulation of mitochondrial capacity, with 60% inhibition of respiratory rates after 12 h of hypoxia. RCR remained stable, while ADP/O increased after 12 h of hypoxia and 1 h of re-oxygenation, suggesting increased phosphorylation efficiency. CCO showed a fast and remarkable increase of its catalytic activity only after 3 h of hypoxia. AOX mRNA levels showed similar patterns in gills and digestive gland, and were upregulated after 12 and 24 h of hypoxia and during re-oxygenation. Results suggest a set of controls regulating mitochondrial functions in response to oxygen fluctuations, and demonstrate the fast and extreme plasticity of oyster mitochondria in response to oxygen variations.

Molecular and phenotypic effects of early exposure to an environmentally relevant pesticide mixture in the Pacific oyster, Crassostrea gigas.

Early life stages are crucial for organism development, especially for those displaying external fertilization, whose gametes and early stages face environmental stressors such as xenobiotics. The pacific oyster, Crassostrea gigas, is considered a model species in ecotoxicology because of its ecological characteristics (benthic, sessile, filter feeding). So far studies have investigated the impact of xenobiotics at embryotoxic, genotoxic and physiological endpoints, sometimes at the multigenerational scale, highlighting the role of epigenetic mechanisms in transmitting alterations induced by exposure to single xenobiotics. However, to date, little is known about the impact of environmentally-mimicking contaminants cocktails. Thus, we examined the impact of an early exposure to environmentally relevant mixture on the Pacific oyster life history. We studied transcriptomic, epigenetic and physiological alterations induced in oysters exposed to 18 pesticides and metals at environmental concentration (nominal sum concentration: 2.85 µg.L-1, measured sum concentration: 3.74 ± 0.013 µg.L-1) during embryo-larval stage (0-48 h post fertilization, hpf). No significant differences in embryo-larval abnormalities at 24 hpf were observed during larval and spat rearing; the swimming behaviour of exposed individuals was disturbed, while they were longer and heavier at specific time points, and exhibited a lower epinephrine-induced metamorphosis rate as well as a higher survival rate in the field. In addition, RNA-seq analyses of gastrula embryos revealed the differential expression of development-related genes (e.g. Hox orthologues and cell cycle regulators) between control and exposed oysters. Whole-genome DNA methylation analyses demonstrated a significant modification of DNA methylation in exposed larvae marked by a demethylation trend. Those findings suggest that early exposure to an environmentally relevant pesticide mixture induces multi-scale latent effects possibly affecting life history traits in the Pacific oyster.

Differences in chemical contaminants bioaccumulation and ecotoxicology biomarkers in Mytilus edulis and Mytilus galloprovincialis and their hybrids.

The Mytilus mussels are spread all over the world and many related species coexist in several areas and can produce hybrid offspring. Mussels have been used for decades in national and international programs to monitor chemical contamination in the environment. Differences in bioaccumulation and biotransformation abilities between species and their hybrids should be evaluated to assess the comparability of the results obtained within the international biomonitoring programs. The objective of this study was to characterize bioaccumulation abilities and biomarker responses in Mytilus edulis, Mytilus galloprovincialis and their hybrids via an in situ transplantation experimentation on their progenies. Four mussel groups (M. edulis, M. galloprovincialis and two hybrids batches) issued from genetically characterized parents were transplanted for one year in Charente Maritime (France) to ensure their exposure to identical sources of contamination. The bioaccumulation of several families of contaminants (trace metals, polycyclic aromatic hydrocarbons, polybrominated diphenyl ethers, polychlorinated biphenyls), the response of several biomarkers (DNA strand breaks level, lysosomal membrane stability, metallothionein content, acetylcholine esterase activity) and some physiological parameters (growth, mortality, gonadal development), were analyzed. Differences were observed between species, however they were contaminant-specific. Variations in contaminants levels were observed between progenies, with higher levels of Cu, PBDE, PCB in M. edulis, and higher levels of Cd, Hg, Zn in M galloprovincialis. This study demonstrated that variations in contaminant bioaccumulation and different biomarker responses exist between Mytilus species in the field. Data on species or the presence of hybrid individuals (or introgression) is an important additional parameter to add to biomonitoring programs databases.

Patterns of spatial genetic structuring in the endangered limpet Patella ferruginea: implications for the conservation of a Mediterranean endemic.

Patella ferruginea Gmelin, 1791 is an endangered marine gastropod endemic to the Western Mediterranean. Its range is restricted to the Sardinian-Corsican region (SCR), North Africa, a few scattered sites in Southern Spain, and Sicily. Inter-simple sequence repeat (ISSR) markers and three different mitochondrial DNA (mtDNA) regions, Cytochrome c Oxidase subunit I, 12S (small-subunit ribosomal RNA gene) and 16S (large-subunit ribosomal RNA gene), were used to investigate the presence of genetic population structuring. The mtDNA sequences showed very low levels of genetic differentiation. Conversely, ISSRs showed the presence of two main genetic groups, corresponding to Spain, North Africa and Sicily and the SCR. The SCR was further split into two subgroups. The ISSR results suggest that, on a regional scale, the genetic structure of P. ferruginea is mainly determined by the restriction of gene flow by dispersal barriers. On a more local scale human harvesting may play a crucial role in population structuring by increasing the effect of genetic drift.

Realistic environmental exposure to microplastics does not induce biological effects in the Pacific oyster Crassostrea gigas.

The aim of the present study was to evaluate the presence and potential toxic effects of plastic fragments (<400 µm) of polyethylene and polypropylene on the Pacific oyster Crassostrea gigas. Oysters were exposed to environmentally relevant concentrations (0, 0.008, 10, 100 µg of particles/L) during 10 days, followed by a depuration period of 10 days in clean seawater. Effects of microplastics were evaluated on the clearance rate of organisms, tissue alteration, antioxidant defense, immune alteration and DNA damage. Detection and quantification of microplastics in oyster's tissues (digestive gland, gills and other tissues) and biodeposits using infrared microscopy were also conducted. Microplastics were detected in oyster's biodeposits following exposure to all tested concentrations: 0.003, 0.006 and 0.05 particles/mg of biodeposits in oysters exposed to 0.008, 10 and 100 µg of particles/L, respectively. No significant modulation of biological markers was measured in organisms exposed to microplastics in environmentally relevant conditions.

Replicated anthropogenic hybridisations reveal parallel patterns of admixture in marine mussels.

Human-mediated transport creates secondary contacts between genetically differentiated lineages, bringing new opportunities for gene exchange. When similar introductions occur in different places, they provide informally replicated experiments for studying hybridisation. We here examined 4,279 Mytilus mussels, sampled in Europe and genotyped with 77 ancestry-informative markers. We identified a type of introduced mussels, called "dock mussels," associated with port habitats and displaying a particular genetic signal of admixture between M. edulis and the Mediterranean lineage of M. galloprovincialis. These mussels exhibit similarities in their ancestry compositions, regardless of the local native genetic backgrounds and the distance separating colonised ports. We observed fine-scale genetic shifts at the port entrance, at scales below natural dispersal distance. Such sharp clines do not fit with migration-selection tension zone models, and instead suggest habitat choice and early-stage adaptation to the port environment, possibly coupled with connectivity barriers. Variations in the spread and admixture patterns of dock mussels seem to be influenced by the local native genetic backgrounds encountered. We next examined departures from the average admixture rate at different loci, and compared human-mediated admixture events, to naturally admixed populations and experimental crosses. When the same M. galloprovincialis background was involved, positive correlations in the departures of loci across locations were found; but when different backgrounds were involved, no or negative correlations were observed. While some observed positive correlations might be best explained by a shared history and saltatory colonisation, others are likely produced by parallel selective events. Altogether, genome-wide effect of admixture seems repeatable and more dependent on genetic background than environmental context. Our results pave the way towards further genomic analyses of admixture, and monitoring of the spread of dock mussels both at large and at fine spacial scales.

Copper induces expression and methylation changes of early development genes in Crassostrea gigas embryos.

Copper contamination is widespread along coastal areas and exerts adverse effects on marine organisms such as mollusks. In the Pacific oyster, copper induces severe developmental abnormalities during early life stages; however, the underlying molecular mechanisms are largely unknown. This study aims to better understand whether the embryotoxic effects of copper in Crassostrea gigas could be mediated by alterations in gene expression, and the putative role of DNA methylation, which is known to contribute to gene regulation in early embryo development. For that purpose, oyster embryos were exposed to 4 nominal copper concentrations (0.1, 1, 10 and 20 µg L-1 Cu2+) during early development assays. Embryotoxicity was monitored through the oyster embryo-larval bioassay at the D-larva stage 24 h post fertilization (hpf) and genotoxicity at gastrulation 7 hpf. In parallel, the relative expression of 15 genes encoding putative homeotic, biomineralization and DNA methylation proteins was measured at three developmental stages (3 hpf morula stage, 7 hpf gastrula stage, 24 hpf D-larvae stage) using RT-qPCR. Global DNA content in methylcytosine and hydroxymethylcytosine were measured by HPLC and gene-specific DNA methylation levels were monitored using MeDIP-qPCR. A significant increase in larval abnormalities was observed from copper concentrations of 10 µg L-1, while significant genotoxic effects were detected at 1 µg L-1 and above. All the selected genes presented a stage-dependent expression pattern, which was impaired for some homeobox and DNA methylation genes (Notochord, HOXA1, HOX2, Lox5, DNMT3b and CXXC-1) after copper exposure. While global DNA methylation (5-methylcytosine) at gastrula stage didn't show significant changes between experimental conditions, 5-hydroxymethylcytosine, its degradation product, decreased upon copper treatment. The DNA methylation of exons and the transcript levels were correlated in control samples for HOXA1 but such a correlation was diminished following copper exposure. The methylation level of some specific gene regions (HoxA1, Hox2, Engrailed2 and Notochord) displayed changes upon copper exposure. Such changes were gene and exon-specific and no obvious global trends could be identified. Our study suggests that the embryotoxic effects of copper in oysters could involve homeotic gene expression impairment possibly by changing DNA methylation levels.

Effects of a parental exposure to diuron on Pacific oyster spat methylome.

Environmental epigenetic is an emerging field that studies the cause-effect relationship between environmental factors and heritable trait via an alteration in epigenetic marks. This field has received much attentions since the impact of environmental factors on different epigenetic marks have been shown to be associated with a broad range of phenotypic disorders in natural ecosystems. Chemical pollutants have been shown to affect immediate epigenetic information carriers of several aquatic species but the heritability of the chromatin marks and the consequences for long term adaptation remain open questions. In this work, we investigated the impact of the diuron herbicide on the DNA methylation pattern of spat from exposed Crassotrea gigas genitors. This oyster is one of the most important mollusk species produced worldwide and a key coastal economic resource in France. The whole genome bisulfite sequencing (WGBS, BS-Seq) was applied to obtain a methylome at single nucleotide resolution on DNA extracted from spat issued from diuron exposed genitors comparatively to control spat. We showed that the parental diuron exposure has an impact on the DNA methylation pattern of its progeny. Most of the differentially methylated regions occurred within coding sequences and we showed that this change in methylation level correlates with RNA level only in a very small group of genes. Although the DNA methylation profile is variable between individuals, we showed conserved DNA methylation patterns in response to parental diuron exposure. This relevant result opens perspectives for the setting of new markers based on epimutations as early indicators of marine pollutions.

Microalgal sensitivity varies between a diuron-resistant strain and two wild strains when exposed to diuron and irgarol, alone and in mixtures.

A wild strain of Chaetoceros calcitrans and wild and diuron-resistant strains of Tetraselmis suecica, were exposed to the PSII inhibitor herbicides diuron and irgarol, individually and in mixtures. The effects of three concentrations of diuron and irgarol and four binary mixtures were evaluated on doubling time, relative reactive oxygen species and lipid content by flow cytometry, and on photosynthetic efficiency by pulse amplitude modulated fluorescence. In both wild strains, significant effects were observed for each molecule at the highest concentration tested: at irgarol 0.5 µg L(-1), C. calcitrans was shown to be more sensitive than T. suecica (+52% and +19% in doubling time, respectively), whereas at diuron 5 µg L(-1), T. suecica was more affected (+125% in doubling time) than C. calcitrans (+21%). Overall, irgarol had a higher toxicity at a lower concentration than diuron (no effect at diuron 0.5 µg L(-1)) for both wild strains. The strongest mixture (irgarol 0.5 µg L(-1) + diuron 5 µg L(-1)) increased doubling time by 356% for T. suecica, thus showing amplified effects when the two compounds were mixed. Sequencing of the diuron-resistant strain demonstrated a single mutation in the psbA gene coding sequence. Although resistance of this strain to diuron was confirmed with no effect at the highest diuron concentration, no resistance to irgarol was shown. In addition, the mutant strain exposed to the strongest mixture showed a 3.5-fold increase in doubling time compared with irgarol alone, thereby supporting the hypothesis of a biochemical interaction between these two compounds.

Exposure of marine mussels Mytilus spp. to polystyrene microplastics: Toxicity and influence on fluoranthene bioaccumulation.

The effects of polystyrene microbeads (micro-PS; mix of 2 and 6 µm; final concentration: 32 µg L(-1)) alone or in combination with fluoranthene (30 µg L(-1)) on marine mussels Mytilus spp. were investigated after 7 days of exposure and 7 days of depuration under controlled laboratory conditions. Overall, fluoranthene was mostly associated to algae Chaetoceros muelleri (partition coefficient Log Kp = 4.8) used as a food source for mussels during the experiment. When micro-PS were added in the system, a fraction of FLU transferred from the algae to the microbeads as suggested by the higher partition coefficient of micro-PS (Log Kp = 6.6), which confirmed a high affinity of fluoranthene for polystyrene microparticles. However, this did not lead to a modification of fluoranthene bioaccumulation in exposed individuals, suggesting that micro-PS had a minor role in transferring fluoranthene to mussels tissues in comparison with waterborne and foodborne exposures. After depuration, a higher fluoranthene concentration was detected in mussels exposed to micro-PS and fluoranthene, as compared to mussels exposed to fluoranthene alone. This may be related to direct effect of micro-PS on detoxification mechanisms, as suggested by a down regulation of a P-glycoprotein involved in pollutant excretion, but other factors such as an impairment of the filtration activity or presence of remaining beads in the gut cannot be excluded. Micro-PS alone led to an increase in hemocyte mortality and triggered substantial modulation of cellular oxidative balance: increase in reactive oxygen species production in hemocytes and enhancement of anti-oxidant and glutathione-related enzymes in mussel tissues. Highest histopathological damages and levels of anti-oxidant markers were observed in mussels exposed to micro-PS together with fluoranthene. Overall these results suggest that under the experimental conditions of our study micro-PS led to direct toxic effects at tissue, cellular and molecular levels, and modulated fluoranthene kinetics and toxicity in marine mussels.

Regulation of a truncated isoform of AMP-activated protein kinase α (AMPKα) in response to hypoxia in the muscle of Pacific oyster Crassostrea gigas.

AMP-activated protein kinase α (AMPKα) is a key regulator of energy balance in many model species during hypoxia. In a marine bivalve, the Pacific oyster Crassostrea gigas, we analyzed the protein content of adductor muscle in response to hypoxia during 6 h. In both smooth and striated muscles, the amount of full-length AMP-activated protein kinase α (AMPKα) remained unchanged during hypoxia. However, hypoxia induced a rapid and muscle-specific response concerning truncated isoforms of AMPKα. In the smooth muscle, a truncated isoform of AMPKα was increased from 1 to 6 h of hypoxia, and was linked with accumulation of AKT kinase, a key enzyme of the insulin signaling pathway which controls intracellular glucose metabolism. In this muscle, aerobic metabolism was maintained over the 6 h of hypoxia, as mitochondrial citrate synthase activity remained constant. In contrast, in striated muscle, hypoxia did not induce any significant modification of neither truncated AMPKα nor AKT protein content, and citrate synthase activity was altered after 6 h of hypoxia. Together, our results demonstrate that hypoxia response is specific to muscle type in Pacific oyster, and that truncated AMPKα and AKT proteins might be involved in maintaining aerobic metabolism in smooth muscle. Such regulation might occur in vivo during tidal intervals that cause up to 6 h of hypoxia.

Transcriptomic response of the Pacific oyster Crassostrea gigas to hypoxia.

Marine intertidal organisms commonly face hypoxic stress during low tide emersion; moreover, eutrophic conditions and sediment nearness could lead to hypoxic phenomena; it is indeed important to understand the molecular processes involved in the response to hypoxia. In this study the molecular response of the Pacific oyster Crassostrea gigas to prolonged hypoxia (2mg O(2) L(-1) for 20d) was investigated under experimental conditions. A transcriptomic approach was employed using a cDNA microarray of 9058 C. gigas clones to highlight the genetic expression patterns of the Pacific oyster under hypoxic conditions. Lines of oysters resistant (R) and susceptible (S) to summer mortality were used in this study. ANOVA analysis was used to identify the genes involved in the response to hypoxia in comparison to normoxic conditions. The hypoxic response was maximal at day 20. The principal biological processes up-regulated by hypoxic stress were antioxidant defense and the respiratory chain compartment, suggesting oxidative stress caused by hypoxia or an anticipatory response for normoxic recovery. This is the first study employing microarrays to characterize the genetic markers and metabolic pathways responding to hypoxic stress in C. gigas.