Metabolomics and Biomarkers in Retinal and Choroidal Vascular Diseases

The retina is one of the most important structures in the eye, and the vascular health of the retina and choroid is critical to visual function. Metabolomics provides an analytical approach to endogenous small molecule metabolites in organisms, summarizes the results of “gene-environment interactions”, and is an ideal analytical tool to obtain “biomarkers” related to disease information. This study discusses the metabolic changes in neovascular diseases involving the retina and discusses the progress of the study from the perspective of metabolomics design and analysis. This study advocates a comparative strategy based on existing studies, which encompasses optimization of the performance of newly identified biomarkers and the consideration of the basis of existing studies, which facilitates quality control of newly discovered biomarkers and is recommended as an additional reference strategy for new biomarker discovery. Finally, by describing the metabolic mechanisms of retinal and choroidal neovascularization, based on the results of existing studies, this study provides potential opportunities to find new therapeutic approaches.


Introduction
Vasculogenesis and angiogenesis are two processes that are indispensable for the formation of vascular networks, and they both promote neovascularization in a specific microenvironment [1,2]. Neovascularization in the physiological state contributes to wound healing. Pathological neovascularization results when exogenous stimuli, such as hypoxia and ischemia, disrupt the regulation of the dynamic balance between endogenous proangiogenic and anti-angiogenic factors [3]. Abnormal neovascularization is associated with many diseases, including cancer, hypertension, and ocular neovascularization disease.
The retina is one of the most important structures in the eye, converting incident light into neural signals that can be received by the visual cortex of the brain, thus allowing us to receive information from the outside world. The high metabolic demands of the retina are met by a well-organized vascular system. The types of neovascularization affecting the retina include both retinal neovascularization and subchoroidal or choroidal neovascularization (CNV) [4]. Incomplete endothelial cells of neovascularization undergo spontaneous rupture, producing fundus hemorrhage, and exudation, leading to loss of central vision, visual distortion, and reduced contrast sensitivity in patients, eventually progressing to irreversible vision loss and leading to blindness. Neovascularization associated with retinal vascular disease is a comprehensive phenotype of a group of diseases, such as diabetic retinopathy (DR), retinopathy of prematurity (ROP) and neovascular age-related macular degeneration (nAMD), pathological myopia, and idiopathic CNV, and associated risk factors have been reported in different studies [5]. Frequently reported risk factors include age [6], hyperglycemia [7], diet [8], transmissibility, trauma, and inflammation [9].

Application of Metabolomics in Retinal and Choroidal Neovascularization Studies
The eye is a relatively independent organ, and the blood-ocular barrier system regulates intraocular fluid content and maintains a relatively stable state of metabolic environment in the general circulation of the body [14]. The blood-retinal barrier is a physiological barrier to metabolic exchange in the retina, divided into internal and external barriers. Retinal capillary endothelial cells limit the efflux of plasma and its components into the interstitial space of the retina, forming an internal blood-retinal barrier. Retinal pigment epithelial cells located between the choroid and the outer segments of photoreceptor cells form the outer blood-retinal barrier, which separates retinal tissue fluid from choroidal tissue fluid and has selective permeability and active transport functions. The choroidal vessels accommodate up to 85% of the blood volume in the eye, allowing plasma to bathe the retinal pigment epithelium (RPE) and meet the high metabolic demands of the photoreceptors [15].
When exogenous stimuli cause pathological neovascularization, it disrupts the barrier system and the diffusion barrier between the retinal extracellular fluid and the adjacent vitreous disappears, disrupting the local environmental homeostasis of the aqueous humor (AH) and vitreous, which also offers the possibility to study neovascular disease through ocular tissue and blood [13]. Neovascularization is a common pathological manifestation of a group of vascular diseases involving the retina, of which CNV is mainly a proliferation of neovascularization in the choroidal capillary layer that breaks through and grows across Bruch's membrane toward the retina and enters between Bruch's membrane and RPE or between the RPE layer and the retinal neurosensory layer, ultimately affecting vision, of which nAMD is the main representative [16]. Neovascularization of the retinal capillary network is represented by DR and ROP [17]. There is great etiological variation among these eye diseases, but they exhibit the same neovascularization phenomenon, which may arise as a cumulative result of gene-environment interactions. Therefore, metabolomics studies may contribute to the understanding of disease mechanisms. This section summarizes studies on the use of metabolomics platforms to reveal metabolic changes in several classes of diseases presenting with retinal and choroidal neovascularization changes, including AMD, PDR, and ROP. Metabolites that undergo changes in pathological states are summarized and the mechanisms of disease are further elucidated on the basis of pathway analysis. We emphasize the importance of technical and biological validation and clinical utility assessment, advocating a comparative results strategy based on existing studies, which covers both performance optimization of newly identified biomarkers and consideration based on the existing research base as additional optimized reference standards for the discovery of desirable biomarkers.

