Spatialized Metabolomic Annotation Combining MALDI Imaging and Molecular Networks.

MALDI mass spectrometry imaging has gained major interest in the field of chemical imaging. This technique makes it possible to locate tens to hundreds of ionic signals on the sample surface without any a priori. One of the current challenges is still the limited ability to annotate signals in order to convert m/z values into probable chemical structures. At the same time, data obtained by LC-MS/MS have benefited from the development of numerous chemoinformatics tools, in particular molecular networks, for their efficient annotation. For the first time, we present here the combination of MALDI-FT-ICR imaging with molecular networks from MALDI-MS/MS data directly acquired on plant tissue sections. Annotation improvements are demonstrated, paving the way for new annotation pipelines for MALDI imaging.

Calcium influx and spermatogenesis in the testis and liver enzyme activities in the zebrafish are rapidly modulated by the calcium content of the water.

This study investigated the effects of varying environmental Ca2+ concentrations on the influx of Ca2+ to the testis, testicular morphology, and liver enzymes in the zebrafish. Adult zebrafish (Danio rerio) were held in water containing low (0.02 mM), control (0. 7 mM) or high (2 mM) Ca2+ concentrations for 12 h. Testes were then incubated in vitro with 0.1 µCi/mL 45Ca2+ to measure Ca2+ influx at 30 and 60 min and qualitative and quantitative testicular histological analyses were conducted. In addition, activity of alanine aminotransferase (ALT), aspartate aminotransferase (AST) and gamma-glutamyl transpeptidase (GGT), enzymes that indicate tissue damage, were evaluated in the liver. The testes from zebrafish exposed in vivo to low (0.02 mM) and high (2 mM) Ca2+ content water had a higher Ca2+ influx than the control group after 30 min of incubation, and at 60 min (high Ca2+ group only). There were morphological changes in the testes from the low and high Ca2+ groups including spermatozoa distributed in dense agglomerates and apoptotic cells. Furthermore, zebrafish exposed to high Ca2+ containing water had an increased density of haploid cells (spermatids and spermatozoa). In addition, both low and high Ca2+ water affected liver function by increasing ALT and GGT activities. Collectively, these studies show that alterations in calcium homeostasis in the testis, stimulation of the spermatogenic wave and hepatic injury were rapid responses to changes in the concentration of Ca2+ in the water.

Antiparasitic Ovalicin Derivatives from Pseudallescheria boydii, a Mutualistic Fungus of French Guiana Termites.

Social insects are in mutualism with microorganisms, contributing to their resistance against infectious diseases. The fungus Pseudallescheria boydii SNB-CN85 isolated from termites produces ovalicin derivatives resulting from the esterification of the less hindered site of the ovalicin epoxide by long-chain fatty acids. Their structures were elucidated using spectroscopic analysis and semisynthesis from ovalicin. For ovalicin, these compounds displayed antiprotozoal activities against Plasmodium falciparum and Trypanosoma brucei, with IC50 values of 19.8 and 1.1 µM, respectively, for the most active compound, i.e., ovalicin linoleate. In parallel, metabolomic profiling of a collection of P. boydii strains associated with termites made it possible to highlight this class of compounds together with tyroscherin derivatives in all strains. Finally, the complete genome of P. boydii strains was obtained by sequencing, and the cluster of potential ovalicin and ovalicin biosynthesis genes was annotated. Through these metabolomic and genomic analyses, a new ovalicin derivative named boyden C, in which the 6-membered ring of ovalicin was opened by oxidative cleavage, was isolated and structurally characterized.

How the variability between computer-assisted analysis procedures evaluating immune markers can influence patients' outcome prediction.

