Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 15 September 2025

Sec. Parasite and Host

Volume 15 - 2025 | https://doi.org/10.3389/fcimb.2025.1646537

This article is part of the Research TopicHigh throughput approaches to understanding the interactions between host and gastrointestinal parasitesView all articles

Transcriptomic profile of the zoonotic parasite Anisakis pegreffii upon in vitro exposure to human dendritic cells

  • 1Department of Ecological and Biological Sciences, “Tuscia University”, Viterbo, Italy
  • 2Department of Experimental Medicine, Sapienza University of Rome, Rome, Italy
  • 3Department of Biology, Health and Environment, Section of Parasitology, Faculty of Pharmacy and Food Sciences, University of Barcelona, Barcelona, Spain
  • 4Department of Public Health and Infectious Diseases, University Hospital “Policlinico Umberto I”, Sapienza-University of Rome, Rome, Italy

Anisakis pegreffii is a zoonotic marine nematode whose life-cycle involves marine organisms–small crustaceans, fish and squids as intermediate/paratenic hosts, and marine mammals, mainly cetaceans–as definitive ones. When its third-stage larvae (L3) are accidentally ingested by humans with the consumption of raw or undercooked parasitized fish and/or squids, the parasite fails to complete its life cycle, leading to human anisakiasis. Despite increasing interest in understanding the molecular basis of pathogenic effects in human anisakiasis, the transcriptomic response of A. pegreffii L3 to interaction with human immune cells, remains poorly understood. Thus, in this study, the transcriptomic profile of A. pegreffii L3 larvae under in vitro exposure to human dendritic cells (DCs) was performed for the first time. A total of 3914 differentially expressed genes (DEGs) were identified in A. pegreffii L3 after exposure to immature DCs (iDCs), by RNA-seq, allowing to detect 1868 upregulated and 2046 downregulated transcripts. Upregulated genes were significantly enriched in pathways related to energy metabolism, oxidative stress response and structural maintenance, suggesting active metabolic and structural adaptation to immune-induced stress. Conversely, genes involved in cytoskeletal organization and intracellular trafficking were downregulated, potentially reflecting the parasite’s developmental arrest in an unsuitable host such as humans. These findings provide novel insights into the molecular response pathways activated by this zoonotic parasite during the early stages of interaction with the human immune system.

Introduction

Anisakis pegreffii is a zoonotic nematode within the A. simplex (s. l.) complex, with an indirect life cycle involving cetaceans as definitive hosts, and small crustaceans, teleost fish, and cephalopods as intermediate/paratenic hosts (Mattiucci et al., 2018; 2022). Humans become accidental hosts by consuming raw or undercooked fish or squids infected with live third stage larvae (L3) of A. pegreffii (Mattiucci et al., 2018; Guardone et al., 2018). In the human, however, the parasite is unable to complete its development (Trumbić et al., 2021) but it can trigger inflammation leading to acute gastrointestinal symptoms and, if not promptly removed, marked eosinophilic infiltration, with abscess and granuloma formation (Nieuwenhuizen and Lopata, 2013; Mattiucci et al., 2013; Trumbić et al., 2021). Exposure to Anisakis antigens may also elicit IgE hypersensitization, provoking allergic reactions (Nieuwenhuizen and Lopata, 2013; Mattiucci et al., 2013; Brusca et al., 2023; EFSA, 2024). The pathogenicity of A. pegreffii is likely driven by a combination of factors, involving both mechanical invasive capacity and the release of excretory/secretory products (ESPs) (Mehrdana et al., 2017)—some of which are transported via extracellular vesicles (EVs) (Cavallero et al., 2022; Palomba et al., 2023) and miRNA (Cavallero et al., 2022; Bellini et al., 2025) —as well as complex interactions with host immune cells (Maizels and McSorley, 2016). Recent transcriptomic studies in a rodent model have shown that A. pegreffii L3 larvae upregulate genes associated with ribosomal function and oxidative phosphorylation, potentially indicating a failed attempt to molt (Trumbić et al., 2021). In contrast, in natural heterothermic fish hosts, autophagy-related pathways are activated, suggesting a dormancy-like state, highlighting the parasite’s ability to modulate its gene expression in response to diverse host environments (Trumbić et al., 2021).

However, very little is known about how A. pegreffii responds to specific components of the human immune system. Dendritic cells (DCs), as antigen-presenting cells (APCs), play a pivotal role in initiating cellular immunity (Maizels and McSorley, 2016; Lundie et al., 2016; Peng et al., 2022). Generally speaking, upon exposure to helminth-derived products, DCs undergo phenotypic changes that promote the triggering of a Th2-polarized immune response, contributing to immune tolerance; this response is frequently associated with elevated IgE levels, eosinophilia, and mast cell activation (Sher et al., 2003). Previous in vitro studies have shown that co-cultures of DCs with A. pegreffii L3 can impair DC function by reducing their maturation and immunogenic potential (Napoletano et al., 2018). Similarly, antigens from A. simplex (s. s.) larvae have been shown to suppress the expression of costimulatory molecules on DCs while promoting Treg expansion and modulating IL - 10 and IFN-γ production in a host-genotype-dependent manner (Zamora et al., 2019; 2021). However, the reciprocal interaction between the parasite and DCs remains poorly understood—both in terms of how DCs are modulated by the parasite and how the parasite’s gene expression responds to this key immune cell type during their interaction.

Therefore, this study aims to characterize the transcriptomic response of A. pegreffii L3 larvae during in vitro exposure to human DCs, in order to better understand early parasites’ molecular interactions with components of the human immune system.

Materials and methods

Anisakis L3 sampling

Anisakis L3 larvae were carefully removed from the body cavity of European hake (Merluccius merluccius) specimens, caught approximately 12h before, from the Adriatic Sea (off San Benedetto del Tronto coast)–a fishing area with a known high prevalence of Anisakis infection (Cipriani et al., 2018). To control for variability and minimize possible host-related effects, larvae were collected from three different M. merluccius individuals. After their removal, L3 were checked for their integrity and viability under a dissecting microscope. A viable larva in this context was considered physically intact and motile, as measured in terms of its ability to move spontaneously or by stimulation with tweezers and a needle (EFSA, 2024). Alive and not disrupted larvae were washed in a sterile 1X phosphate-buffered saline solution (PBS, Sigma, St Louis, MO) several times and then treated for 1 min with 4% acetic acid (Carlo Erba, Cornaredo, Italy) to inhibit bacterial contamination.

Molecular identification of Anisakis spp.

