
==== Front
mSystems
mSystems
msystems
mSystems
2379-5077
American Society for Microbiology 1752 N St., N.W., Washington, DC

39191377
msystems00545-24
10.1128/msystems.00545-24
msystems.00545-24
Research Article
editors-pickEditor’s Picksystems-biologySystems BiologyGut microbiota metabolically mediate intestinal helminth infection in zebrafish
https://orcid.org/0009-0003-4753-6347
Hammer Austin J. 1 Formal analysis Investigation Methodology Software Visualization Writing – original draft Writing – review and editing
https://orcid.org/0000-0001-7741-7822
Gaulke Christopher A. 2 Conceptualization Data curation Formal analysis Investigation Methodology Project administration Writing – original draft Writing – review and editing
https://orcid.org/0000-0003-2987-2913
Garcia-Jaramillo Manuel 3 Conceptualization Data curation Formal analysis Investigation Methodology Validation Writing – review and editing
https://orcid.org/0009-0003-4753-6347
Leong Connor 1 4 Formal analysis Investigation Methodology Visualization Writing – original draft Writing – review and editing
Morre Jeffrey 5 Data curation Formal analysis Investigation
https://orcid.org/0000-0002-8332-3408
Sieler Jr. Michael J. 1 Visualization Writing – review and editing
https://orcid.org/0000-0002-6348-9347
Stevens Jan F. 6 7 Writing – review and editing
Jiang Yuan 8 Funding acquisition Investigation Methodology Writing – review and editing
https://orcid.org/0000-0003-0743-8956
Maier Claudia S. 5 Conceptualization Formal analysis Funding acquisition Investigation Project administration Supervision Writing – original draft Writing – review and editing
https://orcid.org/0000-0001-7968-5142
Kent Michael L. 1 Conceptualization Formal analysis Funding acquisition Investigation Methodology Project administration Resources Supervision Writing – original draft Writing – review and editing
https://orcid.org/0000-0002-5996-7764
Sharpton Thomas J. 1 8 Conceptualization Data curation Formal analysis Funding acquisition Investigation Methodology Project administration Resources Supervision Writing – original draft Writing – review and editing Thomas.Sharpton@oregonstate.edu

1 Department of Microbiology, Oregon State University , Oregon, USA
2 Department of Pathobiology, University of Illinois Urbana Champaign , Illinois, USA
3 Department of Environmental and Molecular Toxicology, Oregon State University , Oregon, USA
4 Department of Biomedical Sciences, Oregon State University , Corvallis, Oregon, USA
5 Department of Chemistry, Oregon State University , Oregon, USA
6 Department of Pharmaceutical Sciences, Oregon State University , Oregon, USA
7 Linus Pauling Institute, Oregon State University , Oregon, USA
8 Department of Statistics, Oregon State University , Oregon, USA
Editor Rawls John F. Duke University School of Medicine , Durham, North Carolina, USA

Address correspondence to Thomas J. Sharpton, Thomas.Sharpton@oregonstate.edu
Austin J. Hammer and Christopher A. Gaulke contributed equally to this article. Author order was determined in order of reverse seniority.

The authors declare no conflict of interest.

9 2024
27 8 2024
27 8 2024
9 9 e00545-2419 4 2024
29 7 2024
Copyright © 2024 Hammer et al.
2024
Hammer et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

Intestinal helminth parasite (IHP) infection induces alterations in the composition of microbial communities across vertebrates, although how gut microbiota may facilitate or hinder parasite infection remains poorly defined. In this work, we utilized a zebrafish model to investigate the relationship between gut microbiota, gut metabolites, and IHP infection. We found that extreme disparity in zebrafish parasite infection burden is linked to the composition of the gut microbiome and that changes in the gut microbiome are associated with variation in a class of endogenously produced signaling compounds, N-acylethanolamines, that are known to be involved in parasite infection. Using a statistical mediation analysis, we uncovered a set of gut microbes whose relative abundance explains the association between gut metabolites and infection outcomes. Experimental investigation of one of the compounds in this analysis reveals salicylaldehyde, which is putatively produced by the gut microbe Pelomonas, as a potent anthelmintic with activity against Pseudocapillaria tomentosa egg hatching, both in vitro and in vivo. Collectively, our findings underscore the importance of the gut microbiome as a mediating agent in parasitic infection and highlight specific gut metabolites as tools for the advancement of novel therapeutic interventions against IHP infection.

IMPORTANCE

Intestinal helminth parasites (IHPs) impact human health globally and interfere with animal health and agricultural productivity. While anthelmintics are critical to controlling parasite infections, their efficacy is increasingly compromised by drug resistance. Recent investigations suggest the gut microbiome might mediate helminth infection dynamics. So, identifying how gut microbes interact with parasites could yield new therapeutic targets for infection prevention and management. We conducted a study using a zebrafish model of parasitic infection to identify routes by which gut microbes might impact helminth infection outcomes. Our research linked the gut microbiome to both parasite infection and to metabolites in the gut to understand how microbes could alter parasite infection. We identified a metabolite in the gut, salicylaldehyde, that is putatively produced by a gut microbe and that inhibits parasitic egg growth. Our results also point to a class of compounds, N-acyl-ethanolamines, which are affected by changes in the gut microbiome and are linked to parasite infection. Collectively, our results indicate the gut microbiome may be a source of novel anthelmintics that can be harnessed to control IHPs.

KEYWORDS

gut microbiome
host-parasite relationship
helminths
multi-omic integration
mediation analysis
National Science Foundation (NSF) 2025457 Sharpton Thomas J. HHS | National Institutes of Health (NIH) R21AI135641 Sharpton Thomas J. HHS | National Institutes of Health (NIH) R01ES030226 Sharpton Thomas J. HHS | National Institutes of Health (NIH) S10RR027878 Stevens Jan F. cover-dateSeptember 2024
==== Body
pmcINTRODUCTION

Intestinal helminth parasitic infections present a significant global health burden, affecting at least one-quarter of the global population (1, 2), and are disproportionately experienced by individuals in impoverished nations, particularly children (3). These infections are also prominent among domestic animals (4, 5), which places tremendous strain on livestock management and veterinary practices (4, 6, 7). Among infected individuals intestinal helminth parasite (IHP) infections can contribute to anemia (8), cognitive impairment (9), physical wasting (10), as well as a host of other conditions that contribute to the equivalent of millions of disability-adjusted life years (11). Unfortunately, the extreme burden of IHP infection may be exacerbated by the emergence of drug-resistant parasites. High levels of broad anthelmintic drug resistance in animal populations have been observed globally for decades (12, 13) posing a tremendous risk to humans because practices such as broad blanket anthelmintic administration (14) and prior widespread prophylactic administration of anthelmintic drugs (15) have provided strong selective forces for virulent drug-resistant organisms (16). While improved helminth management practice may help slow the rate of anthelmintic resistance (16, 17), the future of controlling IHP infection may depend on innovating new methods and resources in anthelmintic discovery to stay ahead in the pugilistic battle with drug-resistance in infectious nematodes.

In the search for new approaches to control IHPs, there is suggestive evidence that the intestinal microbiome can enhance or reduce parasite infection (18). Microorganisms produce a diverse trove of metabolic compounds, including anthelmintic drugs. For example, the avermectin class of compounds, which include the most widely administered anthelmintic drugs on the planet, was originally derived from soil-borne bacteria such as Streptomyces avermectilis (19). Besides the soil, locations where helminths and microbes have evolved to co-locate, such as the gut, may offer a rich resource of microbially derived anthelmintic compounds (18) as their evolution by microbial community members may have been critical to microbial exploitation of the shared ecological niche. However, the complex and variegated metabolic landscape of microbially derived compounds that are found in the gut and relevant to IHP infection remains unexplored.

Little is known about the existence of anthelmintic compounds derived from gut bacteria, but it is known that gut microbiota can drive intestinal helminth infection through a variety of mechanisms. For instance, bacteria can alter the protective integrity of mucosal barriers (20–23), drive peristaltic activity (24, 25), and produce inhibitory antibiotic compounds that limit pathogen and parasite survival (19, 26–28). Moreover, signals from specific gut bacteria can cause egg hatching of helminthic parasites in the gut, such as in the case of Trichuris muris infections in mice (29). In addition, gut microbes engage in extensive interaction with the vertebrate host immune system (30), and intestinal helminths possess a diverse suite of immunomodulatory tools (31, 32). Much of this cross-talk depends on the production of metabolite products by host, microbe, and parasite. However, our current understanding of the complex set of metabolite interactions that may directly or indirectly drive parasite colonization success is limited. Clarifying the set of metabolites that mediate the interaction between the host, parasites, and the gut microbiome may provide a toolkit of compounds to control parasite infection.

Efforts to understand the role of the gut microbiome in health and disease conditions can benefit from the application of analytic techniques that examine possible mediating roles of intestinal microbes and metabolites. To this end, mediation inference techniques seek to examine whether the relationship between two variables depends on the hypothesized mediating effect of some third variable. While early iterations of these methods relied on regression-based structural equation modeling (33) with strict assumptions regarding the model type, new methods are being innovated that account for data-specific assumptions, such as sparsity and compositionality in microbiome data (34–37), and the suite of mediation tools available to researchers is expanding rapidly. Recent applications of mediation analysis in microbiome science have identified lipid compounds produced by Akkermansia muciniphila that modulate murine immunity and metabolism (38), clarified the role of gut microbiota in mediating the relationship between diet and immune inflammation (39), and established the role of the gut microbiome in the development of childhood asthma (40). However, the extension of mediation techniques to high-dimensional multi-omic data where there are a large number of both independent and mediating features remains limited. Large sample sizes are often required to attain sufficient power to detect mediating effects (41), which can impose substantial logistical challenges on researchers seeking to use these techniques with expensive vertebrate models. Despite this challenge, mediation techniques have played an important recent role in mapping the gut microbiome to metabolites that may be involved in disease conditions such as anorexia nervosa (42), and these methods offer a dynamic opportunity to ascertain the intricate mechanisms by which the gut microbiome links to complex diseases, such as parasite infection.

Robustly identifying novel connections between the gut microbiome and parasite infection via mediating metabolites presents a challenge that requires experimental investigation using an organism that displays robust patterns of quantifiable IHP infection, displays a tractable set of gut microbes, and may be scaled in the laboratory to deal with inherent statistical limitations of mediation methods. In line with these requirements, the zebrafish model provides a powerful tool for modeling parasite infections and shedding light on the intricate relationship between host-microbiota interactions (43) and disease outcomes (44, 45). The model possesses a well-characterized taxonomic gut microbial composition (46, 47) with a functional composition that resembles that of humans (48), zebrafish offers a well-established model of intestinal helminth parasite infection (49), and can be experimentally scaled in a cost-effective manner (50–52). Zebrafish have previously been used to discover and assay disease-related natural products (53–55) with broad relevance. Collectively, the zebrafish-IHP model can be highly controlled to investigate intricate routes by which the gut microbiome may mediate parasite infection, and insights gleaned from connections between gut microbiota and gut metabolites may offer translational potential for understanding the metabolite-based interactions between the gut microbiota, intestinal helminth parasites, and the host.