Differences between computer-assisted image analysis (CAI) algorithms may cause discrepancies in the identification of immunohistochemically stained immune biomarkers in biopsies of breast cancer patients. These discrepancies have implications for their association with disease outcome. This study aims to compare three CAI procedures (A, B and C) to measure positive marker areas in post-neoadjuvant chemotherapy biopsies of patients with triple-negative breast cancer (TNBC) and to explore the differences in their performance in determining the potential association with relapse in these patients. A total of 3304 digital images of biopsy tissue obtained from 118 TNBC patients were stained for seven immune markers using immunohistochemistry (CD4, CD8, FOXP3, CD21, CD1a, CD83, HLA-DR) and were analyzed with procedures A, B and C. The three methods measure the positive pixel markers in the total tissue areas. The extent of agreement between paired CAI procedures, a principal component analysis (PCA) and Cox multivariate analysis was assessed. Comparisons of paired procedures showed close agreement for most of the immune markers at low concentration. The probability of differences between the paired procedures B/C and B/A was generally higher than those observed in C/A. The principal component analysis, largely based on data from CD8, CD1a and HLA-DR, identified two groups of patients with a significantly lower probability of relapse than the others. The multivariate regression models showed similarities in the factors associated with relapse for procedures A and C, as opposed to those obtained with procedure B. General agreement among the results of CAI procedures would not guarantee that the same predictive breast cancer markers were consistently identified. These results highlight the importance of developing additional strategies to improve the sensitivity of CAI procedures.

Structural revision of the Mcl-1 inhibitor MIM1: synthesis and biological studies on ovarian cancer cells with evaluation of designed analogues.

In the area of cancer research, the development of new and potent inhibitors of anti-apoptotic proteins is a very active and promising topic. The small molecule MIM1 has been reported earlier as one of the first selective inhibitors of the anti-apoptotic protein Mcl-1. In the present paper, we first revised the structure of this molecule based on extensive physicochemical analyses. Then we designed and synthesized a focused library of analogues for the corrected structure of MIM1. Next, these molecules were subjected to a panel of in cellulo biological studies, allowing the identification of dual Bcl-xL/Mcl-1 inhibitors, as well as selective Mcl-1 inhibitors. These results have been complemented by fluorescence polarization assays with the Mcl-1 protein. Preliminary structure-activity relationships were discussed and extensive molecular modelling studies allowed us to propose a rationale for the biological activity of this series of new inhibitors, in particular for the selectivity of inhibition of Mcl-1 versus Bcl-xL.

Centrosome, the Newly Identified Passenger through Tunneling Nanotubes, Increases Binucleation and Proliferation Marker in Receiving Cells.

Type 1 tunneling nanotubes (TNTs-1) are long, cytoplasmic protrusions containing actin, microtubules and intermediate filaments that provide a bi-directional road for the transport of various components between distant cells. TNT-1 formation is accompanied by dramatic cytoskeletal reorganization offering mechanical support for intercellular communication. Although the centrosome is the major microtubule nucleating center and also a signaling hub, the relationship between the centrosome and TNTs-1 is still unexplored. We provide here the first evidence of centrosome localization and orientation towards the TNTs-1 protrusion site, which is implicated in TNT-1 formation. We also envision a model whereby synchronized reorientation of the Golgi apparatus along with the centrosome towards TNTs-1 ensures effective polarized trafficking through TNTs-1. Furthermore, using immunohistochemistry and live imaging, we observed for the first time the movement of an extra centrosome within TNTs-1. In this regard, we hypothesize a novel role for TNTs-1 as a critical pathway serving to displace extra centrosomes and potentially to either protect malignant cells against aberrant centrosome amplification or contribute to altering cells in the tumor environment. Indeed, we have observed the increase in binucleation and proliferation markers in receiving cells. The fact that the centrosome can be both as the base and the user of TNTs-1 offers new perspectives and new opportunities to follow in order to improve our knowledge of the pathophysiological mechanisms under TNT control.

Multimodal Imaging Mass Spectrometry to Identify Markers of Pulmonary Arterial Hypertension in Human Lung Tissue Using MALDI-ToF, ToF-SIMS, and Hybrid SIMS.

Pulmonary arterial hypertension (PAH) is a rare and deadly disease affecting roughly 15-60 people per million in Europe with a poorly understood pathology. There are currently no diagnostic tools for early detection nor does a curative treatment exist. The lipid composition of arteries in lung tissue samples from human PAH and control patients were investigated using matrix-assisted laser desorption ionization (MALDI) imaging mass spectrometry (IMS) combined with time-of-flight secondary ion mass spectrometry (TOF-SIMS) imaging. Using random forests as an IMS data analysis technique, it was possible to identify the ion at m/z 885.6 as a marker of PAH in human lung tissue. The m/z 885.6 ion intensity was shown to be significantly higher around diseased arteries and was confirmed to be a diacylglycerophosphoinositol PI(C18:0/C20:4) via MS/MS using a novel hybrid SIMS instrument. The discovery of a potential biomarker opens up new research avenues which may finally lead to a better understanding of the PAH pathology and highlights the vital role IMS can play in modern biomedical research.

