Metabolomic analysis of swainsonine poisoning in renal tubular epithelial cells

1 Introduction

Locoweed, a plant species with a global distribution (1), has been documented in various countries including the United States (2), Brazil (3), Argentina (4), Australia (5), and China (6). The earliest instances of poisonous plant poisoning diseases in the western region of the country were documented in the United States in 1873 (7). In China, the distribution of locoweed is mainly concentrated in Qinghai, Xinjiang, Tibet, Gansu, and other regions (8). Swainsonine (SW) is an indolizidine alkaloid, which is the main toxic component of locoweed (9–11). SW-induced toxicity is primarily characterized by the inhibition of α-mannosidase activity and the induction of widespread cellular vacuolar degeneration (12, 13). Previous studies have unequivocally demonstrated that SW can inflict damage on multiple organs, encompassing the cerebrum, cerebellum, liver, kidneys, pancreas, and thyroid gland (1). Notably, renal injury is observed at relatively lower dosages or at earlier time points (4, 14).

Metabolites are important characteristics of an organism’s phenotype, understanding the metabolomics provides important insights for the biological processes and their underlying mechanisms. Metabolomics is the systematic quantification and analysis of all small molecules (15) present in biological samples such as cells, tissues, or biological fluids. It involves the qualitative and quantitative characterization of these substances and their comprehensive analysis. Metabolomics has the ability to directly capture the molecular phenotype of a species (16) as well as the functional state of an individual (17–19). Through this analysis, metabolites of biological significance and statistically significant differences can be identified thus facilitating disease diagnosis and metabolite analysis. Metabolomics technology has been extensively employed in the identification of biomarkers for disease diagnosis and screening. Studies have utilized this technology to screen for specific metabolite associated with depression (20), neurodegenerative diseases (21), venous thromboembolism (22), osteoporosis (23), and diabetes mellitus (24). Currently, the metabolomics research surrounding locoweed and SW remains scant, with the bulk of studies focusing on elucidating the intricate relationship between the secondary metabolites of endophytic fungi and the biosynthesis of SW. Unfortunately, research papers exploring the intoxication mechanisms of SW are highly limited, with only a solitary study identified. In this study, the author cleverly harnessed targeted metabolomics and high-throughput sequencing techniques to reveal that SW has the capacity to profoundly alter bile acid metabolism and disrupt the delicate balance of intestinal microbiota in mice, ultimately triggering inflammatory reactions in the liver (25). Nevertheless, our understanding of the kidney toxicity and associated metabolites triggered by SW remains tenuous, necessitating further investigation in this vital area.

Currently, the prevention and control of livestock poisoning primarily rely on the implementation of preventive measures, as there is currently no specific antidote drug available for affected animals. As the early stage of locoweed poisoning is reversible, timely intervention can help mitigate the symptoms of locoweed poisoning by feeding non-toxic pasture grasses. Therefore, timely identification of the metabolic changes due to the locoweed poisoning is important for effective intervention. At present, studies on SW primarily focuses on its biosynthesis and toxicity mechanisms, with limited reports on the metabolic pathways it participates in after entering cells. Previous research has indicated that locoweed poisoning primarily affects the kidney, causing severe damage to this organ and the renal epithelial cells, especially the proximal tubular epithelial cells, are the most significantly affected (4, 14). In this study, rat primary renal tubular epithelial cells were utilized as an experimental model to investigate the metabolomic perturbation following locoweed poisoning. Untargeted metabolomics technology was employed to analyze the metabolic profile of SW-induced alteration in renal tubular cells in order to aid the diagnosis of locoweed poisoning. The identification of these metabolomic changes is important for the early detection and prevention of locoweed poisoning in livestock.

2 Materials and methods

2.1 Ethical approval

All procedures were conducted in accordance with the Code of Ethics and approved by the Laboratory Animal Management and Ethics Committee of Northwest A&F University, Project number: XN2023-1007. All efforts were made to minimize the suffering of animals.

2.2 Chemical reagents and instrumentation

Swainsonine (provided by the Laboratory of Animal Toxicology Team, Northwest A&F University, purity >98%), methanol (Merck, Germany), acetonitron (Merck, Germany), L-2-chlorophenylalanine (Aladin, China), formic acid (TCI, Japan), ultra-high performance liquid phase (Waters, United States), High resolution mass spectrometry (Waters, United States), chromatographic column (Waters, Unites States), and Cell incubator (Thermo, United States).

2.3 Cell culture and treatment

Primary rat renal tubular epithelial cells were obtained by tissue block culture, the detailed operational steps are as follows.

Euthanize the rats by decapitation, and aseptically extract the kidneys. Delicately remove the capsule and separate the cortical portion. Cut the cortical tissue into 1 mm³ tissue blocks. Thoroughly rinse the tissue blocks with PBS buffer until the supernatant becomes clear and free of turbidity. Centrifuge the tissue blocks at 1,000 rpm for 3 min. Add 0.1% type IV collagenase to the tissue blocks and incubate the mixture at 37°C for 1 h to disperse the tissue. To terminate the enzymatic digestion, add an equal volume of culture medium to the mixture. Pass the content of the centrifuge tube through a 100-mesh screen to collect the filtered solution. Repeat the filtration process using a 400-mesh screen and collect the upper layer material, which represents the renal tubule segments. Resuspend the renal tubule segments in an appropriate amount of culture medium. Transfer the resuspended renal tubule segments to culture dishes and incubate them in a cell culture incubator. Replace the culture medium every 48 h to maintain the optimal growth conditions.