To uncover gut microbial metabolites that mediate IHP infections, we used a zebrafish model of intestinal helminth infection by the nematode Pseudocapillaria tomentosa (49, 56), as it affords access to the large sample sizes needed to disentangle these relationships. In particular, instead of housing fish together we individually housed parasite-infected fish hosts to understand why a small number of zebrafish bear a disproportionate parasite burden in the absence of social or co-housing dynamics, and produced paired microbiome and metabolome data from infected and uninfected fish to investigate links between the microbiome, parasitic infection, and intestinal metabolites. We observed that the gut microbiome explains the variation in infection burden across individuals. We then utilized a mediation inference framework to identify microbe-metabolite interactions that statistically mediate worm burden. This work reveals a potent anthelmintic, salicylaldehyde, whose effect on infection burden is mediated by members of the gut microbiome. Analysis of the paired microbiome-metabolome data also implicates N-acylethanolamines (NAEs) in the association between microbiota and parasite infection burden. Collectively, our work demonstrates that the zebrafish gut microbiome metabolically mediates IHP infection outcomes and reveals novel microbiome-sourced anthelmintic drug leads.

RESULTS

Intestinal helminthic parasite infections are overdispersed among socially isolated zebrafish

IHP infections frequently manifest overdispersed distributions across wildlife populations, agricultural settings, and scientific laboratories (57–59). Prior work has shown that social behavior and interactions can drive differences in parasite infection burden (60–62) but it remains unclear if such behavior and interactions underlie the distribution of burden. For example, in zebrafish, social hierarchies and behaviors (63–65) may impact feeding, which could result in interindividual biases in oral exposure to infectious agents. To understand whether zebrafish parasite burden overdispersion occurs in the absence of co-housing or social dynamics, we individually housed 100 zebrafish in 1.2 L tanks and exposed 50 fish to P. tomentosa eggs. Stool samples were collected from all surviving individuals at several time points: prior to P. tomentosa exposure, immediately before exposure, and 29 days following parasite exposure, which prior work has shown is the peak of infection (49, 66). At this final time, point fish were sacrificed and infection burden was quantified through cytological analysis of dissected intestinal tissue.

Interrogation of the distribution of infection burden among P. tomentosa-exposed fish 29 days post-exposure (dpe) reveals that infection is overdispersed across the population (σ²/μ = 4.719, Fig. S1), indicating that relatively small numbers of exposed individuals carry the bulk of mature worms in their guts. Our investigation of socially isolated individuals reveals that community social dynamics alone are not the sole drivers of overdispersed helminth parasite worm burden among zebrafish populations and indicates that other factors underlie this phenomenon.

Gut microbiome composition associated with parasite exposure and infection burden

Gut microbiomes display highly personalized forms across individuals (67, 68), and variation in gut microbial communities has been consistently linked to parasite infection (66, 69, 70). Thus, it is conceivable that overdispersion in parasite infection burden outcomes results from the intricate interplay between parasite exposure and gut microbiome composition that occurs in each individual, where bacterial consortia could tip the scales toward susceptibility or resilience. We investigated whether zebrafish gut microbiome community structure is related to parasite exposure and subsequent overdispersion of parasite infection burden. To investigate how gut microbiome composition changes before and after parasite exposure, we generated 16S rRNA gene sequence data from stool samples collected at both a pre-exposure baseline and at 29 days post-exposure (dpe). To test whether the initial microbial community state influences the microbiome’s association with infection burden at the peak of infection, we initially altered fish gut microbial communities by administering antibiotics three days prior to P. tomentosa parasite exposure (Fig. 1). This strategy was employed to ensure variability in the gut microbiome compositions across the individually housed fish, which was necessary for analyzing the potential role of the gut microbiome in parasite infection. The resultant microbial diversity enabled us to explore the association between distinct microbiome profiles and the differential parasite burden outcomes. Microbial communities at 29 dpe were significantly stratified by both parasite exposure (PERMANOVA, F = 7.1618, P < 0.0001) as well as parasite burden (PERMANOVA, F = 12.1514, P < 0.0001). These results are consistent with our earlier findings (66), but in this case individually housing fish eliminates possible impacts of co-housing that may drive homogeneity of the gut microbiome among infected versus non-infected individuals. Thus, this design and these results provide particularly compelling evidence that zebrafish gut microbial communities are connected to P. tomentosa and later infection success, and may contribute to the overdispersion of infection burden across individuals.

Fig 1 Schematic of zebrafish husbandry and treatment events and timeline. (1) Briefly, 100 adult fish were placed in individual tanks, (2b) half of the fish were subsequently exposed to antibiotics, (3b) then fish were randomly exposed to the zebrafish parasite Pseudocapillaria tomentosa. Fecal samples were collected (2a) prior to antibiotic exposure, (3a) just prior to parasite exposure, and (4) 29 days post-parasite exposure (dpe) after which fish intestinal histopathology was assessed. Samples were split and processed for untargeted fecal metabolomic analysis as well as fecal 16S rRNA DNA amplicon sequencing.

Study involves 100 fish reared individually, split into antibiotic and control groups, then further divided into parasite and control treatments. Samples and checks include fecal 16S, metabolome, and histopathology at different time points.

Effect of IHP exposure on the gut microbiome composition depends on the pre-exposure microbiome state

Given that parasite colonization and the effects of parasite infection have been linked with gut microbiome composition in this work and elsewhere (71), we reasoned that altering the initial state of the microbiome may shape its relationship to subsequent IHP infection. Analysis of 16S rRNA gene sequences generated from fish stool collected after this antibiotic exposure but prior to P. tomentosa exposure reveals that antibiotic administration successfully altered the composition of the zebrafish gut microbiome (PERMANOVA, F = 27.565, P < 0.0001). A corresponding analysis at 29 dpe shows that the relationship between fecal microbial community composition and IHP exposure depends on this prior antibiotic exposure (PERMANOVA, F = 3.16, P = 0.009). Moreover, we find that parasite infection burden is strongly linked to the composition of the microbiome (PERMANOVA, F = 7.1618, P < 0.0001, Fig. 2) but that this relationship fundamentally depends on whether hosts were exposed to antibiotics first (PERMANOVA, F = 4.2087, P = 0.002). This interaction is particularly noteworthy because at 29 dpe no strong relationship is observed between the composition of the gut microbiome and antibiotic exposure alone (PERMANOVA, F = 1.32, P = 0.2277). These findings collectively suggest that perturbations to the initial state of the microbiome, such as through antibiotic exposure, have a cryptic effect on the successional interplay between IHP infection and the gut microbiome, even when the statistical effects of antibiotic exposure are no longer prominently observed.

Fig 2 Principal coordinates analysis (PCoA) ordination of Bray-Curtis dissimilarity of microbiome communities at 29 days following parasite exposure. Each point represents an individual fish. The halo intensity around points represents the number of quantified parasites in the gut at dissection. Point colors represent antibiotic-exposed (blue) and -unexposed (red) groups. The arrow illustrates an envfit relationship for worm burden, depicting the linear direction of association between parasite burden and Bray-Curtis dissimilarity.

Plot visualizes the PCoA of microbial communities and depicts two principal coordinates explaining 33.51 and 21.46 percent of the variance. Points are colored by antibiotic exposure and point halo transparency by worm burden.

This observation may be of special interest given that geographic locales which have higher levels of IHP infection also tend to be locations where microbiomes may be disrupted by the use of antibiotics to manage bacterial infections (72, 73). Given the varied influential roles of the microbiome on parasite infection (26–29), disruption of the initial microbiome state by antibiotics may interfere with the ability of the microbiome to protect the host from helminth infection. Elucidating how antibiotic use alters the microbes involved in helminth resistance, as well as their specific interactions with parasites, is crucial. This knowledge could guide the development of microbiome-based therapies that supplement those protective elements and reduce the adverse impacts of helminth infection.

Fecal metabolites, including salicylaldehyde and N-acylethanolamines, inversely associate with helminth infection burden

Metabolic products convey information about the presence of pathogens and microbes (70, 74), produce signals that underlie immune control (75), and interact within a complex network of gut microbes, the host immune system, and invading parasites, driving shifts in the intestinal environment. To understand the metabolomic landscape wherein zebrafish gut microbes and parasites co-locate, we performed untargeted metabolomic profiling of the fecal samples collected from the same individuals and at the same time point as our microbiome profiling analysis. Prior research has shown that the zebrafish gut metabolome is composed of a diverse array of lipids and fatty acids (76, 77), as well as amino acids and various biogenic amines (78, 79). Consistent with previous research, our annotated gut metabolomic data are also dominated by a large number of complex lipids, vitamins, and amino acid derivatives, as well as polar metabolites from many compound classes. We first sought to find metabolites from this diverse metabolite set that are statistically associated with parasitic worm burden. Due to the high level of overdispersion frequently observed in parasite infection data (57–59), we utilized negative-binomial generalized linear models (GLMs) to examine the statistical relationship between a set of 303 annotated metabolites and P. tomentosa worm burden. We uncovered 35 metabolites which associated with parasite worm burden among infected hosts (Negative Binomial GLM, FDR < 0.1, Table S1).

Numerous compounds that are associated with helminth parasite worm burden have also been linked to parasite infection in other work. Salicylaldehyde, a compound previously noted as a soil and plant nematicide (80–82), shares an inverse relationship with worm burden. This analysis also identifies compounds that have previously been linked to parasite infection such as a major form of Vitamin E, gamma-tocopherol (83). In addition, at 29 dpe, we find that six of the eight compounds classed as N-acylethanolamines (NAEs) such as oleoyl ethanolamide (OEA), linoleoyl ethanolamide (LEA), 2-linoleoyl glycerol, and related NAE precursor compounds (e.g., glycerophospho-N-oleoyl ethanolamine), manifest an abundance profile that sharply distinguishes infected versus uninfected individuals, where infected individuals display higher metabolite abundances (Wilcoxon-Rank Sum Test, P < 0.05, Fig. 3A). Furthermore, six of the eight NAE-related compounds in these data are strongly inversely associated with parasite worm burden (Negative Binomial GLM, FDR < 0.1, Fig. 3B; Table S1).

Fig 3 (A) The abundance of NAE-related compounds significantly differs between infected and uninfected fish. (B) NAE abundance inversely associated with parasite exposure for six of eight identified compounds. “*” indicates P < 0.05, “**” P < 0.01, and “***” P < 0.001.

Box plots of log metabolite abundance for various compounds compare actively infected and no detected worm burden groups, with significance levels indicated. Heatmap of z-scaled metabolite abundance across samples, with worm count data alongside.

