Next Article in Journal
A Bayesian Network Meta-Analysis and Systematic Review of Guidance Techniques in Botulinum Toxin Injections and Their Hierarchy in the Treatment of Limb Spasticity
Previous Article in Journal
Metabolomics Reveals Strain-Specific Cyanopeptide Profiles and Their Production Dynamics in Microcystis aeruginosa and M. flos-aquae
Previous Article in Special Issue
Fumonisin B Series Mycotoxins’ Dose Dependent Effects on the Porcine Hepatic and Pulmonary Phospholipidome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Metabolomics and Lipidomics Analysis Reveals Lipid Metabolic Disorder in NCM460 Cells Caused by Aflatoxin B1 and Aflatoxin M1 Alone and in Combination

1
Key Laboratory of Quality & Safety Control for Milk and Dairy Products of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
Laboratory of Quality and Safety Risk Assessment for Dairy Products of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
3
Milk and Milk Products Inspection Center of Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
4
State Key Laboratory of Animal Nutrition, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China
5
Research and Development Institute, Heilongjiang Feihe Dairy Co., Ltd., Qiqihar 161000, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxins 2023, 15(4), 255; https://doi.org/10.3390/toxins15040255
Submission received: 16 February 2023 / Revised: 14 March 2023 / Accepted: 27 March 2023 / Published: 31 March 2023
(This article belongs to the Special Issue Toxicology Research on Mycotoxins)

Abstract

:
Aflatoxin B1 (AFB1) and aflatoxin M1 (AFM1) are universally found as environmental pollutants. AFB1 and AFM1 are group 1 human carcinogens. Previous sufficient toxicological data show that they pose a health risk. The intestine is vital for resistance to foreign pollutants. The enterotoxic mechanisms of AFB1 and AFM1 have not been clarified at the metabolism levels. In the present study, cytotoxicity evaluations of AFB1 and AFM1 were conducted in NCM 460 cells by obtaining their half-maximal inhibitory concentration (IC50). The toxic effects of 2.5 μM AFB1 and AFM1 were determined by comprehensive metabolomics and lipidomics analyses on NCM460 cells. A combination of AFB1 and AFM1 induced more extensive metabolic disturbances in NCM460 cells than either aflatoxin alone. AFB1 exerted a greater effect in the combination group. Metabolomics pathway analysis showed that glycerophospholipid metabolism, fatty acid degradation, and propanoate metabolism were dominant pathways that were interfered with by AFB1, AFM1, and AFB1+AFM1. Those results suggest that attention should be paid to lipid metabolism after AFB1 and AFM1 exposure. Further, lipidomics was used to explore the fluctuation of AFB1 and AFM1 in lipid metabolism. The 34 specific lipids that were differentially induced by AFB1 were mainly attributed to 14 species, of which cardiolipin (CL) and triacylglycerol (TAG) accounted for 41%. AFM1 mainly affected CL and phosphatidylglycerol, approximately 70% based on 11 specific lipids, while 30 specific lipids were found in AFB1+AFM1, mainly reflected in TAG up to 77%. This research found for the first time that the lipid metabolism disorder caused by AFB1 and AFM1 was one of the main causes contributing to enterotoxicity, which could provide new insights into the toxic mechanisms of AFB1 and AFM1 in animals and humans.
Key Contribution: 1. Enterotoxicity effects of AFB1 and AFM1 on metabolism level were firstly investigated by metabolomics and lipidomics analysis simultaneously. 2. Metabolomic analysis revealed that both AFB1 and AFM1 exerted enterotoxicity mainly affecting the lipid metabolism. 3. AFB1 induced the enormous perturbations in metabolism level and AFB1 might be the one that wields dominant toxicity compared with AFM1, it also proved by the IC50 value of AFB1 and AFM1 on NCM460 cells.

Graphical Abstract

1. Introduction

Aflatoxins are mainly secondary metabolites produced by Aspergilus flavus and Aspergillus parasiticus that induce toxic or carcinogenic effects. The main direct and indirect routes of human and animal exposure are ingestion of contaminated food or contact with environment aflatoxins [1]. Studies have shown that acute aflatoxin poisoning can cause abdominal pain, vomiting, and even death [2]. Another study by Akinrinmade et al. [3] pointed to AFB1-induced intestinal injuries in rats. Infants or children are negatively affected by aflatoxins, including nutrient maldigestion and malabsorption, chronic immune activation, impaired bone growth, and remodeling gastrointestinal diseases [4,5]. Among these toxins, Aflatoxin B1 (AFB1) is considered the most recurrent and also the most harmful. AFB1 often occurs in agricultural production and cereal crops, as well as in the soil ecosystem and environment [6,7]. A study of animal feed from China during 2018–2020 showed that AFB1 was present in 81.9% of the samples [8]. As a result of their stable physical and chemical properties, aflatoxins are currently recognized as some of the most toxic carcinogens. AFB1 has genotoxic, carcinogenic, teratogenic, and mutagenic activities [9,10,11,12]. Previous toxicological data indicated that AFB1 are group 1 human carcinogens [13]. AFB1 has been the research focus of most studies since its discovery, mainly because of the intermediate metabolite aflatoxin AFB1-exo-8,9 epoxide. As a highly unstable molecule, AFB1-exo-8,9 epoxide reacts with cellular macromolecules to induce metabolic disorders [14]. Another toxin causing great concern is Aflatoxin M1 (AFM1); after ingestion of AFB1 by cows, aflatoxin M1 (AFM1) remained in milk through hydroxylation (hence the M in the name) [1]. AFM1 was verified to have higher enteric carcinogenicity owning to its high polarity [15]. It is frequently found in human and animal milk [16,17] and infant formula [18,19]. Thus, it is possible for humans to come in contact with AFB1 and AFM1 through their daily dietary patterns of cereal and milk. Compared with single aflatoxins, combinations may cause more serious effects, posing a greater risk to animal and human health [20,21]. In particular, AFB1 and AFM1 are known to be closely associated with animal and human health, but the interaction mechanism is not clear, especially from a small molecular level.
The intestines are the first barrier to human contact with foodborne pollutants, and they play critical roles in maintaining human and animal health. The intestines have become one of the important targets of aflatoxins. AFB1 is absorbed rapidly from the small intestines into the mesenteric venous blood. However, many studies have focused on the mechanism of renal and liver toxicity [22,23,24], and there have been few studies of the combined enterotoxicity of AFB1 and AFM1. Our previous research revealed the additive and synergistic enterotoxic effects of a combination of AFB1 and AFM1 in vivo and in vitro [25,26]. Apart from our research, few studies have investigated the enterotoxicity of a combination of AFB1 and AFM1 and its underlying mechanisms, let alone from metabolite levels, especially for lipids metabolites to decode the mechanism. Intestinal metabolism is of considerable importance to health and disease [27,28]. Unlike other nutrients, the absorption way of lipids has its own distinguishing features [29]. The normal function of intestinal lipid metabolism is vital to ensure an adequate energy supply for the various organs in the body. The immune toxic effect of AFB1 mainly by reducing the intestinal production of metabolites, including short-chain fatty acids, deoxycholic bile acid, and indole [30]. 3,4,4’-trichloro-carbanilide exerted enterotoxicity mainly by disrupting intestinal bile acid and amino acid metabolism in mice [31]. In addition, a novel probiotic CP9 was proved to attenuate the enteric infection in IPEC-J2 cells induced by Escherichia coli through metabolic modulation [32]. Research of the enterotoxicity of AFB1 and AFM1 using metabolomics and lipidomics is needed. NCM460 is a normal human colonic mucosal epithelial cell line and has become an important in vitro experimental model for intestinal research owning to its potential characteristics before carcinogenesis when the culture is in vitro [33,34]. Hence, this study aimed to explore the enterotoxicity process of AFB1 and AFM1 in NCM460 cells at the metabolic level.
Metabolomics is a common method to reveal corresponding changes in metabolism derived from external stimuli [35]. It aims to monitor changes in metabolites so as to provide critical targets to study mechanisms at the metabolic level. Previous research has used metabolomics to explore the mechanisms of aflatoxicosis, with a focus on the kidneys and liver [36,37,38]. Lipidomics is a rapidly growing method derived from metabolomics, attempting to construct an all-embracing cellular liposome with lipid biology, technology, and medicine [39]. Recent research has shown that lipidomics can precisely identify markers of aflatoxin cytotoxicity [40]. Metabolomics and lipidomics are becoming emerging approaches for the mechanistic investigation of the targets of pollutants [41,42]. They have also been utilized in numerous research fields, including toxicology, drug development, and environmental science [41,43]. The relationship between polybrominated diphenyl ether-47 and breast cancer in mice was unveiled by using metabolomics and lipidomics [42]. These two omics methods were also performed in Neuro-2a cells to investigate the toxic effect of neonicotinoids imidacloprid and acetamiprid [41]. A combination of both omics methods could provide mechanistic clues about the biological processes of toxic substances, which could help us identify toxic targets more quickly. In addition, it can cover the majority of small molecules and thus offer a biomarker of cellular damage after toxicant exposure. So, it is important to study aflatoxins at a small molecular level.
However, no research has simultaneously revealed perturbations of AFB1 and AFM1 exposure on metabolites and lipids in cells. In the present study, we used metabolomics and lipidomics to investigate the changes in metabolites and lipids induced by AFB1 and AFM1 in NCM460 cells.