When the cells were fused to 80%, SW is dissolved in the culture medium to achieve a concentration of 0.8 mg/mL. Subsequently, the SW-containing culture medium is filtered through a 0.22 μm filter membrane to ensure sterility. The sterilized medium is then applied to the cells and incubated for 12 h.

The cell culture medium used in our study was Dulbecco’s Modified Eagle Medium (Gibco, America, 12800-017), supplemented with 15% Fetal Bovine Serum Standard (Newzerum, New Zealand, FBS-S500) for cultivation, without the addition of any other substances.

2.4 Metabolites extraction

The extraction solution with internal label (Vmethanol:Vacetonitrile = 1:1, internal standard concentration 20 mg/L) was prepared, and the extraction solution was added three times (300, 300, and 400 μL), and the cell samples were completely transferred into the EP tube, and vortex mixed for 30 s. Add steel beads, process for 10 min on a 45 Hz grinder, sonicate for 10 min (ice-water bath); leave at −20°C for 1 h; centrifuge at 12,000 rpm for 15 min at 4°C, and extract the supernatant (500 μL); dry the extract in a vacuum concentrator; after drying, add 160 μL of the extract solution (V water:V acetonitrile = 1:1) to reconstitute the extract solution; vortex for 30 s, and sonicate for 10 min in an ice-water bath. The extract was centrifuged at 12,000 rpm for 15 min at 4°C. 120 μL of the supernatant was taken into the injection bottle, and 10 μL of each sample was mixed to form quality control samples (QC) for the assay.

2.5 LC–MS/MS analysis

The liquid-mass spectrometry system for metabolomics analysis consisted of an Acquity I-Class PLUS (Waters, United States) ultra-high performance liquid chromatography (UHPLC) tandem with a Waters Xevo G2-XS QTOF (Waters, United States) high-resolution mass spectrometer. The column used was Acquity UPLC HSS T3 column (1.8 μm, 2.1 mm × 100 mm) (Waters, United States). The mobile phase of positive ion mode (POS) and negative ion mode (26) were the same: mobile phase A: 0.1% formic acid aqueous solution; mobile compositions B: 0.1% formic acid acetonitrile, and the injection volume was 1 μL. The elution gradients were: 98% mobile phase A, 2% mobile phase B, 0–0.25 min; 2% mobile phase A, 98% mobile phase B, 10–13 min; 98% mobile phase A, 2% mobile phase B, 13.1–15 min; and the flow rate was 400 μL/min. The high-resolution mass spectrometer (Xevo G2-XS QTOF, Waters, United States) is capable of primary and secondary mass spectrometry data acquisition in MSe mode under the control of acquisition software (MassLynx V4.2, Waters, United States). Dual-channel data acquisition for both low and high crash energy is performed in each data acquisition cycle. The low collision energy was 2 V, the high collision energy interval was 10–40 V, and the scanning frequency was 0.2 s. The parameters of the ESI ion source were as follows: capillary voltage: 2,500 V (positive ion mode) or −2,000 V (negative ion mode); cone-well voltage: 30 V; temperature of the ion source: 100°C; temperature of the desolvent gas was 500°C; flow rate of the blowback gas was 50 L/h; flow rate of the desolvent gas was 800 L/h; and mass-to-nucleus ratio (m/z) acquisition range was 50–1,200. Flow rate: 800 L/h; mass-to-core ratio (m/z) acquisition range 50–1,200.

2.6 Characterization and quantification of metabolites

The raw data collected by MassLynx V4.2 were processed by Progenesis QI software for peak extraction and peak alignment, and then identified based on the Progenesis QI software online METLIN database, Human Metabolome Database (HMDB), and Biomark’s self-built library, and the theoretical fragments were also identified, and the deviation of mass number of parent ions is 100 ppm, and that of fragments is 50 ppm or less. The deviation of the parent ion mass number is 100 ppm and the deviation of the fragment ion mass number is 50 ppm or less (27).

2.7 Classification and functional annotation of metabolomes

The metabolites were annotated using the HMDB database¹ to obtain superclass and class information; the pathways in which the identified metabolites were found were annotated using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database² to annotate the pathways in which the identified metabolites were found.

2.8 Data processing and analysis

Principal component analysis, which downscales the dimension of high-throughput metabolic data and categorizes them by the similarity of their principal components, is an unsupervised classification mode that responds to the overall distribution of samples among groups and the magnitude of differences between samples within groups. The supervised discriminant analysis statistical method of partial least squares regression orthogonal projections to latent structures-discriminant analysis (OPLS-DA) was used to visualize within-group differences by filtering out orthogonal variables unrelated to categorical variables. To check the reliability of the OPLS-DA model, a permutation test is required. The groupings of the samples were randomly disrupted and the OPLS-DA modeling was performed and the evaluation parameters coefficient of determination for Y (R²Y) and predictability for Y (Q²Y) were obtained according to the new groupings. Variable importance in projection (VIP) reflects both the loading weight of each metabolite in the model and the variability of the response explained by that metabolite, and can be used for ANOVA. In order to analyze the metabolic patterns of metabolites under different experimental conditions, all the metabolites with differences between the obtained comparison pairs were clustered into classes with the same or similar metabolic patterns for hierarchical clustering analysis.

3 Results

3.1 Quality control of survey data