The NAE class of metabolites represents a broad family of lipid messengers that play a well-established role in energy metabolism and feeding behavior (84–86), as well as inflammation (85), and prior work has established a relationship between the gut microbiome and NAEs (87–89), so we next tested if the abundance of these metabolites is related to microbiome composition and antibiotic exposure. Strikingly, we observed a robust association between the abundance of these NAE compounds and gut microbiome composition at a time point prior to parasite exposure and following antibiotic treatment (PERMANOVA, P < 0.05, Table S2). Notably, at this time, there is evidence that the relationship between the zebrafish gut microbiome and the abundance of six NAEs depends on whether or not fish were exposed to prior antibiotics (PERMANOVA, P < 0.1, Table S2). These results are further underscored by our finding that the composition of the gut microbiome at 29 dpe is still strongly linked to NAE abundance for six different compounds (PERMANOVA, P < 0.0002, Table S2). Furthermore, the relationship between the gut microbiome and NAEs is underscored by our results that show antibiotic exposure interacts with NAE abundance in a manner that is significantly related to gut microbiome composition for five of the eight NAE metabolites (PERMANOVA, P < 0.0002, Table S2).

These findings reinforce the emerging view that the gut microbiome plays a fundamental role in regulating NAE levels, an observation which is especially noteworthy because recent work has shown that intestinal nematodes that infect humans, mice, and even insects have genes that encode functions for degradation of NAEs (90). The collective impact of our NAE analysis shows that their abundance is starkly different in parasite-uninfected versus infected hosts, that NAE abundance is linearly related to parasite infection burden, and that the gut microbiome is a principal driver of NAE abundance in zebrafish hosts. Given the multifaceted role of NAEs in host physiology, immune response, and intestinal microbiome control, plus the ability of IHPs to degrade and produce these compounds, emphasizes the significance of these compounds as a nexus in the battle between helminth parasite and vertebrate host.

Connections between fecal metabolite abundance and P. tomentosa worm burden are putatively mediated by fecal microbiota

The prior analysis linking metabolites to parasite worm burden highlighted several compounds that are also known to drive changes in intestinal microbiome composition. Given the complex interplay between gut microbes and metabolite production, these findings open a line of inquiry into microbiome-metabolite interactions. Zebrafish gut microbial taxa have been linked to parasite infection (66), and given the connection between some of the aforementioned metabolites, such as NAEs, we hypothesized that the relationship between members of the gut microbial community and parasite infection depends on metabolite-related cross-talk. To explore the interconnected role of metabolites, microbiota, and parasite worm burden, we selected parasite burden-linked metabolites and prevalent taxa, then statistically analyzed individual metabolite-microbe pairings to identify relationships that may be relevant to parasite burden. Our workflow applied partial correlation to a set of worm burden-linked metabolites and prevalent amplicon sequence variants (ASVs) to identify microbe-metabolite associations which were robust after accounting for the controlling influence of other microbial taxa. In addition, we used mediation inference methods to quantify whether a metabolite’s statistical relationship to worm burden may be hypothetically mediated by members of the microbiota. The results of this approach provided a set of potential interactions between 25 metabolites and 17 members of the microbiota (Adjusted Causal Mediating Effect FDR < 0.3, Fig. 4), whose microbe-metabolite relationship is uniquely strong in the context of other ASVs, and whose interacting relationship may be relevant to parasite infection.

Fig 4 Network of microbe-metabolite interactions predicted to mediate parasite worm burden. Nodes represent fecal metabolites and gut bacteria. Edges represent statistically significant relationships, with colors indicating the direction of correlation (blue: positive, red: negative).

Bipartite network diagram illustrates the associations between various metabolites and bacterial ASVs. Lines represent significant correlations, with color intensity indicating the strength and direction of the associations.

One particular edge in this possible interaction set includes salicylaldehyde and an ASV from the Pelomonas genera. Salicylaldehyde, also known as 2-hydroxybenzaldehyde, is an organic compound that occurs naturally in some foods such as buckwheat (91) and it is known that some salicylaldehyde derivatives exhibit antibacterial activity (92). As described earlier, salicylaldehyde manifests a robust inverse relationship among parasite exposed fish 29 days following parasite exposure (Fig. 5A) and the microbe to which it is correlated, ASV 4, is one member of a small set of microorganisms in this work that is negatively related to total helminth worm burden (Negative Binomial GLM P = 0.002, Fig. 5B). In addition to this inverse relationship following parasite infection, we also find that it is a particularly important feature for predicting later parasite infection burden. We constructed a random forest regression model considering parasite worm burden as a function of all gut microbial taxa present prior to parasite exposure and found that Pelomonas, specifically ASV 4, shows up as one of the most important taxa for predicting parasite worm burden among a feature set that includes hundreds of different microbial taxa (Fig. S2). This same ASV is also strongly correlated with salicylaldehyde abundance (Spearman’s Correlation, ρ = 0.62, P = 0.002, Fig. 5C) and is predicted to mediate the relationship between salicylaldehyde and worm burden (Average Causal Mediation Effect (ACME) FDR < 0.3). Little is known about the function and metabolic capacity of Pelomonas, but given its relationship with salicylaldehyde and its strong inverse relationship with helminth worm burden we surveyed available Pelomonas genomes available in the NCBI genome repository to understand if the taxonomic group possesses genetic pathways associated with salicylaldehyde metabolism. Several available Pelomonas genomes possess salicylaldehyde dehydrogenase (93), an enzyme responsible for metabolism of salicylaldehyde, typically as part of the naphthalene degradation pathway (94), and it has been long demonstrated that some taxa from the Pelomonas genera are capable of metabolizing diverse polycyclic aromatic hydrocarbons (PAHs) including naphthalene (95). Pelomonas metabolic capacity to interact with salicylaldehyde could very well explain the predicted mediating role of this taxon in the relationship between salicylaldehyde and worm burden.

Fig 5 (A) Scatterplot of log salicylaldehyde abundance against worm burden among fish that are 29 dpe. (B) Scatterplot of Pelomonas ASV 4 relative abundance against worm burden among parasite exposed fish 29 dpe. (C) The linear relationship between log salicylaldehyde abundance and Pelomonas relative abundance. (D) P. tomentosa eggs were exposed to salicylaldehyde at a concentration of 2 mg/L in an in vitro assay. The y-axis depicts the % of eggs that are unlarvated or dead. (E) In vivo salicylaldehyde exposure assay comparing the number of mature female P. tomentosa worms that produced eggs.

Statistical plots and bar charts depict relationships between worm burden, salicylaldehyde abundance, and Pelomonas relative abundance. It also depicts the percentage of dead or unlarvated eggs and the count of mature female worms in samples.

We also found that the relationship between γ-tocopherol, which is a major form of vitamin E, and parasite worm burden is putatively mediated by ASVs from the Mycobacterium and Pseudomonas genera and parasite worm burden. The abundance of these ASVs was positively linked to worm burden in this work (Negative Binomial GLM, FDR < 0.05), and in an earlier investigation, the genera-level abundance of both of these taxa displayed a positive correlation with IHP burden in zebrafish (66). γ-tocopherol has been shown to modulate the composition of the gut microbiome (83, 96) and mitigate colitis caused by LPS-induced inflammation signaling in mice (97), and reduced abundance of this metabolite leaves mouse hosts more susceptible to helminth infections (98, 99). If this compound displays similar colitis mitigating and gut microbiome-altering effects in zebrafish, reduced abundance of this metabolite, either by parasite infection or through metabolism by parasite-promoting taxa, may help clarify the basis through which these members of the microbiota associate with increased worm burden.

Salicylaldehyde administration completely inhibits P. tomentosa egg larvation

To determine the anthelmintic effect of salicylaldehyde, we utilized both in vitro and in vivo drug exposure assays to determine how exposure to the drug impacts P. tomentosa. One mechanism common to many anthelmintic compounds such as albendazole and ivermectin is inhibition of worm egg larvation (100). Therefore, we innovated an in vitro assay to evaluate how salicylaldehyde exposure impacts P. tomentosa egg larvation. Specifically, exposed 131 P. tomentosa eggs to salicylaldehyde at a dose of 2 mg/L, 136 eggs at 15 mg/L, and a group of 242 eggs that were reared without exposure to salicylaldehyde. For each group, we quantified larvation rates after 5 days. In the salicylaldehyde-unexposed group, 174/242 (72%) of the eggs are successfully larvated. Conversely, 0 eggs from either salicylaldehyde-treated groups larvated (Fisher’s Exact Test, P < 2e-16, Fig. 5D). The complete inhibition of egg larvation demonstrates notable inhibition of egg hatching in P. tomentosa.

We followed up on this finding by determining if salicylaldehyde can disrupt active P. tomentosa infection in vivo. To do so, we exposed 60 fish to P. tomentosa parasite eggs, then split the fish equally into 30 fish that were exposed to salicylaldehyde at a concentration of 2 mg/L, and 30 fish that were not exposed to salicylaldehyde, but to a DMSO control. Mortality in salicylaldehyde exposed and unexposed host fish was 10 and 5 fish, respectively, which started at 16 dpe. All examined moribund or fresh dead fish showed infections. The experiment was terminated at 24 dpe, and the intestines of all surviving fish showed 100% infection in both groups. Results indicate a slight reduction in infection burden among salicylaldehyde-exposed fish, with an average of 15 worms/fish in salicylaldehyde-treated fish versus an average of 19 worms/fish among untreated controls (Wilcoxon Rank-Sum Test, P = 0.3; Fig. 5E). Although more strikingly, we found that drug treatment interfered with worm development: 7/25 (28%) the fish that were not exposed to salicylaldehyde contained gravid female worms with eggs in their guts, whereas no female parasites or free eggs were recovered from the salicylaldehyde-exposed treatment group (Fisher’s Exact Test, P = 0.01, Fig. 5E).

Overall, these assays indicate that salicylaldehyde represents an enticing anthelmintic drug lead, although further work is required to clarify the mechanism by which it acts, as well as to explore the specificity of IHP types it is active against. Furthermore, elucidating the potential salicylaldehyde-producing capability of Pelomonas or other microbes may offer a microbiome-oriented route for the control of parasite infection.

DISCUSSION

The rise of anthelmintic-resistant IHPs presents an exigent challenge to identify new drugs and techniques to control IHP infection. To do so, we must expand our understanding of the factors that underlie infection, as this knowledge can be leveraged to design new biocontrol strategies. Based on the historical discovery of novel anthelmintic compounds among bacteria (19) as well as accumulating evidence that the intestinal microbiome affects the colonization success of IHPs (18, 69), we reasoned that exploring the metabolic landscape of the gastrointestinal microbiota in the context of infection could unlock new leads in the quest for novel anthelmintic strategies. To explore this idea, we investigated the relationship between gut microbes, metabolites, and IHP infection in a zebrafish model using a multi-omic approach to statistical mediation. In so doing, we identified a variety of metabolites that associated with infection burden in a manner that is dependent on specific microbial taxa. These relationships are valuable to resolve not only because the metabolites and taxa in question may be utilized to develop novel infection control strategies (e.g., anthelmintic drugs or probiotics) but because they underscore the putative mechanisms by which the microbiome influences IHP infection outcomes.

Experimental validation of one particular lead, salicylaldehyde, reveals that our approach can uncover microbially mediated compounds that elicit anthelmintic effects. Salicylaldehyde is a nematicide used in agricultural settings and our data is among the metabolites that are most strongly inversely associated with worm burden. We performed follow-up in vitro and in vivo tests to confirm that salicylaldehyde elicits effects against P. tomentosa. In particular, this compound completely inhibits in vitro P. tomentosa egg development and maturation. Cessation of egg production following anthelmintic treatment is a common phenomenon with gastrointestinal nematodes in production animals. For example, moxidectin has been demonstrated to inhibit egg production in female worms of Cooperia that survive treatment in the early stages of resistance (101). We observe a similar phenomenon here, where there is a complete absence of parasite worm eggs recovered from hosts that are exposed to salicylaldehyde. Previous work has also shown that salicylaldehyde prevented hatching of the potato cyst nematode, Globodera pallida (81), but this compound is not currently used for the control of helminth infections in animal populations, and our experiments are the first to demonstrate that this compound is also capable of inhibiting egg maturation in a parasite which infects vertebrate organisms. While the precise mechanism of this action is uncertain, P. tomentosa eggs have relatively delicate shells and are quite susceptible to desiccation and chemical agents (102). Thus, it is possible that salicylaldehyde either damages the egg shell directly, or possibly translocates to unlarvated worms where it may impair their larvation physiology. Future work should seek to uncover the specific mechanisms of action and whether salicylaldehyde can elicit broad effects across IHP species that infect other hosts. Regardless, our study suggests that repurposing salicylaldehyde as an anthelmintic drug against vertebrate IHPs may help control infection in the face of rising multi-drug-resistant IHPs.

As zebrafish are not capable of producing salicylaldehyde, we used statistical mediation to identify a microbe from the Pelomonas genus as a potential source of the compound. In particular, our mediation analysis finds that the relative abundance of Pelomonas positively correlates with the abundance of salicylaldehyde and significantly explains the variation in the underlying relationship between the compound and infection burden. While a variety of processes may underlie these associations, prior work supports the hypothesis that these patterns result from Pelomonas-induced metabolism of salicylaldehyde. The Pelomonas genus is known to possess enzymes that may aid in the metabolism of salicylaldehyde, such as salicylaldehyde dehydrogenase (93), and is a member of a diverse set of Betaproteobacteria which are known to generate and utilize this compound during the degradation of naphthalene (94). While the degradation of naphthalene represents one parsimonious explanation for the origin of salicylaldehyde, another important alternative hypothesis regarding the origin and synthesis of salicylaldehyde begins with phenylalanine. The results of our analysis point to phenylalanine as a metabolite which is inversely related to parasite worm burden and whose relationship with parasite worm burden is potentially mediated by Pelomonas (Fig. 4). The conversion of phenylalanine to trans-cinnamic acid is known to be performed by bacteria (103–105). The following in a phenylalanine to salicylaldehyde metabolism might first require the conversion of trans-cinnamic acid to o-coumaric acid. The terminal conversion of o-coumaric acid to salicylaldehyde is known in tobacco plants (106), as well as in a species of fungus (107), although bacterial catalysis of this terminal reaction is not characterized. While this might represent a parsimonious explanation for the biosynthesis of salicylaldehyde, more work is warranted to establish the distribution of microbial participation in these pathways, especially with respect to Pelomonas. Collectively, future research should seek to establish the precise route by which Pelomonas synthesizes salicylaldehyde to affect IHP infections and whether related mechanisms exist among the microbiota of other vertebrate hosts.

Our multi-omic analysis reveals substantial alterations in NAEs and NAE precursors between parasitic infected and uninfected hosts, with a potential route of control by the gut microbiome. Given that NAEs play in a broad range of physiological functions such as immune regulation as well as energy metabolism and feeding behavior (84–86, 108), IHPs are capable of modulating NAEs to enhance infection (90), and their alteration is associated with changes in the gut microbiome (87), these compounds represent a potentially rich area to understand the intersection of parasite infection with the gut microbiome and host immune regulation. Our results demonstrate that alteration of the gut microbiome by antibiotic exposure appears to drive changes in NAE abundance that are sustained in the profile of several NAE compounds several weeks after initial antibiotic exposure. The abundance of these compounds is also linked to IHP infection and is linearly related to IHP infection burden. The NAE axis represents one route by which the intestinal microbiome drives aspects of host physiology. The functions which are regulated by NAE signaling such as feeding behavior and inflammation response are likely relevant in explaining some aspects of parasite infection, especially given that parasites also possess enzymes involved in regulating NAEs and their associated physiological functions. The synthesis of these findings is complex, but elucidating the principal effects of NAE changes on host physiology and parasite infection, in addition to identifying taxa whose abundance is related to NAE changes, may reveal how intestinal microbiota participate in mediating host response to IHPs and could reveal drug or probiotic targets that aid in the control of helminth infection.

While we highlight several compounds, such as salicylaldehyde and NAEs, which help explain the relationship between the gut microbiome and parasite infection, there exists a rich unexplored repository of metabolites in this data that may display similar anthelmintic activity. To better understand the metabolomic landscape relevant to IHP infection, we modeled the relationship between parasite burden and metabolite abundance and elucidated a collection of compounds that may be harnessed and further investigated in efforts to control parasitic colonization and success. For instance, trunkamide A, a compound that is of known bacterial origin (109), and which has been examined for its antibiotic and antitumor activity (110), shows up as being inversely associated with parasite burden in our data. Given the established ability of this compound to be produced by bacteria, understanding its distribution among aquatic and gastrointestinally associated bacteria may offer an enticing route to the discovery of novel probiotic microbes that could be supplemented and stimulated to produce trunkamide A under specific parasite-related conditions. In addition, several compounds that have been explored for anthelmintic activity, including baliospermin (111), and genistin (112), also show up in our results as being significantly linked to parasite burden. Some of these metabolites may be of bacterial origin, or modified by members of the gut microbiota, such as genistin (113), in a manner that enhances their anthelmintic activity. Applying mediation methods to understand the possible relationships between these compounds and microbes involved in parasite infection may help to establish microbe-dependent routes for the control of parasite infection. Collectively, these findings represent an additional set of metabolomic compounds that may be explored by mining the gut microbiome for potential solutions to the urgent challenge posed by increasing anthelmintic drug resistance.

In summary, this work unravels interactions within the zebrafish gut ecosystem, yielding a deeper understanding of the dynamic relationships among microbiota, metabolites, and parasitic infections. These results extend the application of mediation inference methodologies to reveal specific bacterial metabolites that may serve as key mediators of host-parasite interactions and identified novel anthelmintic drug leads. Notably, salicylaldehyde emerges as a compelling anthelmintic compound, and we demonstrate its ability to inhibit parasitic egg maturation in zebrafish. This work also establishes the involvement of other metabolites, like N-acylethanolamines, in host-parasite-microbiome dynamics emphasizing the need for further research to elucidate their influence. Collectively, our findings support the hypothesis that gut microbiota play a role in parasite infection, and understanding the chemical means by which microbiota are involved in helminth colonization may yield tools for infection control.

MATERIALS AND METHODS

Zebrafish husbandry, parasite exposure, and parasite scoring

All zebrafish research was conducted under the approval of IACUC protocol 2022-0280. Tropical 170-day-old 5D zebrafish were obtained from the Sinnhuber Aquatic Research Laboratory (Corvallis, OR). The fish were housed in a flow-through vivarium on a 14 h:10 h light:dark cycle, and fed Gemma Micro 300 once a day, except on the weekends. Water temperature was recorded daily and ranged from 23 to 28ºC. There was also weekly testing of ammonia (0–0.25 ppm), pH (7.6), hardness (0–25 ppm), and conductivity (90–110 uS/cm) to ensure high water quality. Prior to initial antibiotic exposure or parasite exposure, fish were allowed to equilibrate in their tanks for 2 weeks. Each of the 100 fish used in this protocol was randomly assigned to one of four unique exposure groups, no parasite and no antibiotic exposure, no parasite and antibiotic exposure, parasite exposure and no antibiotics, and both parasite and antibiotic exposure. After this period of equilibration, 50 out of the 100 fish were exposed to a combination of antibiotics, 10 mg/L colistin and 10 mg/L vancomycin. Of the 100 fish used in the experiment, 50 were exposed to P. tomentosa eggs at a dose of 88 eggs per fish.

Zebrafish gut metabolite mass spectrometry

Fecal pellets collected from individual zebrafish were promptly lyophilized to minimize the leaching of metabolomic content into water and then split approximately equally for use in metabolomic and 16S rRNA library preparation and sequencing analyses. Due to the small sample size and high water content of the fecal pellets, precise weighing was not feasible. Therefore, total ion abundance normalization was implemented to account for variation in sample input. Ivermectin and two isotopically labeled amino acids were incorporated into the extraction mix to monitor injection accuracy and platform performance throughout the extended batch run time. Zebrafish fecal samples were prepared for untargeted metabolomic analysis using a modified extraction protocol. An extraction solvent consisting of equal parts 100% ethanol (Sigma 1.11727.1000) and methanol (Fisher A456-1) was prepared and chilled overnight at −20°C. Three 1.4 mm zirconium oxide beads (VWR 10144-554) were added to each 2 mL screw-top bead beating tube (Fisher 02-682-558). An extraction mix was then created by diluting isotope standards (Cambridge Isotope Laboratories MSK-A2.1.2; 1:20 dilution) and 100 µM Ivermectin standard (Sigma PHR1380) 1:10 in the chilled ethanol/methanol solution. Twenty-five microliters of this extraction mix was added to each tube containing a fecal sample. The samples were then homogenized using a Precelly’s 24 bead beater (program: 2 × 5,400 rpm for 45 s, 5 s wait interval), centrifuged at 16,000 × g for 10 minutes at 4°C, and the resulting supernatant (15–25 µL) was transferred to glass vials (Microsolv 9532S-3CP-RS). In cases of large sample volumes, an additional centrifugation step was performed.

Extracts were submitted for analysis using untargeted LC-HRMS/MS-based metabolomics. An AB Sciex TripleTOF 5600 mass spectrometer coupled to a Shimadzu Nexera UHPLC system was used as described previously (114). Metabolite extracts were separated using an Inertsil Phenyl-3 column (2.1 × 140 mm, 100 Angstrom, 5 µm; GL Sciences, Torrance, CA, USA). Column was held at 50°C. We used a 50-minute binary gradient system consisting of: Solvent A, water (LC-MS grade) with 0.1% vol/vol formic acid and solvent B, methanol (LC-MS grade) with 0.1% vol/vol formic acid. Metabolites were eluted using the following gradient program: 1 min, 5% B; 11 min, 30% B; 23 min, 100% B; 35 min, 100% B; 37 min, 5% B; 37 min, 5% B; and 50 min, 5%B. The flow rate was 0.5 mL/min. Injection volume was 10 µL. Positive IonSpray voltage was set to 5,200 V, negative IonSpray voltage was set to 4,200 V. Source temperature was 350°C. The q-TOF mass spectrometer was operated in the data-dependent mode using the following settings: period cycle time 950 ms; accumulation time 100 ms; m/z scan range 50–1,200 Da; and collision energy 40 V. Mass calibration of the TOF analyzer was performed automatically after every fifth LC run.