Anisakis L3 were identified by using a multilocus molecular approach. The mitochondrial cytochrome c oxidase 2 (mtDNA cox2) gene locus was amplified using the primers 211F (forward; 5′-TTTTCTAGTTATATAGATTGRTTYAT-3′) and 210R (reverse; 5′- CACCAACTCTTAAAATTATC - 3′) (Valentini et al., 2006; Mattiucci et al., 2014). The successful PCR products were purified, and Sanger sequenced through an Automated Capillary Electrophoresis Sequencer 3730 DNA Analyzer (Applied Biosystems), using the BigDye® Terminator v3.1 Cycle Sequencing Kit (Life Technologies). Additionally, a direct genotyping determination of the nuclear metallopeptidase 10 gene locus (nas10 nDNA) was performed by the amplification-refractory mutations system (ARMS) PCR assay at nas10 nDNA by the combined use of OUT-F1 (forward; 5’- TATGGCAAATATTATTATCGTA - 3’), OUT-R1 (reverse; 5’-TATTTCCGACAGCAAACAA-3’), INN-F1 (forward; 5’-GCATTGTACACTTCGTATATT-3’), INN-R1 (reverse; 5’-ATTTCTYCAGCAATCGTAAG-3’), following the procedures reported in Palomba et al., 2020. PCR products were separated by electrophoresis using agarose gel (1.5%) stained with GelRed. The distinct banding patterns were detected using ultraviolet transillumination.

Exposure of A. pegreffii L3 to iDCs

For exposure of A. pegreffii L3 to DCs, immature DCs (iDCs) were generated by immunoselected CD14+ monocytes cells obtained from healthy donors, in presence of GM-CSF and IL - 4 as previously described (Dionisi et al., 2018; Napoletano et al., 2007). Each iDCs culture underwent quality control assessment before being used with A. pegreffii (File S1, Supplementary Figure S1). iDCs-L3 larvae co-cultures were established, as previously described (Napoletano et al., 2018). iDCs and L3 larvae were co-cultured in transwell plate (C24 Transwell plate, Corning Costar): L3 larvae were placed in the upper chambers (7 larvae/mL), while iDCs were seeded in the lower chamber (6x105 cells/mL). In detail, a total of 63 L3 larvae were used for the co-culture condition, distributed across nine upper chambers, with 7 larvae per each chamber. Each upper chamber was paired with a corresponding lower chamber containing iDCs at a density of 6×105 cells/mL. To account for biological variability, every three chambers shared iDCs from the same healthy donor, resulting in three experimental groups: Chambers 1-3: Paired with iDCs from Donor 1; Chambers 4-6: Paired with iDCs from Donor 2; Chambers 7-9: Paired with iDCs from Donor 3. As a control, an additional set of nine chambers was prepared under the same conditions but without iDCs, maintaining 7 L3 larvae per chamber. After 24 h, the surviving Anisakis L3 were collected, washed in PBS (twice), and stored in RNA later at -80°C for DNA and RNA extraction.

Ethics statement

All research was conducted in accordance with relevant guidelines and regulations. The study protocol was approved by the Ethics Committee of University Hospital “Policlinico Umberto I” ‘Sapienza’ University of Rome (Protocol nr. 4212). Peripheral blood mononuclear cells (PBMCs) were isolated by Ficoll-Hypaque gradient (1.077 g/mL; Pharmacia LKB, Uppsala, Sweden) from buffy coats of healthy donors obtained from the Transfusion Center, Policlinico Umberto I. Formal written informed consent was obtained from all the healthy blood donors, prior to blood sample collection. Research involving human participants was conducted in compliance with the Declaration of Helsinki.

DNA/RNA extraction

Both DNA and RNA were extracted from each Anisakis larva using TRIzol reagent (Invitrogen, Carlsbad, CA, USA), as previously described (Palomba et al., 2019). Specifically, one larva was collected from each chamber for extraction. This resulted in 9 larvae from the co-culture condition with iDCs (3 larvae for each donor, from different chambers), 9 larvae from the control condition without iDCs (1 larva for each chamber). DNA obtained from each larva was used for species-level identification, while RNA was treated with DNase (DNase I, Invitrogen) to remove any genomic DNA contamination. Then, treated RNA from three larvae (each obtained from a separate chamber with iDCs from different donators) was pooled together. Pooling of larvae was performed to reduce the impact of individual-level transcriptional noise and to obtain a representative average expression profile of A. pegreffii larvae under each experimental condition. This process was repeated three times (3L3/pool, 3 pools in total) to provide triplicate samples for RNA sequencing. RNA concentration, purity and integrity were verified and measured on agarose gel (2%) and by a Bioanalyzer 2100 (Agilent Technologies, Waldbronn, Germany).

Library preparation and RNA sequencing

Following the manufacturer’s protocol, the cDNA library was prepared using a TruSeq Stranded mRNA kit (Illumina, San Diego, USA). In brief, polyadenylated (PolyA+) RNA was purified from 10 μg of total RNA of A. pegreffii L3 using Sera-Mag oligo (dT) beads, fragmented to a length of 100 - 500 nucleotides and reverse transcribed to cDNA using random hexamers. The size-fractionated cDNA was end-repaired and adaptor-ligated according to the manufacturer’s protocol (Illumina). Ligated products of 200 bp were excised from agarose gels and PCR amplified. Products were cleaned using a MinElute PCR purification kit (Qiagen Hilden, Germany) and single-end sequenced on an Illumina HiSeq 2000, according to the manufacturer’s protocol.

Transcript quantification