Generation of a Molecular Network from Electron Ionization Mass Spectrometry Data by Combining MZmine2 and MetGem Software.

Molecular networking (MN) allows one to organize tandem mass spectrometry (MS/MS) data by spectral similarities. Cosine-score used as a metric to calculate the distance between two spectra is based on peak lists containing fragments and neutral losses from MS/MS spectra. Until now, the workflow excluded the generation of the molecular network from electron ionization (EI) MS data as no selection of the putative parent ion is achieved when performing classical gas chromatography (GC)-EI-MS analysis. In order to fill this gap, new functionalities on MetGem 1.2.2 software ( https://github.com/metgem/metgem/releases ) have been implemented, and results from a large EI-MS database and GC-EI-MS analysis will be exemplified.

Do nitrogen- and sulphur-remobilization-related parameters measured at the onset of the reproductive stage provide early indicators to adjust N and S fertilization in oilseed rape (Brassica napus L.) grown under N- and/or S-limiting supplies?

MAIN CONCLUSION: Specific combinations of physiological and molecular parameters associated with N and S remobilization measured at the onset of flowering were predictive of final crop performances in oilseed rape. Oilseed rape (Brassica napus L.) is a high nitrogen (N) and sulphur (S) demanding crop. Nitrogen- and S-remobilization processes allow N and S requirements to reproductive organs to be satisfied when natural uptake is reduced, thus ensuring high yield and seed quality. The quantification of physiological and molecular indicators of early N and S remobilization could be used as management tools to correct N and S fertilization. However, the major limit of this corrective strategy is to ensure the correlation between final performances-related variables and early measured parameters. In our study, four genotypes of winter oilseed rape (OSR) were grown until seed maturity under four nutritional modalities combining high and/or low N and S supplies. Plant final performances, i.e., seed production, N- and S-harvest indexes, seed N and S use efficiencies, and early parameters related to N- or S-remobilization processes, i.e., photosynthetic leaf area, N and S leaf concentrations, leaf soluble protein and leaf sulphate concentrations, and leaf RuBisCO abundance at flowering, were measured. We demonstrated that contrasting final performances existed according to the N and S supplies. An optimal N:S ratio supply could explain the treatment-specific crop performances, thus justifying N and S concurrent managements. Specific combinations of early measured plant parameters could be used to predict final performances irrespective of the nutritional supply and the genotype. This work demonstrates the potential of physiological and molecular indicators measured at flowering to reflect the functioning of N- and S-compound remobilization and to predict yield and quality penalties. However, because the predictive models are N and S independent, instant N and S leaf analyses are required to further adjust the adequate fertilization. This study is a proof of a concept which opens prospects regarding instant diagnostic tools in the context of N and S mineral fertilization management.

Morphological and molecular criteria allow the identification of putative germ stem cells in a lophotrochozoan, the Pacific oyster Crassostrea gigas.

While our knowledge of bivalve gametogenesis recently progressed, data on early stages of gametogenesis remain to be developed, especially when dealing with germinal stem cells (GSC) and their niche in these organisms. Here, we wish to develop a strategy to identify putative GSC in Pacific oyster Crassostrea gigas based on morphological criteria combined with vasa marker expression. A histological quantitative approach, based on stereology, allowed us to identify two types of early germ cells in the germinal epithelium, one presenting round nuclei and the other irregular ones. Both early germ cell types present slightly condensed chromatin in nucleus, are vasa-positive and the Oyvlg (oyster vasa-like gene) expression in these cells is recorded throughout the whole gametogenesis process. The microenvironment of an early germ cell in oyster includes an associated somatic cell presenting an immunolabeling for BMP2/4 and a close myoid cell. In agreement with the GSC characteristics in other species, we postulate that putative germ stem cells in C. gigas correspond to the early germ cell type with irregular nucleus shape; those early germ cells with a round nucleus may consist in progenitors.

Structural Identification of Antibacterial Lipids from Amazonian Palm Tree Endophytes through the Molecular Network Approach.

A library of 197 endophytic fungi and bacteria isolated from the Amazonian palm tree Astrocaryum sciophilum was extracted and screened for antibacterial activity against methicillin-resistant Staphylococcus aureus (MRSA). Four out of five antibacterial ethyl acetate extracts were also cytotoxic for the MRC-5 cells line. Liquid chromatography coupled to tandem mass spectrometry (UPHLC-HRMS/MS) analyses combined with molecular networking data processing were carried out to allow the identification of depsipeptides and cyclopeptides responsible for the cytotoxicity in the dataset. Specific ion clusters from the active Luteibacter sp. extract were also highlighted using an MRSA activity filter. A chemical study of Luteibacter sp. was conducted leading to the structural characterization of eight fatty acid exhibiting antimicrobial activity against MRSA in the tens of µg/mL range.

MetGem Software for the Generation of Molecular Networks Based on the t-SNE Algorithm.

Molecular networking (MN) is becoming a standard bioinformatics tool in the metabolomic community. Its paradigm is based on the observation that compounds with a high degree of chemical similarity share comparable MS2 fragmentation pathways. To afford a clear separation between MS2 spectral clusters, only the most relevant similarity scores are selected using dedicated filtering steps requiring time-consuming parameter optimization. Depending on the filtering values selected, some scores are arbitrarily deleted and a part of the information is ignored. The problem of creating a reliable representation of MS2 spectra data sets can be solved using algorithms developed for dimensionality reduction and pattern recognition purposes, such as t-distributed stochastic neighbor embedding (t-SNE). This multivariate embedding method pays particular attention to local details by using nonlinear outputs to represent the entire data space. To overcome the limitations inherent to the GNPS workflow and the networking architecture, we developed MetGem. Our software allows the parallel investigation of two complementary representations of the raw data set, one based on a classic GNPS-style MN and another based on the t-SNE algorithm. The t-SNE graph preserves the interactions between related groups of spectra, while the MN output allows an unambiguous separation of clusters. Additionally, almost all parameters can be tuned in real time, and new networks can be generated within a few seconds for small data sets. With the development of this unified interface ( https://metgem.github.io ), we fulfilled the need for a dedicated, user-friendly, local software for MS2 comparison and spectral network generation.

Insights into the MALDI Process after Matrix Deposition by Sublimation Using 3D ToF-SIMS Imaging.

Imaging mass spectrometry (IMS) has become a powerful tool to characterize the spatial distribution of biomolecules in thin tissue sections. In the case of matrix-assisted laser desorption ionization (MALDI) IMS, homogeneous matrix deposition is critical to produce high-quality ion images, and sublimation in particular has shown to be an excellent matrix deposition method for the imaging of lipids. Matrix deposition by sublimation is, however, a completely solvent-free system, which ought to prevent the mixing of matrix and analytes thought to be necessary for successful MALDI. Using 3D time-of-flight secondary ion imaging mass spectrometry, we have studied the matrix-tissue interface in 3D with high resolution to understand the MALDI process of lipids after matrix deposition by sublimation. There is a strong indication that diffusion is the process by which lipids migrate from the tissue to the matrix layer. We show that triacylglycerols and phospholipids have a delayed migratory trend as compared to diacylglycerols and monoacylglycerols, which is dependent on time and matrix thickness. Additional experiments show that a pure lipid's capacity to migrate into the matrix is dependent on its fluidity at room temperature. Furthermore, it is shown that cholesterol can only migrate in the presence of a (fluid) lipid and appears to fluidize lipids, which could explain its colocalization with the diacylglycerols and monoacylglycerols in the matrix.

A role for RASSF1A in tunneling nanotube formation between cells through GEFH1/Rab11 pathway control.

BACKGROUND: By allowing intercellular communication between cells, tunneling nanotubes (TNTs) could play critical role in cancer progression. If TNT formation is known to require cytoskeleton remodeling, key mechanism controlling their formation remains poorly understood. METHODS: The cells of human bronchial (HBEC-3, A549) or mesothelial (H2452, H28) lines are transfected with different siRNAs (inactive, anti-RASSF1A, anti-GEFH1 and / or anti-Rab11). At 48 h post-transfection, i) the number and length of the nanotubes per cell are quantified, ii) the organelles, previously labeled with specific tracers, exchanged via these structures are monitored in real time between cells cultured in 2D or 3D and in normoxia, hypoxia or in serum deprivation condition. RESULTS: We report that RASSF1A, a key-regulator of cytoskeleton encoded by a tumor-suppressor gene on 3p chromosome, is involved in TNTs formation in bronchial and pleural cells since controlling proper activity of RhoB guanine nucleotide exchange factor, GEF-H1. Indeed, the GEF-H1 inactivation induced by RASSF1A silencing, leads to Rab11 accumulation and subsequent exosome releasing, which in turn contribute to TNTs formation. Finally, we provide evidence involving TNT formation in bronchial carcinogenesis, by reporting that hypoxia or nutriment privation, two almost universal conditions in human cancers, fail to prevent TNTs induced by the oncogenic RASSF1A loss of expression. CONCLUSIONS: This finding suggests for the first time that loss of RASSF1A expression could be a potential biomarker for TNTs formation, such TNTs facilitating intercellular communication favoring multistep progression of bronchial epithelial cells toward overt malignancy.

