Abstract
Introduction
Intestinal helminth parasites trigger the host immune response through epithelial sensory tuft cells, but helminth-derived molecules that may activate tuft cells are poorly characterized.
Objectives
The study aimed to identify small molecules released in vitro by two nematode parasites, that infect rodents (Nippostrongylus brasiliensis) and ruminants (Haemonchus contortus), and to test candidate ligands in an in vivo model of tuft cell differentiation.
Methods
Small molecules were analyzed by hydrophilic interaction liquid chromatography (HILIC) of material released by adult parasites incubated in serum-free media, followed by mass spectrometry; selected molecules were administered to mice and tuft cell expansion enumerated after 5 days.
Results
A range of different conditions (culture media, timing, oxygenation) were tested, and comparisons made between the conditions, and between the two nematode species at selected points. Common products across the conditions and species included carboxylic acids (malate, succinate), medium chain fatty acids (such as decanoic and undecanoic acids), purines (guanine, xanthine and their derivatives), and phosphocholine compounds. We selected 19 of the prominent molecules for in vivo testing by oral administration, including succinate, a known activator of tuft cell differentiation. Malate elicited a low but significant level of tuft cell expansion, while undecanoic acids with or without a bromine substitution were also able to induce significant differentiation comparable to succinate. Other molecules including phosphorylcholine had no effect.
Conclusion
Multiple molecular species including decanoic and undecanoic acids released by helminths may contribute to activation of tuft cells in vivo.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Helminth parasites and the mammalian immune system have co-adapted through long evolutionary time, with elaborate parasite immune evasion strategies matched by a fine-tuned host immune system to detect and eliminate parasites (Douglas et al., 2021; Girgis et al., 2013; Maizels & Gause, 2023). Nonetheless, helminth infections present a major health burden across the tropics with an estimated 1.45 billion infections with soil-transmitted helminths alone (Pullan et al., 2014), while livestock infestations inflict enormous economic losses across the world (Charlier et al., 2020; Roeber et al., 2013).
The mammalian immune system has established multiple mechanisms in response to helminth parasites, broadly classed as type 2 immunity (Hammad et al., 2022; Lloyd & Snelgrove, 2018; Maizels & Gause, 2023). Responder cells at barrier surfaces release alarmins and cytokines that signal the infection to other immune cells, including type-2 innate lymphoid cells (ILC2s), alternatively activated macrophages (M2s) and eosinophils (Grencis, 2015; McDaniel et al., 2023; Zaiss et al., 2024). Once adaptive immunity has been induced, antigen-specific Th2 and B cells mobilise to reduce parasite burdens or expel worms completely (Harris & Loke, 2017; Maizels & Gause, 2023).
Notably, expulsion of gastrointestinal (GI) parasites is dependent on a fully-responsive epithelial population. Thus, goblet cells, instructed by IL-13-producing ILC2s and Th2 cells, execute anti-helminth effects by secreting mucus and resistin-like molecule β (RELM-β) (Finkelman et al., 2004; Herbert et al., 2009). More recently, another epithelial cell type, the tuft cell, has been recognized as instrumental in activating the immune response against helminths and protists, with the intestinal nematode Nippostrongylus brasiliensis being a particularly potent stimulus for tuft cell expansion (Gerbe et al., 2016; Howitt et al., 2016; Kotas et al., 2023; von Moltke et al., 2016). Tuft cells are also endowed with critical effector functions (Billipp et al., 2020; Campillo Poveda et al., 2023), including the production of acetylcholine which acts in the intestinal tract to increase fluid egress (Billipp et al., 2024) and directly compromise helminth neurological function (Ndjim et al., 2024).
Tuft cells, characterised by an apical tuft-like structure of microvilli reaching into the intestinal lumen, are chemosensory cells expressing a range of G-protein coupled receptors (GPCRs) including small molecule taste sensors, signalling through the cation channel Trpm5 (Billipp et al., 2020; Campillo Poveda et al., 2023; O’Leary et al., 2019). In the intestinal tract, tuft cells are crucial components of the anti-helminth immune response as tuft cell-deficient Pou2f3–/– mice are unable to clear an acute Nippostrongylus brasiliensis roundworm infection (Gerbe et al., 2016). In response to infection, tuft cells release IL-25 that activates IL-13 production by ILC2s, creating a positive feed-forward loop that drives goblet and tuft cell hyperplasia for parasite expulsion (Gerbe et al., 2016; von Moltke et al., 2016). Similar expansion of tuft cells is found in ruminants infected with GI nematodes (Hildersley et al., 2021), demonstrating the evolutionary conservation of this pathway.
Infection with the protist Tritrichomonas muris also induces a tuft cell response, through Trpm5 activation (Howitt et al., 2016) involving the succinate GPCR (Sucnr1/GPR91) (Nadjsombati et al., 2018; Schneider et al., 2018). Succinate administration via drinking water directly drives Sucnr1-dependent tuft cell and subsequent ILC2 activation in distal parts of the small intestine (Lei et al., 2018; Nadjsombati et al., 2018). Although both the protist T. muris and the helminth N. brasiliensis have been shown to produce succinate (Nadjsombati et al., 2018), Sucnr1-deficient mice were incapable of clearing the protist infection, but were able to expel N. brasiliensis worms even in the absence of Sucnr1 (Lei et al., 2018; Nadjsombati et al., 2018), suggesting redundancy and involvement of additional receptors.
Tuft cells express bitter taste receptors of the Tas2r family, some of which are upregulated in helminth infection, and stimulation with bitter compounds induces IL-25 release (Luo et al., 2019), although it is not known whether Tas2r receptors are required for parasite immunity. Tuft cells also express the sweet/umami taste receptor Tas1r3, although this receptor is not required for tuft cell responses to another intestinal nematode, Heligmosomoides polygyrus (Howitt et al., 2020). The receptor(s) activating tuft cells in helminth infection therefore remain unclear and any helminth-derived metabolites driving tuft cell activation have yet to be identified.
Some indications of possible ligands arise from analyses of other micro-organisms, and from studies on nematode biological interactions. In the case of Shigella bacterial infection, the metabolite N-undecanoylglycine was found to activate tuft cells, and to do so via the vomeronasal GPCR Vmn2r26 that is associated with olfactory sensing (Xiong et al., 2022). In the helminth setting, a widely represented small molecule family is the ascarosides, small dideoxy sugars with short chain fatty acid moieties, which function as pheromone signals between individual nematodes (Choe et al., 2012). Among these, ascarosides released from N. brasiliensis adult worms mediate immunosuppressive effects in mouse models of inflammation (Manosalva et al., 2015; Shinoda et al., 2021).
For tuft cells to function as sentinels of infection, they must sense products released by incoming parasites; proteomics has identified excretory/secretory proteins from helminths, released during in vitro culture (Harnett, 2014; Hewitson et al., 2009), but only recently has attention turned to the small molecule repertoire exported from helminth parasites (Chen et al., 2021; Wangchuk et al., 2019a, 2019b, 2019c; Whitman et al., 2021; Yeshi et al., 2020), recently reviewed by Wangchuk (Wangchuk et al., 2023). Analyses of small molecules released both by adult worms (Wangchuk et al., 2019b) and infective L3 larvae of N. brasiliensis (Yeshi et al., 2020), as well as the intestinal contents of infected mice (Chen et al., 2021), identified 45–55 polar metabolites and several hundred lipid species.
In searching for possible tuft cell ligands, we surmised that the host system evolved to sense the presence of molecules common to multiple helminth species, possibly resulting from an obligatory metabolic pathway that is physiologically essential to the parasites. Such a product may be, like succinate, a familiar molecule which activates at high concentrations, or a novel component foreign to the mammalian system that can be more sensitively detected. To pursue this inquiry, we conducted parallel metabolomic analyses on two major GI nematode parasites, the laboratory model mentioned above, N. brasiliensis, and the ruminant “barber’s pole” roundworm Haemonchus contortus. While the metabolic biochemistry of these species has been well described (Harder, 2016), this has been at the whole-organism level rather than, as we report here, a focus on released molecules which may interact with the host.
2 Methods
2.1 Animals and parasites
N. brasiliensis was maintained in rats as previously described (Camberis et al., 2003). H. contortus (MHco3(ISE) strain) was recovered from infected sheep as previously described (Laing et al., 2013). Sheep were aged four to nine months and were reared and housed indoors under parasite-free conditions prior to H. contortus infection. Parasite maintenance and animal experiments in which mice received metabolites in drinking water were conducted under a UK Home Office licence and approved by the University of Glasgow Animal Welfare and Ethical Review Board (AWERB). All experimental procedures in sheep were examined and approved by the Moredun Research Institute AWERB and were conducted under approved UK Home Office licence in accordance with the 1986 Animal (Scientific Procedures) Act, UK.
2.2 Parasite cultures
N. brasiliensis adult worms were incubated at a concentration of 90–200 worms/mL in Earle’s Balanced Salt Solution (EBSS, GIBCO 24010-043), Hanks’ BSS (HBSS, Sigma H8264), or RPMI1640 (Sigma R7509), each with gentamycin (1/1000 dilution; Sigma G1397, 50 mg/mL); each culture condition was set up as 1 mL volumes in 3 separate wells of a 48-well plate, with additional control wells containing media without worms. Worms were then incubated at 37C° for 4 or 24 h. Under normoxic conditions, cultures were in 5% CO2 in air; under hypoxic conditions a mixture of 1% O2 5% CO2 94% N2 was used for 4 or 24 h inside an airtight Billups-Rothenberg hypoxic chamber.
H. contortus adult worms were washed 6 times in warm HBSS followed by 3 washes in HBSS + penicillin/streptomycin (1/1000 dilution; Invitrogen 15140-122, 10,000 units/mL). Worms were then washed in the respective media (HBSS, EBSS or RPMI1640) and male and female worms cultured separately at a concentration of 10 worms/mL as described for N. brasiliensis. For both parasites, media were collected at 4 and 24 h and immediately quenched in a dry ice/ethanol bath before centrifuging in a chilled microfuge at 1900 g for 10 min.
2.3 Metabolomic analyses
Analyses were performed on 10 µL of culture supernatant, to which was added 200 μL chloroform:methanol:water (CMW, 1:3:1) extraction solvent at 4 °C. Following mixing, samples were centrifuged at 13,000 g, 4 °C for 5 min, and the supernatants (~ 195 μL) transferred to fresh 2 mL screw-top Eppendorf tubes, prior to displacing air with argon gas. Samples were stored at -80 °C until analysis.
Samples were subjected to hydrophilic interaction liquid chromatography (HILIC) on a Dionex UltiMate 3000 RSLC system (Thermo Fisher Scientific) using a ZIC-pHILIC column (150 mm x 4.6 mm, 5 μm column, Merck Sequant) maintained at 40 °C and eluted with a linear gradient of 20 mM ammonium carbonate (20% to 95%) versus acetonitrile (80% to 5%) at a flow rate of 0.3 mL/min for 20 min.
Fractions were then injected into an Orbitrap Fusion (Thermo Fisher Scientific), operated in polarity switching mode with the following settings: Resolution 120,000; AGC 2e5; m/z range 70–1000; Sheath gas 40; Auxiliary gas 5; Sweep gas 1; Probe temperature 150 °C; Capillary temperature 325 °C. For positive mode ionization: source voltage + 4.3 kV. For negative mode ionization: source voltage − 3.2 kV. S-Lens RF Level 60.00%.
Annotated data files were analyzed using the online service of Metabolanalyst 5.0 (https://www.metaboanalyst.ca/), R studio and Prism for comprehensive data normalization, statistical and functional analyses.
2.4 In vivo administration of metabolites
Metabolites administered were: L-alanine (A7627, Sigma-Aldrich); 7-aminomethyl-7-carbaguanine (PreQ1 dihydrochloride, SML0807, Sigma-Aldrich); L-carnitine hydrochloride (C0283, Sigma-Aldrich), creatine (C0780, Sigma-Aldrich); NG-NG-dimethyl-L-arginine (D0390, Sigma-Aldrich); 1,2-dioleoyl-sn-glycero-3-phosphocholine (P6354, Sigma-Aldrich); hypoxanthine (H9377, Sigma-Aldrich), 7- and 9-methylguanine (67073, 67074, Sigma-Aldrich); malic acid (M1000, Sigma-Aldrich); phosphocholine chloride (P0378, Sigma-Aldrich); pyruvic acid (107360, Sigma-Aldrich); uridine (U6381, Sigma-Aldrich), xanthine (X7375, Sigma-Aldrich), 3-methylguanine (sc-483803, Santa Cruz Biotechnology), undecanoic acid (171476, Sigma-Aldrich), 10-bromodecanoic acid (541397, Sigma-Aldrich) and 11-bromoundecanoic acid (B82804, Sigma-Aldrich). Details of concentration, delivery route and timing are detailed in the relevant Figure Legends.
2.5 Immunohistochemistry
C57BL/6 female mouse small intestines were embedded in paraffin using the Swiss-rolling technique (Bialkowska et al., 2016), which enables multiple villi to be investigated in a single cut section. Transverse sections were made using a microtome through the gut rolls at a thickness of 5 µm before mounting on glass slides. Sections were deparaffinized by immersing slides in xylene, then hydrated through 100%, 90%, and 70% ethanol successively. Heat-induced epitope retrieval was performed in 10 mM citrate buffer pH 6 (Thermo Fisher Scientific), and then sections were blocked using 2.5% normal horse serum blocking solution (Vector Laboratories) for 1 h at room temperature. Slides were then incubated overnight at 4 °C with rabbit anti-mouse Dclk1 antibody (Abcam) diluted 1/1000 in 2.5% normal horse serum blocking solution. Polyclonal rabbit IgG (Abcam) was used as an isotype control. After washing, sections were incubated with swine anti-rabbit-IgG/FITC (Agilent Dako), washed, and mounted using Vectashield Vibrance antifade mounting medium with DAPI (Vector Laboratories). Slides were imaged using a Leica DMi8 inverted microscope and Leica Application Suite (LAS) X software (Leica Microsystems). The resulting image files were analyzed using ImageJ/Fiji. For each sample, 5 areas that included 5 villus-crypt units were randomly selected, and the number of stained tuft cells counted; from this the number of tuft cells per villus-crypt axis for each mouse was recorded.
2.6 Statistics
Functional annotation of distinguishing metabolites was performed using GraphPad Software (Boston, Massachusetts USA) and RStudio Version 2023. Statistical analysis generated q-values, which provide a more accurate estimate than p values of the false discovery rate (FDR) when dealing with multiple comparisons. Tuft cell numbers from the in vivo experiments were analyzed using a one-way ANOVA test with multiple comparisons with normality and lognormality tests of tuft cell counts.
3 Results
Analyses of small molecules was undertaken by hydrophilic interaction liquid chromatography (HILIC) of material released by adult parasites incubated in serum-free media, followed by mass spectrometry.. A range of different conditions (culture media, timing, oxygenation) were tested, and comparisons were made between the conditions, and between the two nematode species, as described below. In addition, the major metabolites were administered in vivo to experimentally determine effects on tuft cell induction. An overview of the study is shown in Fig. 1.
Methodology of small molecule analysis. Intestinal nematodes (H. contortus and N. brasiliensis adult worms) were cultured under various conditions for 4–24 h. Conditioned media were collected; small molecules extracted into chloroform/methanol and subject to mass spectrometry; identified peaks were compared and statistically verified; individual molecules of interest were administered to mice in drinking water to test for induction of tuft cells, enumerated by histological staining for the DCLK1 marker
3.1 Culture media and metabolite profiles for adult N. brasiliensis
We first explored optimal culture conditions for small molecule identification. In collecting excretory-secretory proteins from helminths, most protocols have used nutrient-enriched media that allows survival for days or even weeks in vitro (Harnett, 2014; Hewitson et al., 2009). For the purposes of collecting small molecules released over shorter time-windows, simpler conditions were selected. We compared incubation of adult worms at 37 °C in EBSS, HBSS or RPMI1640 tissue culture media. While BSS media are simply isotonic buffered salt solutions, RPMI1640 contains all 20 amino acids and 11 vitamins required for longer-term cell survival in tissue culture.
For each media type, three replicate parasite cultures were incubated and supernatants recovered after 4 h. Parasites cultured in either EBSS or HBSS released a greater number of metabolites identified above the significance threshold than parasites cultured in RPMI1640. This was in part due to the richer set of constituents in RPMI1640, which resulted in a high baseline for those metabolites. For example, the inclusion of choline chloride in RPMI1640 evidently confounded identification of this and related molecules. The greatest number of identifications were made in HBSS (Fig. 2 A, B), most of which were also found in EBSS, and many in RPMI1640-derived samples (Table 1).
Small molecules released by N. brasiliensis adult worms. In each case, 200 adults were incubated in 1 mL media; after 4 h, supernatants were recovered and extracted in chloroform:methanol:water for MS analysis. A Graphical plot of metabolites released in HBSS at 4 h showing peak intensities (with medium alone subtracted). Amino acids Met, Phe and Val are annotated, and together with all others (Arg, Glu, His, Leu, Lys, Pro, Trp, Tyr) coloured yellow. Named molecules coloured green. G-PC, sn-Gycero-3-Phosphocholine; MEP, 2-C-Methyl-D-Erythritol 4-Phosphate. Stearic acid is octadecanoic acid. C As A but expanded x axis from 0–2.5 × 106, and coloured as in A ACG, 7-aminomethyl-7-carbaguanine,. ADMA asymmetric dimethylarginine. Stearic acid is octadecanoic acid
A number of amino acids were amongst the most abundant molecules, most notably phenylalanine and methionine (indicated with yellow symbols in Fig. 2 A, B). Other than amino acids, a number of purine precursors and derivatives, and amino acid metabolites, were noted as previously reported for N. brasiliensis (Chen et al., 2021). Among the purine derivatives was methylguanine, that may be related to the 5′ cap of mRNA that is often trans-spliced with a trimethyl (2′2′7′-trimethylated) guanine in nematodes (Thomas et al., 1988). A few citric acid and tricarboxylic acid cycle metabolites were found, most conspicuously malate, and succinate. Both a short-chain (valeric acid) and long-chain (stearic acid) fatty acids were identified, the former having been reported as interfering with histone deacetylation required for CD4+ T cell activation (Luu et al., 2019). We also identified high levels of choline phosphate (phosphocholine) which is widely expressed across a spectrum of nematode parasites (Harnett et al., 1993; Maizels et al., 1987; Péry et al., 1979).
Identifying metabolites in parasite cultures depended on comparison with a standardised set of controls for mammalian and bacterial metabolism, therefore we were not able to conclusively identify ascarosides that are known to be released by nematodes (Choe et al., 2012; Shinoda et al., 2021). However, a prominent peak with mass 276 Daltons was noted in both EBSS- and HBSS-cultured media that is likely to represent ascr#1 (Supplementary Fig. 1).
3.2 Metabolite profiles of Haemonchus contortus culture media
We next analysed another intestinal nematode species, H. contortus that colonises the abomasum (true stomach) of ruminant livestock, and is, like N. brasiliensis, a member of the Trichostrongyloid superfamily. Adult worms were incubated in EBSS, HBSS or RPMI1640 under similar conditions, for 4 h. As observed for N. brasiliensis, fewer molecules surpassed the significance threshold in RPMI1640. A plot of the individual metabolites identified in HBSS is presented in Fig. 3 A.
Metabolites from H. contortus adult worms and comparison with those from N. brasiliensis. Plots of molecules identified from H. contortus incubated in HBSS. Amino acids are coloured yellow (only Ala is annotated), selected identified metabolites coloured green. NB, palmitic acid is hexadecanoic acid, stearic acid is octadecanoic acid. A. Plots of molecules identified from H. contortus incubated in HBSS. Amino acids are coloured yellow (only Ala is annotated), selected identified metabolites coloured green. NB, palmitic acid is hexadecanoic acid, stearic acid is octadecanoic acid. B Comparative intensity of 15 major metabolites released by the two species. NB Caproic acid is hexanoic acid. C Volcano plot comparison of molecules released by the two species with FDR at 0.5
Many of the identified molecules from H. contortus, including choline compounds, malate and succinate, mirrored those found to be released from N. brasiliensis as described in more detail below. However, longer chain fatty acids (palmitic and stearic) were more prominent albeit at a low level of statistical significance in supernatants of H. contortus. We also noted the presence of bromo-undecanoic acid in cultures of this parasite.
3.3 Comparative analysis of two helminth species
Metabolomic data generated from analyses of conditioned media from cultures with N. brasiliensis and H. contortus were combined to identify products released by both species. These analyses identified many common metabolites released into HBSS medium (Fig. 3 B) with similar profiles of small molecules in EBSS and RPMI1640 cultures (Suppl. Figure 3 A, B) Quantitatively, most common products were of similar intensity in the two species, however, some exceptions were noted. While malate and choline compounds (indicated in yellow) were present at high levels for both species, some molecules such as creatine were more prominent in N. brasiliensis conditioned media. A volcano plot comparing metabolomic data from the two species, showing the ratios for each molecule, reveals some key findings such as the presence of succinate and the relative abundance of medium chain undecanoic acids in the supernatants from H. contortus parasites (Fig. 3 C).
3.4 Time-dependent release of metabolites
We then compared metabolomic data from N. brasiliensis conditioned media collected after 4 h and 24 h incubation; products that showed a significant increase over this time period are less likely to have been carried over by the parasites from their intestinal environment. In this setting we employed both HBSS and RPMI160, as the latter medium is more conducive to maintaining parasite viability. As shown in Fig. 4 A (for HBSS) and B (for RPMI1640), most of the key molecules identified at 4 h were significantly increased after 24 h; these included malate, choline phosphate and guanine derivatives; in the case of RPMI1640, N. brasiliensis parasites also elaborated asymmetric dimethylarginine (ADMA), a reported inhibitor of nitric oxide synthase.
Analysis of small molecules released by N. brasiliensis adult worms between 4 and 24 h, and under different oxygenation conditions. In each case, 90 adults were incubated in 1 mL media; after 4 h, supernatants were replaced with fresh media, and conditioned media collected at 24 h; both 4 h and 24 h media were analysed and compared. A Fold-change increases in peak intensity from 4 to 24 h in HBSS. B As A, in RPMI1640. C Comparison of products released under normoxia or hypoxia (1% O2) in HBSS
3.5 Effects of hypoxia on the release of metabolites
As intestinal parasites are typically found in low-oxygen environments, we then asked whether the profile of small molecule release was affected by manipulation of oxygen levels, with parallel cultures set up under normoxia or hypoxia (1% O2). The principal effect of hypoxia was reduced production of numerous metabolites, with the exception of succinate and propanoate, indicating that while N. brasiliensis responds to differential oxygenation conditions, it shows a more active general release of metabolites in normoxia than hypoxia (Fig. 4 C).
3.6 In vivo effects on tuft cell induction
The major small molecule metabolites that had been identified in this study were then tested by administration in drinking water (Fig. 5 A), using protocols that had previously been shown to induce intestinal tuft cell differentiation in mice (Drurey et al., 2022; Nadjsombati et al., 2018).
Effect of small molecules on tuft cell differentiation in vivo. A Protocol for administration in drinking water. B Tuft cell counts, data pooled from two biological replicate experiments (n = 5 mice per experiment, each mouse count based on 5 randomly selected sections). Tuft cell numbers were analysed using one-way ANOVA test with multiple comparisons and Normality and Lognormality Tests for data distribution. C Staining for tuft cells in sections from mice treated with succinate (100 mM), 7-amino-methyl-7-carbaguanine (ACG, 10 mM), choline phosphate (20 mM), malate (100 mM) and pyruvate (100 mM). Sections stained at day 7 for DCLK1 expression. Yellow = DCLK1 and blue = DAPI. Scale bars represent 100 µm in all panels. *, p < 0.05, **** p < 0.0001
The concentrations chosen for each of the metabolites were based on their known safe levels regarding toxicity in rats, as reported on www.drugbank.com, and their solubility in water at 20 °C. Of 14 metabolites listed in Table 1 (3 of which were isomers of methylguanine) and tested in vivo, only succinate was able to drive marked expansion of tuft cell numbers (p < 0.0001), although malate did induce a small but significant increase in DCLK-1+ tuft cell numbers (p < 0.05, Fig. 5B), together with generalized disruption of the epithelial layer (Fig. 5 C). In addition, we tested uridine and pyruvate as less prominent but ubiquitous metabolites, and similarly found no induction of tuft cell expansion.
Furthermore, the medium-chain (10- and 11-carbon) fatty acids decanoic and undecanoic acid were also present in the identified molecules, in some cases brominated, although the type of substitution cannot be determined from mass values. As these compounds are not soluble in water, decanoic acid and two brominated variants conjugated to C-10 or C11, were administered by oral gavage as previously described by Shoji et al. (Shoji et al., 2022), applying compounds daily for 5 days (Fig. 6 A). This shorter treatment, and its mode of delivery, is clearly not as effective as drinking water administration, as shown by the lower numbers of tuft cells induced with succinate (Fig. 6 B). However, both 10- and 11-carbon compounds stimulated tuft cell expansion when given to mice under these conditions, with undecanoic acid activating the same response as succinate (Fig. 6 B, C), implicating medium-chain fatty acids in tuft cell hyperplasia.
Effect of medium chain fatty acids on tuft cell differentiation in vivo. A Protocol for administration by oral gavage (OG). B Staining for tuft cells in sections from mice treated with indicated molecules, data pooled from two biological replicate experiments (n = 5 mice per experiment). Tuft cell numbers were analysed using one-way ANOVA test with multiple comparisons and Normality and Lognormality Tests for data distribution. C Sections stained for DCLK1 in from mice treated for 5 days with succinate, undecanoic acid and negative control. **** p < 0.0001
4 Discussion
The molecular “arms race” between host and parasite has primarily been conceived at the level of macromolecules, in particular protein components (Maizels et al., 2018). However, small molecules from parasites are increasing in prominence through systematic analyses on metabolite production (Wangchuk et al., 2023), and identification of key pheromones such as ascarosides (Choe et al., 2012; Whitman et al., 2021), and small peptides (Chen et al., 2022) in parasite development. Such small molecules may also be detected by the host immune system to flag the presence of parasites. In this setting, host epithelial tuft cells can act as sentinels of helminth infection by detecting molecules released by invading parasites in the gut (O’Leary et al., 2019; Schneider et al., 2018). As tuft cells express sensory receptors associated with taste perception, and signal through a conserved GPCR-associated pathway, we and others have searched for candidate small molecules that may act as a helminth signature at the initiation of infection.
In this study we took a comparative approach to identify small metabolites common to two intestinal nematode species, N. brasiliensis which is a widely used rodent model that stimulates an untrammeled type 2 response, and H. contortus that is a parasite of sheep and goats globally. Although these parasite species employ different entry strategies, with larvae penetrating through cutaneous and gastric surfaces respectively, in both cases adult worms occupy gastrointestinal niches with similar mucosal immune challenges.
To chart and compare parasite products, we used HILIC for high-sensitivity analysis of small metabolites. Our data compare with studies using alternative analytical platforms, such as gas chromatography-mass spectrometry (GC–MS) or liquid chromatography-MS (LC–MS). For example, all three approaches identified stearic acid as one of the major fatty acids released by parasitic nematodes including hookworms, Trichuris muris and N. brasiliensis (Wangchuk et al., 2023). Similar to previous studies (Wangchuk et al., 2019b, 2023) common and species-specific molecules were identified.
From our comparative data, we selected 19 molecular species to test in a mouse model of tuft cell induction, administered by oral gavage as previously conducted with succinate that effectively drives expansion of this target cell population (Drurey et al., 2022; Schneider et al., 2018), through the succinate receptor GPR91. While malate showed a low level of induction, we found 10- and 11-carbon fatty acids to be equal to succinate; a previous report had identified bacterial N-undecanoylglycine to activate tuft cells, and to do so via the olfactory Vmn2r26 receptor expressed by these cells (Xiong et al., 2022). Medium chain fatty acids (9–14 carbons) bind GPR84, but as far as is known, this receptor is restricted to myeloid lineages such as macrophages and neutrophils (Wang et al., 2006). Future studies will aim to test the role of these receptors in helminth-induced tuft cell hyperplasia, and in guarding the host against helminth invasion. Notably, succinate and decanoic (capric) acid were also reported in secreted products of other GI nematodes, including hookworms and T. muris (Wangchuk et al., 2023).
Our HILIC and analysis technique focused on metabolites from conventional pathways, and hence would not have detected substances unique to nematodes. We conducted a directed search for ascarosides, small pheromone molecules associated with nematode communication (Choe et al., 2012; Whitman et al., 2021) that also show interactions with immune responses (Shinoda et al., 2021); we noted ascr#1 from N. brasiliensis but not H. contortus, perhaps reflecting the requirement for species specificity within the pheromone compartment. Future work on this class of small molecules is likely to identify further key host interactions. Within the standard metabolites, only malate was found to show a low level of activation. However, malate may be converted to succinate via fumarate, which could account for stimulating a low level of tuft cell activation by this compound.
The changes in amino acid levels, both consumption and production, are likely to have physiological effects on host cells. For example, group 2 innate lymphoid cells (ILC2s) are highly dependent on large neutral amino acids such as phenylalanine (Hodge et al., 2023). Release of this amino acid, most notably by N. brasiliensis, could potentially enhance ILC2 responses to infection. The differential effect of hypoxic conditions included some interesting results; it is known from microbial and C. elegans nematode metabolism that an anaerobic environment promotes succinate production (Del Borrello et al., 2019).
Our approach of oral delivery does have some limitations, in that we could not evaluate molecules with poor water solubility, and had little control over whether substances reached the small intestine in sufficient concentrations, or were converted or absorbed earlier in their journey. Targeted studies to determine the prevailing concentrations of, for example, medium chain fatty acids in the small intestinal lumen would be informative in this regard. Although we based the concentrations of administered molecules to be below any potential toxic levels, they are nevertheless given at relatively high molar concentrations, and nonspecific perturbations of the epithelium may result; in this setting, damage repair pathways may include IL-13 production, giving rise to tuft cell expansion in the absence of any sentinel receptor activation. More broadly, effects from oral administration of test substances does not necessarily show a direct mode of action, for example if molecules are metabolised or converted in vivo, or if they act on microbial populations in the intestinal tract, resulting in gain or loss of microbiome mediators. Medium chain fatty acids such as decanoic (capric) acid have indeed been reported to have anti-bacterial activity that could impact on tuft cell induction. Future studies in organoid culture systems will minimize a number of these concerns and facilitate demonstration of direct activities.
In conclusion, we have identified a second class of small molecules, the 10/11-carbon fatty acids, that can induce tuft cell differentiation in vivo, although their direct action through a specific receptor has yet to be established. Importantly, these results indicate that multiple ligands may have the capacity to act on intestinal epithelial tissues to promote development of secretory lineages including tuft cells; these findings raise the possibility that redundant mechanisms have evolved to detect the presence of helminths, and that additive or combinatorial signaling may be employed to maximise sensitivity and specificity for helminth infection.
Data availability
Data is provided within the manuscript or supplementary information files.
References
Bialkowska, A. B., Ghaleb, A. M., Nandan, M. O., & Yang, V. W. (2016). Improved Swiss-rolling technique for intestinal tissue preparation for immunohistochemical and immunofluorescent analyses. Journal of Visualized Experiments. https://doi.org/10.3791/54161
Billipp, T. E., Fung, C., Webeck, L. M., Sargent, D. B., Gologorsky, M. B., Chen, Z., McDaniel, M. M., Kasal, D. N., McGinty, J. W., Barrow, K. A., et al. (2024). Tuft cell-derived acetylcholine promotes epithelial chloride secretion and intestinal helminth clearance. Immunity,57(1243–1259), e1248.
Billipp, T. E., Nadjsombati, M. S., & von Moltke, J. (2020). Tuning tuft cells: New ligands and effector functions reveal tissue-specific function. Current Opinion in Immunology,68, 98–106.
Camberis, M., Le Gros, G., and Urban, J., Jr. (2003). Animal model of Nippostrongylus brasiliensis and Heligmosomoides polygyrus. In: R. Coico (ed.) Current Protocols in Immunology. Wiley. pp. 19.12.11–19.12.27.
Campillo Poveda, M., Britton, C., Devaney, E., McNeilly, T. N., Gerbe, F., Jay, P., & Maizels, R. M. (2023). Tuft cells: detectors, amplifiers, effectors and targets in parasite infection. Cells,12, 2477.
Charlier, J., Rinaldi, L., Musella, V., Ploeger, H. W., Chartier, C., Vineer, H. R., Hinney, B., von Samson-Himmelstjerna, G., Bacescu, B., Mickiewicz, M., et al. (2020). Initial assessment of the economic burden of major parasitic helminth infections to the ruminant livestock industry in Europe. Preventive Veterinary Medicine,182, 105103.
Chen, R., Wang, J., Gradinaru, I., Vu, H. S., Geboers, S., Naidoo, J., Ready, J. M., Williams, N. S., DeBerardinis, R. J., Ross, E. M., & Collins, J. J., 3rd. (2022). A male-derived nonribosomal peptide pheromone controls female schistosome development. Cell,185(1506–1520), e1517.
Chen, Y., Zhang, M., Ding, X., Yang, Y., Chen, Y., Zhang, Q., Fan, Y., Dai, Y., & Wang, J. (2021). Mining anti-inflammation molecules from Nippostrongylus brasiliensis-derived products through the metabolomics approach. Frontiers in Cellular and Infection Microbiology,11, 781132.
Choe, A., von Reuss, S. H., Kogan, D., Gasser, R. B., Platzer, E. G., Schroeder, F. C., & Sternberg, P. W. (2012). Ascaroside signaling is widely conserved among nematodes. Current Biology,22, 772–780.
Del Borrello, S., Lautens, M., Dolan, K., Tan, J. H., Davie, T., Schertzberg, M. R., Spensley, M. A., Caudy, A. A., & Fraser, A. G. (2019). Rhodoquinone biosynthesis in C. elegans requires precursors generated by the kynurenine pathway. eLife. https://doi.org/10.7554/eLife.48165
Douglas, B., Oyesola, O., Cooper, M. M., Posey, A., Tait Wojno, E., Giacomin, P. R., & Herbert, D. R. (2021). Immune system investigation using parasitic helminths. Annual Review of Immunology,39, 639–665.
Drurey, C., Lindholm, H. T., Coakley, G., Poveda, M. C., Loser, S., Doolan, R., Gerbe, F., Jay, P., Harris, N., Oudhoff, M. J., & Maizels, R. M. (2022). Intestinal epithelial tuft cell induction is negated by a murine helminth and its secreted products. Journal of Experimental Medicine,219, e20211140.
Finkelman, F. D., Shea-Donohue, T., Morris, S. C., Gildea, L., Strait, R., Madden, K. B., Schopf, L., & Urban, J. F., Jr. (2004). Interleukin-4- and interleukin-13-mediated host protection against intestinal nematode parasites. Immunological Reviews,201, 139–155.
Gerbe, F., Sidot, E., Smyth, D. J., Ohmoto, M., Matsumoto, I., Dardalhon, V., Cesses, P., Garnier, L., Pouzolles, M., Brulin, B., et al. (2016). Intestinal epithelial tuft cells initiate type 2 mucosal immunity to helminth parasites. Nature,529, 226–230.
Girgis, N. M., Gundra, U. M., & Loke, P. (2013). Immune regulation during helminth infections. PLoS Pathogens,9, e1003250.
Grencis, R. K. (2015). Immunity to helminths: Resistance, regulation, and susceptibility to gastrointestinal nematodes. Annual Review of Immunology,33, 201–225.
Hammad, H., Debeuf, N., Aegerter, H., Brown, A. S., & Lambrecht, B. N. (2022). Emerging paradigms in Type 2 immunity. Annual Review of Immunology,40, 443–467.
Harder, A. (2016). The biochemistry of Haemonchus contortus and other parasitic nematodes. Advances in Parasitology,93, 69–94.
Harnett, W. (2014). Secretory products of helminth parasites as immunomodulators. Molecular and Biochemical Parasitology,195, 130–136.
Harnett, W., Houston, K. M., Amess, R., & Worms, M. J. (1993). Acanthocheilonema viteae: Phosphorylcholine is attached to the major excretory-secretory product via an N-linked glycan. Experimental Parasitology,77, 498–502.
Harris, N. L., & Loke, P. (2017). Recent advances in Type-2-cell-mediated immunity: Insights from helminth infection. Immunity,47, 1024–1036.
Herbert, D. R., Yang, J. Q., Hogan, S. P., Groschwitz, K., Khodoun, M., Munitz, A., Orekov, T., Perkins, C., Wang, Q., Brombacher, F., et al. (2009). Intestinal epithelial cell secretion of RELM-beta protects against gastrointestinal worm infection. Journal of Experimental Medicine,206, 2947–2957.
Hewitson, J. P., Grainger, J. R., & Maizels, R. M. (2009). Helminth immunoregulation: The role of parasite secreted proteins in modulating host immunity. Molecular and Biochemical Parasitology,167, 1–11.
Hildersley, K. A., McNeilly, T. N., Gillan, V., Otto, T. D., Loser, S., Gerbe, F., Jay, P., Maizels, R. M., Devaney, E., & Britton, C. (2021). Tuft cells increase following ovine intestinal parasite infections and define evolutionarily conserved and divergent responses. Frontiers in Immunology,12, 781108.
Hodge, S. H., Krauss, M. Z., Kaymak, I., King, J. I., Howden, A. J. M., Panic, G., Grencis, R. K., Swann, J. R., Sinclair, L. V., & Hepworth, M. R. (2023). Amino acid availability acts as a metabolic rheostat to determine the magnitude of ILC2 responses. Journal of Experimental Medicine,220, e20221073.
Howitt, M. R., Cao, Y. G., Gologorsky, M. B., Li, J. A., Haber, A. L., Biton, M., Lang, J., Michaud, M., Regev, A., & Garrett, W. S. (2020). The taste receptor TAS1R3 regulates small intestinal tuft cell homeostasis. Immunohorizons,4, 23–32.
Howitt, M. R., Lavoie, S., Michaud, M., Blum, A. M., Tran, S. V., Weinstock, J. V., Gallini, C. A., Redding, K., Margolskee, R. F., Osborne, L. C., et al. (2016). Tuft cells, taste-chemosensory cells, orchestrate parasite type 2 immunity in the gut. Science,351, 1329–1333.
Kotas, M. E., O’Leary, C. E., & Locksley, R. M. (2023). Tuft cells: Context- and tissue-specific programming for a conserved cell lineage. Annual Review of Pathology: Mechanisms of Disease,18, 311–335.
Laing, R., Kikuchi, T., Martinelli, A., Tsai, I. J., Beech, R. N., Redman, E., Holroyd, N., Bartley, D. J., Beasley, H., Britton, C., et al. (2013). The genome and transcriptome of Haemonchus contortus, a key model parasite for drug and vaccine discovery. Genome Biology,14, R88.
Lei, W., Ren, W., Ohmoto, M., Urban, J. F., Jr., Matsumoto, I., Margolskee, R. F., & Jiang, P. (2018). Activation of intestinal tuft cell-expressed Sucnr1 triggers type 2 immunity in the mouse small intestine. Proc Natl Acad Sci U S A,115, 5552–5557.
Lloyd, C. M., & Snelgrove, R. J. (2018). Type 2 immunity: Expanding our view. Science Immunology,3, eaat1604.
Luo, X. C., Chen, Z. H., Xue, J. B., Zhao, D. X., Lu, C., Li, Y. H., Li, S. M., Du, Y. W., Liu, Q., Wang, P., et al. (2019). Infection by the parasitic helminth Trichinella spiralis activates a Tas2r-mediated signaling pathway in intestinal tuft cells. Proceedings of the National Academy of Sciences of the United States of America,116, 5564–5569.
Luu, M., Pautz, S., Kohl, V., Singh, R., Romero, R., Lucas, S., Hofmann, J., Raifer, H., Vachharajani, N., Carrascosa, L. C., et al. (2019). The short-chain fatty acid pentanoate suppresses autoimmunity by modulating the metabolic-epigenetic crosstalk in lymphocytes. Nature Communications,10, 760.
Maizels, R. M., Burke, J., & Denham, D. A. (1987). Phosphorylcholine-bearing antigens in filarial nematode parasites: Analysis of somatic extracts and in vitro secretions of Brugia malayi and B.pahangi and infection sera. Parasite Immunology,9, 49–66.
Maizels, R. M., & Gause, W. C. (2023). Targeting helminths: The expanding world of type 2 immune effector mechanisms. Journal of Experimental Medicine,220, e20221381.
Maizels, R. M., Smits, H. H., & McSorley, H. J. (2018). Modulation of host immunity by helminths: The expanding repertoire of parasite effector molecules. Immunity,49, 801–818.
Manosalva, P., Manohar, M., von Reuss, S. H., Chen, S., Koch, A., Kaplan, F., Choe, A., Micikas, R. J., Wang, X., Kogel, K. H., et al. (2015). Conserved nematode signalling molecules elicit plant defenses and pathogen resistance. Nature Communications,6, 7795.
McDaniel, M., Heber, L., & von Moltke, J. (2023). Initiation of type 2 immunity at barrier surfaces. Muc Immunol,16, 86–97.
Nadjsombati, M. S., McGinty, J. W., Lyons-Cohen, M. R., Jaffe, J. B., DiPeso, L., Schneider, C., Miller, C. N., Pollack, J. L., Nagana Gowda, G. A., Fontana, M. F., et al. (2018). Detection of succinate by intestinal tuft cells triggers a Type 2 innate immune circuit. Immunity,49(33–41), e37.
Ndjim, M., Gasmi, I., Herbert, F., Josephine, C., Bas, J., Lamrani, A., Coutry, N., Henry, S., Zimmermann, V. S., Dardalhon, V., et al. (2024). Tuft cell acetylcholine is released into the gut lumen to promote anti-helminth immunity. Immunity,57(1260–1273), e1267.
O’Leary, C. E., Schneider, C., & Locksley, R. M. (2019). Tuft cells-systemically dispersed sensory epithelia integrating immune and neural circuitry. Annual Review of Immunology,37, 47–72.
Péry, P., Luffau, G., Charley, J., Petit, A., Rouze, P., & Bernard, S. (1979). Phosphorylcholine antigens from Nippostrongylus brasiliensis. II. Isolation and partial characterization of phosphorylcholine antigens from adult worms. Annales D’immunologie (Institut Pasteur),130C, 889–900.
Pullan, R. L., Smith, J. L., Jasrasaria, R., & Brooker, S. J. (2014). Global numbers of infection and disease burden of soil transmitted helminth infections in 2010. Parasites & Vectors,7, 37.
Roeber, F., Jex, A. R., & Gasser, R. B. (2013). Impact of gastrointestinal parasitic nematodes of sheep, and the role of advanced molecular tools for exploring epidemiology and drug resistance - an Australian perspective. Parasites & Vectors,6, 153.
Schneider, C., O’Leary, C. E., von Moltke, J., Liang, H. E., Ang, Q. Y., Turnbaugh, P. J., Radhakrishnan, S., Pellizzon, M., Ma, A., & Locksley, R. M. (2018). A metabolite-triggered tuft cell-ILC2 circuit drives small intestinal remodeling. Cell,174(271–284), e214.
Shinoda, K., Choe, A., Hirahara, K., Kiuchi, M., Kokubo, K., Ichikawa, T., Hoki, J. S., Suzuki, A. S., Appleton, J. A., Aroian, R. V., et al. (2021). Nematode ascarosides attenuate mammalian type 2 inflammatory responses. Proceedings of the National Academy of Sciences of the United States of America,119, e2108686119.
Shoji, H., Kunugi, H., & Miyakawa, T. (2022). Acute and chronic effects of oral administration of a medium-chain fatty acid, capric acid, on locomotor activity and anxiety-like and depression-related behaviors in adult male C57BL/6J mice. Neuropsychopharmacology Reports,42, 59–69.
Thomas, J. D., Conrad, R. C., & Blumenthal, T. (1988). The C.elegans trans -spliced leader RNA is bound to Sm and has a trimethylguanosine cap. Cell,54, 533–539.
von Moltke, J., Ji, M., Liang, H. E., & Locksley, R. M. (2016). Tuft-cell-derived IL-25 regulates an intestinal ILC2-epithelial response circuit. Nature,529, 221–225.
Wang, J., Wu, X., Simonavicius, N., Tian, H., & Ling, L. (2006). Medium-chain fatty acids as ligands for orphan G protein-coupled receptor GPR84. Journal of Biological Chemistry,281, 34457–34464.
Wangchuk, P., Constantinoiu, C., Eichenberger, R. M., Field, M., & Loukas, A. (2019). Characterization of tapeworm metabolites and their reported biological activities. Molecules,24, 1480.
Wangchuk, P., Kouremenos, K., Eichenberger, R. M., Pearson, M., Susianto, A., Wishart, D. S., McConville, M. J., & Loukas, A. (2019b). Metabolomic profiling of the excretory-secretory products of hookworm and whipworm. Metabolomics,15, 101.
Wangchuk, P., Shepherd, C., Constantinoiu, C., Ryan, R. Y. M., Kouremenos, K. A., Becker, L., Jones, L., Buitrago, G., Giacomin, P., Wilson, D., et al. (2019). Hookworm-derived metabolites suppress pathology in a mouse model of colitis and inhibit secretion of key inflammatory cytokines in primary human leukocytes. Infection and Immunity,87, e00851.
Wangchuk, P., Yeshi, K., & Loukas, A. (2023). Metabolomics and lipidomics studies of parasitic helminths: Molecular diversity and identification levels achieved by using different characterisation tools. Metabolomics,19, 63.
Whitman, J. D., Sakanari, J. A., & Mitreva, M. (2021). Areas of metabolomic exploration for helminth infections. ACS Infectious Diseases,7, 206.
Xiong, Z., Zhu, X., Geng, J., Xu, Y., Wu, R., Li, C., Fan, D., Qin, X., Du, Y., Tian, Y., & Fan, Z. (2022). Intestinal Tuft-2 cells exert antimicrobial immunity via sensing bacterial metabolite N-undecanoylglycine. Immunity,55(686–700), e687.
Yeshi, K., Creek, D. J., Anderson, D., Ritmejeryte, E., Becker, L., Loukas, A., & Wangchuk, P. (2020). Metabolomes and lipidomes of the infective stages of the gastrointestinal nematodes. Nippostrongylus brasiliensis and Trichuris muris. Metabolites,10, 446.
Zaiss, D. M. W., Pearce, E. J., Artis, D., McKenzie, A. N. J., & Klose, C. S. N. (2024). Cooperation of ILC2s and T(H)2 cells in the expulsion of intestinal helminth parasites. Nature Reviews Immunology,24, 294–302.
Acknowledgements
The work was supported by the Wellcome Trust through a Collaborative Grant to RMM, CB, ED, TMcN and PJ (Ref 211814), and the Wellcome Trust core-funded Wellcome Centre for Integrative Parasitology (Ref: 104111). We gratefully acknowledge assistance and expertise from the Flow Core Facility and the Wolfson Research Facility at the University of Glasgow. For the purpose of open access, the author(s) has applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising from this submission. We gratefully acknowledge assistance and expertise from the Flow Core Facility and the Wolfson Research Facility at the University of Glasgow.
Author information
Authors and Affiliations
Contributions
MCP, SL, VVG, JR, CC, GB, EM and KH performed experiments, prepared and analysed samples; MCP, SL, VG, GB, MB, PJ, ED, TM, CB and RMM designed assays, advised and interpreted results; MCP and RMM drafted the main manuscript text, and all authors contributed edits for final version.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
11306_2025_2248_MOESM1_ESM.pdf
Supplementary file1 (PDF 133 KB)— Suppl Fig. 1: Ascarosides released by N. brasiliensis adult worms. 200 90 adults were incubated in 1 mL media; supernatants were recovered and extracted in chloroform:methanol:water for MS analysis. Data for mass 276.1276, corresponding to ascaroside #1 (276.1573). Suppl Fig. 2: Comparative intensity of 15 major metabolites released by N. brasiliensis and H. contortus, in same format as Figure 3 B, shown for EBSS (A) and RPMI1640 (B).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Poveda, M.C., Löser, S., Gillan, V. et al. Metabolomic and functional analyses of small molecules secreted by intestinal nematodes in the activation of epithelial tuft cells. Metabolomics 21, 55 (2025). https://doi.org/10.1007/s11306-025-02248-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s11306-025-02248-w