The raw sequences obtained were processed using Salmon software (v. 1.5.1) (Patro et al., 2017) for mapping against the A. pegreffii de novo reference transcriptome, available online at figshare (https://figshare.com/articles/online_resource/AP_-_Unigenes/18301772) (Palomba et al., 2022). This process was conducted in two distinct phases: 1) the index command of Salmon was used to generate an index of the reference transcriptome; 2) the quant command was used to quantify the abundance of transcripts in the analyzed samples. Both steps were carried out using the software’s default parameters. The transcript quantification phase was performed on a High-Performance Computing Cluster (Castrignanò et al., 2020; Flati et al., 2020). We applied the TransDecoder tool to the Corset-assembled transcriptome, to predict Open Reading Frames, generating peptide (AnisakisPegr:Longest_ORFs_PEP.fasta) and coding sequence (AnisakisPegr:Longest_ORFs_CDS.fasta) files.

Differential gene expression analysis

Differentially expressed genes (DEGs) were identified using IGUANER (Pinna et al., 2024), a software based on DESeq2algorithm (Love et al., 2014). DEGs were defined based on adjusted p-values (padj< 0.05) and log2 fold change thresholds. For the functional annotation of DEGs identified with IGUANER, we adopted two complementary approaches. Initially, we used DIAMOND software (v. 2.0.11) (Buchfink et al., 2014) to conduct sequence comparisons against the Nr, SwissProt, and TrEMBL databases. This step allowed us to obtain highly accurate and relevant homology-based annotations for each analyzed sequence. As second approach, we used Eggnog Mapper (v. 2) (Cantalapiedra et al., 2021) to perform comparisons of sequences with a broad range of homologous groups catalogued in the EggNOG database. This process enabled the acquisition of detailed functional annotations on the basis of the information contained in Gene Ontology (GO) database. A heatmap was generated to visualize the expression patterns of genes based on GO enrichment analysis. Briefly, DEGs from the GO enrichment analysis were filtered to remove duplicates, categorized into up- and down-regulated genes, and merged into a single dataset. The expression data were extracted from the annotated transcriptome file, filtered for significant DEGs (padj< 0.05, |log2FC| > 1), and visualized after eliminating redundant entries.

The complete bioinformatics pipeline is reported in Supplementary Figure S2. All results obtained from the bioinformatics analyses were deposited on Figshare (Table 1). From the complete list of DEGs (padj< 0.05, |log2FC| > 1), genes were selected based on their upregulation and relevance to detoxificant/antioxidant processes as well as chaperone and allergenic properties. This selection was performed manually, guided by existing literature and known gene functions.

Table 1
www.frontiersin.org

Table 1. Overview of produced data files and their access on Figshare.

Quantitative real-time PCR validation

To validate the results of RNA-seq, six differentially expressed genes were randomly selected for quantitative real-time PCR (RT-qPCR). The RNA samples used for the RT-qPCR assay were the same as those used for RNA-seq. The cDNA was synthesized for the RT-qPCR using the high-capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific, Wilmington, DE, USA) following the manufacturer’s instructions. The RT-qPCR was performed using SYBR Green PCR Master Mix, according to the manufacturer’s instructions. Amplifications were conducted for p-glycoprotein 2, tetraspanin, carboxypeptidase, ras-related protein Rab, galectin, aspartic protease 6, and elongation factor (EF) as reference gene. Specific primer pairs were designed (Supplementary Table S1), and standard curves were generated. All reactions were performed in triplicate in the StepOnePlus Real-Time PCR Detection System (Applied Biosystems) and relative quantification was carried out with the ΔΔCT method (Livak and Schmittgen, 2001) using the abundance of EF mRNA as endogenous housekeeping control. The relative transcription levels as obtained by RT-qPCR analyses were compared with abundance levels detected by RNA-seq. Values from replicate experiments were averaged. Finally, to obtain values suitable for statistical comparisons, the fold change (FC) value was calculated. These values (plotted after conversion in log2 numbers) were used to evaluate the correlation between RNA-seq and RT-qPCR methods, applying statistical evaluation using the Pearson test (in Prism GraphPad software).

Results

Molecular identification of A. pegreffii L3

A multilocus genotyping approach, combining mitochondrial and nuclear markers, was applied to ensure species identification (Mattiucci et al., 2018, 2025). BLAST analysis of the 21 sequences obtained from the Anisakis L3 at the mtDNA cox2 gene locus (~600 bp) retrieved a percentage of identity of 99 - 100% with the sequences of A. pegreffii previously deposited (accession numbers, KY565564-KY565562). Additionally, ARMS-PCR analysis targeting the nas10 locus generated a species-specific 117 bp band, diagnostic of the C-allele, further confirming the identification of A. pegreffii (Palomba et al., 2020).

Sequencing statistics

The quality of extracted RNA from L3 reached the standard of sequencing with RIN values >7, and the concentration of each sample was higher than 50 ng/μL. The total number of raw reads obtained by sequencing was 154047207, while the total number of clean reads after filtration was 153572782 (Table 2). The total number of clean read bases was 9580666822. Q30 was greater than 97.13%, and GC content was 43.02-46.81% (Table 2).

Table 2
www.frontiersin.org

Table 2. Summary of data output quality of various libraries.

The raw data of all samples reported in this study were deposited in the NCBI (National Center for Biotechnology Information) under the accession number PRJNA752284.

Differential gene expression in A. pegreffii following iDC exposure

A total of 3914 differentially expressed genes (DEGs) with|log2FC| > 1 and p.adj< 0.05 were detected in the A. pegreffii L3 larvae treatment group (L3 co-cultured with iDCs) compared to the control group (L3 alone) (Supplementary Table S2). Results from DEG counts are reported in Table 3. Among DEGs, 1868 genes were up-regulated in L3 plus iDC compared to L3 alone, while 2046 genes were down-regulated (Figure 1). The most significantly modulated genes, characterized by the highest log2 fold change values, are highlighted in the plot (Figure 1). Notably, these included: glutathione peroxidase, SH3 domain-containing protein, calponin-like protein, and collagen cuticle N domain-containing protein, which exhibit strong differential expression, i.e., 22.00, 22.04, 28.55, 27.39 log2 fold change, respectively.

Table 3
www.frontiersin.org

Table 3. Overview of differentially expressed gene counts across log2FC thresholds, padj ≤ 0.05.

Figure 1
Volcano plot showing gene expression data with log2 fold change on the x-axis and negative log10 of the adjusted p-value on the y-axis. Points represent genes, with significant ones labeled, and dashed lines indicating thresholds for significance.

Figure 1. Volcano plot of DEGs between the treatment group (A. pegreffii L3 in the presence of iDCs) and control group (A. pegreffii in the absence of iDCs). The plot represents gene expression levels, with log2 Fold Change displayed on the x-axis and log10 (p.adjust) on the y-axis. The vertical and horizontal dashed lines indicate the fold-change cut-off = ± 1 and the p-value = 0.05, respectively. Pink dots indicate up and down-regulated genes, and light blue dots represent genes without significant changes in expression.

Functional enrichment and biological relevance of DEGs

To better understand the functional implications of the transcriptional changes, we conducted GO enrichment analysis on up- and down-regulated DEGs. GO enrichment results revealed significantly enriched categories (p.adj< 0.05) (Figure 2, Supplementary Table S3). Up-regulated genes were mostly involved in collagen and cuticulin-based cuticle development (BP), dipeptidyl-peptidase activity (MF), and anchored component of membrane (CC) (Figure 2A, Supplementary Table S3). Conversely, down-regulated genes were enriched in actin filament bundle organization (BP), nuclear receptor activity in (MF), and stereocilium bundle (CC) (Figure 2B, Supplementary Table S3). The heatmap (Figure 3) visualizes the expression patterns of the GO-annotated genes identified in the enrichment analysis. In addition, we specifically examined the list of DEGs (Supplementary Table S2) to identify genes of relevance. Further analysis identified the upregulation of genes associated with detoxification and antioxidant processes (e.g., glutathione peroxidase, thioredoxin reductase 1, peroxidase mlt-7, and superoxide dismutase), allergenic responses (e.g., heat shock protein 70, anis12, anis11, and anis10), and chaperone-related functions (e.g., calponin, SH3 domain-containing protein, HSP70, and HSP60) (Figure 4). Comparison of the DEGs identified with those previously reported in the extracellular vesicles (EVs) of A. pegreffii (Palomba et al., 2022) revealed 21 overlapping genes (Supplementary Table S4).

Figure 2
Scatter plot with two sections labeled A and B, displaying enriched gene sets. Both plots show descriptions of biological processes, cellular components, and molecular functions on the y-axis and GeneRatio on the x-axis. Dot sizes represent count, and colors indicate p.adjust values on a gradient scale. Panel A features blue shades while panel B uses orange and red shades, representing different p.adjust ranges.