Dual Beam Depth Profiling and Imaging with Argon and Bismuth Clusters of Prenylated Stilbenes on Glandular Trichomes of Macaranga vedeliana.

Using a time-of-flight secondary ion mass spectrometer equipped with an argon cluster ion for sputtering and a bismuth liquid metal ion source for analysis, both surfaces of leaves and fruits of Macaranga vedeliana, an endemic New Caledonian species, have been for the first time analyzed by a dual beam depth profiling. To prevent in-vacuum evaporation of the liquid content of the small glandular trichomes covering fruits and leaves surfaces and also to be able to analyze their liquid content while preventing any sublimation of the latter, the samples were kept frozen during the whole experiment using a nitrogen cooled sample holder. Thus, it was possible to demonstrate that vedelianin, an active metabolite of the family of prenylated stilbenes named schweinfurthins, is only located in these glandular trichomes.

Skin Penetration Enhancement by Natural Oils for Dihydroquercetin Delivery.

Natural oils are commonly used in topical pharmaceutical formulations as emulsifiers, stabilizers or solubility enhancers. They are presented as safe and inert components, mainly used for formulation purposes. It is confirmed that natural oils can affect the skin penetration of various substances. Fatty acids are mainly responsible for this effect. Current understanding lacks reliable scientific data on penetration of natural oils into the skin and their skin penetration enhancement potential. In the current study, fatty acid content analysis was used to determine the principal fatty acids in soybean, olive, avocado, sea-buckthorn pulp, raspberry seed and coconut oils. Time of flight secondary ion mass spectrometry bioimaging was used to determine the distribution of these fatty acids in human skin ex vivo after application of the oils. Skin penetration enhancement ratios were determined for a perspective antioxidant compound dihydroquercetin. The results demonstrated skin penetration of fatty acids from all oils tested. Only soybean and olive oils significantly increased the skin distribution of dihydroquercetin and can be used as skin penetration enhancers. However, no correlation can be determined between the fatty acids' composition and skin penetration enhancement using currently available methodological approaches. This indicates that potential chemical penetration enhancement should be evaluated during formulation of topically applied products containing natural oils.

Design, synthesis of novel pyranotriazolopyrimidines and evaluation of their anti-soybean lipoxygenase, anti-xanthine oxidase, and cytotoxic activities.

The synthesis of 14-(aryl)-14H-naphto[2,1-b]pyrano[3,2-e][1,2,4]triazolo[1,5-c]pyrimidine-2-yl) acetamidoximes 2a-e has been accomplished by reaction of 2-acetonitrile derivatives 1a-e with hydroxylamine. Cyclocondensation reaction of precursors 2a-e with some elctrophilic species such as ethylorthoformate, acetic anhydride, and methyl-acetoacetate provided the new oxadiazole derivatives 3a-e, 4a-e, and 5a-e, respectively. On the other hand, the reaction of precursors 2a-e with 2-chloropropanoyl chloride afforded the new acetimidamides 6a-e which evolve under reflux of toluene to the new oxadiazoles 7a-e. The synthetic compounds were screened for their anti-xanthine oxidase, anti-soybean lipoxygenase, and cytotoxic activities. Moderate to weak xanthine oxidase and soybean lipoxygenase inhibitions were obtained but significant cytotoxic activities were noted. The most cytotoxic activities were recorded mainly (i) 5a was the most active (IC50 = 4.0 µM) and selective against MCF-7 and (ii) 2a was cytotoxic against the four cell lines with selectivity for MCF-7 and OVCAR-3 (IC50 = 17 and 12 µM, respectively) while 2e is highly selective against OVCAR-3 (IC50 = 10 µM).