2. Results

2.1. Effect of AFB1 and AFM1 Treatment on NCM460 Cells Viability

The toxic effects of AFB1 and AFM1 (1.25, 2.5, 5, 10, 15 or 20 μM) on NCM460 cells were examined. Aflatoxins reduced cell viability in a dose−dependent manner (Figure 1). The half−inhibitory concentration (IC50) was calculated according to the dose−response of nonlinear regression. The IC50 of AFB1 was 7.287 μM, while that of AFM1 was 9.720 μM, whereas that of the combination IC50 was 5.813 μM. AFB1 induced stronger cytotoxicity than AFM1, and AFB1 + AFM1 was the most toxic for NCM460 cells. The IC10 and IC20 of AFB1 were 2.218 and 2.632 μM, respectively, compared with 2.193 and 2.675 μM for AFM1 based on the dose−response curve. According to the ISO standard (2009), a substance that reduces cell viability below 70% is regarded as cytotoxic after incubation with pre-cultured cells for at least 24 h. Among the concentrations above, 2.5 μM was the maximum for which the inhibition rate was no more than 30% for both AFB1 and AFM1. Thus, AFB1 and AFM1 alone and in combination at 2.5 μM were used for further metabolomics and lipidomics analysis.

2.2. Multivariate Analysis of NCM460 Cell Metabolic Profiles

Orthogonal projection to latent structure−discriminant analysis (OPLS−DA) was used to examine sample variance and distribution of metabolomics data within and between groups. With an orthogonal signal correction filter, OPLS−DA had greater discriminant ability compared with principal component analysis (PCA). It showed no significant deviation in all samples, therefore, no sample was discarded in subsequent analyses (Figure 2). The ellipse showed 95% of the distribution range, which indicated that the four groups of samples had no obvious outliers. Additionally, all aflatoxin groups were completely separated from the control group (Figure 2A−C), indicating the different expression profiles of metabolites between the control and treated cells. This implies that AFB1 and AFM1 induced metabolic disorders in NCM460 cells.

2.3. Differential Metabolites in NCM460 Cells Induced by AFB1 and AFM1

Differential metabolites were screened out by different pairwise comparisons, according to Dunn’s test (p < 0.05 and FC > 1.5). Figure 3A compared the results between the aflatoxin and control groups. Compared with the control group, AFB1 had 15 upregulated and 19 downregulated metabolites, while AFM1 had 14 upregulated and 18 downregulated metabolites. In the AFB1 + AFM1 group, there were 23 upregulated and 23 downregulated metabolites. Venn analysis was performed to find the specific and common target metabolites of single and combined aflatoxins. There were 10, 6, and 17 specific metabolites in the AFB1, AFM1, and AFB1 + AFM1 groups, respectively (Figure 3B). Boxplots depicted the level of these metabolites in the different groups (Figure S1). In addition, 19 common metabolites were found among the aflatoxin groups.
To distinguish the differences and similarities in the toxic effects of different aflatoxins, 52 significantly changed metabolites were used for heatmap analysis: 33 specific (10 for AFB1, 6 for AFM1, and 17 for AFB1 + AFM1) and 19 common metabolites. AFB1 affected the metabolites in orange font, such as choline, lysophosphatidylcholine (LysoPC) (16:1), LysoPC (18:1), and LysoPC (20:4) (Figure 3C). AFM1 affected the metabolites in pale blue/green font, including IMP and glutaric acid. AFB1 + AFM1 affected the metabolites in sky blue font, such as upregulated taurine and downregulated succinic acid. The metabolites in black font were commonly affected by AFB1, AFM1, and AFB1 + AFM1, including 11 downregulated and 8 upregulated metabolites.

2.4. Metabolic Pathway Analysis Affected by AFB1 and AFM1