Figure 2. Bubble plot of GO terms enrichment analysis of DEGs. “Count” indicates the number of DEGs enriched in the pathway; “GeneRatio” indicates the ratio of enriched DEGs to background genes; “p.adj”, shown by color, indicates a significant level of enrichment results. (A) Up-regulated genes, (B) Down-regulated genes.

Figure 3
Hierarchical clustering heatmap depicting gene expression data. Rows represent various genes, and columns are labeled A1, B2, C3, D2, E3, with subclusters. Color gradient from green to red indicates expression levels from low (-2) to high (2). Top annotation bar highlights two conditions: L3 without IDCs (in red) and L3 with IDCs (in blue).

Figure 3. Heat map showing the log2 fold change of genes enriched in GO ontology categories in A. pegreffii L3.

Figure 4
The image displays three categories of genes: Detoxificant/Antioxidant, Chaperone, and Allergenic genes. Each category lists specific genes with corresponding numerical values shown in color-coded boxes. Detoxificant/Antioxidant genes include Glutathione peroxidase (22), Peroxidase mlt-7 (13.44), Thioredoxin reductase 1 (11.77), and Superoxide dismutase [Cu-Zn] (11.55). Chaperone genes feature Calponin-like protein OV9M (28.55), SH3 domain-containing protein (22.04), Heat shock 70 kDa protein (11.70), and Chaperonin-like protein Hsp-60 (1.83). Allergenic genes display Heat shock 70 kDa (11.77), Ani s 12 allergen (6.80), Ani s 11 allergen (3.49), and Ani s 10 allergen (1.44).

Figure 4. Heat map showing the log2 fold change of upregulated genes in A. pegreffii L3 within selected specific categories.

Correlation between RNA-seq and RT-qPCR results

The abundance of transcripts emerging from transcriptomic analyses was confirmed by RT-qPCR by examining six transcripts differentially modulated in A. pegreffii L3 in the presence/absence of iDCs i.e. p-glycoprotein 2, tetraspanin, carboxypeptidase, ras-related protein Rab, galectin, and aspartic protease 6. The validation analysis was performed using triplicate RNA samples. Correlation analysis indicated significant linear relationships (R2 = 0.70; p< 0.03) between RNA-seq and RT-qPCR results (Table 4, Supplementary Figure S3).

Table 4
www.frontiersin.org

Table 4. Gene expression levels obtained by RNAseq and RT-qPCR.

Discussion

Throughout their life cycle, helminth parasites are exposed to diverse microhabitat conditions within their hosts, including fluctuations in temperature, pH, nutrient availability, and host immune pressures. Among these, the interaction with the host immune system is particularly critical for parasite survival, especially in accidental hosts, such as humans, which represent an evolutionary dead-end and an unsuitable ecological niche for A. pegreffii (Trumbić et al., 2021), placing the parasite at a disadvantageous position. In this case, the host immune system rapidly recognizes the parasite as a threat, activating local effector mechanisms that may lead to sustained inflammation (Nieuwenhuizen, 2016). The transcriptomic analysis, here described, revealed that the A. pegreffii larvae respond to the presence of iDCs primarily through an initial transcriptional pattern, likely reflecting an early attempt to cope with human derived immune-mediated stress. Notably, A. pegreffii is able to regulate a complex network of metabolic pathways, strategically investing energy in activating and upregulating genes associated with energy production (i.e., cGMP-dependent protein kinase, sodium/calcium exchangers, trehalase), as well as nitrogen and amino acid metabolism (e.g., aspartate aminotransferase, histidine ammonia-lyase, amino acid transporters). This metabolic shift likely provides the necessary resources to sustain vital functions under host-induced stress (Stryiński et al., 2024). However, this response is energetically costly and may not be sustainable in the long term. The sustained metabolic demand may deplete internal energy reserves leading to the accumulation of toxic metabolic waste products, and ultimately resulting in the death of the invading Anisakis larva.

Interestingly, this upregulation of metabolic genes appears closely linked to the enhanced expression of genes involved in cuticle maintenance. The synthesis of its components is an energy-intensive process that requires a steady supply of amino acids and efficient protein biosynthesis mechanisms (Johnstone, 1994). The modulation of structural genes (e.g., cuticle collagen dpy-7, forming homology 2, myosin heavy chain 6) is crucial for cuticle formation (Johnstone, 2000). This serves multiple functions, including acting as a protective barrier between the nematode and its environment, maintaining its body shape (Johnstone et al., 1992), and facilitating movement through connections with muscle tissue (Francis and Waterston, 1991; Hresko et al., 1999). These processes are pivotal for nematode survival, enabling interactions with hostile environments, when trying to penetrate the gastric and intestinal submucosa (Johnstone, 1994; Page et al., 2014; Sandhu et al., 2021). Over time, the expression of genes encoding these structural proteins has likely evolved to support both developmental processes and defense mechanisms (Kim et al., 2018). For instance, Kim et al. (2018) reported that collagen-related genes were abundantly transcribed in the Anisakis L3, reflecting the active molting process. Similarly, Trumbić et al. (2021) observed that the genes involved in molting, such as structural constituents of the cuticle, collagen, and cuticulin, were upregulated in A. pegreffii L3 infecting rats, suggesting that molting-related gene expression can be influenced by host-specific factors. In this context, in the presence of iDCs, A. pegreffii L3 stage larvae may perceive them as potential threats, triggering a defense response characterized by the coordinated activation of metabolic, antioxidant, and structural pathways. This upregulation may represent a protective strategy aimed at shielding the nematode from host immune effectors, potentially enhancing its survival in an accidental host. However, the direct involvement of structural gene modulation —particularly collagen and cuticulin genes— remains speculative. Further functional validation is required to confirm their role in human host–parasite interactions.

Interestingly, alongside the upregulation of structural genes involved in cuticle formation, we observed the downregulation of genes associated with procollagen-proline 4-dioxygenase activity, critical for the hydroxylation of proline residues during collagen maturation (Page and Johnstone, 2007). The simultaneous upregulation of cuticle-related genes and downregulation of procollagen-proline 4-dioxygenase activity might suggest a complex regulatory mechanism in which A. pegreffii actively synthesizes structural precursors but limits their full post-translational maturation. However, this interpretation remains speculative, and further molecular and functional analyses are required to confirm the regulatory dynamics involved.

Moreover, the downregulation of genes involved in cytoskeletal organization (e.g., actin filament, bundle assembly) and membrane dynamics (e.g., membrane-to-membrane docking) suggests a strategic suppression of processes related to cellular motility and intracellular trafficking. This downregulation may reflect an attempt to conserve energy by limiting metabolically expensive activities that are not immediately critical for the parasite development.

Alongside this metabolic shift, genes involved in oxidative stress responses and those linked to ion homeostasis (sodium/calcium exchangers, anion exchange proteins), may play a role in modulating the host immune response (Nawaratna et al., 2018). This mechanism could contribute to mitigating the effects of reactive oxygen species (ROS) generated during the immune response, potentially enhancing parasite survival within the host. This suggests that A. pegreffii not only adjusts its metabolic activity, but also engages calcium-dependent pathways and antioxidant-related mechanisms to counteract host-induced oxidative stress. The modulation of antioxidant enzymes is a key strategy for parasites to defend themselves against the host immune system, specifically targeting ROS generated by immune cells such as macrophages, neutrophils and eosinophils. Indeed, our previous studies have shown that iDCs strongly increase ROS production as part of their immune defense against A. pegreffii (Napoletano et al., 2018). In addition, the increased transcription of genes implicated in the detoxification of ROS and bolstering antioxidant defense, the upregulation of glutathione peroxidase, thioredoxin reductase 1, peroxidase mlt-7 and superoxide dismutase suggest a metabolic pathway of A. pegreffii to shield itself from ROS produced by host tissue-resident cells, besides iDCs. For example, glutathione peroxidase plays a critical role in removing harmful ROS by catalyzing the reduction of hydrogen peroxide and organic hydroperoxides in A. pegreffii. Instead, superoxidase dismutase has been shown to be exploited in resistance to radiation (Seo et al., 2006) and xenobiotics (Stryiński et al., 2022, 2024), suggesting a versatile role of antioxidant enzymes in the parasite’s defense strategy against various environmental stressors. In the case of thioredoxin reductase, recent studies in filariasis have shown that inhibiting this enzyme induces apoptotic death of the parasite (Sen et al., 2023), highlighting its crucial role in the parasite’s survival and, consequently, in the disease pathogenesis and progression. Additionally, it is interesting to note that the thioredoxin reductase exerts an immunosuppressive role by downregulating the “inflammasome” pathways in macrophages exposed to the parasite (Joardar et al., 2021).

In addition to antioxidant enzymes, the modulation of molecular chaperones represents another critical mechanism employed by A. pegreffii to manage iDC-induced stress. Among these, HSP70 plays a key role in stabilizing proteins during stress (Hartl, 1996). Beyond its role in protein maintenance, HSP70 is involved in differentiation, protection against host-induced damage (including oxidative stress from free radicals), and promoting virulence (Polla, 1991; Cruz-Laufer et al., 2025). These functions allow parasites to adapt to various host niches, endure immune responses, and respond to environmental changes. In this context, HSP70 also emerges as a potential virulence factor of A. pegreffii. The up-regulation of the Hsp70 gene in A. pegreffii co-cultured with iDCs can be understood as an adaptation to changing environments and host-parasite interactions over time. Additionally, it is noteworthy that HSP70 was detected as being carried by EVs released by L3 larvae, along with other proteins involved in metabolism and immune modulation (Palomba et al., 2023). Other specific genes were also found in EVs of A. pegreffii to be differentially regulated following exposure to iDCs (Supplementary Table S4). Beyond their role in stress responses and parasite adaptation, HSP70s are also recognized allergens in various arthropods, including mites, midges, flies, and cockroaches. Similarly, HSPs (including HSP β-1, HSP - 12.2, HSP70, HSP 75KDa, and HSP 90-α), along with superoxide dismutase and galectin, have been identified in the larval extract of the ascarid Toxocara canis (de Silva et al., 2018). In the larva migrans of T. canis, HSPs serve as key immunomodulatory molecules, playing essential roles in the parasite survival, even during accidental human infection (Pérez-Morales and Espinoza, 2015).

Finally, Anis10, Anis11, and Anis12—previously described as allergens—were detected among the differentially expressed genes (Audicana and Kennedy, 2008; Mattiucci et al., 2017). However, their specific roles remain unclear, although they are likely closely interconnected with the host’s immune response, and future investigations are essential to elucidate the precise functions of these genes and their coded proteins (Castrignanò et al., 2007).

Conclusions

This study provides the first comprehensive transcriptomic overview of A. pegreffii L3 larvae during in vitro exposure to human iDCs, revealing for the first time a coordinated activation of antioxidant enzymes, molecular chaperones, and structural genes in response to early immune recognition. The responses triggered by the parasite likely reflect a short-term survival strategy against immune-mediated stress in a non-permissive host environment. Notably, the upregulation of detoxification enzymes and molecular chaperones, along with the modulation of cuticle-related genes, highlights the parasite’s attempt to mitigate host-induced oxidative stress and maintain structural integrity. However, the concurrent downregulation of genes necessary for collagen maturation and cellular motility suggests that while A. pegreffii initiates defense mechanisms, it may be unable to sustain the energetic cost of adaptation in the human host.

Overall, the molecular characterization of the A. pegreffii L3 stage larvae, interacting with iDCs here described, appear to be consistent with the parasite’s inability to develop. At the same time, our data show that the larva activates a transcriptional response likely aimed at counteracting the immune pressure exerted by iDCs — an early, but ultimately unsuccessful, attempt to ensure survival.

While previous work by Napoletano et al. (2018) has shown that A. pegreffii L3 can impair DC function by inducing apoptosis and blocking their differentiation, maturation, and migration to lymph nodes, our study focuses on the parasite’s point of view — specifically, how the larva itself responds to contact with human iDCs at the transcriptomic level. The expression patterns observed suggest that the parasite attempts to resist immune stress by activating genes involved in energy metabolism, oxidative stress response, and structural maintenance. This may allow temporary adaptation to the host environment and could contribute indirectly to shaping the local immune response. In particular, the modulation of DC activity by the parasite — as reported in our previous studies (Napoletano et al., 2018) — may influence Th1/Th2 polarization, potentially favoring a shift toward a Th2-biased environment, typical of chronic helminth infections (Sher et al., 2003).

Further studies are needed to explore whether the upregulated genes identified in this analysis influence the host immune response. For instance, specific antigens derived from these genes could be tested in vitro to assess their immunomodulatory potential, as already done in other parasite models, such as Trichinella spiralis (Yu et al., 2025). In this context, human organoid systems represent a promising experimental platform to investigate the host–parasite interaction in a physiologically relevant setting. These 3D models can mimic the structure and function of the human gastrointestinal epithelium and allow co-culture with immune cells (Perez et al., 2025). They could be used to assess barrier integrity, local cytokine responses, and immune cell activation, as well as to explore whether A. pegreffii antigens contribute to chronic inflammation or tumor-promoting conditions, particularly in the context of prolonged or repeated exposure (Bellini et al., 2024).

Data availability statement

The datasets generated for this study can be found in the NCBI (National Center for Biotechnology Information) under the accession number PRJNA752284 and in Figshare (see Table 3).

Ethics statement