The base peak intensity (BPI) charts monitor the strongest peaks in each chromatogram. The strongest peak intensity at each point in the analysis is shown (Figure 1). The BPI from the QC sample has a good overlap of peak retention time and peak area, indicating good instrument stability. The positive ion mode 17,312 peaks were identified and 2,563 metabolites were annotated, of which 94.25% had relative standard deviation (RSD) ≥ 0.7; in the negative ion mode 7,825 peaks were identified and 1,990 metabolites were annotated, of which 86.81% had RSD ≥ 0.7 (Table 1). If the RSD is greater than 60%, it is an indication that the sample is reliable and can be used for follow-up testing.

Figure 1

Figure 1. Base peak intensity (BPI) overlapping diagram of QC samples (A: positive ion mode; B: negative ion mode).

Table 1

Table 1. List of quality control indicators.

3.2 Principal component analysis

By performing principal component analysis (PCA) on the samples, we can preliminarily understand the overall metabolic differences between the samples of each group and the degree of variation between the samples within the group. After PCA analysis of the total samples, the results (Figure 2) showed that the two groups of samples were divided into two different groups, and each sample was within the ellipse of the 95% confidence interval, indicating that the samples in the group had good repeatability. The two ellipses were independent of each other and there was no overlapping area, indicating that there was a large difference between the two groups of samples.

Figure 2

Figure 2. Principal component analysis (PCA) of all samples (A: positive ion mode; B: negative ion mode). Where the x-axis indicates the first principal component, the y-axis indicates the second principal component, and the percentage axis indicates the percentage contribution of that principal component to the sample variance. Each point on the graph represents a sample, samples in the same group are shown in the same color and samples in different subgroups are shown in different colors. Ellipses represent 95% confidence intervals.

3.3 Orthogonal partial least squares discriminant analysis

Orthogonal partial least squares discriminant analysis (OPLS-DA) is a supervised analysis that can better screen for differential metabolites by excluding influences that are not relevant to the study. OPLS-DA analysis was performed on the control and SW treatment groups and the results are shown in Figure 3. Under the two modes, the samples of the two groups were clearly grouped, indicating that the metabolites of the two groups were significantly different.

Figure 3

Figure 3. Orthogonal projections to latent structures–discriminant analysis (OPLS-DA) score chart (A: positive ion mode; B: negative ion mode). The blue point A in the image is the control group and the red point B is the SW treatment group. The x-axis (t1) represents the prediction component (the between-group variance component), the y-axis (t2) represents the orthogonal component (the within-group variance component), and the cross y-axis percentage represents the component’s proportion of the total variance. The model parameters are listed below, including R²X, R²Y, Q²Y, root mean square error (RMSEE), pre (number of predicted components), and ort (number of orthogonal components).

To check the reliability of the OPLS-DA model, a replacement test is required. This means that the grouping of the samples is randomly scrambled, and the OPLS-DA is modeled according to the scrambled grouping and its R²Y and Q²Y are calculated. After many iterations, the results of multiple modeling are plotted on a scatterplot, as shown in Figure 4, where Q²Y > 0.9 indicates that the model is excellent. The intercept of the Q²Y fitting regression line is negative, indicating that the model does not overfit; the slope of the Q²Y fitting regression line is positive, indicating that the model is meaningful, and the blue point is mostly above the red point, indicating that the independence of the modeling training set and test set is good. In the negative ion mode, the blue point overlaps the red point. This shows that the independence of the negative ion mode is not as good as that of the positive ion mode.

Figure 4

Figure 4. Replacement test diagram of OPLS-DA model (A: positive ion model; B: negative ion mode). In this image, the X-axis represents the correlation between the replacement group and the original model group, the Y-axis represents the value of R²Y or Q²Y (where Q²Y and R²Y of 1 are the values of the original model), the blue point and the red point represent the R²Y and Q²Y of the model after replacement, respectively, and the two dotted lines represent the regression line fitted by R²Y and Q²Y.

3.4 Metabolite classification and functional annotation

Metabolites in positive and negative ion mode are combined and annotated in HMDB as shown in Figure 5. The metabolites mainly include benzenoids, lipids and lipid-like molecules, nucleosides, nucleotides and analogs, organic acids and derivatives, organic nitrogen compounds, organic oxygen compounds, organoheterocyclic compounds, phenylpropanoids and polyketides, eterocyclic compounds, and pids and lipid-like molecules.

Figure 5

Figure 5. Human Metabolome Database (HMDB) classification summary. Items in the same box in the image represent HMDB level classification information, corresponding to the super class and class information of the HMDB database. Column length represents the number of metabolites annotated by the classification.

The KEGG database is a major public database related to metabolic pathways, which can be used for metabolic analysis in vivo. KEGG database is used to annotate all identified metabolites. Figure 6 shows the most annotated information (top 20 displayed) of Knock Out (KO) pathway level 2. The main metabolic pathways affected include amino acid metabolism, cancer: overview, digestive system, lipid metabolism, membrane transport, metabolism of cofactors and vitamins, overview, digestive system, lipid metabolism, membrane transport, metabolism of cofactors and vitamins, metabolism of other amino acids, metabolism of terpenoids and polyketides, nervous system, metabolism of other amino acids, metabolism of terpenoids and polyketides, nervous system, Nucleotide metabolism, and Xenobiotics biodegradation and metabolism.

Figure 6

