Metabolomic analysis reveals a differential adaptation process of the larval stages of Anisakis simplex to the host environment
Iwona Polak1 
Robert Stryiński
Marta Majewska
Elżbieta Łopieńska-Biernat- 1Department of Biochemistry, Faculty of Biology and Biotechnology, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
- 2Department of Human Physiology and Pathophysiology, School of Medicine, Collegium Medicum, University of Warmia and Mazury in Olsztyn, Olsztyn, Poland
Introduction: Anisakis simplex are parasitic nematodes that cause anisakiasis. The possibility of infection with this parasite is through consumption of raw or undercooked fish products. A. simplex infections are often misdiagnosed, especially in subclinical cases that do not present with typical symptoms such as urticaria, angioedema, and gastrointestinal allergy. The resulting allergic reactions range from rapid-onset and potentially fatal anaphylactic reactions to chronic, debilitating conditions. While there have been numerous published studies on the genomes and proteomes of A. simplex, less attention has been paid to the metabolomes. Metabolomics is concerned with the composition of metabolites in biological systems. Dynamic responses to endogenous and exogenous stimuli are particularly well suited for the study of holistic metabolic responses. In addition, metabolomics can be used to determine metabolic activity at different stages of development or during growth.
Materials and methods: In this study, we reveal for the first time the metabolomes of infectious stages (L3 and L4) of A. simplex using untargeted metabolomics by ultra-performance liquid chromatography-mass spectrometry.
Results: In the negative ionization mode (ESI-), we identified 172 different compounds, whereas in the positive ionization mode (ESI+), 186 metabolites were found. Statistical analysis showed that 60 metabolites were found in the ESI- mode with different concentration in each group, of which 21 were more enriched in the L3 larvae and 39 in the L4 stage of A. simplex. Comparison of the individual developmental stages in the ESI + mode also revealed a total of 60 differential metabolites, but 32 metabolites were more enriched in the L3 stage larvae, and 28 metabolites were more concentrated in the L4 stage.
Discussion: The metabolomics study revealed that the developmental stages of A. simplex differed in a number of metabolic pathways, including nicotinate and nicotinamide metabolism. In addition, molecules responsible for successful migration within their host, such as pyridoxine and prostaglandins (E1, E2, F1a) were present in the L4 stage. In contrast, metabolic pathways for amino acids, starch, and sucrose were mainly activated in the L3 stage. Our results provide new insights into the comparative metabolome profiles of two different developmental stages of A. simplex.
1 Introduction
Anisakis simplex is one of the most important emerging parasitic nematodes in Europe, according to the Risk Management Ranking of Foodborne Parasites, prepared for recommendations by the Food and Agriculture Organization of the United Nations (FAO) and the World Health Organization (WHO) (Stryiński et al., 2020). The disease caused by the genus Anisakis is called anisakiasis and frequently described (Audicana et al., 2002; Aibinu et al., 2019). The life cycle of A. simplex is complex and involves four larval stages parasitizing several intermediate and paratenic hosts (fish, cephalopods, and crustaceans) and the adult stage parasitizing marine mammals (seals, dolphins, and whales). The first case of human infection by a species of the Anisakidae family occurred in the 1960s, when Van Thiel of the Institute of Tropical Medicine in Leiden, the Netherlands, identified patients who suffered severe abdominal pain after eating fish (van Thiel et al., 1960). Transmission of the parasite is clearly associated with the consumption of raw or undercooked fish. In particular, Japanese sushi and sashimi, Dutch salted or smoked herring, Nordic gravlax (dry, cured salmon), Hawaiian lomi-lomi (raw salmon), German rollmops (rolled fillet of marinated/pickled herring), South American cebiche, and Spanish boquerones en vinagre (pickled anchovies) are regular routes of infection (Audicana et al., 2002; Ivanovic et al., 2015). The ingestion of viable larvae might lead to gastrointestinal symptoms (abdominal pain, nausea, vomiting, diarrhea), which may be associated with mild to severe allergic reactions, and the clinical symptoms most often are such as rhinitis, urticaria, and, in worst cases, anaphylactic shock (Aibinu et al., 2019). Although, cooking (or freezing) is expected to kill the parasites, it might not decrease its allergenicity, because A. simplex allergens have high heat and frost resistance; and sensitization may occur after consumption (Audicana and Kennedy, 2008). It is estimated that the incidence of the disease is 0.32 cases per 100,000 individuals/year worldwide (Orpha.net, 2022). The globalization, development of diagnostic tools and better analytical methods, have led to significantly even more anisakiasis being reported. Before 2010, over 20,000 cases of anisakiasis were reported worldwide, with the highest prevalence (over 90%) in Japan (EFSA, 2010; Baird et al., 2014), where now 7,000 cases of the disease are reported annually (Yorimitsu et al., 2013; Suzuki et al., 2021). European countries where cases of anisakiasis have been reported include among others, Spain (Herrador et al., 2019), Italy (Guardone et al., 2018; Mattiucci et al., 2018), France (Yera et al., 2018), Croatia (Mladineo et al., 2015), and Poland (Kołodziejczyk et al., 2020). The ingested L3 larvae rarely develop to L4 stage in humans and consequently die, but both of these larval stages are considered dangerous to humans (Kagei and Isogaki, 1992; Sohn et al., 2015).
The infective L3 of parasitic nematodes require specific stimuli for resumption of development and completion of shedding of the outer cuticle; thus, CO2 appears to be the most important stimulus for hatching or molting of Anisakis nematodes. In addition, temperature, pH, and pepsin are known to trigger the development of infectious L3 to adults in vitro (Iglesias et al., 2001, 2005). All these stimuli are thought to reflect the conditions that prevail in vivo when they reach the gastrointestinal tract of the definitive hosts, where they develop into adults. During this molting process, they adapt to the new environment of their hosts’ digestive tracts and also exhibit pathogenicity toward their hosts. These differences may reflect metabolic adaptations of A. simplex larvae to the host switch from fish (L3) to mammals (L4), i.e., adaptations to a new habitat (Baird et al., 2014; Kuhn et al., 2016).
Metabolomics is a distinct ‘omics’ method that offers a more direct assessment of physiology compared to others. It responds promptly to nutrients, stress, and disease unlike transcriptomic or proteomic approaches. Given this advantage, metabolomics has gained considerable interest in various fields such as environmental toxicology, evolutionary biology and developmental studies, medical diagnosis and treatment responses as well as drug development (Vernocchi et al., 2016). Synthetic biologists also use metabolomic flux analyses for insights into the effect of genetic modifications on metabolic pathways and products. As it pertains directly to molecular response mechanisms resulting from genetic alterations or environmental changes at the ultimate level of biological systems’ metabolism function regarding molecule abundance predictions offer better accuracy than gene expression or protein-level information provides (Liu et al., 2021). Similarly, parasitic nematode research increasingly relies on metabolomics techniques due to its immense potential applications in studying these organisms effectively. The gastrointestinal tract is a dynamic metabolic and immunologically active ecosystem, and its complete set of metabolites reflects both the enzymatic pathways of host and gut inhabitants and the complex network that connects them (Palevich et al., 2021, 2022; Whitman et al., 2021). The use of metabolomics is particularly suited to understanding nematode metabolism, including the identification of novel drug targets, differences between developmental stages, and mechanisms involved in host-parasite interactions (Jasmer et al., 2019; Molenaars et al., 2021).
A major problem with untargeted metabolomics in general is the lack of comprehensive measurements of the whole metabolome using a single technique because of its great complexity. As things stand, it is not possible to measure the entire metabolome of an organism in a single experiment. In addition to some known compounds, a large portion of the metabolome has not even been identified (“metabolic dark matter”) (da Silva et al., 2015). Metabolomics is capable of providing a static snapshot of the current metabolic state. However, in most cases, it remains unclear how this state was reached, and which metabolic pathway was active (Salzer and Witting, 2021). With respect to parasitic nematodes, the field of metabolomics faces similar challenges to untargeted metabolomics in general, as its complex nature makes it difficult to obtain complete measurements of the entire metabolome using a single approach. Nevertheless, the use of metabolomics in the study of parasitic nematodes is a promising area that could provide important insights into their biology and pathogenesis.
Such studies on the metabolome of nematodes from the family Anisakidae are lacking. Therefore, to fill the described gap, it was decided to analyze and characterize the metabolome of A. simplex, a parasite of great public health importance. In the central part of this work, we focused on the identification and characterization of the metabolites of A. simplex in two developmental stages to increase our knowledge of the biology of this organism and to find a way to understand how this parasite was able to adapt to different host environments by developing a unique metabolism.
2 Materials and methods
2.1 Parasites and in vitro culture
The study was performed on the L3 and L4 larval stages of Anisakis simplex Nematodes were collected from the Biobank platform implemented for the PARASITE project (www.parasite-project.eu) at Institute of Marine Research, Spanish National Research Council (IIM-CSIC), Vigo, Spain. The larval stages of A. simplex were isolated: L3 from hake (Merluccius) and L4 from striped dolphin (Stenella coeruleoalba). All used larvae were taxonomically identified using conventional PCR to amplify the ITS region and Cox2 gene as described before by Levsen et al. (2018). Eight samples of L3 larvae (24 larvae in total) and four samples of L4 larvae (8 larvae in total) were preserved in −80°C until the time of next step of the analysis.
2.2 Samples preparation
Samples were thawed and add with 800 μL of 80% methanol. Then all samples were extracted at 4°C with ultrasound for 30 min, kept at −40°C for 1 h. After that, samples were vortexed for 30 s, and centrifuged at 12,000 rpm and 4°C for 15 min. Finally, 200 μL of supernatant and 5 μL of DL-o-chlorophenylalanine (140 μg/ml) was transferred to vial for LC-MS/MS analysis. Quality control (QC) samples were used to evaluate the methodology. The same amount of extract was obtained from each sample and mixed as QC samples. The QC samples (3 in total) were prepared using the same sample preparation procedure.
2.3 UPLC-TOF-MS/MS
The separation of compounds was performed by ultra-performance liquid chromatography coupled with tandem mass spectrometry using Ultimate 3000LC combined with Q Exactive MS (Thermo Fisher Scientific, Waltham, MA, United States) in both polarities (ESI −/+). The LC system was comprised of a Hyper gold C18 column (100 × 2.1 mm 1.9 μm) (Thermo Fisher Scientific, Waltham, MA, United States). The mobile phase was composed of solvent A (0.1% formic acid, 5% acetonitrile, HPLC-grade water) and solvent B (0.1% formic acid, acetonitrile) with a gradient elution (0–1.5 min, 0%–20% B; 1.5–9.5 min, 20%–100% B; 9.5–14.5 min, 100% B; 14.5–14.6 min, 100%–0% B; 14.6–18.0 min, 0% B). The flow rate of the mobile phase was 0.3 ml/min. The column temperature was maintained at 40°C, and the sample manager temperature was set at 4°C. Mass spectrometry parameters in ESI+ and ESI- mode are listed as follows: the sheath gas flow rate was set to 45 a. u. (arbitrary units), aux gas flow rate to 15 a. u., and sweep gas flow rate to 1 a. u., capillary temperature was 350°C, and probe heater temperature was 300°C. Electrospray ionization source was set to 3.0 kV in ESI+ and 3.2 kV in ESI-.
2.4 Data processing and analysis
The raw LC-MS/MS data were acquired and aligned using the Compound Discover (v. 3.0, Thermo Fisher Scientific, Waltham, MA, United States) based on the m/z value and the retention time (RT) of the ion signals. The spectral data were normalized, and auto scaled before statistical analysis. The data was introduced into the SIMCA-P software (version 14.1, Santorius, Goettingen, Germany) for multivariate analysis. A Principal Components Analysis (PCA) was first used as an unsupervised method for data visualization and outlier identification (Yeshi et al., 2020). Supervised regression modeling was performed on the data set by use of Partial Least Squares Discriminant Analysis (PLS-DA). The quality of the models was evaluated with the relevant R2 and Q2. The importance of each ion in the PLS-DA was evaluated by variable importance in the projection (VIP) score (Tambellini et al., 2013; Gao et al., 2016). The VIP score positively reflects the metabolite’s influence on the classification, and metabolites with VIP >1 were considered important in the study.
The chemical structures and IDs for metabolites were identified according to Human Metabolome Database (Wishart et al., 2012), KEGG database (Kanehisa and Goto, 2000; Kanehisa et al., 2016), PubChem Compound ID database (Kim et al., 2022) and ChemSpider databse (https://www.chemspider.com/Default.aspx). When necessary, further identification was performed through comparisons of the retention times and MS/MS fragmentation patterns in other databases: Metabolite and Chemical Entity Database (METLIN) (Smith et al., 2005) or MassBank (https://massbank.eu/MassBank/) (Salek et al., 2013).
The metabolites showing different concentrations among two given groups (DMs) were filtered and confirmed by combining the results of the multivariate analysis (VIP values >1.0) and the results of univariate analyses: t-test (p-value ≤0.05) and fold-change (FC) method (-1.5 ⩽ Log2FC ⩽ 1.5) for both ESI modes. Data were visualized on volcano plots.
In multivariate analyses, the hierarchical clustering analysis (HCA) with Euclidean measured distance, and the average clustering algorithm was used to visualize the differences in the concentration of each statistically significant metabolites between groups in two different ESI modes. Subsequently, the metabolites assignment to the compounds groups and pathway enrichment analysis was performed using Metaboanalyst 5.0 (Xia et al., 2009; Pang et al., 2021). The enrichment overview was based on the KEGG database with Caenorhabditis elegans as a reference (Kanehisa and Goto, 2000; Kanehisa et al., 2016).
Furthermore, the distribution of common and unique metabolites identified in the two ESI modes was analyzed and visualized using the Venn Diagrams tool (https://bioinformatics.psb.ugent.be/webtools/Venn/).
3 Results
3.1 Metabolites identification and statistical analysis
The QC samples were used to demonstrate the stability of the LC-MS system. The QC samples run in positive and negative mode at regular intervals throughout the entire sequence. The ion features of the QC samples were used to calculate the relative standard deviation (RSD). The %RSD distribution for negative and positive modes is presented in Supplementary Figures S1A, B, respectively; an overwhelming majority of the RSD is less than 30%. The base peak intensity (BPI) chromatograms of the QC samples in ESI−/+ modes are presented in Supplementary Figures S1C, D. The normalization after alignment was also performed and the line plot (Supplementary Figures S1E, F) was used to evaluate the methodology. The X axis indicates the number of samples, the Y axis indicates the 95% confidence interval. The line plot demonstrated that the system is relatively stable during sample analysis. Due to all of that, the analysis procedure was robust and could be used for subsequent sample analysis.
To investigate the global metabolism variations, the PCA analysis was used to analyze all observations acquired in both ion modes. As show in PCA plot (Figures 1A, B), the QC samples were successfully separated from the tested samples and clustered together (QC, L3 and L4 developmental stages of A. simplex). The parameters R2 and Q2 confirmed the validity of the PCA model as follows: ESI- mode, R2X = 0.676, Q2 = 0.514; ESI + mode, R2X = 0.661, Q2 = 0.499 (Figures 1A, B). To eliminate any non-specific effects of the operative technique and confirm the presence of DMs, the PLS-DA was performed. The PLS-DA score scatter plots showed that there was significant separation between the L3 developmental stage group and the L4 developmental stage group in the ESI- and ESI + modes (Figures 1C, D). The values of the cumulative R2 and Q2 parameters confirmed the validity of the PLS-DA model: ESI- mode, R2X = 0.691, R2Y = 0.995, and Q2 = 0.987; ESI + mode, R2X = 0.683, R2Y = 0.997, and Q2 = 0.990 (Figures 1C, D). Additionally, the CV-ANOVA analysis assessing the reliability of the PLS-DA model was performed (Supplementary Figure S2). The p-value obtained for ESI- and ESI + modes was less than 0.05: 7.3415 × 10−6, and 3.998 × 10−6, respectively. The CV-score plots of the samples in ESI- and ESI + modes showed, as well, clear separation of the two analyzed groups with F factor = 77.64 and 91.61 for negative and positive ionization modes, respectively (Supplementary Figures S2A, B). According to the permutation test results, the PLS-DA model was proved to have good robustness without over fitting (Supplementary Figures S2C, D).
FIGURE 1. The results of multivariate data analysis. The score scatter plots of PCA model (A, B) comparing LC-MS/MS metabolomic profiles for the L3 and L4 developmental stages of A. simplex in negative and positive ionization mode (ESI−/+), respectively, and QC group. The score scatter plots of PLS-DA model (C, D) comparing LC-MS/MS metabolomic profiles for the L3 and L4 developmental stages of A. simplex in negative and positive ionization modes (ESI−/+), respectively. The R2 and Q2 values are indicated in the figures.
As a result of LC-MS/MS analysis, we identified a total of 3603 and 3877 compounds (Supplementary Tables S1, 2) in ESI- and ESI+, respectively. The peak intensities after normalization against QC samples (Supplementary Tables S3, 4) were further processed and the chemical structures and IDs for metabolites were identified (Supplementary Tables S5, 6). In the ESI- we identified 172 different compounds (Supplementary Table S5), when in ESI +, 186 metabolites were found (Supplementary Table S6). It was checked whether the identified compounds overlap and occur similarly in both used polarization modes (Figure 2A). Analysis showed that 94 metabolites are common for both ESI modes (Supplementary Table S6).
FIGURE 2. Visualization of statistically significant data results. Distribution of common and unique metabolites identified in the study between negative and positives ionization modes (ESI−/+) (A). The metabolites IDs are listed in Supplementary Table S7. Distribution of common and unique metabolites for each of A. simplex developmental stage (L3 and L4) identified in the study in negative and positives ionization modes and (ESI−/+) (B). The metabolites IDs are listed in Supplementary Table S10. Volcano plots of the untargeted metabolomics analysis in negative and positive ionization modes (ESI−/+), respectively (C, D), of two developmental stages of A. simplex. The metabolites showing different concentrations among the given groups (DMs) were filtered and confirmed by combining the results of the multivariate analysis (VIP values >1.0) and the results of univariate analyses (-1.5 ⩽ Log2FC ⩽ 1.5 in normalized ratios of L3 vs L4, p-value ≤0.05) for both ESI modes. Metabolites more enriched in the L3 stage are coloured in red, and in the L4 stage in blue. The metabolites IDs with calculated p-values, FC and VIP are listed in Supplementary Tables S8, 9. Hierarchical clustering analysis of the two sample groups (L3 and L4 stages of A. simplex) in negative (ESI-) and positive (ESI+) ionization modes, respectively (E, F). The HCA shows discrimination between the samples groups and differential abundances of DMs in ESI- (E) and ESI+ (F). The scale bars represent the normalized intensity of metabolites, where blue indicates a decrease/low and red an increase/high.
These data were further processed, and several filters were next applied to obtain the final list of DMs among the given groups: a) -1.5 ⩽ Log2FC ⩽ 1.5 in normalized ratios (L3 vs L4), b) t-test (p-value ≤0.05), c) VIP >1.0. Afterwards, the significantly changed metabolites between the two groups were filtered out and listed in Supplementary Tables S8, 9. Statistical analysis showed that in ESI- mode was 60 DMs, of which 21 were more detected in the L3 larvae and 39 in the L4 stage of A. simplex (Supplementary Table S8). Comparison of the individual developmental stages in ESI + mode revealed, as well, a total of 60 DMs, however with 32 metabolites more enriched in the L3 stage larvae and 28 metabolites more enriched in the L4 stage (Supplementary Table S9). The distribution of common and unique metabolites identified in the two ESI modes specific for each developmental stage was also analyzed and visualized with the use of Venn diagram (Figure 2B, Supplementary Table S10). It was found that there are 11 common metabolites between the ESI −/+ modes for the L3 stage (l-arginine, tryptophan, threonylphenylalanine, homoserine, leucyl-tryptophan, glutamylmethionine, gamma-glutamylleucine, glutamic acid, l-methionine, leucine, and glycylleucine), and eight for the L4 stage (13-HOTE, LysoPE (18:1 (9Z)/0:0), prostaglandin E2, N-lauroylglycine, glyceric acid, LysoPE (16:0/0:0), p-cresol, and 19,20-DiHDPA) (Figure 2B, Supplementary Table S10). The volcano plot representations of DMs are shown in Figures 2C, D. Among the DMs more enriched for L3 stage larvae, found in ESI-, we identified gallic acid and trehalose, as those with the highest fold change. Moreover, in the group of metabolites more enriched in L3 we found, e.g., nicotinic acid, leucine, tryptophan or gamma-glutamylleucine (Figure 2C; Table 1). In the studied polarization (ESI-), such compounds as pimelic acid, sebacic acid, benzoic acid, or prostaglandins E1, E2, F1a were found in increased concentration in L4 larvae compared to L3 (Figure 2C; Table 1). The analysis of the A. simplex metabolome in ESI + mode showed higher concentration in L3 compared to L4 of such metabolites as: glutamylmethionine, serine, piperidine or 3-hydroxyanthranilic acid (Figure 2D; Table 2). Moreover, in ESI + mode, increased amounts of e.g., pyridoxine, psi-pelletierine, styrene, citrulline or N-laurolylglicine were observed in L4 larvae in relation to the L3 stage (Figure 2D; Table 2). The HCA was used to visualize the differences in the concentration of each statistically significant metabolite between L3 and L4 developmental stages of A. simplex in two different ESI modes (Figures 2E,F). An overview on the DMs concentrations from all samples can be observed and clear clustering and grouping trend between L3 and L4 developmental stages of A. simplex is showed.
TABLE 1. Metabolites of developmental stages of A. simplex (L3 and L4) identified in negative ionization mode (ESI-).
TABLE 2. Metabolites of developmental stages of A. simplex (L3 and L4) identified in positive ionization mode (ESI+).
3.2 Metabolite classes and enriched pathways
The DMs identified in ESI- and ESI+ (60 in each mode) were assigned to different metabolite classes/groups (Supplementary Tables S11, 12). The highest number of metabolites identified in ESI- mode were categorized into the fatty acyls (e.g., prostaglandin E2, prostaglandin E1, prostaglandin F1a, 8-HETE, 13-HOTE), organic acids (e.g., sebacic acid, pimelic acid, suberic acid, methylsuccinic acid), and benzenoids (p-cresol, benzoic acid, alpha-hydroxyhippuric acid) (Figure 3A, Supplementary Table S11). Of 60 DMs identified in ESI + mode, 23 of them were organic acids (e.g., ureidopropionic acid, citric acid, l-malic acid), nine were organoheterocyclic compounds (dihydrothymine, indoleacetic acid, pyridoxine, serotonin, pyroglutamic acid, indoleacrylic acid, nicotinic acid, tryptophanol, piperidine), 10 were fatty acyls and benzenoids, and four were glycerophospholipids (e.g., LysoPC(18:3 (9Z,12Z, 15Z)); LysoPE (16:0/0:0); LysoPE (18:1 (9Z)/0:0); LysoPC(P-18:0)) (Figure 3B, Supplementary Table S12).
FIGURE 3. Metabolite classes and enriched pathways. Distribution of DMs by different metabolite classes identified in negative and positive ionization modes, respectively (A, B). The metabolites names assigned to different classes are listed in Supplementary Tables S11, 12. Dot plots of pathway enrichment analysis of DMs identified in negative and positive ionization modes, respectively (C, D). The scale bars represent p-value and enrichment ratio. The pathways names with calculated p-value are listed in Supplementary Tables S13, 14.
Based on the DMs identified in both developmental stages of A. simplex in ESI−/+ ionization modes, the metabolic pathways were mapped, and enrichment analysis was performed (p-value ≤0.05) (Figures 3C, D; Supplementary Tables S13, 14). In negative ionization mode (ESI-) the highest number of identified DMs were products of ubiquinone and other terpenoid-quinone biosynthesis (KEGG 00130), followed by phenylalanine, tyrosine and tryptophan biosynthesis (KEGG 00400), fatty acid biosynthesis (KEGG 00061), valine, leucine and isoleucine biosynthesis and degradation (KEGG 00290, 00280), tyrosine metabolism (KEGG 00350), starch and sucrose metabolism (KEGG 00500), and nicotinate and nicotinamide metabolism (KEGG 00760) (Kanehisa and Goto, 2000; Kanehisa et al., 2016). Mapping and enrichment analysis of pathways for identified DMs in positive ionization mode (ESI+) against known pathways in KEGG database, showed that highest number of DMs were products of lysine degradation (KEGG 00310), and nicotinate and nicotinamide metabolism (KEGG 00760) (Kanehisa and Goto, 2000; Kanehisa et al., 2016).
4 Discussion
One of the most recent additions to the Caenorhabditis elegans toolbox is metabolomics and lipidomics, which allow new and deeper investigations of nematode metabolism. Combining a genetically defensible model organism such as C. elegans with the functional evaluation of metabolomics and/or lipidomics holds great promise for expanding our knowledge of metabolism and metabolic regulation (Salzer and Witting, 2021). While there are numerous published studies on the genomes (Mattiucci et al., 2016; Coghlan et al., 2019; Łopieńska-Biernat et al., 2019), transcriptomes (Cavallero et al., 2018; Kim et al., 2018; Llorens et al., 2018; Łopieńska-Biernat et al., 2018; Cavallero et al., 2020), and proteomes (Stryiński et al., 2019; Stryiński et al., 2022; Polak et al., 2020; Mierzejewski et al., 2021; Kochanowski et al., 2022) of A. simplex, less attention has been paid to its metabolome. Metabolic changes between L3 and L4 larval stages of A. simplex have been studied previously, however in a narrow range, e.g., of one metabolic pathway (Łopieńska-Biernat et al., 2008; Łopieńska-Biernat et al., 2018; Łopieńska-Biernat et al., 2019), and almost nothing is known about the complement of an intermediate or end product of metabolism in Anisakis nematodes, especially in the stages found in humans (L3 and L4).
In the first stage of infection the larvae are invading the host intestinal tissues and establishing in the gastrointestinal tract of the host. Such tissue migration requires larvae to adapt to constantly changing external environments, in terms of temperature, oxygen supply, redox potential and host immune reactions (Huang et al., 2020; Yeshi et al., 2020). In this metabolomic study, we observed substantial differences in the composition and abundance of trehalose, nicotinic acid, prostaglandins, and acyl acids, likely associated with key biological functions, including energy metabolism and nicotinate and nicotinamide metabolism during the parasite’s growth and development. Alterations in metabolomic profiles are likely reflected in the adaptation of A. simplex nematode to changing environments within the host (poikilotherm and homeotherm organisms). Nevertheless, comparative analyses of levels of metabolites in the L3 and L4 of these two larval stages led to the identification of groups of molecules with putative roles in mechanisms of parasite pathogenicity.
Our results revealed that specific to L3 are nicotinic acid and branched-chain and aromatic amino acids (BCAA). The recent metabolomics studies showed that BCAA are positively related to longevity common to Dauer in C. elegans and Haemonhus contortus during mitochondrial biogenesis in ensheathment eggs (Palevich et al., 2022). In this study, BCAA were only identified in L3 stage, which like a Dauer stage, uses endogenous energy sources. Starvation-induced stress may invoke a series of regulations of metabolic pathways (Fuchs et al., 2010; Łopieńska-Biernat et al., 2019). In this study, some pathways, including fatty acid metabolism, nicotinate/nicotinamide metabolism, terpenoid backbone biosynthesis were only identified in L3 stage.
Fatty acids serve as energy storage and structural components in biomembranes. The main chemical composition of the cuticular lipids in A. simplex from Atlantic cod (Godus morhua) was recognized as three fatty acids, 15 triacylglycerols, five sterols and 12 sphingolipids (Mika et al., 2010) what is consisted with our current results. Moreover, we have identified fatty acyls and glycerophospholipids in L4 stage larvae.
Terpenoids are the precursor in the production of steroids and sterols. Therefore, the constant activity of the terpenoid backbone biosynthesis pathway is important for the regulation of steroids metabolism, as they regulate a variety of developmental and biological processes (molting, larval development, and innate immunity) (Hernández-Bello et al., 2011). The ubiquinone/terpenoid quinone biosynthesis (energy metabolism) involved in biosynthesis of other secondary metabolites, and metabolism of cofactors (nicotinate/nicotinamide metabolism) and vitamins is specific to two larval stages. Based on current results and previous studies, energy metabolism of A. simplex is based on carbohydrates metabolism, where L3 larvae is metabolizing trehalose, and L4 larvae glycogen (Łopieńska-Biernat et al., 2006, Łopieńska-Biernat et al., 2007; Łopieńska-Biernat et al., 2008; Łopieńska-Biernat et al., 2018; Łopieńska-Biernat et al., 2019). However, based on the current results we can suppose that L4 larvae energy metabolism is also related to fatty acids (Ma et al., 2019). This could be an adaptation due to different host environments and usage of energy sources.
In addition, we have identified significant amounts of the prostaglandin PGE2 in L4 larvae of A. simplex, which helps larvae successfully migrate into their host and acts as a mediator of host immunity. These molecules mediate and have an immunosuppressive effect on T cells (Tebbey and Buttke, 1990). Genes of the prostaglandin synthesis cascade (Fitzpatrick, 2004), including cyclooxygenase (WormBase ParaSite ASIM_0001219401), alpha-tubulin N-acetyltransferaseare (WormBase ParaSite ASIM_0001506401), and phospholipase A2 (WormBase ParaSite ASIM_0001246201) are present in the genome of the nematode studied (Coghlan et al., 2019), demonstrating that A. simplex larvae can synthesize prostaglandins. Moreover, studies on Trichuris suis may provide insight into the mechanisms by which the worm suppresses inflammatory host responses, with one active component identified as prostaglandin PGE2 (Laan et al., 2017).
Lipidome analysis of Schistosoma mansoni has detected eicosanoids, signaling molecules made by the enzymatic or non-enzymatic oxidation of polyunsaturated fatty acids, that promote a Th2 response (Shepherd et al., 2015; Whitman et al., 2021). Metabolites analysis in this study specific to L4 stage also revealed eicosanoids. However, studies show that some helminth-derived small molecules have distinct immunomodulatory components, suggesting that further studies are needed at this host-pathogen interface (Shepherd et al., 2015; Whitman et al., 2021).
Our metabolomics study showed that A. simplex induced metabolic changes in a variety of metabolic pathways in both larval stages, including amino acid metabolism, phospholipid metabolism, energy metabolism, nicotinate/nicotinamide metabolism and ubiquinone/terpenoid quinone biosynthesis. In addition, several stage-specific metabolites were identified, providing potential clues for understanding the molecular mechanisms of this parasite biology: its pathogenicity and adaptation to the host environment.
By undertaking the first small-scale analysis of the metabolome molecules present in the L3 and L4 of A. simplex that are putatively involved in the host environment, we provide here a ready-to-use molecular groundwork for in-depth studies of the biological pathways specifically involved in parasite growth and development. Further studies under experimental, controlled in vitro conditions, as well as between other developmental stages of the life cycle, are needed to reliably assess the role of these molecules in the pathogenesis of anisakiasis. A snapshot of the metabolic status in specific age of the larvae and/or tissues of larval stages of A. simplex could provide new information on the expression of target metabolites potentially involved in host invasion. Furthermore, 13C flux experiments can provide adequate information on the route of formation of key metabolites.
Gastrointestinal nematodes strategy mediated through the carbohydrates and lipid metabolites and employed for the manipulation of the host immune response towards successful propagation and parasitism, opens a wide perspective that should be compounded by the contemporary use of multi-omics approaches.
Data availability statement
The raw data generated for this study can be found in the EMBL-EBI MetaboLights database (www.ebi.ac.uk/metabolights/MTBLS6458) with the identifier MTBLS6458.
Author contributions
Conceptualization: RS and EŁ-B. Methodology: IP and RS. Formal analysis: IP. Data curation: RS. Software: RS. Writing—original draft preparation: IP, RS, and EŁ-B. Writing—review and editing: RS, MM, and EŁ-B. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by National Science Centre of Poland, grant no. 2020/37/N/NZ9/03312.
Acknowledgments
We thank PARASITE Biobank at the Institute of Marine Research (IIM-CSIC), Vigo, Spain for providing the L4 larvae of Anisakis simplex used in this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1233586/full#supplementary-material
References
Aibinu, I. E., Smooker, P. M., and Lopata, A. L. (2019). Anisakis nematodes in fish and shellfish-from infection to allergies. Int. J. Parasitol. Parasites Wildl. 9, 384–393. doi:10.1016/j.ijppaw.2019.04.007
Audicana, M. T., and Kennedy, M. W. (2008). Anisakis simplex: From obscure infectious worm to inducer of immune hypersensitivity. Clin. Microbiol. Rev. 21, 360–379. doi:10.1128/CMR.00012-07
Audicana, M. T., Ansotegui, I. J., de Corres, L. F., and Kennedy, M. W. (2002). Anisakis simplex: Dangerous — Dead and alive? Trends Parasitol. 18, 20–25. doi:10.1016/S1471-4922(01)02152-3
Baird, F. J., Gasser, R. B., Jabbar, A., and Lopata, A. L. (2014). Foodborne anisakiasis and allergy. Mol. Cell Probes 28, 167–174. doi:10.1016/j.mcp.2014.02.003
Cavallero, S., Lombardo, F., Su, X., Salvemini, M., Cantacessi, C., and D’Amelio, S. (2018). Tissue-specific transcriptomes of Anisakis simplex (sensu stricto) and Anisakis pegreffii reveal potential molecular mechanisms involved in pathogenicity. Parasit. Vectors 11, 31. doi:10.1186/S13071-017-2585-7
Cavallero, S., Lombardo, F., Salvemini, M., Pizzarelli, A., Cantacessi, C., and D’Amelio, S. (2020). Comparative transcriptomics reveals clues for differences in pathogenicity between hysterothylacium aduncum, Anisakis simplex sensu stricto and Anisakis pegreffii. Genes (Basel) 11, 321. doi:10.3390/genes11030321
Coghlan, A., Tyagi, R., Cotton, J. A., Holroyd, N., Rosa, B. A., Tsai, I. J., et al. (2019). Comparative genomics of the major parasitic worms. Nat. Genet. 51, 163–174. doi:10.1038/s41588-018-0262-1
da Silva, R. R., Dorrestein, P. C., and Quinn, R. A. (2015). Illuminating the dark matter in metabolomics. Proc. Natl. Acad. Sci. 112, 12549–12550. doi:10.1073/pnas.1516878112
EFSA (2010). Scientific Opinion on risk assessment of parasites in fishery products. EFSA J. 8, 1543. doi:10.2903/j.efsa.2010.1543
Fitzpatrick, F. (2004). Cyclooxygenase enzymes: Regulation and function. Curr. Pharm. Des. 10, 577–588. doi:10.2174/1381612043453144
Fuchs, S., Bundy, J. G., Davies, S. K., Viney, J. M., Swire, J. S., and Leroi, A. M. (2010). A metabolic signature of long life in Caenorhabditis elegans. BMC Biol. 8, 14. doi:10.1186/1741-7007-8-14
Gao, X., Zhang, W., Wang, Y., Pedram, P., Cahill, F., Zhai, G., et al. (2016). Serum metabolic biomarkers distinguish metabolically healthy peripherally obese from unhealthy centrally obese individuals. Nutr. Metab. (Lond) 13, 33. doi:10.1186/s12986-016-0095-9
Guardone, L., Armani, A., Nucera, D., Costanzo, F., Mattiucci, S., and Bruschi, F. (2018). Human anisakiasis in Italy: A retrospective epidemiological study over two decades. Parasite 25, 41. doi:10.1051/parasite/2018034
Hernández-Bello, R., Ramirez-Nieto, R., Muñiz-Hernández, S., Nava-Castro, K., Pavón, L., Sánchez-Acosta, A. G., et al. (2011). Sex steroids effects on the molting process of the helminth human parasite Trichinella spiralis. J. Biomed. Biotechnol. 2011, 625380–625410. doi:10.1155/2011/625380
Herrador, Z., Daschner, Á., Perteguer, M. J., and Benito, A. (2019). Epidemiological scenario of anisakidosis in Spain based on associated hospitalizations: The tip of the iceberg. Clin. Infect. Dis. 69, 69–76. doi:10.1093/cid/ciy853
Huang, Y., Wu, Q., Zhao, L., Xiong, C., Xu, Y., Dong, X., et al. (2020). UHPLC-MS-Based metabolomics analysis reveals the process of schistosomiasis in mice. Front. Microbiol. 11, 1517. doi:10.3389/fmicb.2020.01517
Iglesias, L., Valero, A., Benitez, R., and Adroher, F. J. (2001). In vitro cultivation of Anisakis simplex: Pepsin increases survival and moulting from fourth larval to adult stage. Parasitology 123, 285–291. doi:10.1017/S0031182001008423
Iglesias, L., Malagón, D., Valero, A., Benítez, R., and Javier Adroher, F. (2005). CO2-fixing enzymes during moulting from third larval to fourth larval stage of Anisakis simplex and Hysterothylacium aduncum (Nematoda: Anisakidae). Parasitol. Res. 96, 212–215. doi:10.1007/s00436-005-1342-6
Ivanovic, J., Baltic, M. Z., Boskovic, M., Kilibarda, N., Dokmanovic, M., Markovic, R., et al. (2015). Anisakis infection and allergy in humans. Procedia Food Sci. 5, 101–104. doi:10.1016/j.profoo.2015.09.028
Jasmer, D. P., Rosa, B. A., Tyagi, R., and Mitreva, M. (2019). Omics driven understanding of the intestines of parasitic nematodes. Front. Genet. 10, 652. doi:10.3389/fgene.2019.00652
Kagei, N., and Isogaki, H. (1992). A case of abdominal syndrome caused by the presence of a large number of Anisakis larvae. Int. J. Parasitol. 22, 251–253. doi:10.1016/0020-7519(92)90111-W
Kanehisa, M., and Goto, S. (2000). Kegg: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi:10.1093/nar/28.1.27
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462. doi:10.1093/nar/gkv1070
Kim, J. H., Kim, J. O., Jeon, C. H., Nam, U. H., Subramaniyam, S., Yoo, S. I., et al. (2018). Comparative transcriptome analyses of the third and fourth stage larvae of Anisakis simplex (Nematoda: Anisakidae). Mol. Biochem. Parasitol. 226, 24–33. doi:10.1016/j.molbiopara.2018.10.002
Kim, S., Cheng, T., He, S., Thiessen, P. A., Li, Q., Gindulyte, A., et al. (2022). PubChem protein, gene, pathway, and taxonomy data collections: Bridging biology and chemistry through target-centric views of PubChem data. J. Mol. Biol. 434, 167514. doi:10.1016/j.jmb.2022.167514
Kochanowski, M., Dąbrowska, J., Różycki, M., Sroka, J., Karamon, J., Bełcik, A., et al. (2022). Proteomic profiling and in silico characterization of the secretome of Anisakis simplex sensu stricto L3 larvae. Pathogens 11, 246. doi:10.3390/pathogens11020246
Kołodziejczyk, L., Szostakowska, B., Sobecka, E., Szczucki, K., and Stankiewicz, K. (2020). First case of human anisakiasis in Poland. Parasitol. Int. 76, 102073. doi:10.1016/j.parint.2020.102073
Kuhn, T., Cunze, S., Kochmann, J., and Klimpel, S. (2016). Environmental variables and definitive host distribution: A habitat suitability modelling for endohelminth parasites in the marine realm. Sci. Rep. 6, 30246. doi:10.1038/srep30246
Laan, L. C., Williams, A. R., Stavenhagen, K., Giera, M., Kooij, G., Vlasakov, I., et al. (2017). The whipworm (Trichuris suis) secretes prostaglandin E2 to suppress proinflammatory properties in human dendritic cells. FASEB J. 31, 719–731. doi:10.1096/fj.201600841R
Levsen, A., Svanevik, C. S., Cipriani, P., Mattiucci, S., Gay, M., Hastie, L. C., et al. (2018). A survey of zoonotic nematodes of commercial key fish species from major European fishing grounds—introducing the FP7 PARASITE exposure assessment study. Fish. Res. 202, 4–21. doi:10.1016/j.fishres.2017.09.009
Liu, R., Bao, Z.-X., Zhao, P.-J., and Li, G.-H. (2021). Advances in the study of metabolomics and metabolites in some species interactions. Molecules 26, 3311. doi:10.3390/molecules26113311
Llorens, C., Arcos, S. C., Robertson, L., Ramos, R., Futami, R., Soriano, B., et al. (2018). Functional insights into the infective larval stage of Anisakis simplex s.s., Anisakis pegreffii and their hybrids based on gene expression patterns. BMC Genomics 19, 592. doi:10.1186/s12864-018-4970-9
Łopieńska-Biernat, E., Żółtowska, K., and Rokicki, J. (2006). The content of carbohydrates in larval stages of Anisakis simplex (Nematoda, Anisakidae). Helminthologia 43, 125–129. doi:10.2478/s11687-006-0024-y
Łopieńska-Biernat, E., Żółtowska, K., and Rokicki, J. (2007). Trehalose catabolism enzymes in L3 and L4 larvae of Anisakis simplex. J. Parasitol. 93, 1291–1294. doi:10.1645/GE-906.1
Łopieńska-Biernat, E., Żółtowska, K., and Rokicki, J. (2008). Glycogen catabolism enzymes and protein fractions in the third and fourth larval stages of Anisakis simplex. J. Helminthol. 82, 45–51. doi:10.1017/S0022149X0787355X
Łopieńska-Biernat, E., Molcan, T., Paukszto, Ł., Jastrzębski, J. P., and Myszczyński, K. (2018). Modelling studies determing the mode of action of anthelmintics inhibiting in vitro trehalose-6-phosphate phosphatase (TPP) of Anisakis simplex s.l. Exp. Parasitol. 184, 46–56. doi:10.1016/j.exppara.2017.11.006
Łopieńska-Biernat, E., Stryiński, R., Dmitryjuk, M., and Wasilewska, B. (2019). Infective larvae of Anisakis simplex (Nematoda) accumulate trehalose and glycogen in response to starvation and temperature stress. Biol. Open 8, bio040014. doi:10.1242/bio.040014
Ma, G., Wang, T., Korhonen, P. K., Young, N. D., Nie, S., Ang, C.-S., et al. (2019). Dafachronic acid promotes larval development in Haemonchus contortus by modulating dauer signalling and lipid metabolism. PLoS Pathog. 15, e1007960. doi:10.1371/journal.ppat.1007960
Mattiucci, S., Acerra, V., Paoletti, M., Cipriani, P., Levsen, A., Webb, S. C., et al. (2016). No more time to stay ‘single’ in the detection of Anisakis pegreffii, A. simplex (s. S.) and hybridization events between them: A multi-marker nuclear genotyping approach Parasitology 143, 998–1011. doi:10.1017/S0031182016000330
Mattiucci, S., Cipriani, P., Levsen, A., Paoletti, M., and Nascetti, G. (2018). Molecular epidemiology of Anisakis and anisakiasis: An ecological and evolutionary road map. Adv. Parasitol. 99, 93–263. doi:10.1016/bs.apar.2017.12.001
Mierzejewski, K., Stryiński, R., Łopieńska-Biernat, E., Mateos, J., Bogacka, I., and Carrera, M. (2021). A complex proteomic response of the parasitic nematode Anisakis simplex s.s. to Escherichia coli lipopolysaccharide. Mol. Cell. Proteomics 20, 100166. doi:10.1016/j.mcpro.2021.100166
Mika, A., Gołębiowski, M., Szafranek, J., Rokicki, J., and Stepnowski, P. (2010). Identification of lipids in the cuticle of the parasitic nematode Anisakis simplex and the somatic tissues of the Atlantic cod Gadus morhua. Exp. Parasitol. 124, 334–340. doi:10.1016/j.exppara.2009.11.013
Mladineo, I., Popović, M., Drmić-Hofman, I., and Poljak, V. (2015). A case report of Anisakis pegreffii (Nematoda, Anisakidae) identified from archival paraffin sections of a Croatian patient. BMC Infect. Dis. 16, 42. doi:10.1186/s12879-016-1401-x
Molenaars, M., Schomakers, B. V., Elfrink, H. L., Gao, A. W., Vervaart, M. A. T., Pras-Raves, M. L., et al. (2021). Metabolomics and lipidomics in Caenorhabditis elegans using a single-sample preparation. Dis. Model Mech. 14, dmm047746. doi:10.1242/dmm.047746
Orpha.net (2022). Prevalence of rare diseases: Bibliographic data, orphanet report series, rare diseases collection, january 2022, number 1: Diseases listed in alphabetical order. Available at: https://www.orpha.net/orphacom/cahiers/docs/GB/Prevalence_of_rare_diseases_by_alphabetical_list.pdf.
Palevich, N., Maclean, P. H., Subbaraj, A. K., and Cao, M. (2021). A multimodal metabolomics approach to elucidate the trigger components of rumen fluid for larval exsheathment in a model gastrointestinal parasitic nematode, Haemonchus contortus. bioRxiv, 1–44. doi:10.1101/2021.01.10.426113
Palevich, N., Maclean, P. H., Candy, P. M., Taylor, W., Mladineo, I., and Cao, M. (2022). Untargeted multimodal metabolomics investigation of the Haemonchus contortus exsheathment secretome. Cells 11, 2525. doi:10.3390/cells11162525
Pang, Z., Chong, J., Zhou, G., de Lima Morais, D. A., Chang, L., Barrette, M., et al. (2021). MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 49, W388–W396. doi:10.1093/nar/gkab382
Polak, I., Łopieńska-Biernat, E., Stryiński, R., Mateos, J., and Carrera, M. (2020). Comparative proteomics analysis of Anisakis simplex s.s.—evaluation of the response of invasive larvae to ivermectin. Genes (Basel) 11, 710. doi:10.3390/genes11060710
Salek, R. M., Steinbeck, C., Viant, M. R., Goodacre, R., and Dunn, W. B. (2013). The role of reporting standards for metabolite annotation and identification in metabolomic studies. Gigascience 2, 13. doi:10.1186/2047-217X-2-13
Salzer, L., and Witting, M. (2021). Quo vadis Caenorhabditis elegans metabolomics—a review of current methods and applications to explore metabolism in the nematode. Metabolites 11, 284. doi:10.3390/metabo11050284
Shepherd, C., Navarro, S., Wangchuk, P., Wilson, D., Daly, N. L., and Loukas, A. (2015). Identifying the immunomodulatory components of helminths. Parasite Immunol. 37, 293–303. doi:10.1111/pim.12192
Smith, C. A., Maille, G. O., Want, E. J., Qin, C., Trauger, S. A., Brandon, T. R., et al. (2005). Metlin: A metabolite mass spectral database. Ther. Drug Monit. 27, 747–751. doi:10.1097/01.ftd.0000179845.53213.39
Sohn, W.-M., Na, B.-K., Kim, T. H., and Park, T.-J. (2015). Anisakiasis: Report of 15 gastric cases caused by Anisakis type I larvae and a brief review of Korean anisakiasis cases. Korean J. Parasitol. 53, 465–470. doi:10.3347/kjp.2015.53.4.465
Stryiński, R., Mateos, J., Pascual, S., González, Á. F., Gallardo, J. M., Łopieńska-Biernat, E., et al. (2019). Proteome profiling of L3 and L4 Anisakis simplex development stages by TMT-based quantitative proteomics. J. Proteomics 201, 1–11. doi:10.1016/j.jprot.2019.04.006
Stryiński, R., Łopieńska-Biernat, E., and Carrera, M. (2020). Proteomic insights into the biology of the most important Foodborne parasites in Europe. Foods 9, 1403. doi:10.3390/foods9101403
Stryiński, R., Mateos, J., Carrera, M., Jastrzębski, J. P., Bogacka, I., and Łopieńska-Biernat, E. (2022). Tandem mass tagging (TMT) reveals tissue-specific proteome of L4 larvae of Anisakis simplex s. s.: Enzymes of energy and/or carbohydrate metabolism as potential drug targets in anisakiasis. Int. J. Mol. Sci. 23, 4336. doi:10.3390/ijms23084336
Suzuki, J., Murata, R., and Kodo, Y. (2021). Current status of anisakiasis and Anisakis larvae in tokyo, Japan. Food Saf. 9, 89–100. doi:10.14252/foodsafetyfscj.D-21-00004
Tambellini, N., Zaremberg, V., Turner, R., and Weljie, A. (2013). Evaluation of extraction protocols for simultaneous polar and non-polar yeast metabolite analysis using multivariate projection methods. Metabolites 3, 592–605. doi:10.3390/metabo3030592
Tebbey, P. W., and Buttke, T. M. (1990). Molecular basis for the immunosuppressive action of stearic acid on T cells. Immunology 70, 379–384.
van Thiel, P., Kuipers, F. C., and Roskam, R. T. (1960). A nematode parasitic to herring, causing acute abdominal syndromes in man. Trop. Geogr. Med. 12, 97–113.
Vernocchi, P., Del Chierico, F., and Putignani, L. (2016). Gut microbiota profiling: Metabolomics based approach to unravel compounds affecting human health. Front. Microbiol. 7, 1144. doi:10.3389/fmicb.2016.01144
Whitman, J. D., Sakanari, J. A., and Mitreva, M. (2021). Areas of metabolomic exploration for helminth infections. ACS Infect. Dis. 7, 206–214. doi:10.1021/acsinfecdis.0c00658
Wishart, D. S., Jewison, T., Guo, A. C., Wilson, M., Knox, C., Liu, Y., et al. (2012). HMDB 3.0—the human metabolome database in 2013. Nucleic Acids Res. 41, D801–D807. doi:10.1093/nar/gks1065
Xia, J., Psychogios, N., Young, N., and Wishart, D. S. (2009). MetaboAnalyst: A web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 37, W652–W660. doi:10.1093/nar/gkp356
Yera, H., Fréalle, É., Dutoit, E., and Dupouy-Camet, J. (2018). A national retrospective survey of anisakidosis in France (2010-2014): Decreasing incidence, female predominance, and emerging allergic potential. Parasite 25, 23. doi:10.1051/parasite/2018016
Yeshi, K., Creek, D. J., Anderson, D., Ritmejerytė, E., Becker, L., Loukas, A., et al. (2020). Metabolomes and lipidomes of the infective stages of the gastrointestinal nematodes, nippostrongylus brasiliensis and Trichuris muris. Metabolites 10, 446. doi:10.3390/metabo10110446
Keywords: Anisakis simplex, larval development, metabolome, trehalose, fatty acids
Citation: Polak I, Stryiński R, Majewska M and Łopieńska-Biernat E (2023) Metabolomic analysis reveals a differential adaptation process of the larval stages of Anisakis simplex to the host environment. Front. Mol. Biosci. 10:1233586. doi: 10.3389/fmolb.2023.1233586
Received: 02 June 2023; Accepted: 04 July 2023;
Published: 13 July 2023.
Edited by:
Elda Dervishi, University of Alberta, CanadaReviewed by:
Paniz Jasbi, Arizona State University, United StatesDhanasekaran Shanmugam, National Chemical Laboratory (CSIR), India
Copyright © 2023 Polak, Stryiński, Majewska and Łopieńska-Biernat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elżbieta Łopieńska-Biernat, ela.lopienska@uwm.edu.pl