The studies involving humans were approved by Ethics Committee of University Hospital “Policlinico Umberto I” ‘Sapienza’ University of Rome (Protocol nr. 4212). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MP: Writing - review & editing, Conceptualization, Methodology, Writing - original draft, Investigation, Resources. AR: Writing - original draft, Resources, Methodology, Writing - review & editing, Investigation, Conceptualization. TC: Investigation, Writing - original draft, Software, Validation, Supervision, Methodology, Writing - review & editing. CN: Writing - review & editing, Writing - original draft, Methodology. XR-G: Writing - review & editing, Writing - original draft, Methodology. VP: Validation, Writing - original draft, Software, Writing - review & editing. FL: Writing - original draft, Software, Writing - review & editing, Validation. DC: Writing - review & editing, Writing - original draft. SM: Writing - original draft, Conceptualization, Funding acquisition, Resources, Project administration, Methodology, Writing - review & editing, Supervision, Investigation.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This study was supported by the Italian Ministry of Health (RF) 2018 - 12367986, title “Innovative approaches and parameters in the diagnosis and epidemiological surveillance of the Anisakisrelated human diseases in Italy”.

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.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fcimb.2025.1646537/full#supplementary-material

References

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

PubMed Abstract | Crossref Full Text | Google Scholar

Bellini, I., Scribano, D., Ambrosi, C., Chiovoloni, C., Rodon, S., Pronio, A., et al. (2025). Anisakis extracellular vesicles elicit immunomodulatory and potentially tumorigenic outcomes on human intestinal organoids. Parasitol. Vectors. 17, 393. doi: 10.1186/s13071-024-06471-7

PubMed Abstract | Crossref Full Text | Google Scholar

Brusca, I., Barrale, M., Zarcone, M., Fruscione, S., Onida, R., De Bella, D. D., et al. (2023). Basophil activation test in the diagnosis of Anisakis allergy: An observational study from an area of high seafood consumption in Italy. Pathogens 12, 777. doi: 10.3390/pathogens12060777

PubMed Abstract | Crossref Full Text | Google Scholar

Buchfink, B., Xie, C., and Huson, D. H. (2014). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | Crossref Full Text | Google Scholar

Cantalapiedra, C. P., Hernandez-Plaza, A., Letunic, I., Bork, P., and Huerta-Cepas, J. (2021). eggNOG-mapper v2: Functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol. Biol. Evol. 38, 5825–5829. doi: 10.1093/molbev/msab293

PubMed Abstract | Crossref Full Text | Google Scholar

Castrignanò, T., De Meo, P. D. O., Carrabino, D., Orsini, M., Floris, M., and Tramontano, A. (2007). The MEPS server for identifying protein conformational epitopes. BMC Bioinform. 8, S6. doi: 10.1186/1471-2105-8-S1-S6

PubMed Abstract | Crossref Full Text | Google Scholar

Castrignanò, T., Gioiosa, S., Flati, T., Cestari, M., Picardi, E., Chiara, M., et al. (2020). ELIXIR-IT HPC @ CINECA: High performance computing resources for the bioinformatics community. BMC Bioinform. 21, 352. doi: 10.1186/s12859-020-03565-8

PubMed Abstract | Crossref Full Text | Google Scholar

Cavallero, S., Bellini, I., Pizzarelli, A., Arca, B., and D’Amelio, S. (2022). A miRNAs catalogue from third-stage larvae and extracellular vesicles of Anisakis pegreffii provides new clues for host-parasite interplay. Sci. Rep. 12, 9667. doi: 10.1038/s41598-022-13594-3

PubMed Abstract | Crossref Full Text | Google Scholar

Cipriani, P., Sbaraglia, G. L., Paoletti, M., Giulietti, L., Bellisario, B., Palomba, M., et al. (2018). The Mediterranean European hake, Merluccius merluccius: Detecting drivers influencing the Anisakis spp. larvae distribution. Fish. Res. 202, 79–89. doi: 10.1016/j.fishres.2017.07.010

Crossref Full Text | Google Scholar

Cruz-Laufer, A. J., Vanhove, M. P. M., Bachmann, L., Barson, M., Bassirou, H., Bitja Nyom, A. R., et al. (2025). Adaptive evolution of stress response genes in parasites aligns with host niche diversity. BMC Biol. 23, 10. doi: 10.1186/s12915-024-02091-w

PubMed Abstract | Crossref Full Text | Google Scholar

da Silva, M. B., Oviedo, Y., Cooper, P. J., Pacheco, L. G., Pinheiro, C. S., Ferreira, F., et al. (2018). The somatic proteins of Toxocara canis larvae and excretory-secretory products revealed by proteomics. Vet. Parasitol. 259, 25–34. doi: 10.1016/j.vetpar.2018.06.015

PubMed Abstract | Crossref Full Text | Google Scholar

Dionisi, M., De Archangelis, C., Battisti, F., Rahimi Koshkaki, H., Belleudi, F., Zizzari, I. G., et al. (2018). Tumor-derived microvesicles enhance cross-processing ability of clinical grade dendritic cells. Front. Immunol. 9. doi: 10.3389/fimmu.2018.02481

PubMed Abstract | Crossref Full Text | Google Scholar

EFSA Panel on Biological Hazards, Koutsoumanis, K., Allende, A., Alvarez-ordóñez, A., Bover-cid, S., Chemaly, M., et al. (2024). Re-evaluation of certain aspects of the EFSA Scientific Opinion of April 2010 on risk assessment of parasites in fishery products, based on new scientific data. EFSA J22, e9090. doi: 10.2903/j.efsa.2024.9090

PubMed Abstract | Crossref Full Text | Google Scholar

Flati, T., Gioiosa, S., Spallanzani, N., Tagliaferri, I., Diroma, M. A., Pesole, G., et al. (2020). HPC-REDItools: A novel HPC-aware tool for improved large-scale RNA-editing analysis. BMC Bioinform. 21, 353. doi: 10.1186/s12859-020-03562-x

PubMed Abstract | Crossref Full Text | Google Scholar

Francis, R. and Waterston, R. H. (1991). Muscle cell attachment in Caenorhabditis elegans. J. Cell. Biol. 14, 465–479. doi: 10.1083/jcb.114.3.465

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Hartl, F. U. (1996). Molecular chaperones in cellular protein folding. Nature 381, 571–579. doi: 10.1038/381571a0

PubMed Abstract | Crossref Full Text | Google Scholar

Hresko, M. C., Schriefer, L. A., Shrimankar, P., and Waterston, R. H. (1999). Myotactin, a novel hypodermal protein involved in muscle-cell adhesion in Caenorhabditis elegans. J. Cell. Biol. 146, 659–672. doi: 10.1083/jcb.146.3.659

PubMed Abstract | Crossref Full Text | Google Scholar

Joardar, N., Bhattacharya, R., Halder, S., Sen, A., Biswas, S. R., Jana, K., et al. (2021). Filarial thioredoxin reductase exerts anti-inflammatory effects upon lipopolysaccharide induced inflammation in macrophages. Int. J. Biol. Macromol. 193, 1379–1390. doi: 10.1016/j.ijbiomac.2021.10.200

PubMed Abstract | Crossref Full Text | Google Scholar