Figure 6. Kyoto Encyclopedia of Genes and Genomes (KEGG) database classification summary. The entries under the same box in the image represent the level classification notes of the KEGG pathway, corresponding to KO pathway Level 1 and KO pathway Level 2. The length of the column indicates the number of metabolites annotated by the pathway.

3.5 Differential metabolite analysis

Based on the results of OPLS-DA, the VIP of the OPLS-DA model can be obtained by multivariate analysis, and the metabolites of different varieties or tissues can be initially screened out. The screening criteria are ① FC ≥ 1; ② VIP ≥ 1; and ③ p < 0.05.

By combining metabolites in positive and negative ion mode, a total of 2,170 differential metabolites were identified, of which 1,342 were upregulated and 828 were downregulated. Details of the differential metabolites are shown in Supplementary Table 1.

3.5.1 Difference multiple analysis

The Fold Changes were compared after qualitative and quantitative analysis of the detected metabolites Table 2 shows the results of upregulated and downregulated log₂FC of the top 20 metabolites in the test group compared with the control group after log conversion treatment for differential metabolite multiple. Among them, the top 5 metabolites with the highest upregative difference multiple were Ttricosanoylglycine, 1-Tricosanol, Sambacolignoside, Kanzonol T, and Nemonoxacin, respectively. The top 5 metabolites with the largest down-variance were 2-(Ethylamino)-4,5-dihydroxybenzamide, Cinobufotalin, Undecylprodigiosin, Hypochoeroside A, and Muzanzagenin.

Table 2

Table 2. Results of differential metabolite multiple analysis.

3.5.2 Volcanic map of differential metabolites

The volcano plot can directly show the overall distribution of the difference in metabolite content in the two groups, and the statistical significance of the difference in metabolite content. As shown in Figure 7, blue dots represent downregulated differentially expressed metabolites, red dots represent upregulated differentially expressed metabolites, and gray dots represent metabolites with insignificant differences. After sorting by p value, the first five metabolite names identified are shown in the Figure 7.

Figure 7

Figure 7. Volcano map of differential metabolites. In the volcano map, each point represents a metabolite. The x-axis represents the fold change of different substances in the comparison group (log₂ transformed), and the y-axis represents the p value (log₁₀ transformed). The size of the scatter point indicates the VIP value of the OPLS-DA model. Larger scatter points correspond to higher VIP values, indicating more reliable selection of differentiated expression metabolites.

3.5.3 Differential metabolite cluster analysis

Cluster Analysis is a common multivariate statistical analysis method, the quantification of metabolites is expressed by color, the redder the color, the higher its content in a sample, and conversely, the greener the color, the lower the content, which is usually simple and intuitive to observe the overall characteristics of the data. Hierarchical cluster analysis was performed on the differential metabolites, and the results were shown in Figure 8. As can be seen from the figure, the color difference of the selected metabolites in the group is small, indicating that the content is similar, and the color difference between different treatment groups is large and the change is obvious, indicating that these different metabolites can clearly distinguish the control group and the SW group.

Figure 8

Figure 8. Differential metabolism clustering heat map. Group A served as the control group, while Group B was treated with SW. On the x-axis, each sample is represented, and the y-axis represents the quantitative value of metabolites standardized by Z-score after hierarchical clustering.

3.6 Enrichment analysis of differential metabolite metabolic pathways

The cluster profiler was used to analyze the annotation results of the selected differential metabolites in KEGG by means of hypergeometric test, and the differential metabolite enrichment diagram was drawn (only the top 20 results were displayed). The results were shown in Figure 9, where the size of the points represented the number of enriched differential metabolites. The color of the dots represents p value, and the redder the color, the more significant the enrichment. The top 6 channels enriched to Bile secretion, Primary bile acid biosynthesis, Riboflavin metabolism, Ferroptosis, Drug metabolism-cytochrome P450, and pyrimidine metabolism. The five pathways with the most significant differences are Primary bile acid biosynthesis, Ferroptosis, Riboflavin metabolism and Valine, leucine, and isoleucine degradation.

Figure 9

Figure 9. Kyoto Encyclopedia of Genes and Genomes enrichment map of differential metabolites. The X-axis represents the ratio of the number of differentially expressed metabolites in the corresponding pathway to the total number of metabolites detected and annotated by the pathway. The Y-axis corresponds to the pathway names. The color depth of the dots represents the log p value, with a deeper shade of red indicating a more significant enrichment. The size of the dots indicates the number of enriched differentiated metabolites.

3.7 Screening for specific metabolites associated with the SW toxicity

Based on metabolic specificity and the size of VIP value and log-fold change (LFC) of the difference multiple was taken into account. Among them, we have identified five metabolites associated with SW toxicity which were shown in Table 3 and diagnostic accuracy was evaluated by the area under the curve.

Table 3

Table 3. Differential metabolites associated with renal tubular epithelial cells treated by SW.

4 Discussion

This study utilized untargeted metabolomics to analyze the differential metabolites in renal tubular epithelial cells treated with SW compared to the control group. A total of 2,170 significantly different metabolites were screened, most of which were not annotated due to limited related research. The annotated differential metabolites were mainly enriched in the following pathways: Bile secretion, Primary bile acid biosynthesis, Riboflavin metabolism, Ferroptosis, Drug metabolism-cytochrome P450, and Primidine metabolism. Among these differential metabolites, SW is the most distinctive one. However, due to its short half-life of only 20 h (9, 12), it may not be detectable, thus requiring corroboration with other indicators. Therefore, combining VIP values and Log₂FC, we have screened out substances such as 3alpha,7alpha-Dihydroxy-5beta-cholestanate, 2-Hydroxyiminostilbene, and glycochenodeoxycholate as potential diagnostic values for SW poisoning.

