Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 07 March 2024
Sec. Metabolomics

A GC–MS-based untargeted metabolomics approach for comprehensive metabolic profiling of mycophenolate mofetil-induced toxicity in mice

Tongfeng Zhao&#x;Tongfeng Zhao1Yaxin Zhao&#x;Yaxin Zhao2Haotian ChenHaotian Chen1Wenxue Sun,
Wenxue Sun3,4*Yun Guan
Yun Guan1*
  • 1Department of Hematology, Jining No.1 People’s Hospital, Jining, China
  • 2Department of Pharmacy, Jining No.1 People’s Hospital, Jining, China
  • 3Translational Pharmaceutical Laboratory, Jining No.1 People’s Hospital, Jining, China
  • 4Postdoctoral of Shandong University of Traditional Chinese Medicine, Jinan, China

Background: Mycophenolate mofetil (MMF), the morpholinoethyl ester of mycophenolic acid, is widely used for maintenance immunosuppression in transplantation. The gastrointestinal toxicity of MMF has been widely uncovered. However, the comprehensive metabolic analysis of MMF-induced toxicity is lacking. This study is aimed to ascertain the metabolic changes after MMF administration in mice.

Methods: A total of 700 mg MMF was dissolved in 7 mL dimethyl sulfoxide (DMSO), and then 0.5 mL of mixture was diluted with 4.5 mL of saline (100 mg/kg). Mice in the treatment group (n = 9) were given MMF (0.1 mL/10 g) each day via intraperitoneal injection lasting for 2 weeks, while those in the control group (n = 9) received the same amount of blank solvent (DMSO: saline = 1:9). Gas chromatography-mass spectrometry was utilized to identify the metabolic profiling in serum samples and multiple organ tissues of mice. The potential metabolites were identified using orthogonal partial least squares discrimination analysis. Meanwhile, we used the MetaboAnalyst 5.0 (http://www.metaboanalyst.ca) and Kyoto Encyclopedia of Genes and Genomes database (http://www.kegg.jp) to depict the metabolic pathways. The percentages of lymphocytes in spleens were assessed by multiparameter flow cytometry analysis.

Results: Compared to the control group, we observed that MMF treatment induced differential expression of metabolites in the intestine, hippocampus, lung, liver, kidney, heart, serum, and cortex tissues. Subsequently, we demonstrated that multiple amino acids metabolism and fatty acids biosynthesis were disrupted following MMF treatment. Additionally, MMF challenge dramatically increased CD4+ T cell percentages but had no significant influences on other types of lymphocytes.

Conclusion: MMF can affect the metabolism in various organs and serum in mice. These data may provide preliminary judgement for MMF-induced toxicity and understand the metabolic mechanism of MMF more comprehensively.

Introduction

Suppression of the immune system is crucial after organ transplantation and is usually administrated for the treatment of a variety of autoimmune diseases. Mycophenolate mofetil (MMF), the morpholinoethyl ester of mycophenolic acid (MPA), is widely used for maintenance immunosuppression in solid organ, bone marrow, and stem cell transplantation (Srinivas et al., 2003a; Moiseev et al., 2016a; Carlone et al., 2021a). MMF also serves as a substitution for cyclosporine A therapy-induced severe nephrotoxicity or hemolytic uremic syndrome in renal transplant recipients (Ulinski et al., 2005). MMF exerts inhibitory effect on the proliferation of both B- and T-lymphocytes via non-competitively and reversibly blocks the de novo synthesis of guanine nucleotides required for DNA and RNA synthesis during lymphocyte proliferation (Wu, 1994). Unfortunately, it has been proved that the application of immunosuppressant is associated with the occurrence of hematologic toxicity or other organ toxicity.

The gastrointestinal (GI) toxicity is the most common adverse effects of MMF including diarrhea, abdominal pain, nausea, and vomiting, which result in dose reduction or drug withdrawal, and the increased risk of rejection and death after transplantation (Calmet et al., 2015). It has been reported that friability on endoscopy was associated with severe disease; and nausea and erythema were related to poor prognosis (Bhattacharya et al., 2022). Recently, a study demonstrated an intact intestinal microbiota was essential to initiate and sustain the MMF-induced GI toxicity, and MMF exposure was associated with alterations of intestinal flora composition, companied the increase in genes involved in lipopolysaccharide (LPS) biosynthesis (Flannigan et al., 2018). However, the underlying mechanism of MMF-induced GI toxicity remains unclarified. Meanwhile, the administration of MMF is more frequently associated with hematologic toxicity, such as anemia, due to bone marrow suppression or hemolysis, leukopenia, and thrombocytopenia (Danesi and Del Tacca, 2004). Neutropenia and leukopenia have been reported as severe hematologic toxicities associated with MMF treatment in transplant recipients (Nogueras et al., 2005; Varnell et al., 2017). In addition, uncommon side effect, such as hepatotoxicity, is also observed in renal transplant recipients (Balal et al., 2005), and a recent study has proved that MMF-induced hepatotoxicity is associated with mitochondrial abnormality in liver transplant recipients and mice (Warren et al., 2021). Therefore, it is paramount to investigate the mechanisms responsible for MMF-associated toxicity.

Metabolomics approach is widely implemented in systems biology research and it can be systematically used to identify potential metabolic biomarkers under the condition of stimulation or impact. In this study, a GC-MS-based untargeted metabolomics approach was used to investigate the underlying mechanism of MMF-associated toxicity. Alterations in metabolites of serum and other organs were identified, and disrupted metabolic pathways were analyzed to elucidate the metabolic profiling after MMF exposure. This study represents the first comprehensive evaluation of metabolic profiling in MMF-treated mice.

Materials and methods

Chemicals and reagents

MMF with purity of 99%, heptadecanoic acid (purity ≥ 98%), methanol (chromatographic grade), and pyridine were purchased from Macklin Biochemical (Shanghai, China). O-methylhydroxylamine hydrochloride (purity ≥ 98%) was obtained from J&K Scientific Ltd. (Beijing, China). N, O-bis(trimethylsilyl)trifluoroacetamide (containing 1% trimethylchlorosilane) was obtained from Sigma Aldrich (St. Louis, MO, United States). The purified water was purchased from Wahaha (Hangzhou, China). RPMI medium was obtained from (Thermo Fisher Scientific, Waltham, MA, United States).

Animals and treatment

A total of 18 C57BL/6 mice weighting 35 ± 5 g (6-week-old) were purchased from Pengyue Experimental Animal Breeding Co., Ltd. (Jinan, China). Animals were housed under a 12 h light/dark cycle at 20°C–22°C. All experimental procedures conformed to the Guidelines for the Use of Laboratory Animals, and approved by the Ethical Committee for Animal Experimentation of Jining First People’s Hospital (Approval No. JNRM-2023-DW-021).

A total of 700 mg MMF was dissolved in 7 mL dimethyl sulfoxide (DMSO), and then 0.5 mL of mixture was diluted with 4.5 mL of saline (100 mg/kg). After 1-week of acclimation, mice in MMF group (n = 9) were intraperitoneally injected with MMF (0.1 mL/10g) every day for 2 weeks, and the dose was adjusted with the increase in body weight. Mice received the same amount of blank solvent (DMSO: saline = 1:9) served as control group.

Sample collection

Mice were anesthetized with 1% sodium pentobarbital (50 mg/kg) after the last administration in mice. Blood samples were collected from each mouse after eyeball enucleation, centrifuged at 4,000 rpm for 10 min at 4°C, and serum samples were then stored at −80°C. Subsequently, mice were sacrificed by cervical dislocation and at the same time, the intestine, hippocampus, lung, liver, kidney, heart, cortex and spleen tissues were immediately collected on ice, frozen in liquid nitrogen, and then stored at −80°C for succeeding experiments.

Sample preparation

A 100-μL aliquot of sample was mixed with 350 μL of heptadecanoic acid (100 μg/mL in methanol). After centrifugation at 14,000 rpm for 15 min at 4°C, the supernatant liquid was dried with liquid nitrogen at 37°C. Subsequently, 80 μL of O-methylhydroxylamine hydrochloride (15 mg/mL in pyridine) was added and incubated at 70°C for 90 min. A total of 100 μL of N, O-bis(trimethylsilyl)trifluoroacetamide containing 1% trimethyl chlorosilane was added to each sample, followed by incubation at 70°C for 60 min.

Tissue sample (50 mg for each) was homogenized with 1 mL of methanol, mixed with 50 μL of heptadecanoic acid (1 mg/mL in methanol), and centrifuged at 14,000 rpm for 15 min at 4°C. O-methylhydroxylamine hydrochloride (80 μL; 15 mg/mL in pyridine) was then added at 70°C for 90 min, mixed with 100 μL of N, O-bis(trimethylsilyl)trifluoroacetamide (containing 1% trimethyl chlorosilane), and incubated at 70°C for 60 min. Pooling 10 µL of each sample of control group and MMF group served as quality control (QC). A 0.22-μm filter was used to purify samples for further GC-MS analysis.

Lymphocytes isolation

The excised spleens were placed in RPMI medium. For the preparation of single cell suspensions, spleen tissues were firstly ground, passed through a 70-mm-pore mesh, and washed with RPMI. After that, the samples were centrifuged (350 × g) at 4°C for 10 min, followed by abandoning the supernatants. In order to lyse the erythrocytes, the cells were resuspended in 3 mL of TRIS-ammonium chloride solution (0.144 M of NH4Cl plus 0.017 M of TRIS; pH 7.2), incubated at 4°C for 2 min, and then washed with RPMI for twice. After centrifugation (350 × g at 4°C for 10 min), the cells were re-suspended in RPMI containing 10% FBS and lymphocytes from spleen were obtained.

Multiparameter flow cytometry analysis

Multiparameter flow cytometry was performed according to a standard protocol. Single-cell suspensions were incubated with anti-CD16/32 monoclonal antibody (Biolegend, San Diego, CA, United States) for 15 min at room temperature to block Fc receptors before staining with the specific antibodies. Spleen mononuclear cells were then stained with the indicated fluorescent monoclonal antibodies for surface molecules. Intracellular staining was performed using fixation and permeabilization buffers (eBioscience) according to the manufacturers’ instructions.

GC-MS based metabolomics analysis

An Agilent 7890B GC system coupled to a 7000C GC/MS Triple Quad Mass Detector (Agilent Technologies, United States), equipped with an HP-5MS fused silica capillary column (30 m × 0.25 mm × 0.25 μm), was implemented for metabolomics analysis. Helium was used as the carrier gas with a flow rate at 1 mL/min. Sample (1 μL) was injected into GC-MS with a split ratio of 50:1. The injection temperature was set to 280°C, transfer line temperature was 250°C, and ion source temperature was 230°C, respectively. Electron collision ionization was set to −70 EV, and the frequency of acquisition was 20 spectra/s. The ionization mode of mass spectrometry is electrospray ionization with a mass/charge (m/z) full scan range of 50–800.

We have successfully uploaded the source data from our GC-MS analysis to MetaboLights, a move aimed at enhancing transparency and simplifying access to the raw data for the scientific community. Our research project is exclusively identified by the code MTBLS9315. To delve into our study, you can follow this link: https://www.ebi.ac.uk/metabolights/MTBLS9315.

Data processing and multivariate analysis

Raw data from GC-MS was analyzed using the Agilent Mass Hunter (Version B.07.00, Agilent Technologies, CA, United States). The preprocessing included alignment, retention time correction, baseline filtration, and deconvolution. For metabolite identification, a library containing all QC samples was established, and the U.S. National Institute of Standards and Technology (NIST 14) GC-MS library was applied to identify the unknown metabolites of QC samples. The metabolites with similarity > 80% were considered as structurally identified (Caterino et al., 2020; Zhang et al., 2021). Subsequently, a new spectrum library named “New Library” was constructed, and used for spectrum matching of metabolites of experimental samples. Manually validation was conducted to reduce deconvolution errors in the process of automatic data processing and to eliminate false identifications.

An integrated data matrix composed of the peak index (RT-m/z pair), sample name, and corresponding peak area was obtained. Peak area was normalized using Microsoft Excel™ (Microsoft, Redmond, WA, United States). SIMCA-P 14.0 (Umetrics, Sartorius Stedim Biotech) was implemented to further analyze the data by multivariate analysis including principal components analysis (PCA) and orthogonal partial least squares discrimination analysis (OPLS-DA). The two-tailed Student’s t-tests were carried out to compare differences of the two groups. Compounds with variable importance in projection (VIP) values > 1.0 and two-tailed Student’s t-test p values < 0.05 were considered as potential differential expressed metabolites. Clustered heatmap diagram and bubble diagram (pathway enrichment) were constructed using MetaboAnalyst 5.0 (http://www.metaboanalyst.ca).

Results

GC–MS total ion chromatograms (TICs) of samples

As illustrated in Figure 1, representative GC-MS TICs of QC from tissue and serum samples are presented. Supplementary Tables S1–S8 includes retention time, metabolites, % similarity, molecular formula, CAS number, peak area, and molecular weight for different tissues and serum. The QC in our study were prepared by pooling 10 µL of each sample from the control group and MMF group. We indicated that there were significant differences in TICs among different QC samples. Supplementary Figure S1 showed that there was a linear positive correlation between retention time and retention index curve.

FIGURE 1
www.frontiersin.org

FIGURE 1. Representative gas chromatography–mass spectrometry (GC–MS) total ion chromatograms (TICs) of quality control (QC). (A) Intestine (B) Hippocampus (C) Lung (D) Liver (E) Kidney (F) Heart (G) Serum (H) Cortex.

Metabolomic data analysis

PCA was used to distinguish the metabolic profiles. As illustrated in Figures 2A–H, clear differences between the control and the MMF groups were observed in intestine, hippocampus, lung, liver, kidney, heart, serum and cortex samples. In PCA model, the total variance explained was 76.1%, 57.3%, 64.3%, 41.5%, 59.3%, 34.6%, 36.0%, and 43.9% respectively, indicating satisfactory construct validity. The parameter scores of OPLS-DA were shown in Table 1. Meanwhile, a ranking test was conducted to verify the validity of OPLS-DA model. The abscissa represents the predicted principal component score of the first principal component, the ordinate represents the orthogonal principal component score, and the color represents different experimental groups. In different tissues and serum, the proportion of variance explained by the OPLS-DA model (R2X) was 72.8%, 51.5%, 53.9%, 42.3%, 60.1%, 42.6%, 36.1%, and 43.7%, respectively. We found that the intersection points between blue regression line (Q2-point) and vertical axis were all negative values (Figures 2A–H). These results confirmed the validity of OPLS-DA model and at the same time, the significantly differences across the MMF and control groups were validated.

FIGURE 2
www.frontiersin.org

FIGURE 2. PCA analysis, OPLS-DA score plots and 200 permutation tests. (A) Intestine (B) Hippocampus (C) Lung (D) Liver (E) Kidney (F) Heart (G) Serum (H) Cortex.

TABLE 1
www.frontiersin.org

TABLE 1. OPLS-DA parameter scores.

Identification of differential metabolites

The differential metabolites in different tissues and serum samples after MMF treatment were exhibited in Table 2. VIP value > 1 and at the same time p-value < 0.05 were considered as the important criteria for potential metabolites. Additionally, the values of fold-change (<1) represented that the metabolite levels appeared decreasing trend, whilst those (>1) showed a rising trend. As a results, we found that there was a total of 11 differential metabolites existed in intestine tissues. All of them were upregulated. In hippocampus tissues, 15 differential metabolites were identified. Among these metabolites, 7 metabolites (L-Isoleucine, Serine, L-Threonine, L-Aspartic acid, L-Phenylalanine, L-Alanine and L-Glutamic acid) were observed to be downregulated, while 8 metabolites [O-Phosphoethanolamine, Myristic acid, Palmitic acid, MG(0:0/18:0/0:0), 4-Hydroxybenzyl alcohol, Octadecane, 2,5-Furandicarboxylic acid and MG(16:0/0:0/0:0)] were upregulated. Moreover, we observed MMF treatment altered a total of 8 metabolites in lung tissues [L-Threonine, Glycine, Uracil, MG(0:0/18:0/0:0), N-Dodecane, Ellagic acid, Azelaic acid and MG(16:0/0:0/0:0)], with a downward trend. Only 5 upward metabolites [Palmitic acid, Stearic acid, Uridine, Inosine and MG(16:0/0:0/0:0)] were observed in liver tissues. In the kidney, 12 altered metabolites including 6 downregulated (L-Alanine, Glycine, Glyoxylic acid, Maleic acid, Inosine and Cholesterol) and 6 upregulated metabolites [Myristic acid, L-Tyrosine, O-Phosphoethanolamine, Palmitic acid, MG(16:0/0:0/0:0) and MG(0:0/18:0/0:0)] were identified following MMF treatment. Additionally, in heart samples, a total of 14 altered metabolites were identified, including 4 downregulated [myo-Inositol, MG(16:0/0:0/0:0), Adenosine and MG(0:0/18:0/0:0)] and 10 upregulated metabolites (Niacinamide, L-Glutamic acid, Octadecane, Myristic acid, Palmitic acid, Stearic acid, Arachidic acid, Uridine, L-Phenylalanine and Inosine). In serum samples, 4 downregulated [D-Lactic acid, L-Aspartic acid, MG(16:0/0:0/0:0) and MG(0:0/18:0/0:0)] and 5 upregulated (Serine, Stearic acid, L-Alanine, myo-Inositol and 4-Hydroxybenzyl alcohol) differential metabolites were identified. In cortex samples, a total of 19 downregulated differential metabolites [L-Lactic acid, Ethanolamine, Glycine, Serine, L-Threonine, L-Aspartic acid, L-Phenylalanine, myo-Inositol, L-Alanine, N-Dodecane, Malic acid, Pyroglutamic acid, 1-Hexadecanol, L-Glutamic acid, scyllo-Inositol, MG(16:0/0:0/0:0), Adenosine, MG(0:0/18:0/0:0) and Cholesterol] were identified. As manifested in Figures 3A–H, based on the cluster analyses of differential metabolites between the MMF and control groups, the above results were further validated.

TABLE 2
www.frontiersin.org

TABLE 2. Differential metabolites in different tissues and serum samples after mycophenolate mofetil treatment.

FIGURE 3
www.frontiersin.org

FIGURE 3. Heatmap of differential metabolites in (A) Intestine (B) Hippocampus (C) Lung (D) Liver (E) Kidney (F) Heart (G) Serum (H) Cortex in the mycophenolate mofetil (MMF) groups compared with controls. The color of each part represents the importance of metabolite changes (blue, downregulated; red, upregulated). Rows represent samples, and columns represent metabolites.

Metabolic pathway enrichment analysis

In order to comprehensively understand the dynamic adaptation of metabolism to MMF treatment, pathway enrichment analysis was performed based on differential metabolites at different tissues. We used Metaboanalyst 5.0 (http://www.metaboanalyst.ca) and KEGG database (http://www.kegg.jp) to analyze the altered metabolic pathways following MMF treatment. Metabolic pathways that the data of Raw p less than 0.05 and Impact more than 0 were considered as potential disturbed pathways (Figure 4; Table 3). A detailed metabolic network was presented in Figure 5. We found that MMF challenge not only affected several kinds of amino acid metabolic pathways, but also Phenylalanine, tyrosine and tryptophan biosynthesis; Valine, leucine and isoleucine biosynthesis; Galactose metabolism.

FIGURE 4
www.frontiersin.org

FIGURE 4. Summary of pathway analysis performed using MetaboAnalyst 5.0. (A) Intestine: (a) Aminoacyl-tRNA biosynthesis; (b) Galactose metabolism; (c) Phenylalanine, tyrosine and tryptophan biosynthesis; (d) Fatty acid biosynthesis. (B) Hippocampus: (c) Phenylalanine, tyrosine and tryptophan biosynthesis; (e) Aminoacyl-tRNA biosynthesis; (f) Alanine, aspartate and glutamate metabolism; (g) Valine, leucine and isoleucine biosynthesis; (h) Arginine biosynthesis; (i) Histidine metabolism; (j) Nitrogen metabolism; (k) D-Glutamine and D-glutamate metabolism; (C) Lung: (e) Aminoacyl-tRNA biosynthesis; (g) Valine, leucine and isoleucine biosynthesis; (l) Glycine, serine and threonine metabolism. (D) Liver: (a) Biosynthesis of unsaturated fatty acids. (E) Kidney: (c) Phenylalanine, tyrosine and tryptophan biosynthesis; (d) Fatty acid biosynthesis; (e) Aminoacyl-tRNA biosynthesis; (l) Glycine, serine and threonine metabolism; (m) Glyoxylate and dicarboxylate metabolism; (n) Primary bile acid biosynthesis. (F) Heart: (a) Biosynthesis of unsaturated fatty acids; (c) Phenylalanine, tyrosine and tryptophan biosynthesis; (d) Fatty acid biosynthesis; (e) Aminoacyl-tRNA biosynthesis; (j) Nitrogen metabolism; (k) D-Glutamine and D-glutamate metabolism. (G) Serum: (e) Aminoacyl-tRNA biosynthesis; (f) Alanine, aspartate and glutamate metabolism; (o) Ascorbate and aldarate metabolism. (H) Cortex: (c) Phenylalanine, tyrosine and tryptophan biosynthesis; (e) Aminoacyl-tRNA biosynthesis; (f) Alanine, aspartate and glutamate metabolism; (h) Arginine biosynthesis; (i) Histidine metabolism; (m) Glyoxylate and dicarboxylate metabolism; (p) Glutathione metabolism; (q) Porphyrin and chlorophyll metabolism. Black font represents pathways with p-value < 0.05 and impact = 0. Red font represents pathways with p-value < 0.05 and impact > 0. No labeling font represents pathways with p-value > 0.05.

TABLE 3
www.frontiersin.org

TABLE 3. Pathway analysis performed using MetaboAnalyst 5.0 software.

FIGURE 5
www.frontiersin.org

FIGURE 5. Schematic diagram of related metabolic pathways affected by MMF in serum and major tissues. The activation metabolic pathways were marked in red box. Solid arrows represent a single process, while dashed arrows represent multiple processes. Differential metabolites enriched in pathways were marked in bold.

Multiparameter flow cytometry analysis for lymphocytes

Spleen contains a large number of lymphocytes, which plays an important role in the immune processes. Therefore, we isolated lymphocytes from the spleen tissues of mice in different groups to validate the immunosuppressive mechanism of MMF. The percentages of CD3+ T cells, CD4+ T cells, CD8+ T cells, Natural Killer (NK) cells, B cells, peripheral granulocytic-myeloid-derived suppressor cells (G-MDSCs), monocytic-myeloid-derived suppressor cells (M-MDSCs) and Treg were calculated. The results indicated that MMF administration dramatically increased CD4+ T cell percentages (p-value < 0.05) but had no significant influences on other types of lymphocytes (Figure 6).

FIGURE 6
www.frontiersin.org

FIGURE 6. Multiparameter flow cytometry analysis for the percentages of lymphocytes from spleen. (A) Representative images of multiparameter flow cytometry analysis. (B) The percentages of CD3+ T cells, CD4+ T cells, CD8+ T cells, Natural Killer (NK) cells, B cells, Treg cells, peripheral granulocytic-myeloid-derived suppressor cells (G-MDSCs) and monocytic-myeloid-derived suppressor cells (M-MDSCs) in the control and MMF groups. *p < 0.05. ns, no significance.

Discussion

MMF is a kind of new immunosuppressive agents, which can prevent the replication of T and B lymphocytes by inhibiting purine synthesis (Allison and Eugui, 1993). The usage of MMF is generally associated with gastrointestinal side effects and among these incidental consequences, treatment of MMF has been confirmed to cause the disorder of several metabolic pathways such as glucuronide/glucoside, plasma bioelement contents, and nucleotide and lipid metabolism (Shipkova et al., 2005; Kaminska et al., 2012; Heischmann et al., 2017). However, the metabolic profiling of each organ following MMF treatment is still poorly understood.

In the current work, a GC-MS approach was utilized to ascertain the alteration of metabolite in serum and tissue samples (intestine, hippocampus, lung, liver, kidney, heart, cortex) of MMF-challenged mice. MMF can influence various organs through diverse pathways, encompassing variations in drug metabolism, organ-specific physiological responses, and heterogeneous drug distribution. Distinct metabolic product profiles may emerge in different organs due to differences in drug metabolism pathways. Moreover, the unique physiological functions and metabolic characteristics of each organ may yield diverse metabolic effects. Additionally, the non-uniform distribution of the drug in the body can result in varying organ exposures to different drug concentrations, contributing to differences in metabolic levels (Liu et al., 2022). This can clarify why MMF has a completely different metabolic output on various organs. Additionally, based on the pathway analysis results via MetaboAnalyst 5.0 software, we found that several amino acid metabolism-related pathways including alanine, aspartate and glutamate metabolism, arginine biosynthesis, histidine metabolism, and phenylalanine, tyrosine and tryptophan biosynthesis were remarkably disrupted in hippocampus and cortex. Specially, the alteration of nitrogen metabolism and D-Glutamine and D-glutamate metabolism in hippocampus, and glutathione metabolism, porphyrin and chlorophyll metabolism, and glyoxylate and dicarboxylate metabolism in cortex were observed. Some previous researches have demonstrated that the onset of seizures and encephalopathy may be related to the use of MMF in transplant patients (Derle et al., 2015; Pellerin et al., 2018). Generally, lipophilic antibiotics can reach the infection compartment through blood-brain barrier and it is well-known that MMF is indeed lipophilic (Iqbal et al., 2020). Therefore, we speculated that MMF may possess neurotoxicity via affecting the metabolism disorders in hippocampus and cortex.

Amino acids, as important substrates, play a regulatory role in many metabolic pathways and function as diagnostic markers of many diseases (Razak et al., 2017). Several kinds of amino acids including glycine, serine, L-isoleucine, L-threonine, L-aspartic acid, L-phenylalanine, L-alanine, and L-glutamic acid were downregulated in MMF group compared to those of control group. Glycine has very vital roles in cytoprotection, immune response, growth, development, metabolism, and survival of many mammals and humans. In central nervous system, glycine is considered as an important neurotransmitter to control food intake and body homeostasis (Rajendra et al., 1997) and at the same time, it can regulate the immune function, superoxide generation and cytokines synthesis by modulating the levels of intracellular Ca2+ (Zhong et al., 2003). Additionally, we also found that excessive fatty acids such as palmitic acid, myristic acid, stearic acid were existed in the intestine, hippocampus, liver, kidney, heart tissues and serum samples of rats following MMF treatment, while high levels of fatty acids have been confirmed to inhibit the amount of glycine (Dasarathy et al., 2009). We speculated that MMF administration may not only induce neurotoxicity but also cause toxicity in other organs.

To further explore the mechanism of MMF-induced GI toxicity, we found the disorder of phenylalanine, tyrosine and tryptophan biosynthesis, and fatty acid biosynthesis in intestine tissues. Apart from the effects of the disruption of these two metabolic pathways on intestine metabolism, we further indicated that glycerol and sorbitol that involved in galactose metabolism were two main differential metabolites in intestine tissues, indicating that MMF treatment also disrupted galactose metabolism, as shown in Figure 5. Lactose is a disaccharide of galactose and glucose that is hydrolyzed to monosaccharides by enterocytes. Galactose can be absorbed by the mature enterocytes at the tips of the villi (Wright, 2013). Generally, galactose metabolism can be catalyzed sequentially by three enzymatic steps, with the aid of enzymes galactokinase, UTP-hexose-1-phosphate uridylyltransferase and UDP-galactose 4′-epimerase, respectively (Coelho et al., 2015). At the same time, immunosuppression can universally repress the activity of UDP-galactose 4′-epimerase (Lee et al., 2014), while deficiency of each one of the three galactose-metabolism-related enzymes can lead to galactosemias, such as Escherichia coli sepsis (Novelli and Reichardt, 2000). Because MMF is widely used for maintenance immunosuppression in solid organ, bone marrow, and stem cell transplantation (Srinivas et al., 2003b; Moiseev et al., 2016b; Carlone et al., 2021b), we therefore speculated that MMF treatment may induce intestine toxicity via immunosuppression. Additionally, it is reported that after efflux from the enterocyte, approximately 88% of the ingested galactose is retained in the liver, while the remaining amounts are transported into brain for amino acid biosynthesis (Roser et al., 2009; Augustin, 2010; Chichlowski et al., 2011). Therefore, the disorder of galactose metabolism may also affect the biosynthesis of amino acids in brains, which may further validate the neurotoxicity of MMF administration.

Glutathione (GSH) is a crucial endogenous antioxidant and glutamate serving as a precursor to exert important functions for GSH generation (Mothet et al., 2005). The disorder of D-Glutamine and D-glutamate metabolism was not only observed in hippocampus, but also in the heart tissues. Meanwhile, we also demonstrated that glutathione metabolism was disrupted in cortex. These results implied that following treatment of MMF, the imbalance of redox may be occurred in hippocampus, heart and cortex. Nitrogen is essential for the synthesis of a variety of biomolecules including amino acids, porphyrins, nucleotides, glutathione and other critical biological compounds. Glutamine and glutamate metabolism are central to nitrogen metabolism, and low levels of glutamine and glutamate can limit the flow of nitrogen within cells and delay biomass production and growth (Kurmi and Haigis, 2020). Therefore, it is speculated that MMF administration may first disrupt the metabolism of glutamine and glutamate and then nitrogen metabolism, eventually disturbing the balance of amino acid biosynthesis, glutathione metabolism and porphyrin and chlorophyll metabolism in multiple organs. Additionally, some researches have uncovered that the absence of tryptophan can promote immune suppression and at the same time, arginine is confirmed to exert vital role in T cell activation, especially the survival of CD4+ T cells (O’Sullivan et al., 2019; Geiger et al., 2016). In our study, we found that phenylalanine, tyrosine and tryptophan biosynthesis was disrupted in multiple organs including intestine, hippocampus, kidney, heart and cortex, while the disorder of arginine biosynthesis was observed in hippocampus and cortex. We hypothesized that MMF treatment may disrupt phenylalanine, tyrosine and tryptophan biosynthesis and arginine biosynthesis to exert immunosuppression functions. Additionally, the percentages of CD3+ T cells, CD4+ T cells, CD8+ T cells, NK cells, B cells, MDSCs and Treg were also calculated. A report conducted by Sousa et al. has demonstrated that the depletion of CD4+ T cells is directly associated with the activation of immune system (Sousa et al., 2002). In the current study, the results indicated that MMF administration dramatically increased CD4+ T cell percentages but had no significant influences on other types of lymphocytes, indicating that the increased percentages of CD4+ T cells may prohibit immune system. These results also elucidated the immunosuppressive mechanism of MMF, to some extent.

Some limitations in this work should not be ignored. First, only a simplex metabonomic analysis was employed in this study and more analytical approaches may be needed. Second, performing GC-MS metabolomics on CD4+ T cells and CD8+ T cells is beneficial to further shed light on the potential mechanism of action of MMF on immune cells.

Conclusion

In a word, the current study analyzed the metabolic profiling of serum samples and organ tissues in mice following MMF treatment, indicating that MMF challenge may not only induce neurotoxicity in hippocampus and cortex, but also cause toxicity in other organs. More importantly, MMF administration may exert immunosuppressive functions via decreasing CD4+ T cell percentages. These findings provide preliminary insights into metabolomic identification after MMF administration.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Ethics statement

All experimental procedures conformed to the Guidelines for the Use of Laboratory Animals, and approved by the Ethical Committee for Animal Experimentation of Jining First People’s Hospital (Approval No. JNRM-2023-DW-021). The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

TZ: Funding acquisition, Investigation, Visualization, Writing–original draft, Writing–review and editing. YZ: Investigation, Visualization, Writing–original draft, Writing–review and editing. HC: Investigation, Writing–review and editing. WS: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing–review and editing. YG: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing–review and editing.

Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by Key R&D Program of Jining (No. 2021YXNS042), Doctoral Research Startup Foundation of Jining No.1 People’s Hospital (No. 2019003), Shandong Medicine and Health Science Technology Development Plan (No. 202203040532), the Projects of medical and health technology development program in Shandong province (No. 202113050502), the Key R&D Program of Jining (No. 2022YXNS118) and Shandong Medicine and Health Science Technology Development Plan (No. 2018WS475).

Acknowledgments

We acknowledge the strong support of Jining First People’s Hospital and Shandong University of Traditional Chinese Medicine.

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/fmolb.2024.1332090/full#supplementary-material

Abbreviations

MMF, Mycophenolate mofetil; GC-MS, Gas chromatography-mass spectrometry; DMSO, dimethyl sulfoxide; KEGG, Kyoto Encyclopedia of Genes and Genomes; MPA, mycophenolic acid; GI, gastrointestinal; LPS, lipopolysaccharide; QC, quality control; TICs, total ion chromatograms; PCA, principal components analysis; OPLS-DA, orthogonal partial least squares discrimination analysis; NK, Natural Killer; G-MDSCs, granulocytic-myeloid-derived suppressor cells; M-MDSCs, monocytic-myeloid-derived suppressor cells; GSH, Glutathione.

References

Allison, A. C., and Eugui, E. M. (1993). The design and development of an immunosuppressive drug, mycophenolate mofetil. Springer Semin. Immunopathol. 14 (4), 353–380. doi:10.1007/BF00192309

PubMed Abstract | CrossRef Full Text | Google Scholar

Augustin, R. (2010). The protein family of glucose transport facilitators: it's not only about glucose after all. IUBMB Life 62 (5), 315–333. doi:10.1002/iub.315

PubMed Abstract | CrossRef Full Text | Google Scholar

Balal, M., Demir, E., Paydas, S., Sertdemir, Y., and Erken, U. (2005). Uncommon side effect of MMF in renal transplant recipients. Ren. Fail. 27 (5), 591–594. doi:10.1080/08860220500200171

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhattacharya, S., Stoleru, G., Patel, P., Abutaleb, A., Stashek, K., and Cross, R. K. (2022). Characterization of mycophenolate mofetil gastrointestinal toxicity and risk factors for severe disease and poor prognosis. Inflamm. Bowel Dis. 28 (5), 811–814. doi:10.1093/ibd/izab254

PubMed Abstract | CrossRef Full Text | Google Scholar

Calmet, F. H., Yarur, A. J., Pukazhendhi, G., Ahmad, J., and Bhamidimarri, K. R. (2015). Endoscopic and histological features of mycophenolate mofetil colitis in patients after solid organ transplantation. Ann. gastroenterology Q. Publ. Hellenic Soc. Gastroenterology 28 (3), 366–373.

PubMed Abstract | Google Scholar

Carlone, G., Simeone, R., Baraldo, M., Maestro, A., Zanon, D., Barbi, E., et al. (2021a). Area-under-the-curve-based mycophenolate mofetil dosage may contribute to decrease the incidence of graft-versus-host disease after allogeneic hematopoietic cell transplantation in pediatric patients. J. Clin. Med. 10 (3), 406. doi:10.3390/jcm10030406

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlone, G., Simeone, R., Baraldo, M., Maestro, A., Zanon, D., Barbi, E., et al. (2021b). Area-under-the-Curve-Based mycophenolate mofetil dosage may contribute to decrease the incidence of graft-versus-host disease after allogeneic hematopoietic cell transplantation in pediatric patients. J. Clin. Med. 10 (3), 406. doi:10.3390/jcm10030406

PubMed Abstract | CrossRef Full Text | Google Scholar

Caterino, M., Ruoppolo, M., Villani, G. R. D., Marchese, E., Costanzo, M., Sotgiu, G., et al. (2020). Influence of sex on urinary organic acids: a cross-sectional study in children. Int. J. Mol. Sci. 21 (2), 582. doi:10.3390/ijms21020582

PubMed Abstract | CrossRef Full Text | Google Scholar

Chichlowski, M., German, J. B., Lebrilla, C. B., and Mills, D. A. (2011). The influence of milk oligosaccharides on microbiota of infants: opportunities for formulas. Annu. Rev. Food Sci. Technol. 2, 331–351. doi:10.1146/annurev-food-022510-133743

PubMed Abstract | CrossRef Full Text | Google Scholar

Coelho, A. I., Berry, G. T., and Rubio-Gozalbo, M. E. (2015). Galactose metabolism and health. Curr. Opin. Clin. Nutr. Metab. Care 18 (4), 422–427. doi:10.1097/MCO.0000000000000189

PubMed Abstract | CrossRef Full Text | Google Scholar

Danesi, R., and Del Tacca, M. (2004). Hematologic toxicity of immunosuppressive treatment. Transplant. Proc. 36, 703–704. doi:10.1016/j.transproceed.2004.03.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Dasarathy, S., Kasumov, T., Edmison, J. M., Gruca, L. L., Bennett, C., Duenas, C., et al. (2009). Glycine and urea kinetics in nonalcoholic steatohepatitis in human: effect of intralipid infusion. Am. J. Physiol. Gastrointest. Liver Physiol. 297 (3), G567–G575. doi:10.1152/ajpgi.00042.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Derle, E., Kibaroglu, S., Ocal, R., Kirnap, M., Can, U., Benli, S., et al. (2015). Neurologic complications after liver transplant: experience at a single center. Exp. Clin. Transpl. 13 (1), 327–330. doi:10.6002/ect.mesot2014.p177

PubMed Abstract | CrossRef Full Text | Google Scholar

Flannigan, K. L., Taylor, M. R., Pereira, S. K., Rodriguez-Arguello, J., Moffat, A. W., Alston, L., et al. (2018). An intact microbiota is required for the gastrointestinal toxicity of the immunosuppressant mycophenolate mofetil. J. Heart Lung Transplant. 37 (9), 1047–1059. doi:10.1016/j.healun.2018.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Geiger, R., Rieckmann, J. C., Wolf, T., Basso, C., Feng, Y., Fuhrer, T., et al. (2016). L-arginine modulates T cell metabolism and enhances survival and anti-tumor activity. Cell 167 (3), 829–842. doi:10.1016/j.cell.2016.09.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Heischmann, S., Dzieciatkowska, M., Hansen, K., Leibfritz, D., and Christians, U. (2017). The immunosuppressant mycophenolic acid alters nucleotide and lipid metabolism in an intestinal cell model. Sci. Rep. 7, 45088. doi:10.1038/srep45088

PubMed Abstract | CrossRef Full Text | Google Scholar

Iqbal, A., Zaman, M., Wahab Amjad, M., Adnan, S., Abdul Ghafoor Raja, M., Haider Rizvi, S. F., et al. (2020). Solid lipid nanoparticles of mycophenolate mofetil: an attempt to control the release of an immunosuppressant. Int. J. Nanomedicine 15, 5603–5612. doi:10.2147/IJN.S255636

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaminska, J., Sobiak, J., Suliburska, J. M., Duda, G., Glyda, M., Krejpcio, Z., et al. (2012). Effect of mycophenolate mofetil on plasma bioelements in renal transplant recipients. Biol. Trace Elem. Res. 145 (2), 136–143. doi:10.1007/s12011-011-9178-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurmi, K., and Haigis, M. C. (2020). Nitrogen metabolism in cancer and immunity. Trends Cell Biol. 30 (5), 408–424. doi:10.1016/j.tcb.2020.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. J., Gravelat, F. N., Cerone, R. P., Baptista, S. D., Campoli, P. V., Choe, S. I., et al. (2014). Overlapping and distinct roles of Aspergillus fumigatus UDP-glucose 4-epimerases in galactose metabolism and the synthesis of galactose-containing cell wall polysaccharides. J. Biol. Chem. 289 (3), 1243–1256. doi:10.1074/jbc.M113.522516

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D., Huang, J., Gao, S., Jin, H., and He, J. (2022). A temporo-spatial pharmacometabolomics method to characterize pharmacokinetics and pharmacodynamics in the brain microregions by using ambient mass spectrometry imaging. Acta Pharm. Sin. B 12 (8), 3341–3353. doi:10.1016/j.apsb.2022.03.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Moiseev, I. S., Pirogova, O. V., Alyanski, A. L., Babenko, E. V., Gindina, T. L., Darskaya, E. I., et al. (2016a). Graft-versus-host disease prophylaxis in unrelated peripheral blood stem cell transplantation with post-transplantation cyclophosphamide, tacrolimus, and mycophenolate mofetil. Biol. Blood Marrow Transplant. 22 (6), 1037–1042. doi:10.1016/j.bbmt.2016.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Moiseev, I. S., Pirogova, O. V., Alyanski, A. L., Babenko, E. V., Gindina, T. L., Darskaya, E. I., et al. (2016b). Graft-versus-Host disease prophylaxis in unrelated peripheral blood stem cell transplantation with post-transplantation cyclophosphamide, tacrolimus, and mycophenolate mofetil. Biol. Blood Marrow Transpl. 22 (6), 1037–1042. doi:10.1016/j.bbmt.2016.03.004

CrossRef Full Text | Google Scholar

Mothet, J. P., Pollegioni, L., Ouanounou, G., Martineau, M., Fossier, P., and Baux, G. (2005). Glutamate receptor activation triggers a calcium-dependent and SNARE protein-dependent release of the gliotransmitter D-serine. Proc. Natl. Acad. Sci. U. S. A. 102 (15), 5606–5611. doi:10.1073/pnas.0408483102

PubMed Abstract | CrossRef Full Text | Google Scholar

Nogueras, F., Espinosa, M., Mansilla, A., Torres, J., Cabrera, M., and Martin-Vivaldi, R. (2005). Mycophenolate mofetil–induced neutropenia in liver transplantation. Transplant. Proc. 37, 1509–1511. doi:10.1016/j.transproceed.2005.02.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Novelli, G., and Reichardt, J. K. (2000). Molecular basis of disorders of human galactose metabolism: past, present, and future. Mol. Genet. Metab. 71 (1-2), 62–65. doi:10.1006/mgme.2000.3073

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Sullivan, D., Sanin, D. E., Pearce, E. J., and Pearce, E. L. (2019). Metabolic interventions in the immune response to cancer. Nat. Rev. Immunol. 19 (5), 324–335. doi:10.1038/s41577-019-0140-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pellerin, D., Singh, K., Maniatis, T., Chalk, C. H., and Green, L. (2018). Mycophenolate mofetil-induced status epilepticus. Can. J. Neurol. Sci. 45 (5), 585–587. doi:10.1017/cjn.2018.326

PubMed Abstract | CrossRef Full Text | Google Scholar

Rajendra, S., Lynch, J. W., and Schofield, P. R. (1997). The glycine receptor. Pharmacol. Ther. 73 (2), 121–146. doi:10.1016/s0163-7258(96)00163-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Razak, M. A., Begum, P. S., Viswanath, B., and Rajagopal, S. (2017). Multifarious beneficial effect of nonessential amino acid, Glycine: a review. Oxid. Med. Cell Longev. 2017, 1716701. doi:10.1155/2017/1716701

PubMed Abstract | CrossRef Full Text | Google Scholar

Roser, M., Josic, D., Kontou, M., Mosetter, K., Maurer, P., and Reutter, W. (2009). Metabolism of galactose in the brain and liver of rats and its conversion into glutamate and other amino acids. J. Neural Transm. (Vienna) 116 (2), 131–139. doi:10.1007/s00702-008-0166-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Shipkova, M., Armstrong, V. W., Oellerich, M., and Wieland, E. (2005). Mycophenolate mofetil in organ transplantation: focus on metabolism, safety and tolerability. Expert Opin. Drug Metab. Toxicol. 1 (3), 505–526. doi:10.1517/17425255.1.3.505

PubMed Abstract | CrossRef Full Text | Google Scholar

Sousa, A. E., Carneiro, J., Meier-Schellersheim, M., Grossman, Z., and Victorino, R. M. (2002). CD4 T cell depletion is linked directly to immune activation in the pathogenesis of HIV-1 and HIV-2 but only indirectly to the viral load. J. Immunol. 169 (6), 3400–3406. doi:10.4049/jimmunol.169.6.3400

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivas, T. R., Kaplan, B., and Meier-Kriesche, H.-U. (2003a). Mycophenolate mofetil in solid-organ transplantation. Expert Opin. Pharmacother. 4 (12), 2325–2345. doi:10.1517/14656566.4.12.2325

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivas, T. R., Kaplan, B., and Meier-Kriesche, H. U. (2003b). Mycophenolate mofetil in solid-organ transplantation. Expert Opin. Pharmacother. 4 (12), 2325–2345. doi:10.1517/14656566.4.12.2325

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulinski, T., Dubourg, L., Saïd, M. H., Parchoux, B., Ranchin, B., and Cochat, P. (2005). Switch from cyclosporine A to mycophenolate mofetil in nephrotic children. Pediatr. Nephrol. 20 (4), 482–485. doi:10.1007/s00467-004-1778-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Varnell, C. D., Fukuda, T., Kirby, C. L., Martin, L. J., Warshaw, B. L., Patel, H. P., et al. (2017). Mycophenolate mofetil-related leukopenia in children and young adults following kidney transplantation: influence of genes and drugs. Pediatr. Transplant. 21 (7), e13033. doi:10.1111/petr.13033

PubMed Abstract | CrossRef Full Text | Google Scholar

Warren, M., Mitsinikos, T., Yanni, G., Sasaki, M., Sasaki, A. T., and Thomas, D. (2021). Mycophenolate mofetil hepatotoxicity associated with mitochondrial abnormality in liver transplant recipients and mice. J. Pediatr. Gastroenterology Nutr. 73 (4), 463–470. doi:10.1097/MPG.0000000000003171

CrossRef Full Text | Google Scholar

Wright, E. M. (2013). Glucose transport families SLC5 and SLC50. Mol. Asp. Med. 34 (2-3), 183–196. doi:10.1016/j.mam.2012.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J. C. (1994). Mycophenolate mofetil: molecular mechanisms of action. Perspect. drug Discov. Des. 2 (1), 185–204. doi:10.1007/bf02171743

CrossRef Full Text | Google Scholar

Zhang, J., Kaneko, G., Sun, J., Wang, G., Xie, J., Tian, J., et al. (2021). Key factors affecting the flesh flavor quality and the nutritional value of grass carp in four culture modes. Foods 10 (9), 2075. doi:10.3390/foods10092075

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, Z., Wheeler, M. D., Li, X., Froh, M., Schemmer, P., Yin, M., et al. (2003). L-Glycine: a novel antiinflammatory, immunomodulatory, and cytoprotective agent. Curr. Opin. Clin. Nutr. Metab. Care 6 (2), 229–240. doi:10.1097/00075197-200303000-00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: mycophenolate mofetil, toxicity, gas chromatography-mass spectrometry, multiparameter flow cytometry analysis, lymphocytes

Citation: Zhao T, Zhao Y, Chen H, Sun W and Guan Y (2024) A GC–MS-based untargeted metabolomics approach for comprehensive metabolic profiling of mycophenolate mofetil-induced toxicity in mice. Front. Mol. Biosci. 11:1332090. doi: 10.3389/fmolb.2024.1332090

Received: 02 November 2023; Accepted: 08 February 2024;
Published: 07 March 2024.

Edited by:

Xavier Domingo-Almenara, Eurecat, Spain

Reviewed by:

Umakant Sahu, Northwestern University, United States
Alma L. Saucedo, National Council of Science and Technology (CONACYT), Mexico

Copyright © 2024 Zhao, Zhao, Chen, Sun and Guan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Wenxue Sun, sunwenxue121@126.com; Yun Guan, gyuanun1989@163.com

These authors have contributed equally to this work

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