Johnstone, I. L. (1994). The cuticle of the nematode Caenorhabditis elegans: a complex collagen structure. Bioessays 16, 171–178. doi: 10.1002/bies.950160307

PubMed Abstract | Crossref Full Text | Google Scholar

Johnstone, I. L. (2000). Cuticle collagen genes expression in Caenorhabditis elegans. Trends Genet. 16, 21–27. doi: 10.1016/s0168-9525(99)01857-0

PubMed Abstract | Crossref Full Text | Google Scholar

Johnstone, I. L., Shafi, Y., and Barry, J. D. (1992). Molecular analysis of mutations in the Caenorhabditis elegans collagen gene dpy-7. EMBO J. 1, 3857–3863. doi: 10.1002/j.1460-2075.1992.tb05478.x

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Livak, K. J. and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2– ΔΔCT method. methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | Crossref Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

Lundie, R. J., Webb, L. M., Marley, A. K., Phythian-Adams, A. T., Cook, P. C., Jackson-Jones, L. H., et al. (2016). A central role for hepatic conventional dendritic cells in supporting Th2 responses during helminth infection. Immunol. Cell. Biol. 94, 400–410. doi: 10.1038/icb.2015.114

PubMed Abstract | Crossref Full Text | Google Scholar

Maizels, R. M. and McSorley, H. J. (2016). Regulation of the host immune system by helminth parasites. J. Allergy Clin. Immunol. 138, 666–675. doi: 10.1016/j.jaci.2016.07.007

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Mattiucci, S., Colantoni, A., Crisafi, B., Mori-Ubaldini, F., Caponi, L., Fazii, P., et al. (2017). IgE sensitization to Anisakis pegreffii in Italy: Comparison of two methods for the diagnosis of allergic anisakiasis. Parasite Immunol. 39, 1–10. doi: 10.1111/pim.12440

PubMed Abstract | Crossref Full Text | Google Scholar

Mattiucci, S., Fazii, P., De Rosa, A., Paoletti, M., Megna, A. S., Glielmo, A., et al. (2013). Anisakiasis and gastroallergic reactions associated with Anisakis pegreffii infection, Italy. Emerg. Infect. Dis. 19, 496–499. doi: 10.3201/eid1903.121017

PubMed Abstract | Crossref Full Text | Google Scholar

Mattiucci, S., Palomba, M., Belli, B., Aco-Alburqueque, R., Cipriani, P., Roca-Gerones, X., et al. (2025). Hybridization and introgression of mitochondrial genome between the two species Anisakis pegreffii and A. simplex (s. s.) using a wide genotyping approach: evolutionary and ecological implications. Parasitology. 152, 293–313 doi: 10.1017/S0031182025000228

PubMed Abstract | Crossref Full Text | Google Scholar

Mattiucci, S., Palomba, M., and Nascetti, G. (2022). “Anisakis.” Encyclopedia of infection and immunity. 2, 408–423. doi: 10.1016/B978-0-12-818731-9.00075-6

Crossref Full Text | Google Scholar

Mehrdana, F., Kania, P. W., Nazemi, S., and Buchmann, K. (2017). Immunomodulatory effects of excretory/secretory compounds from Contracaecum osculatum larvae in a zebrafish inflammation model. PloS One 12, e0181277. doi: 10.1371/journal.pone.0181277

PubMed Abstract | Crossref Full Text | Google Scholar

Napoletano, C., Mattiucci, S., Colantoni, A., Battisti, F., Zizzari, I. G., Rahimi, H., et al. (2018). Anisakis pegreffii impacts differentiation and function of human dendritic cells. Parasite Immunol. 40, e12527. doi: 10.1111/pim.12527

PubMed Abstract | Crossref Full Text | Google Scholar

Napoletano, C., Pinto, D., Bellati, F., Taurino, F., Rahimi, H., Tomao, F., et al. (2007). A comparative analysis of serum and serum-free media for generation of clinical grade DCs. J. Immunother. 30, 567–576. doi: 10.1097/CJI.0b013e318046f396

PubMed Abstract | Crossref Full Text | Google Scholar

Nawaratna, S. S. K., You, H., Jones, M. K., McManus, D. P., and Gobert, G. N. (2018). Calcium and Ca 2+/Calmodulin-dependent kinase II as targets for helminth parasite control. Biochem. Soc Trans. 46, 1743-1751. doi: 10.1042/BST20180480

PubMed Abstract | Crossref Full Text | Google Scholar

Nieuwenhuizen, N. E. (2016). Anisakis–immunology of a foodborne parasitosis. Parasite Immunol. 38, 548–557. doi: 10.1111/pim.12349

PubMed Abstract | Crossref Full Text | Google Scholar

Nieuwenhuizen, N. E. and Lopata, A. L. (2013). Anisakis-a food-borne parasite that triggers allergic host defences. Int. J. Parasitol. 43, 1047–1057. doi: 10.1016/j.ijpara.2013.08.001

PubMed Abstract | Crossref Full Text | Google Scholar

Page, A. P. and Johnstone, I. L. (2007). The cuticle. WormBook, 1–15. doi: 10.1895/wormbook.1.138.1

PubMed Abstract | Crossref Full Text | Google Scholar

Page, A. P., Stepek, G., Winter, A. D., and Pertab, D. (2014). Enzymology of the nematode cuticle: a potential drug target? Int. J. Parasitol. Drugs Drug Resist. 4, 133–141. doi: 10.1016/j.ijpddr.2014.05.003

PubMed Abstract | Crossref Full Text | Google Scholar

Palomba, M., Libro, P., Di Martino, J., Rughetti, A., Santoro, M., Mattiucci, S., et al. (2022). De novo transcriptome assembly and annotation of the third stage larvae of the zoonotic parasite Anisakis pegreffii. BMC Res. Notes 15, 223. doi: 10.1186/s13104-022-06099-9

PubMed Abstract | Crossref Full Text | Google Scholar

Palomba, M., Paoletti, M., Colantoni, A., Rughetti, A., Nascetti, G., and Mattiucci, S. (2019). Gene expression profiles of antigenic proteins of third stage larvae of the zoonotic nematode Anisakis pegreffii in response to temperature conditions. Parasite 26, 52. doi: 10.1051/parasite/2019055

PubMed Abstract | Crossref Full Text | Google Scholar

Palomba, M., Paoletti, M., Colantoni, A., Rughetti, A., Nascetti, G., and Mattiucci, S. (2023). Proteomic characterization of extracellular vesicles released by third stage larvae of the zoonotic parasite Anisakis pegreffii (Nematoda: Anisakidae). Front. Cell. Infect. Microbiol. 13. doi: 10.3389/fcimb.2023.1079991

PubMed Abstract | Crossref Full Text | Google Scholar

Palomba, M., Paoletti, M., Webb, S. C., Nascetti, G., and Mattiucci, S. (2020). A novel nuclear marker and development of an ARMS-PCR assay targeting the metallopeptidase 10 (nas 10) locus to identify the species of the Anisakis simplex (s. l.) complex (Nematoda, Anisakidae). Parasite 27, 39. doi: 10.1051/parasite/2020033