Bile secretion and Primary bile acid biosynthesis were the pathways with the highest number of differentially expressed metabolites. Among them, glycochenodeoxycholate, a component of primary bile acids and a raw material for bile acid synthesis, 3α,7α,12α,26-tetrahydroxy-5β-cholestane, and 7a-Hydroxy-cholestene-3-one, intermediates in the bile acid synthesis pathway (28), were significantly elevated. Previous results genomic study showed consistent enrichment in bile metabolism and drug metabolism-cytochrome P450 when rat renal tubular epithelial cells were exposed to the same conditions as in this study. Additionally, Fu et al. (25) demonstrated that SW induced liver inflammation by altering bile acid metabolism and intestinal microbiota in mice. Thus, these results suggest that exposure to SW may significantly alter the bile acid metabolism whether these changes is associated with the SW-induced cellular damage remain to be investigated.

Cytochrome P450 is one of the most common enzyme families involved in metabolism of the xenobiotics, playing a crucial role in the metabolism of drugs, toxins, and endogenous substrates (29). It plays a crucial role in the clearance of compounds such as drugs, toxins, and endogenous substrates. Our results have found significant increases in 2-Hydroxyiminostilbene, carbamazepine iminoquinone, carbamazepine-O-quinone, morphine-6-glucuronide, and 6-Methylmercaptopurine. These metabolites are associated with the metabolism of Carbamazepine, morphine, and capecitabine, and they rely on P450 enzymes for catalysis (30), suggesting that cells may activate drug metabolism-cytochrome P450 to facilitate the clearance of drugs after exposure to SW. The particular cytochrome P450 induced by the SW however need to be further investigated.

Ferroptosis is a new type of cell death discovered in recent years. Many studies have shown that ferroptosis is related to kidney damage (31–33). A notable elevation in glycochenodeoxycholate (GCDA) levels has been observed in patients suffering from chronic kidney disease. Wang et al. (34) demonstrated that that glycochenodeoxycholate GCDA has the capacity to activate the expression of hepcidin, profoundly influencing the iron homeostasis within the body. Disruption of this iron balance can lead to an excessive accumulation of iron ions, particularly Fe²⁺, which can subsequently trigger ferroptosis (34). Nevertheless, the intricate relationship between GCDA and ferroptosis remains to be thoroughly elucidated and validated.

There are several documented research on related toxicants. For example, another indole alkaloid, vincristine, its impact on plasma metabolites in human children was mainly enriched in the pathways of purine metabolism, arginine biosynthesis, sphingolipid metabolism, glutathione metabolism, glycerophospholipid metabolism, and lysine degradation (35) This differs significantly from the metabolic pathways affected by SW, possibly due to differences in species and tissues of the samples. Rao et al. (36) found that exposure to Matrine in mice resulted in significant changes in metabolites such as cholic acid, taurocholic acid, L-tyrosine, flavin mononucleotide (FMN), UDP-glucuronate, urea, sulfate, and inosine monophosphate, where riboflavin metabolism, purine metabolism, and ascorbate and aldarate metabolism were the most significant pathways. In our results, we also observed significant changes in FMN and UDP, which are involved in riboflavin and pyrimidine metabolism, respectively, suggesting the potential importance of these pathways in the metabolism of alkaloid toxins.

Riboflavin is a component of vitamin B complex and plays an important role in maintaining body health. Flavin mononucleotide (FMN), a coenzyme of riboflavin, is involved in reactions, including drug, lipid and heterogeneous metabolism, energy metabolism, and cellular processes such as protein folding and cell signaling (37, 38). Cytochrome P450 is one of the most common enzyme families involved in drug metabolism, most P450 enzymes use NAD(P)H as a coenzyme, which undergoes dehydrogenation to form NAD(P) and release electrons. The electron is then transferred to FMN, converting it to FMNH+2 (39). These results suggested that SW might induce cell death by inducing ferroptosis, riboflavin metabolism, and drug metabolism-cytochrome P450.

Untargeted metabolomic study can provide the broadest coverage of metabolites and help undiscovered specific metabolites associated with SW toxicity, giving us an objective and comprehensive understanding of the link between the body’s metabolism and external stimuli. It has been shown to be useful in studying the complex diseases. For example, metabolomics has now become a good exploration tool in psychiatry (40, 41). In veterinary medicine, these characteristic metabolic changes related to livestock plant poisoning have great practical application. Clinical symptoms in livestock often have certain specificity but only show when poisoning is severe enough, while these specific changes in the metabolites can help diagnosis of earlier and faster (7), so as to respond more quickly taking corresponding measures to minimize losses. Our study shows that after SW treatment, based on metabolic specificity and the VIP values and Log₂FC of differential metabolites, we have identified five metabolites which may be used for diagnosis for SW poisoning.