LC-HRMS/MS data processing was performed with Progenesis QI software V2.0 (NonLinear Dynamics, United Kingdom) and ABSciex Masterview (ABSciex, USA) entailing peak picking, retention time correction, peak alignment, and metabolite identifications/annotations. Metabolite annotations were facilitated by Progenesis QI and Masterview using an in-house spectral library based on the IROA Mass Spectrometry Metabolite Library of Standards (MSMLS) containing retention times, exact mass, and MS/MS information of >650 metabolite standards (IROA Technologies, Bolton, MA, USA). This workflow allows obtaining high confidence annotations (L1). In addition, tentative metabolite annotations were obtained by searching the METLIN MS/MS library in Progenesis QI.

Microbial 16S rRNA library preparation and sequencing

Fecal samples were collected from individual adult fish at Days 0, 3, and 32. Fish mortality and parasite exposure precluded the collection of fecal pellets from every fish at all time points. For samples and time points at which fecal samples could be collected, DNA was isolated using the Qiagen DNeasy PowerSoil kit, in accordance with the manufacturer’s directions. After a 10-minute incubation at 65°C, samples were subjected to bead beating, using 0.7 mm garnet beads, for 20 minutes using the Vortex Genie 2 (Fisher, Hampton, NH, USA). PCR was performed in triplicate using 1 microliter of purified DNA from the lysis solution to amplify the V4 region of the 16S rRNA gene, using the 806 r and 515 f primer set (115). Amplified DNA collections were quantified using a Qubit HS Kit (Carlsbad, CA, USA). An equal quantity of DNA was selected from each of the 252 samples, for a total DNA mass of 200 ng, and the pooled collection of DNA was cleaned using the QIAGEN QIAquick PCR purification kit then diluted to a final concentration of 10 nM. The final pooled DNA collection was sequenced by the Center for Quantitative Life Sciences at Oregon State University, using an Illumina MiSeq instrument with 250 bp paired-end reads.

Zebrafish gut microbiome community diversity analyses

Read quality filtering was performed using DADA2 (116) with R 4.1.2. Alpha and beta-diversity analyses were performed using a relative abundance-normalized sequence count table. Generalized linear mixed-effects models were used to model species richness and Shannon diversity as a function of longitudinal timepoint, fish id, antibiotic exposure, parasite exposure, and an interaction of antibiotic and parasite exposure. Bray-Curtis dissimilarity and subsequent NMDS ordination were performed using vegan (117). To clarify how antibiotic administration, parasite exposure, and other host factors relate to gut community composition, we used permutational multivariate analysis of variance (PERMANOVA, adonis2, vegan) with 10,000 permutations.

Mediation analysis of gut bacteria and metabolites