Metabolomics to Identify Metabolic Changes in Retinal and Choroidal Neovascular Disease
Metabolomics studies have reported several ocular diseases associated with neovascularization, including AMD, DR, and ROP, and we used these three common eye diseases as examples to develop a review of metabolic changes. We conducted a literature search using the PubMed database. Articles were identified by searching for titles using the following search terms: ("metabolomics" or "metabonomics" or "metabolome" or "metabolic profiling") AND "neovascularization-related disease*". Additional articles were identified by searching the reference lists of the included studies. We included metabolomics studies based on mass spectrometry (MS) or nuclear magnetic resonance (NMR), focusing on the generalization of population-based studies presenting with neovascularization pathological states and looking for differential metabolites, excluding drug evaluation reports, review articles, and abstracts without full text, resulting in the inclusion of twenty-six studies, including five studies of ROP, nine studies of nAMD, and twelve studies of DR. A simplified flow diagram of study selection is shown in Figure 1. Important details about the study design and methods were extracted from selected articles and summarized (Table 1). and biological validation and clinical utility assessment, advocating a comparative resu strategy based on existing studies, which covers both performance optimization of new identified biomarkers and consideration based on the existing research base as additio optimized reference standards for the discovery of desirable biomarkers.

Metabolomics to Identify Metabolic Changes in Retinal and Choroidal Neovascular Diseas
Metabolomics studies have reported several ocular diseases associated with neov cularization, including AMD, DR, and ROP, and we used these three common eye d eases as examples to develop a review of metabolic changes. We conducted a literatu search using the PubMed database. Articles were identified by searching for titles usi the following search terms: ("metabolomics" or "metabonomics" or "metabolome" "metabolic profiling") AND "neovascularization-related disease*". Additional artic were identified by searching the reference lists of the included studies. We includ metabolomics studies based on mass spectrometry (MS) or nuclear magnetic resonan (NMR), focusing on the generalization of population-based studies presenting with n vascularization pathological states and looking for differential metabolites, excludi drug evaluation reports, review articles, and abstracts without full text, resulting in inclusion of twenty-six studies, including five studies of ROP, nine studies of nAMD, a twelve studies of DR. A simplified flow diagram of study selection is shown in Figure  Important details about the study design and methods were extracted from selected a cles and summarized (Table 1).    All studies were observational, and all 25 studies were case-control studies, except for a longitudinal multi-time point observational study of serum in preterm infants established by Nilsson et al. [18]. Twelve studies (46.15%) were conducted in China [19][20][21][22][23][24][25][26][27][28][29][30], four (16%) in the United States [31][32][33][34], two (8%) in Japan [35,36], and the remaining studies were conducted in Spain, Korea Switzerland, and elsewhere. [18,[37][38][39][40][41][42][43] The studies were published between 2009 and 2022, with sample sizes ranging from 40 to 431, and nonequilibrium rates (number of cases in the majority category/number of cases in the minority category) between the case and control groups of the sample ranging from 1 to 3.66. Factors such as age, gender, race and lifestyle affect metabolic processes, with the ROP study recruiting newborns as study participants and both the nAMD study and the DR study recruiting middle-aged and older participants. The study by Mitchell et al. had the lowest proportion of males in the nAMD group (35%) [32] and the study by Tomita et al. study had the proliferative DR group with the highest proportion (77.1%) [36]. Smoking was present in some studies, especially involving elderly populations. Reasonable control of the influence of biological and experimental confounding factors is also important to ensure the reliability of the study results. Some studies reported measures such as collection of biological samples after eating, gender and age matching. Finally, the selection of study controls affects the interpretation of the study results, but due to ethical factors and sample availability, biological rock samples of controls are often obtained from other diseases or other procedures in our department, and not healthy controls, such as those in the study by Han et al., whose study control samples were obtained from cataract patients without AMD [22], Zhu, Sumarriva et al. used plasma from NPDR subjects [26,33], and Wang et al. obtained vitreous samples from MH patients and AH from cataract patients [21]. Regarding analytical platforms, 24 articles reported mass spectrometry-based metabolomics studies, including liquid chromatography-mass spectrometry (LC-MS; n = 22), and gas chromatography-mass spectrometry (GC-MS; n = 2), with different detection platforms covering different ranges of metabolic pouches, which also contributed to the differences in detection across studies.
The included studies generally utilized a combination of univariate and multivariate analyses to screen for between-group differential metabolites, with 11 of the total studies utilizing supervised multivariate analysis methods, such as partial least squares/orthogonal partial least squares discriminant analysis, in combination with univariate analysis methods, such as t/U test/univariate ANOVA/Wilcoxon rank sum test [19][20][21][22][25][26][27][28][29][30], with 4 of the 10 studies additionally adding fold change as a screening criteria [19,28,29] and one adding univariate area under curve value as an additional screening criterion [26]. It is worth noting that further screening criteria for difference variables often differed between studies, with criteria set at variable importance in projection (VIP) values greater than 1 or 2 and p values less than 0.05 or 0.01 for univariate analysis, or a Benjamini-Hochberg correction to control for false discovery rates. VIP values were used to screen metabolites by measuring the impact and explanatory power of metabolites on the increase in categorical discrimination of each group of samples in [44], while p-values were the probability of occurrence of more extreme results than those observed in the obtained samples in [45], and ploidy change methods, i.e., calculating the fold of difference in expression of a metabolite between two groups based on the relative quantitative or absolute quantitative results of the metabolite. Different analytical methods observe data from different perspectives, and different methods can be used in combination depending on the actual study purpose. In the included study, we counted significant metabolites and their trends ( Figure 2). There were 38 metabolites reported repeatedly in different biological samples, 12 of which were reported more than, or equal to, 3 times. Among them, Citrulline appeared most frequently (n = 6), followed by Arginine, Lysine, Proline and Pyruvic acid.

Biomarker Interpretation and Optimal Selection
Some of the reasons that prevent biomarkers identified by experimental studies to be used in clinical applications are the difficulty of explaining changes in complex disease mechanisms by individual biomarkers [46] and the difficulty of validating the results in other validation sets due to the strict limitations of sample screening [47]. Metabolism is a dynamic process, and metabolic data respond to changes in endogenous and exogenous metabolites at a particular point in time [48]. Within an organism, metabolites coordinate with each other to form a complete network of biological functions, and it is difficult to explain complex disease mechanisms with a single metabolite. Pathway analysis is a way to interpret and provide more biofunctional information about biomarkers by mapping metabolites that are changed by abnormal disease states onto an interactive network constructed with a priori knowledge, by mapping individual changed biomarkers into a coherent metabolic network [49]. In this review, biochemical interpretation of all altered metabolites was performed using MetaboAnalyst 5.0 (http://www.metaboanalyst.ca/, accessed on 1 August 2020), a software that allows analysis of the impact of particular compounds on biochemical pathways. Pathway analysis was derived from integrating differential metabolites based on the Kyoto Encyclopedia of Genes and Genomes and Human Metabolome Database [50]. We performed separate pathway analysis of the differential metabolites identified in our study in different biological sample sources, and the results of the study showed that 8, 4, 4 and 10 statistically significant metabolic pathways (p < 0.05) were enriched in plasma, serum, AH and vitreous humor, respectively. Table 2 shows the results of pathway analysis in each biological sample (FDR < 0.05) and matching differential metabolites in each pathway. Arginine biosynthesis, Aminoacyl-tRNA biosynthesis and Glyoxylate and dicarboxylate metabolism pathways were repeatedly enriched in different biological samples. When pathway analysis was performed separately, according to disease type, the following four metabolic pathways were repeatedly enriched (p < 0.05): Phenylalanine, tyrosine and tryptophan biosynthesis; Aminoacyl-tRNA biosynthesis; Linoleic acid metabolism; Lysine degradation. Among them, the Aminoacyl-tRNA biosynthesis pathway was meaningfully enriched in the analysis of different biological samples and disease species, respectively. It is worth noting that the existing databases, such as KEGG, include the results of existing studies, and more detailed information on metabolic pathways needs to be further improved. Therefore, the p-values and Q-values of the pathway analysis are only a reference for the study results, and those pathways that are not significant are also worth interpreting.  The focus of this study was to review metabolomics studies of the pathological phenomenon of retinal and choroidal neovascularization, and, therefore, included differential metabolites among subjects with different disease sources, biological samples, ages, and ethnicities, and screened for repeatedly identified differential metabolites. Biomarkers should have the ability to help identify diseases, and nine studies in this review evaluated the intergroup discriminatory performance of identified potential biomarkers [51], as shown in Table 3. Only the study by Mitchell et al. did an internal validation (test set) for identifying a panel of metabolite markers (159 differential features), obtaining an AUC value of 0.83 AUC value [32]. Ideal biomarkers need to be evaluated and compared in multiple ways, and repeatedly identified differential metabolites may become biological representations of disease characteristics. We validated the extent to which repeated identification of differential metabolites contributed to disease classification using the raw metabolic data mentioned in the four studies included in the review, including the study by Han et al. that included both anionic and cationic patterns and the study by Wang et al. that covered both vitreous and AH biological samples [21,22,37,42]. We randomly divided the case and control groups into training and validation sets containing 70% and 30% of these samples [52], respectively, and used repeated identification of differential metabolites as metabolic biomarker panels to construct multiple machine learning models for repeated identification of differential metabolites, including logistic regression, random forest, support vector machine (SVM) and XGBoost, and performed model evaluation. The results are shown in Table 4. The results show that repeated identification of differential metabolites had good classification performance in different studies and models. In addition, the study suggested that among the included metabolomics studies, random forests obtained better performance in small sample studies and n/p large differences. XGBoost showed better performance in unbalanced sample data and SVM and integration methods were more stable in studies with relatively small sample sizes and numbers of features. The results provide an additional reference for researchers to compare with the gold standard when identifying biomarkers in the future, in addition to duplicate identification markers in studies, and it covers both performance optimization of newly identified biomarkers and considerations based on the existing research base.

Biomarker Interpretation and Optimal Selection
Random Forest (RF) is a statistical learning theory, which uses a bootstrap resampling method to draw multiple samples from the original sample, model decision trees for each sample, and then combines the predictions of multiple decision trees to arrive at the final prediction by voting. Support vector machines (SVM) are based on a model that defines "the classifier with the largest interval on the feature space". The extreme gradient boosting (XGBoost) algorithm is based on the idea of boosting, which is an optimization improvement of the gradient boosting decision tree, and is a collection of many CART regression tree models to form a strong classifier. AUC = area under curve.

References Discriminant Models Discriminant Group; Precision
Zhou (2020) Altered metabolites 10 metabolites obtained AUC larger than 0.7 under positive ion mode; Under negative ion mode, 5 metabolites obtained AUC larger than 0.7.

Metabolomics for the Interpretation and Treatment of Retinal and Choroidal Neovascularization
A well-functioning vascular system is necessary to meet the high-energy metabolic demands of the eye and the function of retinal neurons, and pathological states of the vascular system (neovascularization) are accompanied by a range of metabolic changes and occur in a range of ocular diseases. Metabolomics is being used to unravel disease-associated metabolic disorders and also offers the opportunity to find new therapeutic targets, employing angiogenic phenotypes that may reprogram the metabolism of disease [53].

Amino Acid Metabolism and Neovascularization
A growing number of studies have shown that amino acid metabolism affects functions such as vascular tone, redox homeostasis, and immunity and inflammation, and amino acids with different characteristics have different functions [54,55]. In this study, many amino acid metabolism-related pathways were found to be separately enriched in different biological samples and diseases. Glutamine was identified as a differential metabolite in three studies and is the most abundant free amino acid in plasma [21,25,41], participating in the tricarboxylic acid cycle, providing precursors for gluconeogenesis and converting to glutamate and ammonia by the action of glutaminase. Glutamate is a retinal excitatory neurotransmitter, and high concentrations of glutamate overstimulate ionotropic glutamate receptors (N-methyl-D-aspartate receptors), leading to overload of intracellular calcium ions and ultimately cell death [56]. Glutamate is often involved in maintaining important functions of the retina as a downstream product of, or related to, other amino acid metabolic pathways. Arginine induces the gene expression of insulin-like growth factor 1 (IGF-1) [57], and IGF-1 is involved in vascular neogenesis by enhancing the vascular endothelial-derived growth factor (VEGF) signaling pathway [58], which was observed to be upregulated in four studies [24,27,41,43]. The metabolism of arginine is regulated by nitric oxide synthase (NOS) and arginase, and its metabolites could mediate a variety of angiogenic reactions. L-arginine is catalyzed by NOS to produce nitric oxide (NO), one kind of multifunctional messenger molecule. Numerous studies have shown NOS and NOS-produced NO involvement in promoting VEGF-induced vascular permeability and angiogenesis through the regulation of endothelial cell function [59][60][61]. Palmer et al. found that arginine deficiency in endothelial cell cultures was associated with eNOS dysfunction [62]. Arginase competes with NOS to metabolize L-arginine to produce urea and ornithine. Increased activity of arginase could decrease the availability of L-arginine for NOS, thus reducing bioavailability of NO production, increasing formation of reactive oxygen species, and ultimately resulting to endothelial cell dysfunction [63]. Proline is the downstream product of arginine metabolism and plays a key role in maintaining cell redox homeostasis. It was elucidated that Lproline could reduce the level of NO related marker P-asymmetric dimethylarginine and improve the redox state, thus increasing the bioavailability of NO [64]. Proline is converted to glutamate by further metabolism [65]. In addition, pathways related to branched chain amino acids (BCAAs) (valine, leucine and isoleucine biosynthesis) were found to be significantly enriched in tissue biological samples (AH and Vitreous). Of these, valine was identified in four studies [21,23,30,38], and, interestingly, valine levels showed an opposite direction of change in two DR studies, probably due to differences in study populations and analysis platforms, while leucine showed an increasing trend in both studies. Higher BCAAs level are associated with increased neurotoxic levels of glutamate in the retina [66]. BCAAs exert their function through activation of the mammalian rapamycin (mTOR) pathway, which regulates cell growth, proliferation, and survival [67]. MTOR plays a key role in upregulating the VEGF pathway, which leads to augmented expression of caspase-3 and promotes proliferation and migration of endothelial cells [68,69]. Complement factor H is a tyrosine sulfated protein and there is an association between its specific haplotype and an elevated risk of developing nAMD [70,71]. Complement factor H on the RPE and choroid regulates the formation of the local membrane attack complex, which is essential for the development of laser-induced CNV [72]. Inhibitors of amino acid-related metabolic pathways have important therapeutic value for diseases caused by pathological neovascularization. However, further studies are still needed.

Abnormalities in Neovascular Metabolism and Potential Therapeutic Opportunities
The active metabolism of the retina requires large amounts of energy and oxygen consumption. Vascular dysfunction may lead to oxidative stress and excessive production of reactive oxygen species, resulting in oxidative damage [73]. The process of neovascularization also requires significant oxygen consumption, which can further exacerbate retinal hypoxia. Hypoxia-inducible factor-1α (HIF-1α) is a specific mediator in adaptation to hypoxic environments and pathological responses [74]. HIF-1α deficiency in a nude mouse model significantly reduced the growth rate of malignant tumors and decreased neovascularization [75]. Under hypoxic conditions, HIF-1α enters the nucleus and binds to HIF-1β to form dimers, which finally bind to hypoxia response elements; the bound dimers promote angiogenesis by inducing upregulation of the expression of VEGF, erythropoietin, glucose transporter 1 and platelet-derived growth factor (PDGF) [76]. PDGF is a hypoxia-regulated gene product that, together with VEGF-A, promotes ocular neovascularization, and the combination of the two causes an additive effect [77]. PDGF inhibitor (Fovista; Ophthotech, New York, NY, USA) has been proposed in combination with anti-VEGF (ranibizumab) for the treatment of nAMD, and results from a phase 2b clinical trial showed significant improvement, compared to anti-vascular endothelial growth factor monotherapy [78]. However, to the disappointment of the clinical investigators, the next two phase 3 trials did not show the benefits of the combination therapy.
Oxidative stress mainly affects RPE cells [79]. Proline was identified in four studies and was reported to mediate metabolism between RPE cells and the retina, to act as an alternative energy source during stress or hypoxia, and to act as an antioxidant in maintaining redox homeostasis by controlling mitochondrial function [80]. Oxidative stress also disrupts the regulation of the complement system by RPE cells [81], which is involved in innate immune response and, once activated, forms membrane attack complexes [82]. Nozaki et al. showed that somatic induction of endochoroidal neovascularization increased the levels of C3a and C5a, and that C3a and C5a induced an increase in VEGF secretion by primary human RPE in vitro [83]. Knockout of C3 gene protected mice from choroidal neovascularization after laser treatment. Complement pathway-related therapeutic agents have been proposed, such as an inducer targeting C5 (ARC1905) used in clinical trials for the treatment of nAMD and geographic atrophy [84] and the therapeutic complement receptor 2-factor H used to inhibit angiogenesis and RPE injury in a mouse model of CNV [85].
Endothelial cells rely on glycolysis to obtain ATP, and upregulation of pyruvate (an intermediate product of glycolytic metabolism) was observed in both serum and vitreous metabolomics studies in [21,28,36] The end product of glycolysis, lactate, contributes largely to the angiogenic phenotype through inhibition of prolyl hydroxylase 2 and activation of HIF1α and nuclear factor-κB [86]. In this study, the Citrate cycle (TCA cycle) pathway was found to be enriched, and previous studies found that oral supplementation with TCA cycle metabolites helped maintain retinal function and protected photoreceptor cells and nerves on the inner retinal network in mice [87]. A recent study by Kanow et al. highlights the importance of maintaining a gradient of aerobic glycolytic activity between the retina and the RPE [88]. When the glycolytic activity of the RPE is abnormally high, less glucose is available for retinal consumption, which can also lead to photoreceptor death. This is evidenced by the upregulation of lactate observed in three studies. VEGF enhances the expression of the key glycolytic enzyme 6-phosphofructo-2-kinase/fructose-2,6-biphosphatase 3 (PFKFB3), increasing glycolysis to support the high ATP requirement required for angiogenesis [89,90]. Pharmacological blockade of PFKFB3 has been shown to effectively impair pathological angiogenesis and enhance the anti-angiogenic effects of anti-VEGF in CNV models [89]. Extracellular purine metabolites (ATP, ADP and adenosine) are important regulators of vascular homeostasis. It was shown that adenosine A2a receptor (ADORA2A) activation promotes transcriptional induction of glycolytic enzymes through the HIF-1α pathway, whereas knockdown of ADORA2A reduces hypoxia-induced glycolytic enzyme expression and glycolytic flux and inhibits endothelial cell proliferation and angiogenesis [91]. These findings provide directions for the development of therapeutic drugs targeting the glycolytic pathway for retinal neovascularization.
Pathway analysis revealed that Arachidonic acid metabolism, Biosynthesis of unsaturated fatty acids, Glycerolipid metabolism and Linoleic acid metabolism pathways were significantly enriched, suggesting that fatty acid metabolism may be involved in the mechanism of retinal and choroidal neovascularization. Retina is a tissue rich in polyunsaturated fatty acids (PUFA), which are metabolized mainly by cyclooxygenase, and lipoxygenase, and cytochrome P450 enzymes, 12-hydroxyeicosatetraenoic acid (12-HETE) and 20-hydroxyeicosatetraenoic acid , are among the metabolites [92][93][94]. It is thought that 12-HETE is an important regulator of neovascularization that regulates the expression of VEGF and PEDF. Arachidonic acid metabolism via cytochrome P450-4A/F produces 20-HETE, which has been shown to stimulate c-Src and epidermal growth factor receptor-mediated downstream signaling pathways, including MAPK and PI3K/Akt pathways, the eNOS uncoupling, and NOX/ROS system activation, thereby inducing the proliferation, migration, secretion, and angiogenesis of pro-angiogenic molecules (e.g., HIF-1α, VEGF) in endothelial progenitor cells. In vivo and in vitro studies have shown that eicosapentaenoic acid (EPA) inhibits the activation and expression of VEGF-specific tyrosine kinase receptors, while EPA supplementation significantly reduces the mRNA expression and protein levels of ICAM-1, MCP-1, VEGF and IL-6 [95][96][97], providing an opportunity to find new therapeutic targets. The number of studies is currently limited, and some studies of related diseases exist outside of the study to provide complementary findings [98]. Biomarkers and the search for disease-related metabolic pathways contribute to our understanding of disease mechanisms, but are not yet clinically translatable and require continuous research. In addition, metabolomics studies require strict quality control, from sample sources, inclusion criteria, analytical platforms, and data processing and analysis, which can all affect the results.

Conclusions
With the advances in histological technologies, more and more features of diseases at the molecular level of genes, proteins, metabolism, etc. have been obtained and have unveiled a new understanding of disease features. Retinal and choroidal neovascularization disease is a common feature of a group of diseases and is one of the serious pathologies that threaten visual function. Current metabolomics studies have identified metabolic changes in selected retinal and choroidal neovascular diseases, identified disease-related metabolic pathways, and provided new insights into understanding disease mechanisms. However, metabolomics study design is critical to the accuracy and reproducibility of findings, so this study summarized existing studies and advocates a comparative strategy, based on existing studies, as an additional guiding strategy to direct new biomarker discoveries. In addition, the study discusses existing metabolomics findings to promote understanding of disease mechanisms and potential therapeutic opportunities to provide additional therapeutic opportunities to the current treatment strategy based on anti-vascular endothelial factor therapy.