In summary, the findings of this study reveal that the primary metabolites encompass benzenoids, lipids and lipid-like molecules, nucleosides, nucleotides, and analogs, along with organic acids and derivatives. Notably, the differential metabolites are predominantly enriched in specific pathways such as bile secretion, primary bile acid biosynthesis, riboflavin metabolism, ferroptosis, drug metabolism-cytochrome P450, and primidine metabolism. Upon screening for potential biomarkers of SW poisoning, it has been observed that each metabolite exhibits an area under the ROC curve (AUC) value were 1, indicating a significant distinction between the control group and the SW group. Consequently, these metabolites are viable candidates for potential biomarkers of SW poisoning. However, the diagnostic values of these metabolites for SW poisoning need to be further validated in a variety of livestock and to rule out the potential effects of species, feeding environment, grazing practices, and other factors including the coexisting diseases.

Certainly, this study has certain limitations as it was conducted solely as an in vitro experiment, unable to fully simulate the in vivo environment. In the future, in vivo experiments could be conducted on the target animal, sheep, with various time gradients set to more precisely screen for early biomarkers of SW poisoning, thereby improving the detection and prevention of locoweed poisoning.

Data availability statement

The data presented in the study are deposited in the EMBL-EBI MetaboLights database with the identifier MTBLS9951.

Ethics statement

The animal study was approved by Laboratory Animal Management and Ethics Committee of Northwest A&F University. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

SZ: Formal analysis, Methodology, Validation, Writing – original draft. YiZ: Methodology, Writing – original draft. HY: Methodology, Writing – original draft. YL: Formal analysis, Writing – original draft. LT: Software, Writing – original draft. YaZ: Software, Writing – original draft. PS: Formal analysis, Writing – original draft. KW: Software, Writing – original draft. BZ: Conceptualization, Writing – review & editing. HL: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was funded by the grants from the National Natural Science Foundation of China (Nos. 32072929 and 32072928) and National Key R&D Program of China (No.2023YFD1801100).

Acknowledgments

Widely targeted/non-target metabolomics experiment was carried out at Biomarker Technologies (Beijing, China). We extend our sincere gratitude to them for providing invaluable technical support throughout the process. And we deeply appreciate the insightful and invaluable feedback provided by Yanan Tian, which has significantly contributed to the improvement of our manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1387853/full#supplementary-material

Footnotes

1. ^https://hmdb.ca/

2. ^http://www.genome.jp/kegg/

References

1. Marin, RE, Micheloud, JF, Vignale, ND, Gimeno, EJ, O'Toole, D, Gardner, DR, et al. Intoxication by Astragalus garbancillo var. garbancillo in llamas. J Vet Diagn Invest. (2020) 32:467–70. doi: 10.1177/1040638720914338

PubMed Abstract | Crossref Full Text | Google Scholar

2. Cook, D, Gardner, DR, Lee, ST, Pfister, JA, Stonecipher, CA, and Welsh, SL. A swainsonine survey of North American Astragalus and Oxytropis taxa implicated as locoweeds. Toxicon. (2016) 118:104–11. doi: 10.1016/j.toxicon.2016.04.033

PubMed Abstract | Crossref Full Text | Google Scholar

3. Riet-Correa, F, Medeiros, RM, and Schild, AL. A review of poisonous plants that cause reproductive failure and malformations in the ruminants of Brazil. J Appl Toxicol. (2012) 32:245–54. doi: 10.1002/jat.1754

PubMed Abstract | Crossref Full Text | Google Scholar

4. Marin, RE, Gardner, DR, Armien, AG, Fortunato, RH, and Uzal, FA. Intoxication of llamas by Astragalus punae in Argentina. J Vet Diagn Invest. (2022) 34:674–8. doi: 10.1177/10406387221094272

PubMed Abstract | Crossref Full Text | Google Scholar

5. O'Sullivan, BM, and Goodwin, JA. An outbreak of Swainsona poisoning in horses. Aust Vet J. (1997) 53:446–7. doi: 10.1111/j.1751-0813.1977.tb05496.x

PubMed Abstract | Crossref Full Text | Google Scholar

6. Zhao, M, Gao, X, Wang, J, He, X, and Han, B. A review of the most economically important poisonous plants to the livestock industry on temperate grasslands of China. J Appl Toxicol. (2013) 33:9–17. doi: 10.1002/jat.2789

PubMed Abstract | Crossref Full Text | Google Scholar

7. Green, BT, Welch, KD, Lee, ST, Stonecipher, CA, Gardner, DR, Stegelmeier, BL, et al. Biomarkers and their potential for detecting livestock plant poisonings in Western North America. Front Vet Sci. (2003) 10:1104702. doi: 10.3389/fvets.2023.1104702

Crossref Full Text | Google Scholar

8. Zhou, KS, Tian, W, Zhang, Z, and Tan, CJ. Research advances on chemical constituents of Oxytropis ochrocephala Bunge and its bioactivities. Nat Product Res Dev. (2020) 32:161–71. doi: 10.16333/j.1001-6880.2020.1.022

Crossref Full Text | Google Scholar

9. Li, B, Zhang, C, Zhu, Y, Sun, P, Fan, S, Wang, W, et al. Development of novel formulation for sustained release of drug to prevent Swainsonine-containing plants poisoning in livestock. Animals. (2023) 13:2646. doi: 10.3390/ani13162646

PubMed Abstract | Crossref Full Text | Google Scholar

10. Wang, S, Tan, P, Wang, H, Wang, J, Zhang, C, Lu, H, et al. Swainsonine inhibits autophagic degradation and causes cytotoxicity by reducing CTSD O-GlcNAcylation. Chem Biol Interact. (2023) 382:110629. doi: 10.1016/j.cbi.2023.110629