Regression-based mediation analysis was used to investigate the hypothesis that the relationship between fecal metabolites to nematode parasite burden is mediated by the abundance of gut microbiota. Briefly, the standard approach for these techniques employs regression modeling to analyze the association between two variables, then models the possible effect of a mediating variable on the relationship between the variables in the initial model. We constrained the initial feature space by selecting metabolites for which annotation and research-characterized biological identity were available. The sparsity of gut bacterial data challenges statistical investigation of correlation, so only taxa whose relative abundance was greater than 0 in more than 30% of samples were used. Microbiome data were normalized using relative abundance, and the log of metabolite MS abundance was log-transformed by adding 1 to initial values which were equal to 0. Initial feature selection yielded 40 prominent bacterial taxa, and 27 metabolite compounds which were significantly related to worm burden (NB-GLM, FDR < 0.1). We utilized the nptest package (118) to test for partial correlation for pairwise relationships between metabolites and bacterial taxa. We applied partial correlation for each microbe-metabolite pairing and used all remaining taxa as conditioning variables to isolate direct microbe-metabolite links. Pairings with an FDR-adjusted relationship <0.3 were considered. Concurrently, ASV mediation of metabolite-parasite burden was tested using the mediation package (119) in R, where metabolites were coded to represent high and low metabolite abundances, with ASVs as mediators. The average causal mediation effect (ACME), average direct effect (ADE), and proportion of direct effect mediated were calculated for each model of microbe-metabolite pairing. Family-wise error rates were controlled by applying Benjamini-Hochberg FDR correction (FDR < 0.3) to partial correlation and ACME P-values for tested mediating relationships. Visualizations of putative mediating interactions were visualized using ggplot2 (120). Code to recreate the mediation analysis and visualizations is available here (https://github.com/CodingUrsus/Zebrafish_Microbiome_and_Parasites).

Salicylaldehyde toxicity assay

A toxicity assessment identified the highest SA dose which did not result in significant mortality. Adult zebrafish were aqueously exposed to SA in two 48-h periods spaced 5 days apart. The fish were monitored for any adverse health effects including mortality, and fish were evaluated to ensure they were alive and capable of active movement. Toxicity endpoints were evaluated from the initial exposure to 7 days after the last exposure. The chemicals used in this study were salicylaldehyde (SA, cas: 90-02-8) and dimethyl sulfoxide (DMSO, cas: 67-68-5). Salicylaldehyde was obtained from Tokyo Chemical Industries (lot: J052M-EQ) and DMSO was purchased from VWR (lot: 22H2456964). Dilutions were made with 100% DMSO and stored in closed vials in a desiccator at room temperature.

Fish were aqueously exposed to 0,1, 2, 3, 5, and 10 mg/L of SA with 0.01% DMSO. The 0–2 mg/L groups had two replicate tanks containing 16 males and 16 females. Furthermore, the 3–10 mg/L groups had one replicate tank containing 6 males and 6 females. The 3–10 mg/L exposure groups were intended as positive controls, to ensure that null effects were not a result of no chemical exposure.

Exposure groups 0–2 mg/L SAL occurred in 9L tanks first filled with 3L of system water followed by 1 mL of a concentrated SA stock solution. Following the addition of SA, the remaining 3L of fish water was added to the tank to mix the chemical. The 3–10 mg/L exposures were conducted in 2.8L tanks and followed the same chemical additional pattern. However, only a total of 2L of water was added to the tank. Solutions were refreshed every 24 h using the same method previously described. Lastly, each exposure tank was aerated with an air stone and the fish were not fed, to preserve water quality.

Parasite salicylaldehyde exposure assays

In vivo salicylaldehyde exposure

Based on our previous transmission studies (49, 102), we experimentally infected 100 5D line fish by placing them in a large tank from which ~30 infected fish were removed the day before. The recipient fish were exposed in the infection tank for 5 days. P. tomentosa infection was promoted within the tank by reducing the waterflow, not removing the detritus, and keeping fish carcasses in the tank. Over a 5-day period, 1 L of water from another infected tank was added to the exposure tank to further enhance infection. After 5 days of exposure, the remaining fish (~75) were randomly divided into four 9 L tanks, two control tanks and two salicylaldehyde exposure tanks.

Fish were then exposed to salicylaldehyde in the same manner as the toxicity trials in their aquaria with the water turned off at 14–15 and 21–22 days post-initial exposure (dpe). Using the results of the toxicity trial, fish were dosed with either 2 mg/L SA with 0.01% DMSO or 0.01% DMSO (controls). Fish were monitored daily, and fresh dead fish were examined for the presence of the worm. At 24 dpe, the infection status of the fish was evaluated using wet mounts. The fish were first euthanized using a hypodermic shock. Following euthanasia, the gut was dissected out of the fish and placed on a glass slide with a 60 × 24 mm coverslip with about 200 µL of water. The gut was compressed with a cover slip and viewed on a compound microscope at 50 and 100× magnification. While viewing the guts, the number of immature worms and mature female eggs were counted. Each slide was read by two different readers within about 20 minutes, and the average worm counts were used for future analyses.

Egg larvation assay

We tested the ability of SA to inhibit egg larvation through two separate trials in line with previously used protocols (102). Briefly, P. tomentosa eggs were collected by placing fish in a 10 L static tank for 2 days. After 2 days, the water was filtered through a custom 3D-printed filter apparatus fitted with 105, 40, and 25 µm nylon screens. The material retained on the 25 µm screen was collected in 15 mL conical tubes and centrifuged at ~3,000 × g for 5 minutes. The supernatant was decanted until 1 mL of water remained.

In the first trial, collected eggs were exposed to 0, 2, or 15 mg/L of salicylaldehyde with 0.01% DMSO in 15 mL conical tubes. Eggs were transferred into each exposure group by pipetting 200 µL of the egg pellet formed as described above into each conical tube. The solution was homogenized and transferred to a 30°C room for incubation. P. tomentosa egg larvation was quantified after 5 days of SA exposure (~10 days old). A second trial was conducted where eggs were placed in a container on a shaker to aid in the distribution of DMSO control and SA throughout the egg-hatching solution. Eggs were exposed to 0 or 2 mg/L of SA with 0.01% DMSO in 15 mL conical tubes, with the shaker on speed setting 2.5 (Hoefer, San Francisco CA) throughout the exposure. Larvation of eggs was quantified 3 days after SA exposure (~6 days old). In both trials, after exposure, eggs were collected by centrifuging the tubes at 5,000 rpm for 5 minutes. The eggs were subsampled by taking 20 µL of water from the bottom of the tube and placing them on a glass slide and covered with a glass 24 × 24 mm coverslip. The number of larvated and unlarvated/dead eggs was counted using a compound microscope at 50, 100, or 400 × magnification.

ACKNOWLEDGMENTS

We express our gratitude to Kristin Kasschau and Alexandra Alexiev for their insights on microbial growth and metabolism requirements.

This project was supported by NIH NIAID R21 to T.J.S. (R21AI135641), NIH NIEHS R01 to TJS (R01ES030226), NSF grant (2025457) to T.J.S., NIH grant (S10RR027878) to J.F.S., and a Tartar fellowship to A.J.H.

This work used resources in Oregon State University’s mass spectrometry center supported in part by institutional funds.

DATA AVAILABILITY

The nucleotide data underlying the findings of this study are available in the NCBI Sequence Read Archive (SRA) under BioProject ID PRJNA1132310, and annotated metabolomic data from positive and negative ion modes are available here (https://github.com/CodingUrsus/Zebrafish_Microbiome_and_Parasites/).

SUPPLEMENTAL MATERIAL

The following material is available online at https://doi.org/10.1128/msystems.00545-24.

10.1128/msystems.00545-24.SuF1 Figure S1 msystems.00545-24-s0001.tiff

Distribution of mature Pseucapillaria tomentosa worms.

10.1128/msystems.00545-24.SuF2 Figure S2 msystems.00545-24-s0002.tiff

The 10 most important features based on the increase in node purity for regression of helminth worm burden.

10.1128/msystems.00545-24.SuF3 Legends msystems.00545-24-s0003.docx

Legends for supplemental figures.

10.1128/msystems.00545-24.SuF4 Supplemental Tables msystems.00545-24-s0004.docx

Tables S1 and S2.

ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.
==== Refs
REFERENCES

1 Maizels RM, Smits HH, McSorley HJ. 2018. Modulation of host immunity by helminths: the expanding repertoire of parasite effector molecules. Immunity 49 :801–818. doi:10.1016/j.immuni.2018.10.016 30462997
2 Hall A, Hewitt G, Tuffrey V, de Silva N. 2008. A review and meta-analysis of the impact of intestinal worms on child growth and nutrition. Matern Child Nutr 4 Suppl 1 :118–236. doi:10.1111/j.1740-8709.2007.00127.x 18289159
3 Crompton DWT, Nesheim MC. 2002. Nutritional impact of intestinal helminthiasis during the human life cycle. Annu Rev Nutr 22 :35–59. doi:10.1146/annurev.nutr.22.120501.134539 12055337
4 Charlier J, van der Voort M, Kenyon F, Skuce P, Vercruysse J. 2014. Chasing helminths and their economic impact on farmed ruminants. Trends Parasitol 30 :361–367. doi:10.1016/j.pt.2014.04.009 24888669
5 Morgan E, Charlier J, Hendrickx G, Biggeri A, Catalan D, Von Samson-Himmelstjerna G, Demeler J, Müller E, Van Dijk J, Kenyon F, Skuce P, Höglund J, O’Kiely P, Van Ranst B, De Waal T, Rinaldi L, Cringoli G, Hertzberg H, Torgerson P, Wolstenholme A, Vercruysse J. 2013. Global change and helminth infections in grazing ruminants in Europe: impacts, trends and sustainable solutions. Agriculture 3 :484–502. doi:10.3390/agriculture3030484
6 Charlier J, Rinaldi L, Musella V, Ploeger HW, Chartier C, Vineer HR, Hinney B, von Samson-Himmelstjerna G, Băcescu B, Mickiewicz M, Mateus TL, Martinez-Valladares M, Quealy S, Azaizeh H, Sekovska B, Akkari H, Petkevicius S, Hektoen L, Höglund J, Morgan ER, Bartley DJ, Claerebout E. 2020. Initial assessment of the economic burden of major parasitic helminth infections to the ruminant livestock industry in Europe. Prev Vet Med 182 :105103. doi:10.1016/j.prevetmed.2020.105103 32750638
7 Charlier J, Höglund J, von Samson-Himmelstjerna G, Dorny P, Vercruysse J. 2009. Gastrointestinal nematode infections in adult dairy cattle: impact on production, diagnosis and control. Vet Parasitol 164 :70–79. doi:10.1016/j.vetpar.2009.04.012 19414223
8 Roche M, Layrisse M. 1966. The nature and causes of “hookworm anemia.” Am J Trop Med Hyg 15 :1029–1102. doi:10.4269/ajtmh.1966.15.1032 5334793
9 Ezeamama AE, Friedman JF, Acosta LP, Bellinger DC, Langdon GC, Manalo DL, Olveda RM, Kurtis JD, McGarvey ST. 2005. Helminth infection and cognitive impairment among Filipino children. Am J Trop Med Hyg 72 :540–548. doi:10.4269/ajtmh.2005.72.540 15891127
10 Raj E, Calvo-Urbano B, Heffernan C, Halder J, Webster JP. 2022. Systematic review to evaluate a potential association between helminth infection and physical stunting in children. Parasit Vectors 15 :135. doi:10.1186/s13071-022-05235-5 35443698
11 Hotez PJ, Alvarado M, Basáñez M-G, Bolliger I, Bourne R, Boussinesq M, Brooker SJ, Brown AS, Buckle G, Budke CM, et al. . 2014. The global burden of disease study 2010: interpretation and implications for the neglected tropical diseases. PLoS Negl Trop Dis 8 :e2865. doi:10.1371/journal.pntd.0002865 25058013
12 Echevarria F, Borba MFS, Pinheiro AC, Waller PJ, Hansen JW. 1996. The prevalence of anthelmintic resistance in nematode parasites of sheep in Southern Latin America: Brazil. Vet Parasitol 62 :199–206. doi:10.1016/0304-4017(95)00906-x 8686165
13 Van Wyk JA, Stenson MO, Van der Merwe JS, Vorster RJ, Viljoen PG. 1999. Anthelmintic resistance in South Africa: surveys indicate an extremely serious situation in sheep and goat farming. Onderstepoort J Vet Res 66 :273–284.10689700
14 Gasbarre LC. 2014. Anthelmintic resistance in cattle nematodes in the US. Vet Parasitol 204 :3–11. doi:10.1016/j.vetpar.2014.03.017 24798800
15 Geerts S, Gryseels B. 2001. Anthelmintic resistance in human helminths: a review. Trop Med Int Health 6 :915–921. doi:10.1046/j.1365-3156.2001.00774.x 11703846
16 Fissiha W, Kinde MZ. 2021. Anthelmintic resistance and its mechanism: a review. Infect Drug Resist 14 :5403–5410. doi:10.2147/IDR.S332378 34938088
17 Nielsen MK, Kaplan RM, Abbas G, Jabbar A. 2023. Biological implications of long-term anthelmintic treatment: what else besides resistance are we selecting for? Trends Parasitol 39 :945–953. doi:10.1016/j.pt.2023.08.002 37633759
18 Sharpton TJ, Combrink L, Arnold HK, Gaulke CA, Kent M. 2020. Harnessing the gut microbiome in the fight against anthelminthic drug resistance. Curr Opin Microbiol 53 :26–34. doi:10.1016/j.mib.2020.01.017 32114334
19 Burg RW, Miller BM, Baker EE, Birnbaum J, Currie SA, Hartman R, Kong YL, Monaghan RL, Olson G, Putter I, Tunac JB, Wallick H, Stapley EO, Oiwa R, Omura S. 1979. Avermectins, new family of potent anthelmintic agents: producing organism and fermentation. Antimicrob Agents Chemother 15 :361–367. doi:10.1128/AAC.15.3.361 464561
20 Johansson MEV, Jakobsson HE, Holmén-Larsson J, Schütte A, Ermund A, Rodríguez-Piñeiro AM, Arike L, Wising C, Svensson F, Bäckhed F, Hansson GC. 2015. Normalization of host intestinal mucus layers requires long-term microbial colonization. Cell Host Microbe 18 :582–592. doi:10.1016/j.chom.2015.10.007 26526499
21 Petersson J, Schreiber O, Hansson GC, Gendler SJ, Velcich A, Lundberg JO, Roos S, Holm L, Phillipson M. 2011. Importance and regulation of the colonic mucus barrier in a mouse model of colitis. Am J Physiol Gastrointest Liver Physiol 300 :G327–33. doi:10.1152/ajpgi.00422.2010 21109593
22 Jakobsson HE, Rodríguez-Piñeiro AM, Schütte A, Ermund A, Boysen P, Bemark M, Sommer F, Bäckhed F, Hansson GC, Johansson MEV. 2015. The composition of the gut microbiota shapes the colon mucus barrier. EMBO Rep 16 :164–177. doi:10.15252/embr.201439263 25525071
23 Liang L, Liu L, Zhou W, Yang C, Mai G, Li H, Chen Y. 2022. Gut microbiota-derived butyrate regulates gut mucus barrier repair by activating the macrophage/WNT/ERK signaling pathway. Clin Sci 136 :291–307. doi:10.1042/CS20210778
24 Waclawiková B, Codutti A, Alim K, El Aidy S. 2022. Gut microbiota-motility interregulation: insights from in vivo, ex vivo and in silico studies. Gut Microbes 14 :1997296. doi:10.1080/19490976.2021.1997296 34978524
25 Obata Y, Castaño Á, Boeing S, Bon-Frauches AC, Fung C, Fallesen T, de Agüero MG, Yilmaz B, Lopes R, Huseynova A, Horswell S, Maradana MR, Boesmans W, Vanden Berghe P, Murray AJ, Stockinger B, Macpherson AJ, Pachnis V. 2020. Neuronal programming by microbiota regulates intestinal physiology. Nature 578 :284–289. doi:10.1038/s41586-020-1975-8 32025031
26 Kommineni S, Bretl DJ, Lam V, Chakraborty R, Hayward M, Simpson P, Cao Y, Bousounis P, Kristich CJ, Salzman NH. 2015. Bacteriocin production augments niche competition by enterococci in the mammalian gastrointestinal tract. Nature 526 :719–722. doi:10.1038/nature15524 26479034
27 Zipperer A, Konnerth MC, Laux C, Berscheid A, Janek D, Weidenmaier C, Burian M, Schilling NA, Slavetinsky C, Marschal M, Willmann M, Kalbacher H, Schittek B, Brötz-Oesterhelt H, Grond S, Peschel A, Krismer B. 2016. Human commensals producing a novel antibiotic impair pathogen colonization. Nature 535 :511–516. doi:10.1038/nature18634 27466123
28 Sassone-Corsi M, Nuccio S-P, Liu H, Hernandez D, Vu CT, Takahashi AA, Edwards RA, Raffatellu M. 2016. Microcins mediate competition among Enterobacteriaceae in the inflamed gut. Nature 540 :280–283. doi:10.1038/nature20557 27798599
29 Hayes KS, Bancroft AJ, Goldrick M, Portsmouth C, Roberts IS, Grencis RK. 2010. Exploitation of the intestinal microflora by the parasitic nematode Trichuris muris. Science 328 :1391–1394. doi:10.1126/science.1187703 20538949
30 Zheng D, Liwinski T, Elinav E. 2020. Interaction between microbiota and immunity in health and disease. Cell Res 30 :492–506. doi:10.1038/s41422-020-0332-7 32433595
31 Gazzinelli-Guimaraes PH, Nutman TB. 2018. Helminth parasites and immune regulation. F1000Research 7 . doi:10.12688/f1000research.15596.1
32 Maizels RM, McSorley HJ. 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 27476889
33 Baron RM, Kenny DA. 1986. The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol 51 :1173–1182. doi:10.1037//0022-3514.51.6.1173 3806354
34 Sohn MB, Li H. 2019. Compositional mediation analysis for microbiome studies. Ann Appl Stat 13 :661–681. doi:10.1214/18-AOAS1210
35 Carter KM, Lu M, Jiang H, An L. 2020. An information-based approach for mediation analysis on high-dimensional metagenomic data. Front Genet 11 :148. doi:10.3389/fgene.2020.00148 32231681
36 Zhang H, Chen J, Li Z, Liu L. 2021. Testing for mediation effect with application to human microbiome data. Stat Biosci 13 :313–328. doi:10.1007/s12561-019-09253-3 34093887
37 Zhang J, Wei Z, Chen J. 2018. A distance-based approach for testing the mediation effect of the human microbiome. Bioinformatics 34 :1875–1883. doi:10.1093/bioinformatics/bty014 29346509
38 Zhang Q, Linke V, Overmyer KA, Traeger LL, Kasahara K, Miller IJ, Manson DE, Polaske TJ, Kerby RL, Kemis JH, et al. . 2023. Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut. Nat Microbiol 8 :424–440. doi:10.1038/s41564-023-01326-w 36759753
39 Shi H, Ter Horst R, Nielen S, Bloemendaal M, Jaeger M, Joosten I, Koenen H, Joosten LAB, Schweren LJS, Vasquez AA, Netea MG, Buitelaar J. 2022. The gut microbiome as mediator between diet and its impact on immune function. Sci Rep 12 :5149. doi:10.1038/s41598-022-08544-y 35338162
40 Depner M, Taft DH, Kirjavainen PV, Kalanetra KM, Karvonen AM, Peschel S, Schmausser-Hechfellner E, Roduit C, Frei R, Lauener R, Divaret-Chauveau A, Dalphin J-C, Riedler J, Roponen M, Kabesch M, Renz H, Pekkanen J, Farquharson FM, Louis P, Mills DA, von Mutius E, Ege MJ, PASTURE study group. 2020. Maturation of the gut microbiome during the first year of life contributes to the protective farm effect on childhood asthma. Nat Med 26 :1766–1775. doi:10.1038/s41591-020-1095-x 33139948
41 Schoemann AM, Boulton AJ, Short SD. 2017. Determining power and sample size for simple and complex mediation models. Soc Psychol Personal Sci 8 :379–386. doi:10.1177/1948550617715068
42 Fan Y, Støving RK, Berreira Ibraim S, Hyötyläinen T, Thirion F, Arora T, Lyu L, Stankevic E, Hansen TH, Déchelotte P, et al. . 2023. The gut microbiota contributes to the pathogenesis of anorexia nervosa in humans and mice. Nat Microbiol 8 :787–802. doi:10.1038/s41564-023-01355-5 37069399
43 Xia H, Chen H, Cheng X, Yin M, Yao X, Ma J, Huang M, Chen G, Liu H. 2022. Zebrafish: an efficient vertebrate model for understanding role of gut microbiota. Mol Med 28 :161. doi:10.1186/s10020-022-00579-1 36564702
44 Phillips JB, Westerfield M. 2014. Zebrafish models in translational research: tipping the scales toward advancements in human health. Dis Model Mech 7 :739–743. doi:10.1242/dmm.015545 24973743
45 Adhish M, Manjubala I. 2023. Effectiveness of zebrafish models in understanding human diseases-a review of models. Heliyon 9 :e14557. doi:10.1016/j.heliyon.2023.e14557 36950605
46 Sharpton TJ, Stagaman K, Sieler MJ, Arnold HK, Davis EW. 2021. Phylogenetic integration reveals the zebrafish core microbiome and its sensitivity to environmental exposures. Toxics 9 :10. doi:10.3390/toxics9010010 33467528
47 Roeselers G, Mittge EK, Stephens WZ, Parichy DM, Cavanaugh CM, Guillemin K, Rawls JF. 2011. Evidence for a core gut microbiota in the zebrafish. ISME J 5 :1595–1608. doi:10.1038/ismej.2011.38 21472014
48 Gaulke CA, Beaver LM, Armour CR, Humphreys IR, Barton CL, Tanguay RL, Ho E, Sharpton TJ. 2020. An integrated gene catalog of the zebrafish gut microbiome reveals significant homology with mammalian microbiomes. Microbiology. doi:10.1101/2020.06.15.153924
49 Kent ML, Gaulke CA, Watral V, Sharpton TJ. 2018. Pseudocapillaria tomentosa in laboratory zebrafish Danio rerio: patterns of infection and dose response. Dis Aquat Organ 131 :121–131. doi:10.3354/dao03286 30460918
50 Nathan J, Kannan RR. 2020. Antiangiogenic molecules from marine actinomycetes and the importance of using zebrafish model in cancer research. Heliyon 6 :e05662. doi:10.1016/j.heliyon.2020.e05662 33319107
51 Yang L, Ho NY, Alshut R, Legradi J, Weiss C, Reischl M, Mikut R, Liebel U, Müller F, Strähle U. 2009. Zebrafish embryos as models for embryotoxic and teratological effects of chemicals. Reprod Toxicol 28 :245–253. doi:10.1016/j.reprotox.2009.04.013 19406227
52 Jia H-R, Zhu Y-X, Duan Q-Y, Chen Z, Wu F-G. 2019. Nanomaterials meet zebrafish: toxicity evaluation and drug delivery applications. J Control Release 311–312 :301–318. doi:10.1016/j.jconrel.2019.08.022
53 Crawford AD, Liekens S, Kamuhabwa AR, Maes J, Munck S, Busson R, Rozenski J, Esguerra CV, de Witte PAM. 2011. Zebrafish bioassay-guided natural product discovery: isolation of angiogenesis inhibitors from East African medicinal plants. PLOS ONE 6 :e14694. doi:10.1371/journal.pone.0014694 21379387
54 Pitchai A, Rajaretinam RK, Freeman JL. 2019. Zebrafish as an emerging model for bioassay-guided natural product drug discovery for neurological disorders. Medicines (Basel) 6 :61. doi:10.3390/medicines6020061 31151179
55 Delgadillo-Silva LF, Tsakmaki A, Akhtar N, Franklin ZJ, Konantz J, Bewick GA, Ninov N. 2019. Modelling pancreatic β-cell inflammation in zebrafish identifies the natural product wedelolactone for human islet protection. Dis Model Mech 12 :dmm036004. doi:10.1242/dmm.036004 30679186
56 Kent ML, Harper C, Wolf JC. 2012. Documented and potential research impacts of subclinical diseases in zebrafish. ILAR J 53 :126–134. doi:10.1093/ilar.53.2.126 23382344
57 Crofton HD. 1971. A quantitative approach to parasitism. Parasitology 62 :179–193. doi:10.1017/S0031182000071420
58 Shaw DJ, Dobson AP. 1995. Patterns of macroparasite abundance and aggregation in wildlife populations: a quantitative review. Parasitology 111 Suppl :S111–27. doi:10.1017/s0031182000075855 8632918
59 Rabajante JF. 2022. On spatiotemporal overdispersion and macroparasite accumulation in hosts leading to aggregation: a quantitative framework. Diseases 11 :4. doi:10.3390/diseases11010004 36648869
60 MacIntosh AJJ, Jacobs A, Garcia C, Shimizu K, Mouri K, Huffman MA, Hernandez AD. 2012. Monkeys in the middle: parasite transmission through the social network of a wild primate. PLOS ONE 7 :e51144. doi:10.1371/journal.pone.0051144 23227246
61 Ezenwa VO, Worsley-Tonks KEL. 2018. Social living simultaneously increases infection risk and decreases the cost of infection. Proc Biol Sci 285 :20182142. doi:10.1098/rspb.2018.2142 30487314
62 Otterstatter MC, Thomson JD. 2007. Contact networks and transmission of an intestinal pathogen in bumble bee (bombus impatiens) colonies. Oecologia 154 :411–421. doi:10.1007/s00442-007-0834-8 17713789
63 Oliveira RF. 2013. Mind the fish: zebrafish as a model in cognitive social neuroscience. Front Neural Circuits 7 :131. doi:10.3389/fncir.2013.00131 23964204
64 Jones LJ, Norton WHJ. 2015. Using zebrafish to uncover the genetic and neural basis of aggression, a frequent comorbid symptom of psychiatric disorders. Behav Brain Res 276 :171–180. doi:10.1016/j.bbr.2014.05.055 24954772
65 Green J, Collins C, Kyzar EJ, Pham M, Roth A, Gaikwad S, Cachat J, Stewart AM, Landsman S, Grieco F, Tegelenbosch R, Noldus L, Kalueff AV. 2012. Automated high-throughput neurophenotyping of zebrafish social behavior. J Neurosci Methods 210 :266–271. doi:10.1016/j.jneumeth.2012.07.017 22884772
66 Gaulke CA, Martins ML, Watral VG, Humphreys IR, Spagnoli ST, Kent ML, Sharpton TJ. 2019. A longitudinal assessment of host-microbe-parasite interactions resolves the zebrafish gut microbiome’s link to Pseudocapillaria tomentosa infection and pathology. Microbiome 7 :10. doi:10.1186/s40168-019-0622-9 30678738
67 Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI. 2009. A core gut microbiome in obese and lean twins. Nature 457 :480–484. doi:10.1038/nature07540 19043404
68 Leshem A, Segal E, Elinav E. 2020. The gut microbiome and individual-specific responses to diet. mSystems 5 :e00665-20. doi:10.1128/mSystems.00665-20
69 Zaiss MM, Harris NL. 2016. Interactions between the intestinal microbiome and helminth parasites. Parasite Immunol 38 :5–11. doi:10.1111/pim.12274 26345715
70 Leung JM, Graham AL, Knowles SCL. 2018. Parasite-microbiota interactions with the vertebrate gut: synthesis through an ecological lens. Front Microbiol 9 :843. doi:10.3389/fmicb.2018.00843 29867790
71 Pickard JM, Zeng MY, Caruso R, Núñez G. 2017. Gut microbiota: role in pathogen colonization, immune responses, and inflammatory disease. Immunol Rev 279 :70–89. doi:10.1111/imr.12567 28856738
72 Browne AJ, Chipeta MG, Haines-Woodhouse G, Kumaran EPA, Hamadani BHK, Zaraa S, Henry NJ, Deshpande A, Reiner RC, Day NPJ, Lopez AD, Dunachie S, Moore CE, Stergachis A, Hay SI, Dolecek C. 2021. Global antibiotic consumption and usage in humans, 2000-18: a spatial modelling study. Lancet Planet Health 5 :e893–e904. doi:10.1016/S2542-5196(21)00280-1 34774223
73 Pullan RL, Brooker SJ. 2012. The global limits and population at risk of soil-transmitted helminth infections in 2010. Parasit Vectors 5 :81. doi:10.1186/1756-3305-5-81 22537799
74 Wang L, Cao L, Chang Y, Fu Y, Wang Y, Zhang K, Zhang S, Zhang L. 2023. Microbiome-metabolomics analysis of the impacts of Cryptosporidium muris Infection in BALB/C mice. Microbiol Spectr 11 :e0217522. doi:10.1128/spectrum.02175-22 36533947
75 Yang W, Cong Y. 2021. Gut microbiota-derived metabolites in the regulation of host immune responses and immune-related inflammatory diseases. Cell Mol Immunol 18 :866–877. doi:10.1038/s41423-021-00661-4 33707689
76 da Silva KM, Iturrospe E, Bars C, Knapen D, Van Cruchten S, Covaci A, van Nuijs ALN. 2021. Mass spectrometry-based zebrafish toxicometabolomics: a review of analytical and data quality challenges. Metabolites 11 :635. doi:10.3390/metabo11090635 34564451
77 Medriano CA, Bae S. 2022. Acute exposure to microplastics induces metabolic disturbances and gut dysbiosis in adult zebrafish (Danio rerio). Ecotoxicol Environ Saf 245 :114125. doi:10.1016/j.ecoenv.2022.114125 36183426
78 Aguilar A, Mattos H, Carnicero B, Sanhueza N, Muñoz D, Teles M, Tort L, Boltaña S. 2022. Metabolomic profiling reveals changes in amino acid and energy metabolism pathways in liver, intestine and brain of zebrafish exposed to different thermal conditions. Front Mar Sci 9 . doi:10.3389/fmars.2022.835379
79 Sun B, Liu M, Tang L, Zhou X, Hu C, Chen L. 2023. Variability in fecal metabolome depending on age, PFBS pollutant, and fecal transplantation in zebrafish: a non-invasive diagnosis of health. J Environ Sci (China) 127 :530–540. doi:10.1016/j.jes.2022.06.019 36522083
80 Caboni P, Aissani N, Cabras T, Falqui A, Marotta R, Liori B, Ntalli N, Sarais G, Sasanelli N, Tocco G. 2013. Potent nematicidal activity of phthalaldehyde, salicylaldehyde, and cinnamic aldehyde against Meloidogyne incognita. J Agric Food Chem 61 :1794–1803. doi:10.1021/jf305164m 23379671
81 Danquah WB, Back MA, Grove IG, Haydock PPJ. 2011. In vitro nematicidal activity of a garlic extract and salicylaldehyde on the potato cyst nematode, Globodera pallida. Nematol 13 :869–885. doi:10.1163/138855411X560959
82 Bahiri G, Elliott I. 2013. Composition containing salicylaldehyde nematicides and garlic extract. ES2400745T3
83 Choi Y, Lee S, Kim S, Lee J, Ha J, Oh H, Lee Y, Kim Y, Yoon Y. 2020. Vitamin E (α-tocopherol) consumption influences gut microbiota composition. Int J Food Sci Nutr 71 :221–225. doi:10.1080/09637486.2019.1639637 31298050
84 Hansen HS, Diep TA. 2009. N-acylethanolamines, anandamide and food intake. Biochem Pharmacol 78 :553–560. doi:10.1016/j.bcp.2009.04.024 19413995
85 Tsuboi K, Uyama T, Okamoto Y, Ueda N. 2018. Endocannabinoids and related N-acylethanolamines: biological activities and metabolism. Inflamm Regen 38 :28. doi:10.1186/s41232-018-0086-5 30288203
86 Mennella I, Boudry G, Val-Laillet D. 2019. Ethanolamine produced from oleoylethanolamide degradation contributes to acetylcholine/dopamine balance modulating eating behavior. J Nutr 149 :362–365. doi:10.1093/jn/nxy281 30722047
87 Fornelos N, Franzosa EA, Bishai J, Annand JW, Oka A, Lloyd-Price J, Arthur TD, Garner A, Avila-Pacheco J, Haiser HJ, Tolonen AC, Porter JA, Clish CB, Sartor RB, Huttenhower C, Vlamakis H, Xavier RJ. 2020. Growth effects of N-acylethanolamines on gut bacteria reflect altered bacterial abundances in inflammatory bowel disease. Nat Microbiol 5 :486–497. doi:10.1038/s41564-019-0655-7 31959971
88 Hasan AU, Rahman A, Kobori H. 2019. Interactions between host PPARs and gut microbiota in health and disease. Int J Mol Sci 20 :387. doi:10.3390/ijms20020387 30658440
89 Di Paola M, Bonechi E, Provensi G, Costa A, Clarke G, Ballerini C, De Filippo C, Passani MB. 2018. Oleoylethanolamide treatment affects gut microbiota composition and the expression of intestinal cytokines in Peyer’s patches of mice. Sci Rep 8 :14881. doi:10.1038/s41598-018-32925-x 30291258
90 HashiniM.B, Donovan A. 2018. Host- and helminth-derived endocannabinoids that have effects on host immunity are generated during infection - PMC. Infect Immun. doi: 10.1128/IAI.00441-18
91 Janeš D, Kreft S. 2008. Salicylaldehyde is a characteristic aroma component of buckwheat groats. Food Chem 109 :293–298. doi:10.1016/j.foodchem.2007.12.032 26003350
92 Pelttari E, Karhumäki E, Langshaw J, Peräkylä H, Elo H. 2007. Antimicrobial properties of substituted salicylaldehydes and related compounds. Z Naturforsch C 62 :487–497. doi:10.1515/znc-2007-7-806 17913062
93 Salicylaldehyde pelomonas - protein - NCBI. 2024 Available from: https://www.ncbi.nlm.nih.gov/protein/?term=salicylaldehyde%20pelomonas
94 Jia B, Jia X, Hyun Kim K, Ji Pu Z, Kang M-S, Ok Jeon C. 2017. Evolutionary, computational, and biochemical studies of the salicylaldehyde dehydrogenases in the naphthalene degradation pathway. Sci Rep 7 :43489. doi:10.1038/srep43489 28233868
95 Stringfellow WT, Aitken MD. 1995. Competitive metabolism of naphthalene, methylnaphthalenes, and fluorene by phenanthrene-degrading pseudomonads. Appl Environ Microbiol 61 :357–362. doi:10.1128/aem.61.1.357-362.1995 7887615
96 Pham VT, Dold S, Rehman A, Bird JK, Steinert RE. 2021. Vitamins, the gut microbiome and gastrointestinal health in humans. Nutr Res 95 :35–53. doi:10.1016/j.nutres.2021.09.001 34798467
97 Liu KY, Nakatsu CH, Jones-Hall Y, Kozik A, Jiang Q. 2021. Vitamin E alpha- and gamma-tocopherol mitigate colitis, protect intestinal barrier function and modulate the gut microbiota in mice. Free Radic Biol Med 163 :180–189. doi:10.1016/j.freeradbiomed.2020.12.017 33352218
98 Smith A, Madden KB, Yeung KJA, Zhao A, Elfrey J, Finkelman F, Levander O, Shea-Donohue T, Urban JF Jr. 2005. Deficiencies in selenium and/or vitamin E lower the resistance of mice to Heligmosomoides polygyrus infections. J Nutr 135 :830–836. doi:10.1093/jn/135.4.830 15795443
99 Au Yeung KJ, Smith A, Zhao A, Madden KB, Elfrey J, Sullivan C, Levander O, Urban JF, Shea-Donohue T. 2005. Impact of vitamin E or selenium deficiency on nematode-induced alterations in murine intestinal function. Exp Parasitol 109 :201–208. doi:10.1016/j.exppara.2004.12.015 15755417
100 Sutherland IA, Brown AE, Leathwick DM. 2003. The effect of anthelmintic capsules on the egg output and larval viability of drug-resistant parasites. Vet Res Commun 27 :149–157. doi:10.1023/a:1022867322264 12718508
101 Yazwinski TA, Tucker CA, Wray E, Jones L, Reynolds J, Hornsby P, Powell J. 2013. Control trial and fecal egg count reduction test determinations of nematocidal efficacies of moxidectin and generic ivermectin in recently weaned, naturally infected calves. Vet Parasitol (Amst) 195 :95–101. doi:10.1016/j.vetpar.2012.12.061
102 Kent ML, Watral V, Villegas EN, Gaulke CA. 2019. Viability of Pseudocapillaria tomentosa eggs exposed to heat, ultraviolet light, chlorine, iodine, and desiccation. Zebrafish 16 :460–468. doi:10.1089/zeb.2019.1736 31216245
103 Xiang L, Moore BS. 2005. Biochemical characterization of a prokaryotic phenylalanine ammonia lyase. J Bacteriol 187 :4286–4289. doi:10.1128/JB.187.12.4286-4289.2005 15937191
104 Moffitt MC, Louie GV, Bowman ME, Pence J, Noel JP, Moore BS. 2007. Discovery of two cyanobacterial phenylalanine ammonia lyases: kinetic and structural characterization. Biochemistry 46 :1004–1012. doi:10.1021/bi061774g 17240984
105 Lovelock SL, Turner NJ. 2014. Bacterial anabaena variabilis phenylalanine ammonia lyase: a biocatalyst with broad substrate specificity. Bioorg Med Chem 22 :5555–5557. doi:10.1016/j.bmc.2014.06.035 25037641
106 Malinowski J, Krzymowska M, Godoń K, Hennig J, Podstolski A. 2007. A new catalytic activity from tobacco converting 2‐coumaric acid to salicylic aldehyde. Physiol Plant 129 :461–471. doi:10.1111/j.1399-3054.2006.00837.x
107 Sarkate A, Saini SS, Kumar P, Sharma AK, Sircar D. 2018. Salicylaldehyde synthase activity from venturia inaequalis elicitor-treated cell culture of apple. J Plant Physiol 221 :66–73. doi:10.1016/j.jplph.2017.12.002 29247889
108 Kurlyandchik I, Lauche R, Tiralongo E, Warne LN, Schloss J. 2022. Plasma and interstitial levels of endocannabinoids and N-acylethanolamines in patients with chronic widespread pain and fibromyalgia: a systematic review and meta-analysis. Pain Rep 7 :e1045. doi:10.1097/PR9.0000000000001045 36381652
109 Donia MS, Ravel J, Schmidt EW. 2008. A global assembly line for cyanobactins. Nat Chem Biol 4 :341–343. doi:10.1038/nchembio.84 18425112
110 Salvatella X, Caba JM, Albericio F, Giralt E. 2003. Solution structure of the antitumor candidate trunkamide A by 2D NMR and restrained simulated annealing methods. J Org Chem 68 :211–215. doi:10.1021/jo026464s 12530842
111 In vitro anthelmintic activity of Baliospermum montanum Muell. 2008. Arg Roots - Document - Gale Academic OneFile. Available from: https://go.gale.com/ps/i.do?p=AONE&u=googlescholar&id=GALE|A178741632&v=2.1&it=r&sid=AONE&asid=9e18e60e
112 Tandon V, Das B. 2018. Genistein: is the multifarious botanical a natural anthelmintic too? J Parasit Dis 42 :151–161. doi:10.1007/s12639-018-0984-0 29844617
113 Coldham NG, Darby C, Hows M, King LJ, Zhang A-Q, Sauer MJ. 2002. Comparative metabolism of genistin by human and rat gut microflora: detection and identification of the end-products of metabolism. Xenobiot 32 :45–62. doi:10.1080/00498250110085809
114 Kirkwood JS, Maier C, Stevens JF. 2013. Simultaneous, untargeted metabolic profiling of polar and nonpolar metabolites by LC-Q-TOF mass spectrometry. Curr Protoc Toxicol Chapter 4 :Unit4. doi:10.1002/0471140856.tx0439s56
115 Apprill A, McNally S, Parsons R, Weber L. 2015. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol 75 :129–137. doi:10.3354/ame01753
116 Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. 2016. DADA2: high-resolution sample inference from illumina amplicon data. Nat Methods 13 :581–583. doi:10.1038/nmeth.3869 27214047
117 Dixon P. 2003. VEGAN, a package of R functions for community ecology. J Veg Sci 14 :927–930. doi:10.1111/j.1654-1103.2003.tb02228.x
118 Helwig NE. 2019. Robust nonparametric tests of general linear model coefficients: acomparison of permutation methods and test statistics. Neuroimage 201 :116030. doi:10.1016/j.neuroimage.2019.116030 31330243
119 Tingley D, Yamamoto T, Hirose K, Keele L, Imai K. 2014. Mediation: R package for causal mediation analysis. J Stat Softw 59 :1–38. doi:10.18637/jss.v059.i05 26917999
120 Ggplot2 citation info. 2016. Available from: https://cran.r-project.org/web/packages/ggplot2/citation.html
