Abstract Background Visual fatigue, commonly attributed to excessive eye use or dry conditions, is traditionally associated with deficiencies in liver and kidney yin in Chinese medicine. However, its metabolic aspects remain largely unexplored. Methods Levothyroxine sodium combined with tail-clip stimulation induced a rat model of visual fatigue with liver and kidney yin deficiency. At 3 (M1), 7 (M2) and 14 (M3) days after induction, histopathological changes were observed, and metabolic profiling was completed using untargeted UHPLC-Q-Exactive Orbitrap/MS. Results The rats exhibited signs of liver and kidney yin deficiency and visual fatigue on days 3 and 7, respectively. Compared to the control group, we identified 127 and 96 differential metabolites in the serum on days 7 and 14, respectively, primarily lipids and organic nitrogen compounds. Moreover, we observed disruptions in sphingolipid metabolism and signaling pathways. Conclusion This study enriches our understanding of the metabolic profile associated with liver-kidney-yin deficiency type visual fatigue. Keywords: visual fatigue, liver and kidney yin deficiency, UHPLC-Q-Exactive Orbitrap/MS, metabolomics, signaling pathways 1. Introduction Visual fatigue is caused by a variety of factors, including excessive eye use and environmental or pathological factors, which are mainly manifested as frequent eye discomfort, such as ocular acid, dryness, and blurred vision ([37]1, [38]2). It is common in both children and adults. Statistically, the prevalence of visual fatigue in children (< 18 years old) ranges from 12.4–32.2% ([39]3, [40]4), and that in college students worldwide is 46–71% ([41]5). In addition, the number of people suffering from visual fatigue is increasing annually. These people come from all age groups, especially those who frequently use digital devices such as computers ([42]6). The dangers of visual fatigue are obvious, which not only affects the state of study and work but may also cause the progression of diseases such as dry eye and video terminal syndrome ([43]7). Therefore, considering the impact of visual fatigue on physical and mental health, the identification of methods to alleviate visual fatigue is necessitated. To improve clinical prevention and treatment strategies for visual fatigue, we first need to understand the exact mechanism underlying visual fatigue. According to traditional Chinese medicine (TCM) theory, visual fatigue is related to the liver, kidneys, and lungs ([44]8). Its etiology involves multiple aspects such as liver-kidney yin deficiency, qi blood deficiency, and stomach weakness ([45]9). Zeng et al. classified visual fatigue into three types according to the dialectics of TCM–heart-blood deficiency, liver-kidney deficiency, and qi-blood deficiency, of which liver-kidney yin deficiency is the most common ([46]10). Liver and kidney yin deficiency may cause eye discomfort such as dry eyes, tired eyes, and blurry vision ([47]11). However, the etiology and biochemical basis of visual fatigue with liver-kidney yin deficiency remain unknown. Understanding the effects of liver and kidney yin deficiency on visual fatigue is expected to provide valuable information for clinical applications. In traditional and modern medicine, metabolomics is an important technology for studying the networks of biological systems. Metabolomics describes the metabolic profiles of endogenous substances and their dynamic changes to characterize different physiological and pathological states of an organism ([48]12). Recent studies have emphasized the application of metabolomics in understanding the pathophysiological mechanisms of liver diseases ([49]13, [50]14) and kidney diseases ([51]15, [52]16), advancing precision medicine and patient care. However, traditional research methods often fail to effectively clarify the potential of TCM therapeutic interventions, hindering the modern exploration of TCM-derived functional compounds. Metabolomics can accurately capture the molecular interactions among disease-responsive metabolites, thereby characterizing the systemic fluctuations and molecular complexity of TCM-derived functional compounds ([53]17, [54]18). Clinical examinations have revealed that individuals with symptoms of liver and kidney yin deficiency exhibit enhanced energy metabolism and disturbed free radical metabolism ([55]19). Moreover, micronutrients affect the morphology and function of eye tissues by participating in multiple cellular metabolic pathways to maintain the body’s homeostatic balance ([56]20). Therefore, the identification of metabolite changes during the progression of visual fatigue can provide useful information for improving eye health. However, the metabolites and metabolic pathways involved in visual fatigue have not yet been identified. High-resolution mass spectrometry (MS), especially UHPLC coupled with a Q-Exactive orbitrap/MS, is widely used to analyze the chemical composition of biological samples and herbal plants ([57]21). The advantages of this technique include fast polarity switching and high sensitivity. Therefore, in this study, we utilized levothyroxine sodium combined with a tail clamp and light/dark environment manipulation to create a rat model of visual fatigue with liver and kidney yin deficiency. Then, through untargeted UHPLC-Q-Exactive orbitrap/MS metabolomics analysis identified differential metabolites and perturbed pathways closely related to visual fatigue progression, especially in sphingolipid metabolism and sphingolipid signaling pathway. The obtained data can help to better understand the metabolic characteristics associated with visual fatigue and, reveal the link between sphingolipid metabolic disorders and visual fatigue in the context of liver and kidney yin deficiency. This offers a new perspective for modern neuro-ophthalmology, metabolically corroborating the TCM theory that “liver and kidney yin deficiency leads to inadequate eye nutrition” ([58]11, [59]22), shows the complementarity of TCM and Western medicine, thereby providing a new scientific basis for the targeted treatment of this disease. 2. Methods 2.1. Animal preparation Sprague-Dawley (SD) male rats (specific pathogen-free (SPF) grade, license number SCXK (Jing) 2021–0011), weighting 200 ± 20 g, 6–7 weeks old, were acquired from Beijing Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China) and housed in an SPF laboratory (license number SYXK (Shan) 2018–0009). Before performing the experiments, all animals were acclimatized under controlled conditions (20–25°C, 45–65% relative humidity, 12-h dark/light cycle) for five days, with free access to water and food. The research protocol was in accordance with the ethical requirements and authorized by the Animal Ethics Committee (ethics approval number 20220913). All methods were carried out in accordance with relevant guidelines and regulations and are reported in accordance with ARRIVE guidelines ([60]https://arriveguidelines.org) for the reporting of animal experiments. 2.2. Induction of visual fatigue in rats and grouping The rats were randomly assigned to four groups, namely normal control (K1) and treatment for 3 days (M1), 7 days (M2), and 14 days (M3), with 10 rats per group. Visual fatigue was induced using levothyroxine sodium gavage combined with tail-clamping stimulation. Briefly, rats in the model groups (M1, M2, and M3) were administered levothyroxine sodium (90 g/kg) by gavage in the morning at 1 mL/100 g. Afterward, the tails of the animals were clamped with a clip (approximately 3 cm from the tip of the tail) for 5 min once a day for effective stimulation to induce the symptoms of liver and kidney yin deficiency. The cages of the rats in the model group were covered with black cloth to simulate a dark environment, and an LED light was placed on the cages as the only light source. The light source was 10–15 cm away from the rat, with an intensity of 3000–500 lx. During the modeling period, the cage was knocked every 15 min to stimulate the rats to keep their eyes open for 2 h per day, and water and food were forbidden. 2.3. Sample collection and preparation At the end of light exposure on day 15, the rats were anesthetized intraperitoneally with 20% urethane and blood was harvested from the abdominal aorta. After standing for 2–4 min, the blood samples were centrifuged at 3000 r/min for 10 min, and the supernatants were isolated and stored at -80°C. Tissues, including the spleen, liver, lens, ciliary body, and eyelids, were rapidly excised. After washing with normal saline, the tissues were ground into a homogenate with a grinding rod in a cold environment and centrifuged at 4°C and 12,000 rpm for 10 min. Finally, supernatants were harvested and stored at -20°C until assayed. 2.4. Detection of biochemical parameters According to the instructions of each kit, the levels of the following indicators were quantified: cyclic adenosine monophosphate (cAMP), cyclic guanosine monophosphate (cGMP), glutathione peroxidase (GSH-Px), total antioxidant capacity (T-AOC), and Ca(2+)/Mg(2+)-ATPase in the lens and ciliary body; aquaporin 4 (AQP4), aquaporin 5 (AQP5), and lactate dehydrogenase (LDH) in the tarsal gland and palpebral muscle; Na(+)/K(+)-ATPase and Ca(2+)/Mg(2+)-ATPase in the liver and cAMP, cGMP, estradiol (E[2]), testosterone (T), corticosterone (CORT), triiodothyronine (T3), thyroxine (T4), interleukin 2 (IL-2), interleukin 6 (IL-6), superoxide dismutase (SOD), and malondialdehyde (MDA). 2.5. Untargeted LC-MS/MS metabolomics analysis 2.5.1. Serum sample preprocessing Serum samples were thawed at 4°C and cold methanol/acetonitrile/water solution was added in the ratio of 2:2:1 (v/v). The mixed samples were vortexed and sonicated for 30 min, followed by incubation at -20°C for 10 min. After centrifugation at 14,000 × g and 4°C for 20 min, the supernatant was harvested and dried under vacuum. For MS analysis, 100 μL of acetonitrile/water solution (1:1, v/v) was added to the sample for re-dissolution, vortexed for 30 min, and centrifuged at 14,000 × g for 15 min. Finally, the supernatant was collected for further analysis. 2.5.2. UHPLC-Q-Exactive orbitrap MS assay MS separation of the samples was accomplished using a Vanquish LC UHPLC system coupled with a HILIC column. The column temperature was set at 25 °C. Detection was carried out using a binary solvent system with mobile phases A (water containing 25 mM ammonium acetate and 25 mM ammonia) and B (acetonitrile). The specific gradient elution parameters were as follows: 0 to 1.5 min, 98% B; 1.5 to 12 min, 98% B linearly declined to 2%; 12 to 14 min, 2% B; 14 to 14.1 min, 2% B straightly elevated to 98%; 14.1 to 17 min, 98% B. During the analysis period, the temperature of the autoinjector was maintained at 4°C, and the injection volume was 2 μL with a flow rate of 0.3 mL/min. The MS1 and MS2 spectra were captured using a Q Exactive mass spectrometer, and the chromatographic data were analyzed and detected in the positive (pos) and negative (neg) ion modes of electrospray ionization (ESI). The specific parameters of the ESI source and mass spectrum were as follows: gas 1 and gas 2 were 60; CUR was 30 psi; ion source temperature was 600 °C; spray voltage of positive or negative was 5,500 V; the range of primary mass-to-charge ratio detection was 80–1200 Da, the resolution/secondary resolution was 60,000/30,000, and the cumulative scanning time was 50 ms. 2.6. Data preprocessing and analysis The collected raw data were converted to.mzXML format using ProteoWizard software ([61]23), and XCMS software ([62]24) was utilized for data processing operations, such as peak alignment, retention time correction, and peak area extraction. After matching the raw peak area information with the total ion flow profile, data analysis was conducted. First, the number and classification of the metabolites were identified using the in-house database (Shanghai Applied Protein Technology). Next, multivariate statistical analyses were performed to characterize metabolic perturbations and differences among the different groups, including principal component analysis (PCA) and orthogonal partial least squares-discriminant analysis (OPLS-DA). The variable importance for projection (VIP) was calculated based on the OPLS-DA model, which was applied to mine the differential metabolites with biological significance. The selection criterion was VIP > 1 (P < 0.05). Finally, to reveal the biological functions of the metabolites, KEGG pathway enrichment analysis was conducted for differential metabolites. Differences between groups were measured using Fisher’s Exact Test, and entries with P < 0.05 were selected. 2.7. Statistical analysis Results are expressed as mean ± SEM and were analyzed as well as plotted by GraphPad Prism (version 7.0). Outcomes between the groups were compared using a one-way analysis of variance (ANOVA) and Dunnett’s test. Differences were considered statistically significant at P < 0.05. 3. Results 3.1. Rats with visual fatigue show symptoms of liver and kidney yin deficiency Throughout the animal experiments (0–14 d), we monitored the body weight changes in all rats and found that the weight of the rats in each group increased slowly, with no significant differences among the four groups ([63] Figure 1A ). There were no significant differences in the thymus and spleen indices between the K1 and the three model groups ([64] Supplementary Figures S1A, S1B ). Regarding the liver energy metabolism enzyme activities of rats, compared to the K1 group, the model groups displayed increased Ca(2+)/Mg(2+) ATPase and a decline in Na(+)/K(+) ATPase levels, especially in the M2 group ([65] Figures 1B, C ). We measured key parameters related to the crystalline lenses and tarsal glands of the rats. The results showed that in the model groups, the levels of GSH-Px, AQP5, cAMP, and cGMP in the lens were elevated, whereas those of T-AOC and Ca(2+)/Mg(2+) ATP were decreased; significant changes in most indicators were observed in the M2 and M3 groups ([66] Figures 1D–I ). However, the levels of AQP4 and AQP5 in the tarsal gland and LDH in the palpebral muscle did not significantly differ between the normal and model groups ([67] Figures 1J–L ). Figure 1. [68]Line and box plot graphs showcasing various physiological measurements over time, comparing groups K1, M1, M2, and M3. Graphs A to L display metrics like body weight, enzyme activities, and biochemical levels, indicating statistical differences with symbols such as asterisks and hashtags for specific comparisons within groups. [69]Open in a new tab Biochemical and physiological alterations in rats with visual fatigue (n = 10). (A) Changes in body weight from day 1 to 14 after visual fatigue induction. (B) Liver Ca(2+)/Mg(2+) ATP levels in rats of each group. (C) Level of liver Na(+)/K(+) ATP of each group. (D-I) Disordered lens-related indexes of GSH - Px (D), T-AOC (E), AQP5 (F), cAMP (G), cGMP (H), and Ca(2+)/Mg(2+) ATP (I) in rats with visual fatigue. (J–L) Tarsal gland AQP4 (J) and AQP5 (K) as well as LDH (L) levels in palpebral muscles in rats of each group. Data are shown as the mean ± SEM. K1, normal control group; M1, 3-day model group; M2, 7-day model group; M3, 14-day model group. *P < 0.05 vs. K1; ^# P < 0.05 vs. M1; ^& P < 0.05 vs. M2. Moreover, the levels of inflammatory factors, AQP, antioxidants, and hormones in the serum of rats are important indicators of visual fatigue. As shown in [70]Figures 2A–I , the serum of the model groups presented high levels of IL-6, IL-2, AQP4, AQ5, cAMP, and cGMP, while exhibiting low levels of SOD and cAMP/cGMP. However, no significant differences in MDA levels were observed among the four groups. Compared to the normal (K1) group, the levels of serum hormones (E[2], T, CORT, T3, and T4) increased significantly in rats from the three model groups ([71] Figures 3A–E ). Overall, these findings indicate that rats in the M2 and M3 groups exhibited both liver and kidney yin deficiency and visual fatigue symptoms, confirming the success of our model construction. The rats in these two groups were selected for further analysis. Figure 2. [72]Box plot graphs showing serum levels of various biomarkers across groups labeled K1, M1, M2, and M3. Panels A to I display different biomarkers: IL-6, IL-2, AQP4, AQP5, SOD, MDA, cAMP, cGMP, and cAMP/cGMP ratio. In panels A, B, D, G, and H, the M groups (M1, M2, M3) show higher median levels than K1. Panel E shows SOD levels decreasing from K1 to M3, with the most significant decrease in M3. Panel F shows MDA levels similar across all groups. Panel I shows cAMP/cGMP ratio decreasing from K1 to M3. [73]Open in a new tab Abnormal blood inflammatory factors and antioxidant activity are observed in rats with visual fatigue (n = 10). (A) IL-6. (B) IL-2. (C) AQP4. (D) AQP5. (E) SOD. (F) MDA. (G) cAMP. (H) cGMP. (I) cAMP/cGMP. Data are shown as the mean ± SEM. K1, normal control group; M1, 3-day model group; M2, 7-day model group; M3, 14-day model group. *P < 0.05 vs. K1; ^# P < 0.05 vs. M1. Figure 3. [74]Box plots depicting serum levels in different groups: (A) Serum E2 levels, (B) Serum T levels, (C) Serum CORT levels, (D) Serum T3 levels, (E) Serum T4 levels. Each plot compares levels between group K1 and groups M1, M2, and M3. Significant differences are indicated by asterisks and hash symbols. [75]Open in a new tab Elevated serum hormones in rats with visual fatigue (n = 10). (A) E[2]. (B) T. (C) CORT. (D) T3. (E) T4. Data are shown as the mean ± SEM. K1, normal control group; M1, 3-day model group; M2, 7-day model group; M3, 14-day model group. *P < 0.05 vs. K1; ^# P < 0.05 vs. M1. 3.2. Serum metabolic profiles In the positive and negative ion modes, the spectra of the total ion chromatogram showed that the intensities and retention times of the chromatographic peaks overlapped, indicating that instrument error was negligible ([76] Supplementary Figure S2A ). PCA revealed that the QC samples in the two modes were closely clustered together, indicating good replicability ([77] Supplementary Figure S2B ). Overall, the stability, precision, and reproducibility of the UHPLC-Q Exactive Orbitrap MS method were excellent. 3.3. Metabolite determination and multivariate statistical analysis Statistically, 551 (neg) and 695 (pos) metabolites were detected. For M2 vs. K1, an examination of the PCA data revealed that the K1 and M2 samples could be clearly separated in both the pos and neg ion modes ([78] Figure 4A ). An OPLS-DA model was established to further verify the differences between the two groups. The validity of the model was examined using a permutation test. The R2Y/Q2 values for the Pos ions were 0.978/0.522 and 0.989/0.555 for the Neg ion, suggesting that the model was stable and reliable ([79] Figure 4B ). Similarly, in the comparison between the M3 and K1 groups, we observed a significant separation between the two groups ([80] Figure 5A ). The R2Y and Q2 parameters of the OPLS-DA model under the pos and neg ion modes were satisfactory (pos: R2Y, 0.996; Q2, 0.462; neg: R2Y, 0.943; Q2, 0.106), implying that the model was not overfitted and that the results were reliable ([81] Figure 5B ). Figure 4. [82]Panel A shows two 3D scatter plots for ESI positive and negative ion modes, with data points colored purple and green representing K1 and M2. Axes are labeled PC1, PC2, and PC3 with varying percentages. Panel B includes two permutation plots displaying R2 and Q2 values for ESI positive and negative, indicating model validation across 200 permutations. [83]Open in a new tab Multivariate statistical analysis for K1 and M2 groups. (A) PCA score plot in the positive (left panel) and negative (right panel) modes. (B) Permutation test plot of the OPLS-DA model in positive (left panel) and negative (right panel) modes, R^2 is the explained variance, and Q^2 is the predictive ability of the model. PC, principal component; ESI, electrospray ionization. Figure 5. [84]Panel A displays two 3D scatter plots labeled ESI (+) and ESI (-), showing the distribution of variables K1 and M3 with axes PC1, PC2, and PC3. Panel B shows two permutation plots with R2 and Q2 values for ESI (+) and ESI (-), indicating statistical validation of models over 200 permutations. [85]Open in a new tab Multivariate statistical analysis for K1 and M3 groups. (A) PCA score plot in the positive (left panel) and negative (right panel) modes. (B) Permutation test plot of the OPLS-DA model in positive (left panel) and negative (right panel) modes, R^2 is the explained variance, and Q^2 is the predictive ability of the model. PC, principal component; ESI, electrospray ionization. 3.4. Differential metabolite identification At VIP > 1 and P < 0.05, 127 differential metabolites were identified between the K1 and M2 groups, including 62 in Pos and 65 in Neg ([86] Supplementary Tables S1 , [87]S2 ). In the pos mode, these metabolites were clearly assigned to eight categories, most of which belonged to lipids and lipid-like molecules (21/62, 33.8%), followed by organic nitrogen compounds (9/62, 14.5%) and organic acids and derivatives (9/62, 14.5%). The relative abundances of fingolimod and phytosphingosine were significantly upregulated, whereas that of acebutolol, creatine, and 4-hydroxy-3-methoxybenzyl alcohol were markedly downregulated in the M2 group ([88] Figure 6A ). In the neg mode, these metabolites belonged to eight classes, primarily lipid-like molecules (23/65, 35.4%) and organic acids and derivatives (13/65, 20%). Notably, in the M2 group, the abundances of 6-phosphogluconic acid and isopentenyl pyrophosphate increased, whereas those of Psoralidin and Urocanate declined ([89] Figure 6B ). Figure 6. [90]Two clustered heatmaps labeled A and B display metabolite levels with color gradients from blue to red, indicating intensity from low to high. Metabolites are listed on the y-axis, while samples labeled K1 and M2 groups are on the x-axis. The heatmaps show variations across samples, suggesting group-specific metabolic patterns. [91]Open in a new tab Heatmap of differential metabolites between K1 and M2 groups under the positive (A) and negative modes (B). The color red indicates higher relative expression, while purple indicates lower relative expression. K1, normal control group; M2, 7-day model group. Using the same threshold, 96 differential metabolites were identified between the K1 and M3 groups, with 56 in the Pos and 40 in the Neg modes ([92] Supplementary Tables S3 , [93]S4 ). In Pos mode, these metabolites were explicitly divided into nine categories– lipids and lipid-like molecules (15/56, 26.8%), organic nitrogen compounds (10/56, 17.9%), and organoheterocyclic compounds (7/56, 12.5%). Compared to the K1 group, 43 metabolites (N-acetylhistamine and fingolimod) were elevated, and 13 metabolites (4−phenylbutan−2−ol and phenylalanine) were reduced in the M3 group ([94] Figure 7A ). Moreover, in the Neg mode, the 40 metabolites were classified into seven classes– lipids and lipid-like molecules (13/40, 32.5%), organoheterocyclic compounds (6/40, 15%), and organic oxygen compounds (6/40, 15%). We observed that the abundances of 31 metabolites including 6-phosphogluconic acid and 3-ketofusidic acid, increased, whereas those of 9 metabolites, such as shikimate and eicosapentaenoic acid, decreased in the M3 group compared to the K1 group ([95] Figure 7B ). Figure 7. [96]Heatmap showing metabolite profiles labeled A and B. Each row represents a metabolite, and each column is a sample group, either K1 or M3. Color intensity ranges from blue (low values) to red (high values), indicating the relative abundance of metabolites. Hierarchical clustering on the left organizes metabolites by similarity in both panels. [97]Open in a new tab Heatmap of differential metabolites between K1 and M3 groups under the positive (A) and negative modes (B). The color red indicates higher relative expression, while purple indicates lower relative expression. K1, normal control group; M3, 14-day model group. 3.5. Functional pathway enrichment analyses To explore the potential metabolic pathways involved related to the differential metabolites, pathway enrichment analysis was conducted. As shown in [98]Figure 8A , [99]B , the differential metabolites in the K1/M2 and K1/M3 compartments mainly influenced sphingolipid metabolism and the sphingolipid signaling pathway. Figure 8. [100]Scatter plots labeled A and B show enriched KEGG pathways against rich factors. Points vary in size based on metabolite numbers and color based on -log10(p.value). A highlights pathways like sphingolipid signaling and metabolism. B focuses on sphingolipid metabolism. The color scale ranges from green to red, indicating significance levels. [101]Open in a new tab Functional pathway enrichment analyses of differential metabolites between K1 vs. M2 (A) and K1 vs. M3 (B). K1, normal control group; M2, 7-day model group; M3, 14-day model group. 4. Discussion Unlike books and printed materials, great information is transmitted through video display terminals, causing increased complains of eye discomfort and fatigue ([102]25). Common symptoms of visual fatigue include visual disturbances, ocular discomfort, and flickering of words within close range. TCM suggests that liver blood depletion and excessive use of the brain can lead to a lack of nourishment for the eyes, thus causing visual fatigue, making liver-kidney yin deficiency evidence of a common visual fatigue pattern. However, little information is available on the development of visual fatigue in the liver-kidney yin deficiency syndrome. The identification of metabolic markers of visual fatigue provides a reference for elucidating the pathological processes and contributes to the development of precise interventions. In this study, we established animal models of visual fatigue at different time points (3, 7, and 14 d). The results of biochemical indices showed that rats had disorders in the levels of sex hormones and thyroid hormones 3 days after modeling, which confirmed the appearance of liver-kidney yin deficiency. Seven days after modeling, the rats displayed an abnormal metabolism of free radicals, indicating that the visual fatigue model was successfully constructed. These results were confirmed in previous research ([103]22). During the progression of visual fatigue, 127 and 96 significantly abnormal metabolites were identified at 7 and 14 days (compared to controls) after model induction, respectively. Pathway enrichment analysis revealed that these metabolites were involved in sphingolipid metabolism and sphingolipid signaling pathways. Important causes of visual fatigue include damage to the structure or function of the ocular surface and fundus retina ([104]26). Vision originates from the retina. Damage to the structure of the cells distributed in the retina, especially the retinal pigment epithelium (RPE), is affected ([105]27). RPE cells are extremely active and require large amounts of oxygen to consume energy and nutrients; therefore, they are prone to produce reactive oxygen species (ROS) ([106]28). This process can lead to oxidative stress if ROS are not released in a timely manner ([107]29). On the other hand, continued overuse of the eyes or exposure to abnormal environments (such as bright light) increases the load on the retina and aggravates ocular oxidative stress, leading to dry eyes and eye fatigue ([108]30). Decreased SOD activity has been detected in multiple eye-related diseases, such as cataracts and glaucoma ([109]31, [110]32), and is regarded as a potential diagnostic marker for oxidative stress-related diseases ([111]33). Reduced SOD activity leads to intense lipid peroxidation ([112]34). Abnormalities in this pathway have been implicated in the pathogenesis of cataracts. Specifically, lipid peroxidation alters the internal composition and conformation of cells by affecting the permeability of the cell membrane, causing loss of protein function, and ultimately leading to cataract formation ([113]35). Moreover, Mg is a cofactor for several enzymes that maintain ionic homeostasis in the lens, such as Ca(2+)-ATPase and Na(+)/K(+) ATPase ([114]36); Mg deficiency impairs ATPase function by enhancing oxidative stress ([115]37). Thus, inhibition of Ca(2+)-ATPase function leads to the loss of calcium homeostasis and calcium accumulation within the lens, which may eventually contribute to lens opacity ([116]38). Consistent with the findings of the above studies, we observed decreased concentrations of serum SOD and lens Ca(2+)-ATPase in the model rats. Accordingly, we speculated that visual fatigue can be alleviated by suppressing oxidative stress and restoring the Ca/Mg balance. Furthermore, metabolite annotation results showed that the majority of the identified differential metabolites belonged to lipids and lipid-like molecules, organic nitrogen compounds, as well as organic acids and derivatives, such as fingolimod, phytosphingosine, sphingosine, quillaic acid, phenylalanine and creatine. These metabolites also affect sphingolipid metabolism and the sphingolipid signaling pathway during liver and kidney yin deficiency - type visual fatigue progression. With the deepening of our understanding of ophthalmic diseases, numerous studies have revealed that lipid and sphingolipid signaling pathways, as well as organic nitrogen compounds and organic acids, are involved in the pathological mechanisms of ocular inflammatory diseases ([117]39–[118]41). Among the identified differential metabolites, organic acids are a significant proportion. For example, phenylalanine is an essential amino acid and neurotransmitter precursor, and creatine plays a key role in cellular energy metabolism and antioxidant defense. This study observed reduced levels of phenylalanine and creatine in the model group. Combined with previous research reports, phenylalanine modification may reduce empty capsids to lower virus - vector - induced intraocular inflammatory responses ([119]42). Creatine supplementation may delay cone secondary degeneration in retinitis pigmentosa and protect the retina by reducing oxidative stress and inflammation ([120]43). Thus, maintaining or restoring the homeostasis of these organic acids may help delay visual fatigue progression. Lipids, especially sphingolipids, are enriched in the nervous system and retina ([121]44). Sphingolipid metabolism disorders are associated with visual dysfunction such as visual fatigue. For example, Green et al. found that patients with hereditary neuroretinas exhibit altered levels of complex sphingolipids ([122]45). This study detected abnormal expression of the organic nitrogen compound metabolite fingolimod. Fingolimod (also known as FTY720, a sphingosine analog) is used clinically as an S1PR agonist and immunomodulatory drug for the treatment of recurrent multiple sclerosis ([123]46, [124]47). Interestingly, multiple sclerosis presents symptoms such as visual impairment and fatigue ([125]48). Fingolimod/FTY720 is found to inhibit NFκB and pro-inflammatory cytokines in the retina of diabetic rats, preventing disruption of the blood-retinal barrier; in addition, FTY720 reduces neuronal loss in rats with glaucoma ([126]49, [127]50). In TCM theory, liver and kidney yin deficiency may relate to endocrine and immune disorders ([128]51). As an immunomodulatory agent, we speculate that fingolimod may participate in the pathological process of liver and kidney yin deficiency type visual fatigue by regulating ocular immune response or neuroprotective mechanisms. Despite its promising applications, fingolimod has been documented to cause possible adverse ocular effects (such as macular edema) ([129]52). Therefore, further studies are needed to confirm and evaluate the role of fingolimod in improving visual function disorders like visual fatigue. Porter et al. ([130]53) showed that the sphingosine 1-phosphate signaling pathway is implicated in many key processes in the eye, ranging from light stress/apoptotic responses to retinal/vascular development, suggesting that the sphingolipid signaling pathway serves a role in both vision and apoptosis. Another differential metabolite, quillaic acid, exhibits various biological activities, including anti-inflammatory, antiviral, and antineoplastic effects ([131]54, [132]55). It also exhibits the ability to induce apoptosis in various cell types. Although there is currently no direct evidence of quillaic acid’s role in visual fatigue, based on its known anti-inflammatory properties, we speculate that it may be beneficial in alleviating liver and kidney yin deficiency - induced visual fatigue symptoms. Taken together, the key differential metabolites and disrupted sphingolipid metabolic pathways found in this study’s liver and kidney yin deficiency visual fatigue model suggest that they may contribute to visual fatigue development by affecting biological processes such as retinal inflammation and apoptosis. This provides a new perspective on understanding metabolic characteristics of visual fatigue, innovatively establishing a link between TCM syndromes and sphingolipid metabolic disorders, revealing the complementarity between the traditional TCM syndrome classification system and modern neuro - ophthalmic cognition. However, further research is needed to explore the specific molecular mechanisms. Information on the signature metabolites associated with liver and kidney yin deficiency patterns in visual fatigue is limited, and our work fills this gap. The screened metabolites contribute to understanding the underlying causes of visual fatigue. However, this study has several limitations. First, the inherent physiological differences between rats and humans restrict the direct extrapolation of results. Second, the assessment of visual fatigue and liver-kidney yin deficiency in the animal model mainly depends on biochemical indicators. Future research should integrate more objective behavioral evaluations, such as optokinetic responses and pupil reaction, for a comprehensive visual function assessment. Moreover, the untargeted metabolomics analysis used here is exploratory. The identified differential metabolites future need quantitative validation and mechanistic exploration through targeted metabolomics (e.g., LC-MS/MS) or functional experiments. Therefore, the results presented here should be interpreted with caution. Finally, the mechanisms of sphingolipid metabolism and signaling pathways implicated in visual fatigue need to be further explored in large populations. 5. Conclusions Overall, we revealed the metabolic profiles of an animal model of visual fatigue with liver and kidney yin deficiency syndrome and identified metabolites and signaling pathways prominently associated with this disease. Sphingolipid metabolites are expected to be potential targets for alleviating visual fatigue. Further targeted metabolomics or functional studies are needed to quantitative validation our findings, confirm metabolite targets, and address animal model limitations. This study offers fresh insights into visual fatigue-related metabolic signatures, bridges TCM theory with modern neuro-ophthalmological metabolic mechanisms, and lays an initial scientific foundation for precise disease intervention. Glossary MS mass spectrometry TCM traditional Chinese medicine SD Sprague-Dawley SPF specific pathogen-free cAMP cyclic adenosine monophosphate cGMP cyclic guanosine monophosphate GSH-Px glutathione peroxidase T-AOC total antioxidant capacity AQP4 aquaporin 4 AQP5 aquaporin 5 LDH lactate dehydrogenase E2 estradiol T testosterone CORT corticosterone T3 triiodothyronine T4 thyroxine IL-2 interleukin 2 IL-6 interleukin 6 SOD superoxide dismutase MDA malondialdehyde pos positive neg negative ESI electrospray ionization PCA principal component analysis OPLS-DA orthogonal partial least squares-discriminant analysis VIP variable importance for projection ANOVA analysis of variance RPE retinal pigment epithelium ROS reactive oxygen species. Funding Statement The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the National Key R&D Program of China (grant number: 2018YFC1706800). Data availability statement The original contributions presented in the study are included in the article/[133] Supplementary Material . Further inquiries can be directed to the corresponding author/s. Ethics statement The animal study was approved by The Food and Cosmetics Inspection and Authentication Center, Shaanxi Academy of Traditional Chinese Medicine. The study was conducted in accordance with the local legislation and institutional requirements. Author contributions XZ: Conceptualization, Data curation, Formal analysis, Writing – original draft. XL: Data curation, Methodology, Resources, Writing – original draft. ZH: Data curation, Investigation, Writing – review & editing. YC: Methodology, Writing – review & editing. YZ: Conceptualization, Writing – review & editing. JZ: Resources, Writing – review & editing. JD: Formal analysis, Project administration, Resources, Writing – review & editing. LC: Conceptualization, Methodology, Project administration, Supervision, Writing – review & editing. 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. Generative AI statement The author(s) declare that no Generative AI was used in the creation of this manuscript. 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: [134]https://www.frontiersin.org/articles/10.3389/fendo.2025.1586581/fu ll#supplementary-material Supplementary Figure 1 Changes in thymus index (A) and spleen index (B) in rats in each group (n = 10). K1, normal control group; M1, 3-day model group; M2, 7-day model group; M3, 14-day model group. [135]Image1.tif^ (70.6KB, tif) Supplementary Figure 2 Peak intensity chromatograms of serum samples in positive and negative modes. (A): Total ion chromatogram. (B): Metabolic profiles of the serum samples in the PCA plot. PCA, principal component analysis; ESI, electrospray ionization. [136]Image2.tif^ (702.4KB, tif) Supplementary Table 1 List of differential metabolites between K1 and M2 groups under the positive mode. [137]Table1.xlsx^ (27.6KB, xlsx) Supplementary Table 2 List of differential metabolites between K1 and M2 groups under the negative mode. [138]Table2.xlsx^ (28.7KB, xlsx) Supplementary Table 3 List of differential metabolites between K1 and M3 groups under the positive mode. [139]Table3.xlsx^ (26KB, xlsx) Supplementary Table 4 List of differential metabolites between K1 and M3 groups under the negative mode. [140]Table4.xlsx^ (21.6KB, xlsx) References