PubMed Abstract | Crossref Full Text | Google Scholar

11. Wu, C, Han, T, Lu, H, and Zhao, B. The toxicology mechanism of endophytic fungus and swainsonine in locoweed. Environ Toxicol Pharmacol. (2016) 47:38–46. doi: 10.1016/j.etap.2016.08.018

PubMed Abstract | Crossref Full Text | Google Scholar

12. Stegelmeier, BL, James, LF, Panter, KE, Ralphs, MH, Gardner, DR, Molyneux, RJ, et al. The pathogenesis and toxicokinetics of locoweed (Astragalus and Oxytropis spp.) poisoning in livestock. J Nat Toxins. (1999) 8:35–45.

PubMed Abstract | Google Scholar

13. Cholich, LA, Martinez, A, Micheloud, JF, Pistán, ME, García, EN, Robles, CA, et al. Alpha-mannosidosis caused by toxic plants in ruminants of Argentina. An Acad Bras Cienc. (2021) 93:e20191496. doi: 10.1590/0001-3765202120191496

PubMed Abstract | Crossref Full Text | Google Scholar

14. Wu, CC, Wang, W, Liu, XX, Ma, F, Cao, D, Yang, XW, et al. Pathogenesis and preventive treatment for animal disease due to locoweed poisoning. Environ Toxicol Pharmacol. (2014) 37:336–47. doi: 10.1016/j.etap.2013.11.013

PubMed Abstract | Crossref Full Text | Google Scholar

15. Ramírez, J, House, LK, Karrison, TG, Janisch, LA, Turcich, M, Salgia, R, et al. Prolonged pharmacokinetic interaction between capecitabine and a CYP2C9 substrate, celecoxib. J Clin Pharm Ther. (2019) 59:1632–40. doi: 10.1002/jcph.1476

PubMed Abstract | Crossref Full Text | Google Scholar

16. Zhan, Q, Wang, L, Liu, N, Yuan, Y, Deng, L, Ding, Y, et al. Serum metabolomics study of narcolepsy type 1 based on ultra-performance liquid chromatography-tandem mass spectrometry. Amino Acids. (2023) 55:1247–59. doi: 10.1007/s00726-023-03315-z

PubMed Abstract | Crossref Full Text | Google Scholar

17. Witting, M. (Re-)use and (re-)analysis of publicly available metabolomics data. Proteomics. (2023) 23:e2300032. doi: 10.1002/pmic.202300032

PubMed Abstract | Crossref Full Text | Google Scholar

18. Wang, Z, Zhu, H, and Xiong, W. Advances in mass spectrometry-based multi-scale metabolomic methodologies and their applications in biological and clinical investigations. Sci Bull. (2023) 68:2268–84. doi: 10.1016/j.scib.2023.08.047

PubMed Abstract | Crossref Full Text | Google Scholar

19. Liu, ML, Zheng, P, Liu, Z, Xu, Y, Mu, J, Guo, J, et al. GC-MS based metabolomics identification of possible novel biomarkers for schizophrenia in peripheral blood mononuclear cells. Mol BioSyst. (2014) 10:2398–406. doi: 10.1039/c4mb00157e

PubMed Abstract | Crossref Full Text | Google Scholar

20. Rimti, FH, Shahbaz, R, Bhatt, K, and Xiang, A. A review of new insights into existing major depressive disorder biomarkers. Heliyon. (2023) 9:e18909. doi: 10.1016/j.heliyon.2023.e18909

PubMed Abstract | Crossref Full Text | Google Scholar

21. Dong, R, Lu, Q, Kang, H, Suridjan, I, Kollmorgen, G, Wild, N, et al. CSF metabolites associated with biomarkers of Alzheimer's disease pathology. Front Aging Neurosci. (2023) 15:1214932. doi: 10.3389/fnagi.2023.1214932

PubMed Abstract | Crossref Full Text | Google Scholar

22. Della Corte, V, Riolo, R, Scaglione, S, Pecoraro, R, and Tuttolomondo, A. The role of biomarkers, metabolomics, and COVID-19 in venous thromboembolism-a review of literature. Int J Mol Sci. (2023) 24:13411. doi: 10.3390/ijms241713411

PubMed Abstract | Crossref Full Text | Google Scholar

23. Li, Y, Si, Y, Ma, Y, and Yin, H. Application and prospect of metabolomics in the early diagnosis of osteoporosis: a narrative review. Bioanalysis. (2023) 15:1369–79. doi: 10.4155/bio-2023-0131

PubMed Abstract | Crossref Full Text | Google Scholar

24. Shahisavandi, M, Wang, K, Ghanbari, M, and Ahmadizar, F. Exploring metabolomic patterns in type 2 diabetes mellitus and response to glucose-lowering medications-review. Gene. (2023) 14:1464. doi: 10.3390/genes14071464

PubMed Abstract | Crossref Full Text | Google Scholar

25. Fu, K, Chen, X, Shou, N, Wang, Z, Yuan, X, Wu, D, et al. Swainsonine induces liver inflammation in mice via disturbance of gut microbiota and bile acid metabolism. J Agric Food Chem. (2023) 71:1758–67. doi: 10.1021/acs.jafc.2c08519

PubMed Abstract | Crossref Full Text | Google Scholar