PubMed Abstract | Crossref Full Text | Google Scholar

Patro, R., Duggal, G., Love, M. I., Irizarry, R. A., and Kingsford, C. (2017). Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 14, 417–419. doi: 10.1038/nmeth.4197

PubMed Abstract | Crossref Full Text | Google Scholar

Peng, J., Federman, H. G., Hernandez, C. M., and Siracusa, M. C. (2022). Communication is key: Innate immune cells regulate host protection to helminths. Front. Immunol. 13. doi: 10.3389/fimmu.2022.995432

PubMed Abstract | Crossref Full Text | Google Scholar

Perez, M. G., Coman, D., Neves, J. F., and Britton, C. (2025). Organoids in parasitology: a game-changer for studying host-nematode interactions. Parasitology. 1–44 doi: 10.1017/S0031182025100620

PubMed Abstract | Crossref Full Text | Google Scholar

Pérez-Morales, D. and Espinoza, B. (2015). The role of small heat shock proteins in parasites. Cell. Stress Chaperones 20, 767–780. doi: 10.1007/s12192-015-0607-y

PubMed Abstract | Crossref Full Text | Google Scholar

Pinna, V., Di Martino, J., Liberati, F., Bottoni, P., and Castrignanò, T. (2024). “IGUANER - differential gene expression and fUnctionAlaNalyzER,” in Lecture notes in computer science, vol. 14516 . Eds. Sachdeva, S. and Watanobe, Y.. (Switzerland: Springer Nature), 78–93. doi: 10.1007/978-3-031-58502-9_5

Crossref Full Text | Google Scholar

Polla, B. S. (1991). Immunol. Today 12.3, A38-A41. doi: 10.1016/S0167-5699(05)80011-8

PubMed Abstract | Crossref Full Text | Google Scholar

Sandhu, A., Badal, D., Sheokand, R., Tyagi, S., and Singh, V. (2021). Specific collagens maintain the cuticle permeability barrier in Caenorhabditis elegans. Genetics 217, iyaa047. doi: 10.1093/genetics/iya047

PubMed Abstract | Crossref Full Text | Google Scholar

Sen, A., Dewaker, V., Debnath, U., Jana, K., Rath, J., Joardar, N., et al. (2023). In silico exploration and in vitro validation of the filarial thioredoxin reductase inhibitory activity of Scytonemin and its derivatives. J. Biomol. Struct. Dyn 21, 1–13. doi: 10.1080/07391102.2023.2283876

PubMed Abstract | Crossref Full Text | Google Scholar

Seo, M., Kho, B.-M., Guk, S.-M., Lee, S.-H., and Chai, J.-Y. (2006). Radioresistance of Anisakis simplex third-stage larvae and the possible role of superoxide dismutase. J. Parasitol. 92, 416–418. doi: 10.1645/GE-610R.1

PubMed Abstract | Crossref Full Text | Google Scholar

Sher, A., Pearce, E., and Kaye, P. (2003). Shaping the immune response to parasites: role of dendritic cells. Curr. Opin. Immunol. 15, 421–429.

PubMed Abstract | Google Scholar

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 ss: Enzymes of energy and/or carbohydrate metabolism as potential drug targets in anisakiasis. Int. J. Mol. Sci. 23, 4336. doi: 10.3390/ijms23084336

PubMed Abstract | Crossref Full Text | Google Scholar

Stryiński, R., Polak, I., Gawryluk, A., Rosa, P., and Łopieńska-Biernat, E. (2024). The response of Anisakis simplex (s. s.) to anthelmintics - Specific changes in xenobiotic metabolic processes. Exp. Parasitol., 261, 108751. doi: 10.1016/j.exppara.2024.108751

PubMed Abstract | Crossref Full Text | Google Scholar

Trumbić, Ž., Hrabar, J., Palevich, N., Carbone, V., and Mladineo, I. (2021). Molecular and evolutionary basis for survival, its failure, and virulence factors of the zoonotic nematode Anisakis pegreffii. Genomics 113, 2891–2905. doi: 10.1016/j.ygeno.2021.06.032

PubMed Abstract | Crossref Full Text | Google Scholar

Valentini, A., Mattiucci, S., Bondanelli, P., Webb, S. C., Mignucci-Giannone, A. A., Colom-Llavina, M. M., et al. (2006). Genetic relationships among Anisakis species (Nematoda: Anisakidae) inferred from mitochondrial cox2 sequences, and comparison with allozyme data. J. Parasitol. 92, 156–166. doi: 10.1645/GE-3504.1(2006

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, W., Yuan, X., Zhai, P., Li, X., and Han, C. (2025). Effects of Trichinella spiralis excretory-secretory antigens on expression of indoleamine 2, 3-dioxygenase on dendritic cells in vitro. Parasite 32, 26. doi: 10.1051/parasite/2025018

PubMed Abstract | Crossref Full Text | Google Scholar

Zamora, V., Andreu-Ballester, J. C., Rodero, M., and Cuéllar, C. (2021). Anisakis simplex: Immunomodulatory effects of larval antigens on the activation of Toll like Receptors. International. Immunopharmacol. 100, 108120. doi: 10.1016/j.intimp.2021.108120

PubMed Abstract | Crossref Full Text | Google Scholar

Zamora, V., Rodero, M., Andreu-Ballester, J. C., Mendez, S., and Cuéllar, C. (2019). Induction of tolerogenic properties by Anisakis larval antigens on murine dendritic cells. Parasite Immunol. 41, e12616. doi: 10.1111/pim.12616

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: human anisakiasis, Anisakis pegreffii, human dendritic cells, RNA-seq, gene expression, oxidative stress, energy metabolism

Citation: Palomba M, Rughetti A, Castrignanò T, Napoletano C, Roca-Geronès X, Pinna V, Liberati F, Canestrelli D and Mattiucci S (2025) Transcriptomic profile of the zoonotic parasite Anisakis pegreffii upon in vitro exposure to human dendritic cells. Front. Cell. Infect. Microbiol. 15:1646537. doi: 10.3389/fcimb.2025.1646537

Received: 13 June 2025; Accepted: 26 August 2025;
Published: 15 September 2025.

Edited by:

Collette Britton, University of Glasgow, United Kingdom

Reviewed by:

Daniel A. Abugri, Alabama State University, United States
Saurabh Gupta, GLA University, India

Copyright © 2025 Palomba, Rughetti, Castrignanò, Napoletano, Roca-Geronès, Pinna, Liberati, Canestrelli and Mattiucci. 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: Simonetta Mattiucci, c2ltb25ldHRhLm1hdHRpdWNjaUB1bmlyb21hMS5pdA==; Tiziana Castrignanò, dGl6aWFuYS5jYXN0cmlnbmFub0B1bml0dXMuaXQ=

These authors have contributed equally to this work

Disclaimer: 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.