Time-of-flight secondary ion mass spectrometry imaging of biological samples with delayed extraction for high mass and high spatial resolutions.

RATIONALE: In Time-of-Flight Secondary Ion Mass Spectrometry (TOF-SIMS), pulsed and focused primary ion beams enable mass spectrometry imaging, a method which is particularly useful to map various small molecules such as lipids at the surface of biological samples. When using TOF-SIMS instruments, the focusing modes of the primary ion beam delivered by liquid metal ion guns can provide either a mass resolution of several thousand or a sub-µm lateral resolution, but the combination of both is generally not possible. METHODS: With a TOF-SIMS setup, a delayed extraction applied to secondary ions has been studied extensively on rat cerebellum sections in order to compensate for the effect of long primary ion bunches. RESULTS: The use of a delayed extraction has been proven to be an efficient solution leading to unique features, i.e. a mass resolution up to 10000 at m/z 385.4 combined with a lateral resolution of about 400 nm. Simulations of ion trajectories confirm the experimental determination of optimal delayed extraction and allow understanding of the behavior of ions as a function of their mass-to-charge ratio. CONCLUSIONS: Although the use of a delayed extraction has been well known for many years and is very popular in MALDI, it is much less used in TOF-SIMS. Its full characterization now enables secondary ion images to be recorded in a single run with a submicron spatial resolution and with a mass resolution of several thousand. This improvement is very useful when analyzing lipids on tissue sections, or rare, precious, or very small size samples.

A methodology to ensure and improve accuracy of Ki67 labelling index estimation by automated digital image analysis in breast cancer tissue.

INTRODUCTION: Immunohistochemical Ki67 labelling index (Ki67 LI) reflects proliferative activity and is a potential prognostic/predictive marker of breast cancer. However, its clinical utility is hindered by the lack of standardized measurement methodologies. Besides tissue heterogeneity aspects, the key element of methodology remains accurate estimation of Ki67-stained/counterstained tumour cell profiles. We aimed to develop a methodology to ensure and improve accuracy of the digital image analysis (DIA) approach. METHODS: Tissue microarrays (one 1-mm spot per patient, n = 164) from invasive ductal breast carcinoma were stained for Ki67 and scanned. Criterion standard (Ki67-Count) was obtained by counting positive and negative tumour cell profiles using a stereology grid overlaid on a spot image. DIA was performed with Aperio Genie/Nuclear algorithms. A bias was estimated by ANOVA, correlation and regression analyses. Calibration steps of the DIA by adjusting the algorithm settings were performed: first, by subjective DIA quality assessment (DIA-1), and second, to compensate the bias established (DIA-2). Visual estimate (Ki67-VE) on the same images was performed by five pathologists independently. RESULTS: ANOVA revealed significant underestimation bias (P < 0.05) for DIA-0, DIA-1 and two pathologists' VE, while DIA-2, VE-median and three other VEs were within the same range. Regression analyses revealed best accuracy for the DIA-2 (R-square = 0.90) exceeding that of VE-median, individual VEs and other DIA settings. Bidirectional bias for the DIA-2 with overestimation at low, and underestimation at high ends of the scale was detected. Measurement error correction by inverse regression was applied to improve DIA-2-based prediction of the Ki67-Count, in particularfor the clinically relevant interval of Ki67-Count < 40%. Potential clinical impact of the prediction was tested by dichotomising the cases at the cut-off values of 10, 15, and 20%. Misclassification rate of 5-7% was achieved, compared to that of 11-18% for the VE-median-based prediction. CONCLUSIONS: Our experiments provide methodology to achieve accurate Ki67-LI estimation by DIA, based on proper validation, calibration, and measurement error correction procedures, guided by quantified bias from reference values obtained by stereology grid count. This basic validation step is an important prerequisite for high-throughput automated DIA applications to investigate tissue heterogeneity and clinical utility aspects of Ki67 and other immunohistochemistry (IHC) biomarkers.

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.