26. van Zogchel, LMJ, Zappeij-Kannegieter, L, Javadi, A, Lugtigheid, M, Gelineau, NU, Lak, NSM, et al. Specific and sensitive detection of neuroblastoma mRNA markers by multiplex RT-qPCR. Cancer. (2021) 13:150. doi: 10.3390/cancers13010150

PubMed Abstract | Crossref Full Text | Google Scholar

27. Wang, J, Zhang, T, Shen, X, Liu, J, Zhao, D, Sun, Y, et al. Serum metabolomics for early diagnosis of esophageal squamous cell carcinoma by UHPLC-QTOF/MS. Metabolomics. (2016) 12:116. doi: 10.1007/s11306-016-1050-5

Crossref Full Text | Google Scholar

28. Deo, AK, and Bandiera, SM. Identification of human hepatic cytochrome P450 enzymes involved in the biotransformation of cholic and chenodeoxycholic acid. Drug Metab Dispos. (2008) 36:1983–91. doi: 10.1124/dmd.108.022194

PubMed Abstract | Crossref Full Text | Google Scholar

29. Zhao, M, Ma, J, Li, M, Zhang, Y, Jiang, B, Zhao, X, et al. Cytochrome p450 enzymes and drug metabolism in humans. Int J Mol Sci. (2021) 22:12808. doi: 10.3390/ijms222312808

Crossref Full Text | Google Scholar

30. Pearce, RE, Uetrecht, JP, and Leeder, JS. Pathways of carbamazepine bioactivation in vitro: II. The role of human cytochrome P450 enzymes in the formation of 2-hydroxyiminostilbene. Drug Metab Dispos. (2005) 33:1819–26. doi: 10.1124/dmd.105.004861

PubMed Abstract | Crossref Full Text | Google Scholar

31. Li, J, Cao, F, Yin, HL, Huang, ZJ, Lin, ZT, Mao, N, et al. Ferroptosis: past, present, and future. Cell Death Dis. (2020) 11:88. doi: 10.1038/s41419-020-2298-2

PubMed Abstract | Crossref Full Text | Google Scholar

32. Tonnus, W, Meyer, C, Steinebach, C, Belavgeni, A, von Mässenhausen, A, Gonzalez, NZ, et al. Dysfunction of the key ferroptosis-surveilling systems hypersensitizes mice to tubular necrosis during acute kidney injury. Nat Commun. (2021) 12:4402. doi: 10.1038/s41467-021-24712-6

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ni, L, Yuan, C, and Wu, X. Targeting ferroptosis in acute kidney injury. Cell Death Dis. (2022) 13:182. doi: 10.1038/s41419-022-04628-9

PubMed Abstract | Crossref Full Text | Google Scholar

34. Wang, LJ, Zhao, GP, Wang, XF, Liu, XX, Li, YX, Qiu, LL, et al. Glycochenodeoxycholate affects Iron homeostasis via up-regulating hepcidin expression. Nutrients. (2022) 14:3176. doi: 10.3390/nu14153176

PubMed Abstract | Crossref Full Text | Google Scholar

35. Verma, P, Devaraj, J, Skiles, JL, Sajdyk, T, Ho, RH, Hutchinson, R, et al. A metabolomics approach for early prediction of vincristine-induced peripheral neuropathy. Sci Rep. (2020) 10:9659. doi: 10.1038/s41598-020-66815-y

PubMed Abstract | Crossref Full Text | Google Scholar

36. Rao, SW, Duan, YY, Zhao, DS, Liu, CJ, Xu, SH, Liang, D, et al. Integrative analysis of transcriptomic and metabolomic data for identification of pathways related to matrine-induced hepatotoxicity. Chem Res Toxicol. (2022) 35:2271–84. doi: 10.1021/acs.chemrestox.2c00264

PubMed Abstract | Crossref Full Text | Google Scholar

37. Lee, TY, Farah, N, Chin, VK, Lim, CW, Chong, PP, Basir, R, et al. Medicinal benefits, biological, and nanoencapsulation functions of riboflavin with its toxicity profile: a narrative review. Nutr Res. (2023) 119:1–20. doi: 10.1016/j.nutres.2023.08.010

Crossref Full Text | Google Scholar

38. Powers, HJ. Riboflavin (vitamin B-2) and health. Am J Clin Nutr. (2003) 77:1352–60. doi: 10.1093/ajcn/77.6.1352

Crossref Full Text | Google Scholar

39. Wang, L, Xu, Z, Zhang, Y, Wang, R, Wang, J, Yang, S, et al. Recent insights into function, structure, and modification of cytochrome P450 153 a family. Mol Biol Rep. (2023) 50:6955–61. doi: 10.1007/s11033-023-08553-8

PubMed Abstract | Crossref Full Text | Google Scholar

40. Pedrini, M, Cao, B, Nani, JVS, Cerqueira, RO, Mansur, RB, Tasic, L, et al. Advances and challenges in development of precision psychiatry through clinical metabolomics on mood and psychotic disorders. Prog Neuro-Psychopharmacol Biol Psychiatry. (2019) 93:182–8. doi: 10.1016/j.pnpbp.2019.03.010

Crossref Full Text | Google Scholar

41. Quinones, MP, and Kaddurah-Daouk, R. Metabolomics tools for identifying biomarkers for neuropsychiatric diseases. Neurobiol Dis. (2009) 35:165–76. doi: 10.1016/j.nbd.2009.02.019

PubMed Abstract | Crossref Full Text | Google Scholar