Metabolite set enrichment analysis (MSEA) and metabolic pathway analysis (MetPA) were used to focus on the pathways and metabolite sets after aflatoxins treatment. Pathway analysis revealed that the metabolism of purine (p = 0.005) and glycerophospholipid (p = 0.014) were the pathways most significantly affected by AFB1 when compared with the control group (Figure 4A). AFB1+AFM1 treatment significantly affected the glycerophospholipid metabolism pathway of NCM460 cells (p = 0.043) (Figure 4B). For AFM1, there was no significantly affected pathway, but MSEA showed that AFM1 mainly affected fatty acid metabolites, such as propanoate, and β oxidation of very long−chain fatty acids (Figure S2). This may be because AFB1 and AFM1 exerted their toxicity in different ways. Glycerophospholipid metabolism was identified as the key pathway for AFB1 and AFB1 + AFM1, while AFM1 mainly disrupted fatty acid metabolism. AFB1 might be the predominant toxin when combined with AFM1. Thus, the different structural compositions of AFB1 and AFM1 may result in different enterotoxicity for NCM460 cells.
To explore further the different toxic effects of AFB1 and AFM1, KEGG analysis (https://www.omicshare.com/tools/, 21 July 2022) of significantly changing specific metabolites in different groups was carried out. Glycerophospholipid metabolism was the main pathway affected by AFB1 (Figure 5A). AFM1 induced changes in fatty acid degradation (Figure 5B), and AFB1 + AFM1 affected propanoate and sulfur metabolism (Figure 5C). Glycerophospholipid metabolism was affected by AFB1, while AFM1 was inclined to affect fatty acid degradation. AFB1 + AFM1 interfered with more metabolic pathways, and the most affected was propanoate metabolism. Based on the results of metabolomic pathway analysis, aflatoxins mainly affect lipid metabolism. To elucidate further the specific changes in lipids metabolism, lipidomics analysis was performed.

2.5. Lipid Profiles on NCM460 Induced by AFB1 and AFM1

To assess the changes in lipid metabolites, lipid profiles in NCM460 cells were determined after aflatoxin treatment. In total, 27 lipid subclasses were detected, including 441 lipid molecules. Triacylglycerol (TAG) was the most frequently detected, comprising 128 lipid molecules and accounting for 29% of all lipids, followed by cardiolipin (CL), diacylglycerol (DAG), and sphingomyelin (SM) (Figure 6A). PCA of lipids demonstrated clear separation among the three treatment groups and control group, which suggested significant differences in the effects of aflatoxins on lipid metabolism in NCM460 cells (Figure 6B). This also proved that the metabolic effect of AFB1 and AFM1 was mainly through changes in lipid metabolism.

2.6. Analysis of Differential Lipids in NCM460 Cells

According to the Dunn’s test p < 0.05 and FC > 1.5, differentially expressed lipids were filtrated by pairwise comparison of control and aflatoxin−treated groups. Compared with the control group, there were 63 (22 upregulated and 41 downregulated), 23 (6 upregulated and 17 downregulated), and 52 (40 upregulated and 12 downregulated) lipids that satisfied the screening conditions in the AFB1, AFM1 and AFB1 + AFM1 groups (Figure 7A). The changes in lipid metabolism on NCM460 cells induced by AFB1 were more intense than those induced by AFM1. Venn analysis was performed to find the specific differentially expressed lipids between the cells treated with AFB1 and AFM1 alone and in combination. There were only five lipids identified in common among the AFB1, AFM1 and AFB1 + AFM1 groups, which showed that AFB1 and AFM1 had different effects on lipid profile (Figure 7B). There were 34, 11, and 30 specific lipids identified in the AFB1, AFM1 and AFB1 + AFM1 groups, respectively.
Heatmap analysis was used to display the expression profile of the specific lipids identified above. For AFB1, 34 specific lipids were identified, 27 lipids were decreased, and 12 of these, including CL, TAG, phosphatidyl ethanolamine (PE) and SM, accounted for 63% of the total downregulated lipids (Figure 7C). Only seven metabolites were increased by AFB1: TAG46:2 (18:2), TAG54:6 (18:2), TAG60:3 (18:2), TAG 58:3 (24:0), lysophosphatidic acid (LBPA)38:6 and LBPA36:4, and phosphatidylcholine (PC) 36:4. AFM1 tended to exert its toxicity effects through reducing free fatty acids (FFAs), phosphatidylinositol, CL, cholesteryl and phosphatidylglycerol (PG) (Figure 7D), and CL and PG accounted for 70% of the total. Only the level of TAG48:0 (16:0) was increased by AFM1. AFB1 + AFM1 significantly increased all the TAGs (Figure 7E), and PC40:7, DAG32:0 (16:0/16:0), PC36:3p, LysoPC20:4 and PE42:2p were decreased. In short, AFB1 exerted its toxic effects on NCM460 cells through changing the level of CL, TAG, PE, and SM; AFM1 affected CL and PG; while increased metabolic level of TAG played a critical role in AFB1+AFM1−induced enterotoxicity.

3. Discussion

Aflatoxins are regarded as a health issue of increasing international concern, and information on potential ways to mitigate the negative effects of AFB1 and AFM1 exposure remains limited. Ammoniation and adsorption on clays or organic adsorbents are commonly used in feed [44,45]. Recently, studies have paid attention to the AFM1 binding ability of different strains of Lactobacilli around AFM1 in milk [46,47]. The toxicity of AFB1 and AFM1 has been determined in various organs and tissues in vivo and in vitro, and the intestines are key target organs. Aflatoxins have been shown to impair intestinal mRNA and protein levels [26,48], but few studies focused on metabolism and lipid levels. Therefore, the present study investigated the metabolic effects and underlying toxic mechanisms of AFB1 and AFM1 on NCM460 cells, which could provide a theoretical basis for the risk assessment and treatment of enterotoxicity.
The viability of NCM460 cells was reduced by AFB1 and AFM1 in a dose-dependent manner. The order of severity of toxic effects on NCM460 cells was AFB1 + AFM1 > AFB1 > AFM1, as evaluated by IC50 (Figure 1). Only one study to assess IC50 value on the intestinal cell was found, pointing out that the IC50 of AFB1 on Caco-2 cells was 39.1 ± 13.2 μM, while AFB1 had so few cytotoxic effects on the two HCT116 cell lines that IC50 could not be calculated [49]. The IC50 value on the Caco-2 cell was higher than our results, which could be explained by the different cell lines. The cancer cell line Caco-2 tolerance of aflatoxins is stronger than normal intestinal cell line NCM460. This article also confirmed that AFB1 has different toxicity to different intestinal cells, and there is also a great difference in IC50 value. Regardless of the number of specific differential metabolites and enrichment pathways, the effects of AFB1 + AFM1 were stronger than those of AFB1 and AFM1 alone, which was consistent with the phenotypic tests.
Metabolomics analysis pointed out that glycerophospholipid was significantly increased by AFB1. Glycerophospholipids are an essential part of the lipid bilayer and have an important influence on membrane properties and dynamics [50,51]. It also alleviated ulcerative colitis in mice [52]. Thus, the significant changes in glycerophospholipids induced by AFB1 may influence the overall metabolism of NCM460 cells, which infers that AFB1 causes intestinal lesions and even cancer by disrupting glycerophospholipid metabolism [37]. The decreased level of LysoPC (16:1), LysoPC (18:1), and LysoPC (20:4) and increased choline in the AFB1 treatment group indicates that AFB1 can damage the intestinal barrier and immune function, which may even develop into cancer [53,54]. Some studies confirmed that AFB1 exerted a negative effect on intestinal cancer cells [55]. In addition, both specific metabolites LysoPC and xanthosine participated in intestinal inflammation in the AFB1 group [56,57]. Changes in the level of these metabolites are also reported to be related to intestinal inflammation and mucosal homeostasis [58]. Xanthosine was also found in the liver [59], but the opposite trend was presented after AFB1 treatment. Higher levels of choline induced cell swelling because choline, phosphocholine, and PC are essential elements for the structural integrity of cell membranes [60]. Choline was increased by AFB1 in this study, which has become the endogenous biomarker of cancer [53]. Thus, AFB1 exerted enteric toxicity through disruption of intestinal homeostasis, reflected by the metabolic disorder of glycerophospholipids. The specific influence of AFM1 was reflected in fatty acid degradation, and glutaric acid was significantly enrichment in this pathway (Figure 5B). Recent research identified that glutaric acid may be involved in regulating intestinal microbiota [61], with increased glutaric acid levels exerting a negative effect on the intestine. Therefore, the influence of fatty acid metabolism may be the main cause of AFM1-induced enterotoxicity. Propanoate metabolism and sulfur metabolism were specifically enriched in the AFB1 + AFM1 group (Figure 5C), and the downregulation of succinic acid and upregulation of taurine were identified in these two pathways. Prior research has demonstrated that succinic acid plays a critical role in regulating intestinal immune responses to external stimulation [62,63]. AFB1 + AFM1 increased the taurine level, which may induce intestinal proinflammatory and oxidative effects [64,65]. Taurine regulates the production of short-chain fatty acids and the occurrence of colitis [66,67]. Therefore, increased taurine probably promotes intestinal inflammatory response. The effects on glycerophospholipid metabolism, fatty acid degradation, and propanoate metabolism all reflect that AFB1 and AFM1 affect lipid metabolism.
Lipid metabolites were identified as the pivotal type of differential metabolites; thus, it is necessary to define the role of lipid metabolism in intestinal injury induced by AFB1 and AFM1. Aflatoxins are highly lipophilic [68], and lipid metabolic disorders are a common hallmark of aflatoxicosis [59,69]. Each part of the intestine has its own functions and maintains normal metabolism of lipids. Dietary lipids undergo digestion, emulsification, hydrolysis, and micellization in the duodenum. The jejunum promotes lipid assimilation, and the ileum provides the place for this process and enterohepatic circulation [70]. Disorders of lipid metabolism occur when aflatoxins decrease the length of the small intestine induced [71]. A lipidomics study identified 48 lipids among the DAGs, TAGs, FFAs, and phospholipids that have metabolic dysregulation upon exposure to mycotoxin [72]. To understand further the mechanism of aflatoxin-induced enteric toxicity at the metabolic level, we performed lipidomics analysis.
PCA of lipidomics analysis proved the different focus of AFB1 and AFM1 on lipid metabolism (Figure 6B). In addition, the results showed that AFB1 evoked greater changes than AFM1 in lipid metabolism. AFB1 exerted enterotoxic effects on NCM460 cells mainly by changing the levels of CL, TAG, PE, and SM (Figure 7C). The changed level of CL was reported to affect mitochondrial function and inflammatory response [73]. The synthesis of SM participated in cell normal physiology, and dysregulation of SM level may affect the tumor transformation [74]. PC and DAG were also found in the AFB1 group. Synthesis of PC is necessary for lipid absorption in dietary and metabolic homeostasis in the intestines [75]. Fecal lipidomics in F344 rats revealed that AFB1 induced intestinal epithelial toxicity through changes in the level of PE and DAG [55]. The decrease in PC32:1, PE, and DAG may be accompanied by inflammation. Previous research showed that AFB1 disrupted intestinal homeostasis, which was marked by a reduction in TAG and cholesterol [76]. Another study pointed out that the plasma and liver levels of TAG in rats were increased after acute exposure to AFB1 [69]. The reason for the inconsistent results may be that there are many different types of TAG. In the present study, AFB1 simultaneously increased and decreased the different TAG molecules. The specific mechanism of changes at different molecular levels is worth further investigation.
The enterotoxic effects of AFM1 were mainly mediated through decreasing levels of CL and PG (Figure 7D). Different CL species have been identified in cells treated with AFB1 or AFM1, which can block Toll-like receptor 4 response to lipopolysaccharides and maintain intestinal homeostasis [77]. In contrast, decreased CL level may induce intestinal homeostasis disorder. PG was reported to reduce the production of tumor necrosis factor-α and exert anti-inflammatory effects in vitro [78]. The decreased level of PG in our study indicated that AFM1 induced a proinflammatory response in NCM460 cells.
The accumulation of TAG accompanied by inhibition of intestinal fat absorption causes severe adiposity [79]. Ugbaja et al. (2020) reported that the synthesis of TAG occurred via complex pathways that involved PC or DAG or both membrane lipids [80]. Thus, the increase in TAG by AFB1 + AFM1 meant that not only were cellular lipids altered, but also priority interfered with metabolic pathways of lipids even involved in other cellular. Therefore, AFB1 + AFM1 exerted enterotoxicity by raising the level of TAG metabolism. Targeting lipid metabolism or signaling is becoming a trend for disease therapeutics [81]. Regulation of lipid metabolism alters the antitumor function of T cells in cancer mouse models [82,83]. The regulation of lipid metabolic pathways can alter the production and function of Treg cells, which may also be applied to the treatment of autoimmune diseases [84].
How AFB1 and AFM1 altered the level of lipids in NCM460 cells remains unclear. We propose a hypothesis that AFB1 and AFM1 exert their enterotoxicity through lipid dysmetabolism. This needs to be confirmed in future studies.

4. Conclusions

Enterotoxic effects of AFB1 and AFM1 alone and in combination were investigated by simultaneous metabolomics and lipidomics analysis. Metabolomics analysis revealed AFB1 and AFM1, both individually and in combination, disrupt lipid metabolism in different specific pathways: glycerophospholipid metabolism was the dominant interrupted by AFB1; fatty acid degradation was disrupted by AFM1; propanoate metabolism was the main interrupted in AFB1 + AFB1. A combination of AFB1 and AFM1 caused more extensive metabolic disturbances than either aflatoxin alone. AFB1 and AFM1 played an enterotoxic role by disrupting intestinal lipid metabolism. Further, lipidomics analysis revealed that 75 specific lipids with significant changes were screened after AFB1 and AFM1 treatment. Of all the differential lipids, TAGs were the most frequently observed, with a proportion of 43%. Importantly, AFB1 caused major perturbations in lipid metabolism levels in NCM460 cells and might show greater enterotoxicity compared with AFM1. The results pointed the lipid metabolism disorder caused by AFB1 and AFM1 should be emphasized in environmental toxicological studies, which provide a novel direction for the study of aflatoxin-related health issues, as well as novel molecular targets for preventing AFB1 and AFM1 contamination. Lipid metabolism identified in this study also can be a regulatory mechanism for reducing the damage of aflatoxins and may even be for clinical disease treatment.

5. Materials and Methods

5.1. Chemicals and Reagents

AFB1 and AFM1 (95% purity) were provided by Pribolab (Qingdao, Shangdong, China), and a stock liquor of 200 μM was formulated by dissolving in dimethyl sulfoxide (DMSO) (Sigma−Aldrich, St. Louis, MO, USA). RPMI 1640 and fetal bovine serum (FBS) used for cell growth were from Gibco (Grand Island, NY, USA). The antibiotics (100 U/mL penicillin and 100 μg/mL streptomycin) and nonessential amino acids (NEAAs) were obtained from Life Technologies (Carlsbad, CA, USA). L-Glutamine was supplied by Sigma-Aldrich. 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl-2H-tetrazolium bromide (MTT) and phosphate buffer solution were provided by Beyotime Biotechnology (Shanghai, China).

5.2. Cell Culture and Treatment

The normal human colon cell line NCM460 cells were purchased from the American Type Culture Collection (Manassas, VA, USA) and cultured as previously [85]. NCM460 cells were maintained in RPMI160 with 10% (v/v) FBS, 1% antibiotics, L–glutamine, and NEAA in a 37 °C incubator in 5% CO2 (Thermo, Waltham, MA, USA). Water used in the incubation system was purified through a Milli−Q system (Millipore, Bedford, MA, USA). When NCM460 cells were seeded in 6-well and 96−well plates (Corning, Corning, NY, USA) at 2 × 105 cells/well and 5 × 104 cells/well for 24 h and challenged with AFB1/AFM1 (0−20 μM) for 48 h. The cells in 96−well and 6−well were respectively collected for subsequent cell viability, metabolomics, and lipidomics analysis.

5.3. Cell Viability Assay

The viability of NCM460 cells in 96−well plates after exposure to AFB1 and AFM1 was evaluated by using the MTT assay. After the cells grew up to a confluence of 90%, AFB1 and AFM1 were added to 96−well plates (5 × 104 cells/well) at 0, 1.25, 2.5, 5, 10, 15, and 20 μM for 48 h. A total of 5 mg/mL MTT solution was prepared by dissolving 25 mg MTT with 5 mL MTT solvent. The cells were treated with 10 μL 5 mg/mL MTT solution incubated for 4 h. Then the cells were cultured with 100 μL formazan solvent (37 °C for 3−4 h). The plate was gently shaken for 10 min to completely dissolve the crystals, and absorbance was recorded at 570 on the microplate reader (Thermo Scientific, Waltham, MA, USA). GraphPad Prism 8 was used to analyze the IC10, IC20, and IC50 according to the dose–response of nonlinear regression.

5.4. Metabolomics Analysis

Metabolome extraction was performed as described previously [86]. Each treatment group had eight duplicate samples. Cell samples were mixed with 80% ice−cold methanol in water, incubated at 1500 rpm at 4 °C for 30 min, and centrifuged at 12,000 rpm at 4 °C for 10 min. The supernatant obtained above was transferred into a clean 1.5 mL centrifuge tube. Genevac miVac (Tegent Scientific, Shanghai, China) was used to redissolve the dried extracts in water with 5% acetonitrile. Moreover, the upper liquids were collected for LC−MS analysis.
The LC−MS parameters were as described previously [86]. ACQUITY UPLC HSS T3 1.8 μm, 3.0 × 100 mm columns were used for chromatographic analysis. LC−MS experiments of Q−Exactive MS with a HESI probe controlled by Xcalibur 2.3 software were carried out on the Dionex UltiMate3000 HPLC system. The MS detection parameters were sheath gas flow rate 40 psi, auxiliary gas flow rate 11 arb, sweep gas flow rate 0, and full MS scan, and the scan range was 80−1000 m/z, capillary temperature 350 °C, spray voltage 3.5 kV and 3.2 kV, respectively, for positive and negative mode, s-lens RF level 55, and auxiliary gas heater temperature 220 °C. The Q−Orbitrap performed data−dependent scans in full MS/dd−MS2 to obtain product ion information with normalized collision energy (NCE), and the mass resolution was set at 17,500 FWHM (m/z 200) and NCE 35%. All detected ions were compared with HMDB, METLIN, and standard references to annotate ion ID. High−purity nitrogen was used as the atomization gas and the collision gas for higher energy collisional dissociation. Differential metabolites were screened based on the Dunn test p < 0.05 and fold change (FC) >1.5. Two methods were adopted for pathway analysis: metabolite set enrichment analysis (MSEA) and metabolic pathway analysis (MetPA). MSEA analysis can find out the set of metabolites that are enriched in this group of metabolites. The test method used is Over−representation analysis (ORA). MetPA has two important components: ORA and pathway topological analysis. ORA detects whether a metabolite set includes more metabolites with significant differences than would be expected at random. ORA tests can be either a hypergeometric test or Fisher’s exact test.

5.5. Lipidomic Analysis

An improved version of Bligh and Dyer’s method was used to extract lipids [87]. Each group was composed of three duplicate samples, and a total of four groups were applied in the lipidomics study. NCM460 cells were homogenized in 750 µL chloroform: methanol: MilliQ H2O (3:6:1) (v/v/v) and incubated at 1500 rpm at 4 °C for 1 h. After incubation, chloroform (250 µL) and deionized water (350 µL) were used to induce phase separation. The lower organic phase containing lipids was extracted into a clean tube after centrifugation. Chloroform (450 µL) was added to the remaining cells in the aqueous phase to repeatedly extract lipids, and the lipid extracts were pooled into one tube and dried in the OH mode of the SpeedVac. The samples were kept at −80 °C for further testing.
An ExionLC−AD combined with Sciex QTRAP 6500 PLUS was used to perform lipidomic analyses at LipidALL Technologies (Guangzhou, China) [88]. Each lipid of polar was separated by normal-phase HPLC using a TUP−HB silica column (i.d. 150 × 2.1 mm, 3 µm) under the specific conditions: mobile phase A (chloroform: methanol: ammonium hydroxide, 89.5:10:0.5) and mobile phase B (chloroform: methanol: ammonium hydroxide: water, 55:39:0.5:5.5). Multiple reaction monitoring (MRM) transitions were established to comparative analysis of various polar lipids. Individual lipid species were quantified by referencing spiked internal standards. Free fatty acids (FFA) were quantitated using FFA 19:0 (Cayman Chemicals). Glycerol lipids, including diacylglycerols (DAG) and triacylglycerols (TAG), were quantified using a modified version of reverse−phase HPLC/MRM. Separation of neutral lipids was achieved on a Phenomenex Kinetex−C18 column (i.d. 4.6 × 100 mm, 2.6 µm) using an isocratic mobile phase containing chloroform: methanol:0.1 M ammonium acetate 100:100:4 (v/v/v) at a flow rate of 300 µL for 10 min. Levels of TAG were calculated by referencing the spiked internal standard TAG (16:0) 3-d5 obtained from CDN isotopes, respectively. DAG was quantified using d5-DAG16:0/16:0 and d5-DAG18:1/18:1 as internal standards (Avanti Polar Lipids). Free cholesterols and cholesteryl esters (CEs) were analyzed under atmospheric pressure chemical ionization mode on an Agilent 1260 HPLC connected to Sciex QTTRAP 5500, using d6-cholesterol and d6-C18:0 CEs (CDN isotopes) as internal standards [89].
Principal component analysis (PCA) was used to reduce dimensions that identify a small new set of independent variables that best represent the information of the original data so that the transformed data were easier to visualize and interpret. FactoMineR package was used for PCA. Plots of differential metabolites, boxplots, and barplot are all based on hypothetical tests. We used the ggplot2 package for these plots. Original concentrations of individual lipids (MFP) were used in Dunn’s test. Concentrations of individual lipids (MFP) were standardized to z scores for the heatmap plot. All statistical analyses were performed using R 4.0 software. Two-sided p < 0.05 was considered statistically significant. Dunn’s t-test with p < 0.05 and FC > 1.5 were identified as differential lipids.

5.6. Statistical Analysis

Statistical analysis of changes in cell viability and inhibition after AFB1 and AFM1 treatment compared with the control cells was performed by GraphPad Prism 8.0 (San Diego, CA, USA). Analysis of variance combined with Tukey’s multiple analysis was used to analyze the difference after confirming normality by Shapiro−Wilk tests. p < 0.05 represented statistically significant differences.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxins15040255/s1, Figure S1: Boxplot further exhibited the level of each specific metabolites in different groups. Different letters (ac) indicate significant differences (p < 0.05). Figure S2: he metabolites sets enrichment analysis of differential expression metabolites in AFM1/CTL.

Author Contributions

Conceptualization, X.Y.; Methodology, Y.G.; Software, Y.G.; Validation, X.L.; Formal analysis, X.L.; Investigation, X.Y.; Data curation, X.L.; Writing—original draft, X.Y.; Writing—review and editing, J.W.; Supervision, N.Z.; Project administration, J.W.; Funding acquisition, N.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the China Agriculture Research System of MOF and MARA (CARS36), the National Key R&D Program of China (22022YFD1301003), and the Agricultural Science and Technology Innovation Program (ASTIP-IAS12).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available upon request; please contact the contributing authors.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Marchese, S.; Polo, A.; Ariano, A.; Velotto, S.; Costantini, S.; Severino, L. Aflatoxin B1 and M1: Biological Properties and Their Involvement in Cancer Development. Toxins 2018, 10, 214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Mohd-Redzwan, S.; Jamaluddin, R.; Mutalib, A.; Sokhini, M.; Ahmad, Z. A mini review on aflatoxin exposure in Malaysia: Past, present and future. Front. Microbiol. 2013, 4, 334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Akinrinmade, F.J.; Akinrinde, A.S.; Amid, A. Changes in serum cytokine levels, hepatic and intestinal morphology in aflatoxin B1-induced injury: Modulatory roles of melatonin and flavonoid-rich fractions from Chromolena odorata. Mycotoxin Res. 2016, 32, 53–60. [Google Scholar] [CrossRef]
  4. Tessema, M.; De Groote, H.; Brouwer, I.D.; De Boevre, M.; Corominas, A.V.; Stoecker, B.J.; Stoecker, B.J.; Feskens, E.J.; Belachew, T.; Karakitsou, A.; et al. Exposure to aflatoxins and fumonisins and linear growth of children in rural Ethiopia: A longitudinal study. Public Health Nutr. 2021, 24, 3662–3673. [Google Scholar] [CrossRef] [PubMed]
  5. Turner, P.C.; Moore, S.E.; Hall, A.J.; Prentice, A.M.; Wild, C.P. Modification of immune function through exposure to dietary aflatoxin in Gambian children. Environ. Health Perspect. 2003, 111, 217–220. [Google Scholar] [CrossRef] [Green Version]
  6. Fouche, T.; Claassens, S.; Maboeta, M. Aflatoxins in the soil ecosystem: An overview of its occurrence, fate, effects and future perspectives. Mycotoxin Res. 2020, 36, 303–309. [Google Scholar] [CrossRef]
  7. Taniwaki, M.H.; Pitt, J.I.; Magan, N. Aspergillus species and mycotoxins: Occurrence and importance in major food commodities. Curr. Opin. Food Sci. 2018, 23, 38–43. [Google Scholar] [CrossRef] [Green Version]
  8. Zhao, L.; Zhang, L.; Xu, Z.; Liu, X.; Chen, L.; Dai, J.; Karrow, N.A.; Sun, L. Occurrence of Aflatoxin B1, deoxynivalenol and zearalenone in feeds in China during 2018–2020. J. Anim. Sci. Biotechnol. 2021, 12, 74. [Google Scholar] [CrossRef]
  9. Benkerroum, N. Retrospective and Prospective Look at Aflatoxin Research and Development from a Practical Standpoint. Int. J. Environ. Res. Public Health 2019, 16, 3633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Raiola, A.; Tenore, G.C.; Manyes, L.; Meca, G.; Ritieni, A. Risk analysis of main mycotoxins occurring in food for children: An overview. Food Chem. Toxicol. 2015, 84, 169–180. [Google Scholar] [CrossRef]
  11. Chen, X.; Li, C.; Chen, Y.; Ni, C.; Chen, X.; Zhang, L.; Xu, X.; Chen, M.; Ma, X.; Zhan, H. Aflatoxin B1 impairs leydig cells through inhibiting AMPK/mTOR-mediated autophagy flux pathway. Chemosphere 2019, 233, 261–272. [Google Scholar] [CrossRef] [PubMed]
  12. Damiano, S.; Jarriyawattanachaikul, W.; Girolami, F.; Longobardi, C.; Nebbia, C.; Andretta, E.; Lauritano, C.; Dabbou, S.; Avantaggiato, G.; Schiavone, A.; et al. Curcumin Supplementation Protects Broiler Chickens Against the Renal Oxidative Stress Induced by the Dietary Exposure to Low Levels of Aflatoxin B1. Front. Vet. Sci. 2022, 8, 822227. [Google Scholar] [CrossRef] [PubMed]
  13. Ostry, V.; Malir, F.; Toman, J.; Grosse, Y. Mycotoxins as human carcinogens-the IARC Monographs classification. Mycotoxin Res. 2017, 33, 65–73. [Google Scholar] [CrossRef] [PubMed]
  14. Rushing, B.R.; Selim, M.I. Structure and Oxidation of Pyrrole Adducts Formed between Aflatoxin B(2a) and Biological Amines. Chem. Res. Toxicol. 2017, 30, 1275–1285. [Google Scholar] [CrossRef]
  15. Cullen, J.M.; Ruebner, B.H.; Hsieh, L.S.; Hyde, D.M.; Hsieh, D.P. Carcinogenicity of dietary aflatoxin M1 in male Fischer rats compared to aflatoxin B1. Cancer Res. 1987, 47, 1913–1917. [Google Scholar]
  16. Jafari, T.; Fallah, A.A.; Kheiri, S.; Fadaei, A.; Amini, S.A. Aflatoxin M1 in human breast milk in Shahrekord, Iran and association with dietary factors. Food Addit. Contam. Part B Surveil. 2017, 10, 128–136. [Google Scholar] [CrossRef] [PubMed]
  17. Atanda, O.; Oguntubo, A.; Adejumo, O.; Ikeorah, J.; Akpan, I. Aflatoxin M1 contamination of milk and ice cream in Abeokuta and Odeda local governments of Ogun State, Nigeria. Chemosphere 2007, 68, 1455–1458. [Google Scholar] [CrossRef]
  18. Agriopoulou, S.; Stamatelopoulou, E.; Varzakas, T. Advances in Occurrence, Importance, and Mycotoxin Control Strategies: Prevention and Detoxification in Foods. Foods 2020, 9, 137. [Google Scholar] [CrossRef]
  19. Awaisheh, S.S.; Rahahleh, R.J.; Algroom, R.M.; Al-Bakheit, A.A.; Al-Khaza’leh, J.M. Al-Dababseh BA: Contamination level and exposure assessment to Aflatoxin M1 in Jordanian infant milk formulas. Ital. J. Food Saf. 2019, 8, 8263. [Google Scholar]
  20. Dey, D.K.; Kang, J.I.; Bajpai, V.K.; Kim, K.; Lee, H.; Sonwal, S.; Simal-Gandara, J.; Xiao, J.; Ali, S.; Huh, Y.S.; et al. Mycotoxins in food and feed: Toxicity, preventive challenges, and advanced detection techniques for associated diseases. Crit. Rev. Food Sci. Nutr. 2022, 1–22, ahead of print. [Google Scholar] [CrossRef]
  21. Streit, E.; Schatzmayr, G.; Tassis, P.; Tzika, E.; Marin, D.; Taranu, I.; Tabuc, C.; Nicolau, A.; Aprodu, I.; Puel, O. Current situation of mycotoxin contamination and co-occurrence in animal feed--focus on Europe. Toxins 2012, 4, 788–809. [Google Scholar] [CrossRef] [Green Version]
  22. Dlamini, N.Z.; Somboro, A.M.; Amoako, D.G.; Arhin, I.; Khumalo, H.M.; Khan, R.B. Toxicogenicity and mechanistic pathways of aflatoxin B1 induced renal injury. Environ Toxicol. 2021, 36, 1857–1872. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, K.; Jia, S.; Xue, D.; Rajput, S.A.; Liu, M.; Qi, D.; Wang, S. Dual effects of zearalenone on aflatoxin B1-induced liver and mammary gland toxicity in pregnant and lactating rats. Ecotoxicol. Environ. Saf. 2022, 245, 114115. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, L.; He, L.; Zeng, H.; Fu, W.; Wang, J.; Tan, Y.; Zheng, C.; Qiu, Z.; Luo, J.; Lv, C. Low-dose microcystin-LR antagonizes aflatoxin B1 induced hepatocarcinogenesis through decreasing cytochrome P450 1A2 expression and aflatoxin B1-DNA adduct generation. Chemosphere 2020, 248, 126036. [Google Scholar] [CrossRef]
  25. Gao, Y.; Bao, X.; Meng, L.; Liu, H.; Wang, J.; Zheng, N. Aflatoxin B1 and Aflatoxin M1 Induce Compromised Intestinal Integrity through Clathrin-Mediated Endocytosis. Toxins 2021, 13, 184. [Google Scholar] [CrossRef]
  26. Gao, Y.N.; Yang, X.; Wang, J.Q.; Liu, H.M.; Zheng, N. Multi-Omics Reveal Additive Cytotoxicity Effects of Aflatoxin B1 and Aflatoxin M1 toward Intestinal NCM460 Cells. Toxins 2022, 14, 368. [Google Scholar] [CrossRef]
  27. Wu, J.; Wei, Z.; Cheng, P.; Qian, C.; Xu, F.; Yang, Y.; Wang, A.; Chen, W.; Sun, Z.; Lu, Y. Rhein modulates host purine metabolism in intestine through gut microbiota and ameliorates experimental colitis. Theranostics 2020, 10, 10665–10679. [Google Scholar] [CrossRef] [PubMed]
  28. La Vecchia, S.; Sebastian, C. Metabolic pathways regulating colorectal cancer initiation and progression. Semin. Cell Dev. Biol. 2020, 98, 63–70. [Google Scholar] [CrossRef] [PubMed]
  29. Ko, C.W.; Qu, J.; Black, D.D.; Tso, P. Regulation of intestinal lipid metabolism: Current concepts and relevance to disease. Nat. Rev. Gastroenterol. Hepatol. 2020, 17, 169–183. [Google Scholar] [CrossRef]
  30. Zhou, J.; Tang, L.; Wang, J.S. Assessment of the adverse impacts of aflatoxin B1 on gut-microbiota dependent metabolism in F344 rats. Chemosphere 2019, 217, 618–628. [Google Scholar] [CrossRef] [PubMed]
  31. Song, Y.; Zhang, C.; Lei, H.; Qin, M.; Chen, G.; Wu, F.; Chen, C.; Cao, Z.; Zhang, C.; Wu, M. Characterization of triclosan-induced hepatotoxicity and triclocarban-triggered enterotoxicity in mice by multiple omics screening. Sci. Total Environ. 2022, 838 Pt 4, 156570. [Google Scholar] [CrossRef]
  32. Sudan, S.; Zhan, X.; Li, J. A Novel Probiotic Bacillus subtilis Strain Confers Cytoprotection to Host Pig Intestinal Epithelial Cells during Enterotoxic Escherichia coli Infection. Microbiol. Spectr. 2022, 10, e0125721. [Google Scholar] [CrossRef] [PubMed]
  33. Moyer, M.P.; Manzano, L.A.; Merriman, R.L.; Stauffer, J.S.; Tanzer, L.R. NCM460, a normal human colon mucosal epithelial cell line. In Vitro Cell. Dev. Biol.-Anim. 1996, 32, 315–317. [Google Scholar] [CrossRef]
  34. Long, X.; Wong, C.C.; Tong, L.; Chu, E.S.H.; Ho Szeto, C.; Go, M.Y.Y.; Coker, O.O.; Chan, A.W.H.; Chan, F.K.L.; Sung, J.J.Y. Peptostreptococcus anaerobius promotes colorectal carcinogenesis and modulates tumour immunity. Nat. Microbiol. 2019, 4, 2319–2330. [Google Scholar] [CrossRef] [PubMed]
  35. Warth, B.; Raffeiner, P.; Granados, A.; Huan, T.; Fang, M.; Forsberg, E.M.; Benton, H.P.; Goetz, L.; Johnson, C.H.; Siuzdak, G. Metabolomics Reveals that Dietary Xenoestrogens Alter Cellular Metabolism Induced by Palbociclib/Letrozole Combination Cancer Therapy. Cell Chem. Biol. 2018, 25, 291–300.e293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wang, Z.; Gao, Y.; Huang, X.; Huang, S.; Yang, X.; Wang, J.; Zheng, N. Metabolomics analysis underlay mechanisms in the renal impairment of mice caused by combination of aflatoxin M1 and ochratoxin A. Toxicology 2021, 458, 152835. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, S.; Yang, X.; Liu, F.; Wang, X.; Zhang, X.; He, K.; Wang, H. Comprehensive Metabolomic Analysis Reveals Dynamic Metabolic Reprogramming in Hep3B Cells with Aflatoxin B1 Exposure. Toxins 2021, 13, 384. [Google Scholar] [CrossRef] [PubMed]
  38. Marchese, S.; Sorice, A.; Ariano, A.; Florio, S.; Budillon, A.; Costantini, S.; Severino, L. Evaluation of Aflatoxin M1 Effects on the Metabolomic and Cytokinomic Profiling of a Hepatoblastoma Cell Line. Toxins 2018, 10, 436. [Google Scholar] [CrossRef] [Green Version]
  39. Lam, S.M.; Tian, H.; Shui, G. Lipidomics, en route to accurate quantitation. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2017, 1862, 752–761. [Google Scholar] [CrossRef]
  40. Chen, L.; Chen, F.; Liu, T.; Feng, F.; Guo, W.; Zhang, Y.; Feng, X.; Lin, J.M.; Zhang, F. Lipidomics Profiling of HepG2 Cells and Interference by Mycotoxins Based on UPLC-TOF-IMS. Anal. Chem. 2022, 94, 6719–6727. [Google Scholar] [CrossRef]
  41. Wang, X.; Qiu, J.; Xu, Y.; Liao, G.; Jia, Q.; Pan, Y.; Wang, T.; Qian, Y. Integrated non-targeted lipidomics and metabolomics analyses for fluctuations of neonicotinoids imidacloprid and acetamiprid on Neuro-2a cells. Environ. Pollut. 2021, 284, 117327. [Google Scholar] [CrossRef]
  42. Wei, J.; Li, X.; Xiang, L.; Song, Y.; Liu, Y.; Jiang, Y.; Cai, Z. Metabolomics and lipidomics study unveils the impact of polybrominated diphenyl ether-47 on breast cancer mice. J. Hazard. Mater. 2020, 390, 121451. [Google Scholar] [CrossRef] [PubMed]
  43. Rombouts, C.; De Spiegeleer, M.; Van Meulebroek, L.; De Vos, W.H.; Vanhaecke, L. Validated comprehensive metabolomics and lipidomics analysis of colon tissue and cell lines. Anal. Chim. Acta 2019, 1066, 79–92. [Google Scholar] [CrossRef] [PubMed]
  44. Jard, G.; Liboz, T.; Mathieu, F.; Guyonvarc’h, A.; Lebrihi, A. Review of mycotoxin reduction in food and feed: From prevention in the field to detoxification by adsorption or transformation. Food Addit. Contam. Part A 2011, 28, 1590–1609. [Google Scholar] [CrossRef] [PubMed]
  45. Park, D.L.; Price, W.D. Reduction of aflatoxin hazards using ammoniation. Rev. Environ. Contam. Toxicol. 2001, 171, 139–175. [Google Scholar]
  46. Giovati, L.; Magliani, W.; Ciociola, T.; Santinoli, C.; Conti, S.; Polonelli, L. AFM1 in Milk: Physical, Biological, and Prophylactic Methods to Mitigate Contamination. Toxins 2015, 7, 4330–4349. [Google Scholar] [CrossRef] [Green Version]
  47. Jebali, R.; Abbès, S.; Salah-Abbès, J.B.; Ben Younes, R.; Haous, Z.; Oueslati, R. Ability of Lactobacillus plantarum MON03 to mitigate aflatoxins (B1 and M1) immunotoxicities in mice. J. Immunotoxicol. 2015, 12, 290–299. [Google Scholar] [CrossRef] [Green Version]
  48. Bao, X.Y.; Li, S.L.; Gao, Y.N.; Wang, J.Q.; Zheng, N. Transcriptome analysis revealed that aflatoxin M1 could cause cell cycle arrest in differentiated Caco-2 cells. Toxicol. In Vitro 2019, 59, 35–43. [Google Scholar] [CrossRef]
  49. Gauthier, T.; Duarte-Hospital, C.; Vignard, J.; Boutet-Robinet, E.; Sulyok, M.; Snini, S.P.; Alassane-Kpembi, I.; Lippi, Y.; Puel, S.; Oswald, I.P.; et al. Versicolorin A, a precursor in aflatoxins biosynthesis, is a food contaminant toxic for human intestinal cells. Environ. Int. 2020, 137, 105568. [Google Scholar] [CrossRef]
  50. Carrizo, D.; Chevallier, O.P.; Woodside, J.V.; Brennan, S.F.; Cantwell, M.M.; Cuskelly, G.; Elliott, C.T. Untargeted metabolomic analysis of human serum samples associated with exposure levels of Persistent organic pollutants indicate important perturbations in Sphingolipids and Glycerophospholipids levels. Chemosphere 2017, 168, 731–738. [Google Scholar] [CrossRef]
  51. Rodriguez-Cuenca, S.; Pellegrinelli, V.; Campbell, M.; Oresic, M.; Vidal-Puig, A. Sphingolipids and glycerophospholipids—The “ying and yang” of lipotoxicity in metabolic diseases. Prog. Lipid Res. 2017, 66, 14–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Yuan, Z.; Yang, L.; Zhang, X.; Ji, P.; Hua, Y.; Wei, Y. Mechanism of Huang-lian-Jie-du decoction and its effective fraction in alleviating acute ulcerative colitis in mice: Regulating arachidonic acid metabolism and glycerophospholipid metabolism. J. Ethnopharmacol. 2020, 259, 112872. [Google Scholar] [CrossRef]
  53. Glunde, K.; Jacobs, M.A.; Bhujwalla, Z.M. Choline metabolism in cancer: Implications for diagnosis and therapy. Expert Rev. Mol. Diagn. 2006, 6, 821–829. [Google Scholar] [CrossRef]
  54. Zhou, J.; Tang, L.; Wang, J.S. Aflatoxin B1 Induces Gut-Inflammation-Associated Fecal Lipidome Changes in F344 Rats. Toxicol. Sci. 2021, 183, 363–377. [Google Scholar] [CrossRef]
  55. Kim, J.; Park, S.-H.; Do, K.H.; Kim, D.; Moon, Y. Interference with mutagenic aflatoxin B1-induced checkpoints through antagonistic action of ochratoxin A in intestinal cancer cells: A molecular explanation on potential risk of crosstalk between carcinogens. Oncotarget 2016, 7, 39627–39639. [Google Scholar] [CrossRef] [Green Version]
  56. Chattopadhyay, A.; Navab, M.; Hough, G.; Grijalva, V.; Mukherjee, P.; Fogelman, H.R.; Hwang, L.H.; Faull, K.F.; Lusis, A.J.; Reddy, S.T.; et al. Tg6F ameliorates the increase in oxidized phospholipids in the jejunum of mice fed unsaturated LysoPC or WD. J. Lipid Res. 2016, 57, 832–847. [Google Scholar] [CrossRef] [Green Version]
  57. Gao, J.; Xu, K.; Liu, H.; Liu, G.; Bai, M.; Peng, C.; Li, T.; Yin, Y. Impact of the Gut Microbiota on Intestinal Immunity Mediated by Tryptophan Metabolism. Front. Cell. Infect. Microbiol. 2018, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  58. Ju, T.; Kennelly, J.P.; Jacobs, R.L.; Willing, B.P. Insufficient dietary choline aggravates disease severity in a mouse model of Citrobacter rodentium-induced colitis. Br. J. Nutr. 2021, 125, 50–61. [Google Scholar] [CrossRef]
  59. Lu, X.; Hu, B.; Shao, L.; Tian, Y.; Jin, T.; Jin, Y.; Ji, S.; Fan, X. Integrated analysis of transcriptomics and metabonomics profiles in aflatoxin B1-induced hepatotoxicity in rat. Food Chem. Toxicol. 2013, 55, 444–455. [Google Scholar] [CrossRef] [PubMed]
  60. Klein, J. Membrane breakdown in acute and chronic neurodegeneration: Focus on choline-containing phospholipids. J. Neural Transm. 2000, 107, 1027–1063. [Google Scholar] [CrossRef] [PubMed]
  61. Fu, M.H.; Yu, H.Z.; Bao, T.R.; Ba, G.N. Anti-obesity and lipid-lowering mechanism of Corydalis Bungeanae Herba: Based on intestinal microflora and metabolomics. Zhongguo Zhong Yao Za Zhi 2022, 47, 3049–3058. [Google Scholar]
  62. Connors, J.; Dawe, N.; Van Limbergen, J. The Role of Succinate in the Regulation of Intestinal Inflammation. Nutrients 2018, 11, 25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. 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. Detection of Succinate by Intestinal Tuft Cells Triggers a Type 2 Innate Immune Circuit. Immunity 2018, 49, 33–41.e37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Gondo, Y.; Satsu, H.; Ishimoto, Y.; Iwamotom, T.; Shimizu, M. Effect of taurine on mRNA expression of thioredoxin interacting protein in Caco-2 cells. Biochem. Biophys. Res. Commun. 2012, 426, 433–437. [Google Scholar] [CrossRef]
  65. Kim, C.; Cha, Y.N. Taurine chloramine produced from taurine under inflammation provides anti-inflammatory and cytoprotective effects. Amino Acids 2014, 46, 89–100. [Google Scholar] [CrossRef] [PubMed]
  66. Yu, H.; Guo, Z.; Shen, S.; Shan, W. Effects of taurine on gut microbiota and metabolism in mice. Amino Acids 2016, 48, 1601–1617. [Google Scholar] [CrossRef]
  67. Zhao, Z.; Satsu, H.; Fujisawa, M.; Hori, M.; Ishimoto, Y.; Totsuka, M.; Nambu, A.; Kakuta, S.; Ozaki, H.; Shimizu, M. Attenuation by dietary taurine of dextran sulfate sodium-induced colitis in mice and of THP-1-induced damage to intestinal Caco-2 cell monolayers. Amino Acids 2008, 35, 217–224. [Google Scholar] [CrossRef]
  68. Wu, Q.; Jezkova, A.; Yuan, Z.; Pavlikova, L.; Dohnal, V.; Kuca, K. Biological degradation of aflatoxins. Drug Metab. Rev. 2009, 41, 1–7. [Google Scholar] [CrossRef]
  69. Rotimi, O.A.; Rotimi, S.O.; Duru, C.U.; Ebebeinwe, O.J.; Abiodun, A.O.; Oyeniyi, B.O.; Faduyile, F.A. Acute aflatoxin B1—Induced hepatotoxicity alters gene expression and disrupts lipid and lipoprotein metabolism in rats. Toxicol. Rep. 2017, 4, 408–414. [Google Scholar] [CrossRef]
  70. Iqbal, J.; Hussain, M.M. Intestinal lipid absorption. Am. J. Physiol. Endocrinol. Metab. 2009, 296, E1183–E1194. [Google Scholar] [CrossRef] [Green Version]
  71. Siloto, E.V.; Oliveira, E.F.; Sartori, J.R.; Fascina, V.B.; Martins, B.A.; Ledoux, D.R.; Rottinghaus, G.E.; Sartori, D.R. Lipid metabolism of commercial layers fed diets containing aflatoxin, fumonisin, and a binder. Poult. Sci. 2013, 92, 2077–2083. [Google Scholar] [CrossRef] [PubMed]
  72. Bechynska, K.; Kosek, V.; Fenclova, M.; Muchova, L.; Smid, V.; Suk, J.; Chalupsky, K.; Sticova, E.; Hurkova, K.; Hajslova, J.; et al. The Effect of Mycotoxins and Silymarin on Liver Lipidome of Mice with Non-Alcoholic Fatty Liver Disease. Biomolecules 2021, 11, 1723. [Google Scholar] [CrossRef] [PubMed]
  73. Pizzuto, M.; Pelegrin, P. Cardiolipin in Immune Signaling and Cell Death. Trends Cell Biol. 2020, 30, 892–903. [Google Scholar] [CrossRef]
  74. D’Angelo, G.; Moorthi, S.; Luberto, C. Role and Function of Sphingomyelin Biosynthesis in the Development of Cancer. Adv. Cancer Res. 2018, 140, 61–96. [Google Scholar] [PubMed]
  75. Kennelly, J.P.; van der Veen, J.N.; Nelson, R.C.; Leonard, K.A.; Havinga, R.; Buteau, J.; Kuipers, F.; Jacobs, R.L. Intestinal de novo phosphatidylcholine synthesis is required for dietary lipid absorption and metabolic homeostasis. J. Lipid Res. 2018, 59, 1695–1708. [Google Scholar] [CrossRef] [Green Version]
  76. Ugbaja, R.N.; Okedairo, O.M.; Oloyede, A.R.; Ugwor, E.I.; Akinloye, D.I.; Ojo, O.P.; Ademuyiwa, O. Probiotics consortium synergistically ameliorates aflatoxin B1-induced disruptions in lipid metabolism of female albino rats. Toxicon 2020, 186, 109–119. [Google Scholar] [CrossRef] [PubMed]
  77. Coats, S.R.; Hashim, A.; Paramonov, N.A.; To, T.T.; Curtis, M.A.; Darveau, R.P. Cardiolipins Act as a Selective Barrier to Toll-Like Receptor 4 Activation in the Intestine. Appl. Environ. Microbiol. 2016, 82, 4264–4278. [Google Scholar] [CrossRef] [Green Version]
  78. Klein, M.E.; Mauch, S.; Rieckmann, M.; Martinez, D.G.; Hause, G.; Noutsias, M.; Hofmann, U.; Lucas, H.; Meister, A.; Ramos, G. Phosphatidylserine (PS) and phosphatidylglycerol (PG) nanodispersions as potential anti-inflammatory therapeutics: Comparison of in vitro activity and impact of pegylation. Nanomedicine 2020, 23, 102096. [Google Scholar] [CrossRef]
  79. Luo, H.; Jiang, M.; Lian, G.; Liu, Q.; Shi, M.; Li, T.Y.; Song, L.; Ye, J.; He, Y.; Yao, L. AIDA Selectively Mediates Downregulation of Fat Synthesis Enzymes by ERAD to Retard Intestinal Fat Absorption and Prevent Obesity. Cell Metab. 2018, 27, 843–853.e846. [Google Scholar] [CrossRef] [Green Version]
  80. Bates, P.D.; Browse, J. The significance of different diacylgycerol synthesis pathways on plant oil composition and bioengineering. Front. Plant Sci. 2012, 3, 147. [Google Scholar] [CrossRef] [Green Version]
  81. Lim, S.A.; Su, W.; Chapman, N.M.; Chi, H. Lipid metabolism in T cell signaling and function. Nat. Chem. Biol. 2022, 18, 470–481. [Google Scholar] [CrossRef] [PubMed]
  82. Manzo, T.; Prentice, B.M.; Anderson, K.G.; Raman, A.; Schalck, A.; Codreanu, G.S.; Nava Lauson, C.B.; Tiberti, S.; Raimondi, A.; Jones, M.A.; et al. Accumulation of long-chain fatty acids in the tumor microenvironment drives dysfunction in intrapancreatic CD8+ T cells. J. Exp. Med. 2020, 217, e20191920. [Google Scholar] [CrossRef] [PubMed]
  83. Liu, X.; Hartman, C.L.; Li, L.; Albert, C.J.; Si, F.; Gao, A.; Huang, L.; Zhao, Y.; Lin, W.; Hsueh, E.C.; et al. Reprogramming lipid metabolism prevents effector T cell senescence and enhances tumor immunotherapy. Sci. Transl. Med. 2021, 13, eaaz6314. [Google Scholar] [CrossRef]
  84. Mariño, E.; Richards, J.L.; McLeod, K.H.; Stanley, D.; Yap, Y.A.; Knight, J.; McKenzie, C.; Kranich, J.; Oliveira, A.C.; Rossello, F.J.; et al. Gut microbial metabolites limit the frequency of autoimmune T cells and protect against type 1 diabetes. Nat. Immunol. 2017, 18, 552–562. [Google Scholar] [CrossRef] [PubMed]
  85. Yang, X.; Gao, Y.; Huang, S.; Su, C.; Wang, J.; Zheng, N. Whole transcriptome-based ceRNA network analysis revealed ochratoxin A-induced compromised intestinal tight junction proteins through WNT/Ca2+ signaling pathway. Ecotoxicol. Environ. Saf. 2021, 224, 112637. [Google Scholar] [CrossRef] [PubMed]
  86. Jang, C.; Chen, L.; Rabinowitz, J.D. Metabolomics and Isotope Tracing. Cell 2018, 173, 822–837. [Google Scholar] [CrossRef] [Green Version]
  87. Song, J.W.; Lam, S.M.; Fan, X.; Cao, W.J.; Wang, S.Y.; Tian, H.; Chua, G.H.; Zhang, C.; Meng, F.P.; Xu, Z. Omics-Driven Systems Interrogation of Metabolic Dysregulation in COVID-19 Pathogenesis. Cell Metab. 2020, 32, 188–202.e5. [Google Scholar] [CrossRef]
  88. Lam, S.M.; Zhang, C.; Wang, Z.; Ni, Z.; Zhang, S.; Yang, S.; Huang, X.; Mo, L.; Li, J.; Lee, B. A multi-omics investigation of the composition and function of extracellular vesicles along the temporal trajectory of COVID-19. Nat. Metab. 2021, 3, 909–922. [Google Scholar] [CrossRef]
  89. Shui, G.; Cheong, W.F.; Jappar, I.A.; Hoi, A.; Xue, Y.; Fernandis, A.Z.; Tan, B.K.; Wenk, M.R. Derivatization-independent cholesterol analysis in crude lipid extracts by liquid chromatography/mass spectrometry: Applications to a rabbit model for atherosclerosis. J. Chromatogr. A 2011, 1218, 4357–4365. [Google Scholar] [CrossRef]
Figure 1. Dose−response curves of cell viability for effect of aflatoxin B1 (AFB1) (A), aflatoxin M1 AFM1 (B), and combination of AFB1 and AFM1 (AFB1+AFM1) (C) on NCM460 cells (n = 10).
Figure 1. Dose−response curves of cell viability for effect of aflatoxin B1 (AFB1) (A), aflatoxin M1 AFM1 (B), and combination of AFB1 and AFM1 (AFB1+AFM1) (C) on NCM460 cells (n = 10).
Toxins 15 00255 g001
Figure 2. Overview of metabolic profiles on NCM460 cells after aflatoxin treatment. Orthogonal projection to latent structure−discriminant analysis (OPLS−DA) for different treatment groups (n = 8): (A) AFB1, (B) AFM1, and (C) AFB1 + AFM1 compared with control (CTL).
Figure 2. Overview of metabolic profiles on NCM460 cells after aflatoxin treatment. Orthogonal projection to latent structure−discriminant analysis (OPLS−DA) for different treatment groups (n = 8): (A) AFB1, (B) AFM1, and (C) AFB1 + AFM1 compared with control (CTL).
Toxins 15 00255 g002
Figure 3. The analysis of differential metabolites after aflatoxins treatment alone or in combination. (A) The number of differential metabolites. (B) The Venn analysis of different metabolites. (C) Heatmap visualization displaying a relative abundance of specific and common metabolites (n = 8). A total of 52 significantly changed specific common metabolites are listed on the map. Each row in the map represents one metabolite, and each column represents one sample. The gradient color bar from purple to green denotes the concentrations of metabolites from high to low. The metabolites with different colors on the right represent the metabolites of different treatment groups, in which orange is the specific metabolite of the AFB1 group, pale blue-green is AFM1, sky blue is AFB1 + AFM1, and black is the common differential metabolites.
Figure 3. The analysis of differential metabolites after aflatoxins treatment alone or in combination. (A) The number of differential metabolites. (B) The Venn analysis of different metabolites. (C) Heatmap visualization displaying a relative abundance of specific and common metabolites (n = 8). A total of 52 significantly changed specific common metabolites are listed on the map. Each row in the map represents one metabolite, and each column represents one sample. The gradient color bar from purple to green denotes the concentrations of metabolites from high to low. The metabolites with different colors on the right represent the metabolites of different treatment groups, in which orange is the specific metabolite of the AFB1 group, pale blue-green is AFM1, sky blue is AFB1 + AFM1, and black is the common differential metabolites.
Toxins 15 00255 g003
Figure 4. Metabolic pathway analysis (MetPA) of differentially expressed metabolites on NCM460 after aflatoxins treatment, (A) AFB1 and (B) combination of AFB1 + AFM1. The darker the color, the larger the circle, the higher the degree of enrichment.
Figure 4. Metabolic pathway analysis (MetPA) of differentially expressed metabolites on NCM460 after aflatoxins treatment, (A) AFB1 and (B) combination of AFB1 + AFM1. The darker the color, the larger the circle, the higher the degree of enrichment.
Toxins 15 00255 g004aToxins 15 00255 g004b
Figure 5. KEGG pathway enrichment analysis of specific metabolites in (A) AFB1, (B) AFM1, and (C) combination of AFB1 + AFM1. The size of the node indicates the degree of enrichment. The solid line indicates that there is a connection relationship between pathways, and the dotted line indicates the connection line between the complementary pathway and the input pathway.
Figure 5. KEGG pathway enrichment analysis of specific metabolites in (A) AFB1, (B) AFM1, and (C) combination of AFB1 + AFM1. The size of the node indicates the degree of enrichment. The solid line indicates that there is a connection relationship between pathways, and the dotted line indicates the connection line between the complementary pathway and the input pathway.
Toxins 15 00255 g005aToxins 15 00255 g005b
Figure 6. Lipids profiles on NCM460 cells after single and combination exposure of AFB1 and AFM1. (A) Subclasses and quantity of detected lipids. (B) PCA analysis of lipids, each point in the quadrants represents a single sample (n = 3).
Figure 6. Lipids profiles on NCM460 cells after single and combination exposure of AFB1 and AFM1. (A) Subclasses and quantity of detected lipids. (B) PCA analysis of lipids, each point in the quadrants represents a single sample (n = 3).
Toxins 15 00255 g006
Figure 7. The analysis of differential lipids induced by single and combination of aflatoxins. (A) The numbers of the differential lipids after pairwise comparisons. (B) Venn analysis of the differential lipids compared with the control group (CTL). The heatmap analysis of specific metabolites in (C) AFB1, (D) AFM1, (E) AFB1 + AFM1 (n = 3).
Figure 7. The analysis of differential lipids induced by single and combination of aflatoxins. (A) The numbers of the differential lipids after pairwise comparisons. (B) Venn analysis of the differential lipids compared with the control group (CTL). The heatmap analysis of specific metabolites in (C) AFB1, (D) AFM1, (E) AFB1 + AFM1 (n = 3).
Toxins 15 00255 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, X.; Li, X.; Gao, Y.; Wang, J.; Zheng, N. Integrated Metabolomics and Lipidomics Analysis Reveals Lipid Metabolic Disorder in NCM460 Cells Caused by Aflatoxin B1 and Aflatoxin M1 Alone and in Combination. Toxins 2023, 15, 255. https://doi.org/10.3390/toxins15040255

AMA Style

Yang X, Li X, Gao Y, Wang J, Zheng N. Integrated Metabolomics and Lipidomics Analysis Reveals Lipid Metabolic Disorder in NCM460 Cells Caused by Aflatoxin B1 and Aflatoxin M1 Alone and in Combination. Toxins. 2023; 15(4):255. https://doi.org/10.3390/toxins15040255

Chicago/Turabian Style

Yang, Xue, Xue Li, Yanan Gao, Jiaqi Wang, and Nan Zheng. 2023. "Integrated Metabolomics and Lipidomics Analysis Reveals Lipid Metabolic Disorder in NCM460 Cells Caused by Aflatoxin B1 and Aflatoxin M1 Alone and in Combination" Toxins 15, no. 4: 255. https://doi.org/10.3390/toxins15040255

APA Style

Yang, X., Li, X., Gao, Y., Wang, J., & Zheng, N. (2023). Integrated Metabolomics and Lipidomics Analysis Reveals Lipid Metabolic Disorder in NCM460 Cells Caused by Aflatoxin B1 and Aflatoxin M1 Alone and in Combination. Toxins, 15(4), 255. https://doi.org/10.3390/toxins15040255

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop