- 1College of Life and Environment Science, Huangshan University, Huangshan, Anhui, China
- 2Nanxiang Branch of Ruijin Hospital, Shanghai, China
- 3Department of Experimental Management, School of Integrative Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 4Department of Vascular Disease, Shanghai TCM-Integrated Hospital, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 5Traditional Chinese Recovery and Treatment Center, Zhejiang Rehabilitation Medical Center, Hangzhou, China
- 6College of Pharmacy, Hubei University of Chinese Medicine, Wuhan, China
- 7Department of Public Health, International College, Krirk University, Bangkok, Thailand
Introduction: This study investigates the role of hypoxia-related genes in the neuroprotective efficacy of Yang Xue oral liquid (YXKFY) in Alzheimer’s disease (AD) and Parkinson’s disease (PD).
Methods and results: Using differential expression and weighted gene co-expression network analysis (WGCNA), we identified 106 and 9 hypoxia-associated genes in AD and PD, respectively, that are implicated in the transcriptomic and proteomic profiles. An artificial intelligence-driven hypoxia signature (AIDHS), comprising 17 and 3 genes for AD and PD, was developed and validated across nine independent cohorts (n = 1713), integrating 10 machine learning algorithms and 113 algorithmic combinations. Significant associations were observed between AIDHS markers and immune cells in AD and PD, including naive CD4+ T cells, macrophages, and neutrophils. Interactions with miRNAs (hsa-miR-1, hsa-miR-124) and transcription factors (USF1) were also identified. Single-cell RNA sequencing (scRNA-seq) data highlighted distinct expression patterns of AIDHS genes in various cell types, such as high expression of TGM2 in endothelial cells, PDGFRB in endothelial and mesenchymal cells, and SYK in microglia. YXKFY treatment was shown to repair cellular damage and decrease reactive oxygen species (ROS) levels. Notably, genes with previously dysfunctional expression, including FKBPL, TGM2, PPIL1, BLVRB, and PDGFRB, exhibited significant recovery after YXKFY treatment, associated with riboflavin and lysicamine.
Conclusion: The above genes are suggested to be central to hypoxia and neuroinflammation responses in AD and PD, and are potential key mediators of YXKFY’s neuroprotective action.
Introduction
Alzheimer’s disease (AD) and Parkinson’s disease (PD) are the two most common neurodegenerative disorders that have a significant impact on the aging population worldwide (Nowell et al., 2023). AD is estimated to affect approximately 4% of individuals aged 65 and above, whereas PD is estimated to affect around 1% of the total population of individuals aged 60 and above (Macdonald et al., 2018; Aborode et al., 2022). AD is characterized by progressive memory loss, cognitive decline, and the accumulation of amyloid-beta (Aß) plaques and neurofibrillary tangles in the brain (Breijyeh and Karaman, 2020). PD, on the other hand, is characterized by motor symptoms such as tremors, rigidity, and bradykinesia, as well as the loss of dopaminergic neurons in the substantia nigra (Hayes, 2019). While the etiology of AD and PD remains complex and multifactorial, emerging evidence suggests that hypoxia is detrimental to the brain and plays a crucial role in the pathogeneses of these diseases (Burtscher et al., 2021). Hypoxia can result from various factors, including impaired blood flow, mitochondrial dysfunction, and oxidative stress, all of which have been implicated in the development and progression of AD and PD (Grabska-Kobyłecka et al., 2023; Su et al., 2023). Understanding the role of hypoxia-related features in the pathogenesis of AD and PD is crucial for the development of novel therapeutic strategies. Targeting hypoxia signaling pathways, such as HIF-1α, may offer potential therapeutic avenues for the treatment and prevention of these devastating neurodegenerative disorders (Baillieul et al., 2017).
Traditional Chinese medicines (TCMs) have been used in the treatment of neurological disease for thousands of years (Su et al., 2024). Yangxue oral liquid (YXKFY) is composed of four traditional Chinese medicines: Melanteritum, Fructus Crataegi, Fructus Hippophae, and Ziziphus jujuba. Since 7000 B.C., the fruit of Crataegus pinnatifida has been used to produce prehistoric fermented beverages (Moreira et al., 2023). Recent studies have shown that the fruit possesses neuroprotective properties that could be effective against AD (Lee et al., 2019; Moreira et al., 2023). Fructus Hippophae has been traditionally used in the treatment of brain conditions, with the aim of enhancing cognitive abilities such as learning and memory, as well as alleviating pathological damage in mice with AD through modulating oxidative stress and inflammatory processes (Dinkar Gore et al., 2023; Zhao et al., 2023). The fruit of Ziziphus jujuba has been proven to possess anti-inflammatory (Al-Reza et al., 2010) and neuroprotective (Kaeidi et al., 2015) abilities due to the presence of neuroprotective compounds such as triterpenoids, flavonoids, polysaccharides, saponins and alkaloids (Kim et al., 2021). Hence, YXKFY exhibits a promising potential for therapeutic intervention in the treatment of neurodegenerative diseases. Nevertheless, further research is needed to clarify the material basis and targets of YXKFY in addressing neurodegenerative diseases.
Recently, there has been a growing body of evidence supporting the widespread use of transcriptomics and proteomics in enhancing our understanding of pathophysiological mechanisms and facilitating the development of diagnostic tools for various diseases (Li et al., 2020; Mu et al., 2020; Yang et al., 2021; Yang et al., 2023). Moreover, the utilization of weighted gene coexpression network analysis (WGCNA) offers the ability to illustrate the interconnections between diverse genes through the creation of a co-expression network, which not only enables the identification of modules that are associated with specific phenotypes but also proves to be more efficient in the investigation of crucial pathways and genes involved in numerous human disorders (Langfelder and Horvath, 2008; Zhang et al., 2024). In addition, the limitations of expression-based multigene signatures in clinical settings can be attributed to their lack of uniqueness and appropriateness in selected modeling methods, and the absence of strict validation in large multicenter cohorts (Wang et al., 2022; Yang et al., 2024). Hence, our objective was to conduct a comprehensive analysis of the molecular mechanisms, as well as the diagnostic and therapeutic targets of YXKFY on AD/PD, using transcriptome, proteome, and hypoxia-related features.
In this study, we performed an integrated analysis to discovery a consensus artificial intelligence-driven hypoxia signature (AIDHS) from 110 kinds of algorithm combinations across transcriptome and proteome. Then, we employed ultra-high-performance liquid chromatography coupled with Quadrupole Exactive-Orbitrap high-resolution mass spectrometry (UHPLC-Q Exactive-Orbitrap HR-MS) to identify the chemicals present in YXKFY. Finally, the key AIDHS was applied in in vitro experimental validation to assess potential hypoxia-related mechanisms of YXKFY in AD/PD. Figure 1 showed the flowchart of the study design.
Materials and methods
Data acquisition from the GEO database
The mRNA and protein expression data for AD/PD were acquired from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/gds). A total of 36 AD, 48 PD and 47 normal serum samples were collected from GSE339087 (protein) for differential expression analysis. GSE5281 (87 AD cases and 74 controls) and GSE26927 (12 PD cases and 8 controls) datasets were selected for weighted gene co-expression network analysis. Then, we downloaded GSE150696 (9 AD, 12 PD cases and 9 controls), GSE132903 (97 AD cases and 98 controls), GSE63060 (145 AD cases and 104 controls), GSE63061 (140 AD cases and 134 controls), GSE20295 (40 PD cases and 53 controls), GSE99039 (205 PD cases and 233 controls) and GSE57475 (93 PD cases and 49 controls) for further analysis and validation.
Data preprocessing and differential expression analysis
Various intergroup comparisons were conducted using Sangerbox 3.0 (http://vip.sangerbox.com/home.html), which is a freely accessible online platform for data analysis and visualization. Initially, normalization of protein expression in the GSE339087 dataset (AD and PD) was carried out. Following this, volcano plots and heatmaps were generated for the differentially expressed proteins (DEPs). DEPs with an absolute fold change of ≥ 1.2 and a p-value < 0.05 were deemed to be statistically significant in terms of differential expression.
Weighted gene co-expression network analysis (WGCNA)
To identify modules of highly correlated genes and hub genes, we utilized the WGCNA (Langfelder and Horvath, 2008) approach to construct scale-free co-expression networks for two independent datasets: GSE5281 (87 AD cases and 74 controls) and GSE26927 (12 PD cases and 8 controls). In the GSE5281 dataset, the brain regions examined were the entorhinal cortex (EC), hippocampus (HIP), medial temporal gyrus (MTG), posterior cingulate (PC), superior frontal gyrus (SFG), and primary visual cortex (PVC). For the GSE26927 dataset, the substantia nigra (SN) was the focused brain region.
We initially employed Pearson correlation-based hierarchical clustering to group all the genes and samples. Subsequently, we determined the soft threshold power value to establish the co-expression network, ensuring that it conformed to a scale-free co-expression network. The adjacency matrices were then transformed into a topological overlap matrix to identify gene modules. Moreover, similar modules were merged and clustered together. Finally, we generated module-trait relationship diagrams and obtained the gene list associated with each module.
Intersection of hypoxia-related genes involved in AD and PD shared by transcriptome and proteome
The top four modules that have the closest association with AD/PD were identified by analyzing the module-trait relationship diagrams of the GSE5281 and GSE26927 datasets. A total of 8,371 hypoxia-related genes (HRG) were obtained from COREMINE (https://www.coremine.com/medical/#search). To identify the common HRGs between the DEGs and the aforementioned module genes, Venn diagram analysis was conducted. The overlapping genes were then subjected to functional enrichment analysis using the Metascape (https://metascape.org/gp/index.html) (Zhou Y. et al., 2019).
Signature obtained from artificial intelligence-driven integrative approaches
In order to ensure a highly accurate and stable performance for the consensus artificial intelligence-driven hypoxia signature (AIDHS), we integrated a total of 10 machine learning algorithms and explored 113 algorithm combinations. These integrative algorithms include Cox, CoxBoost, elastic network (Enet), generalised boosted regression modelling (GBM), Lasso, stepwise partial least squares regression for Cox (plsRcox), random survival forest (RSF), Ridge, supervised principal components (SuperPC), and survival support vector machine (survival-SVM).
The procedure for generating the signature was as follows: (a) We performed the 113 algorithm combinations on the HRGs, fitting prediction models based on the leave-one-out cross-validation (LOOCV) framework using the GSE39087 cohort; (b) We then validated all these models using nine additional datasets (GSE5281, GSE26927, GSE150696, GSE132903, GSE63060, GSE63061, GSE20295, GSE99039, and GSE57475); (d) For each model, we calculated the area under the curve (AUC) score across all validation datasets. In the case of PD, we considered the model with the highest average AUC as optimal. For AD, the optimal model was determined based on higher average AUC and the least number of genes involved.
Profile of the immune cell infiltration
In order to assess the role of the immune microenvironment in the development of AD and PD, we analyzed and compared the immune infiltration patterns in both AD/PD and normal groups through CIBERSORTx tool (https://cibersortx.stanford.edu/), which specifically designed for accurately quantifying the relative proportions of 22 immune cells subtypes within a complex mixture of gene expressions. Lastly, to assess the correlation between hub genes and differential infiltrated immune cells, a Pearson correlation test analysis was carried out. For further investigation, a correlation threshold of |r| > 0.2 was established.
Signaling information network and TF-miRNA regulatory networks analysis of the hub genes
SIGNOR 2.0 (https://signor.uniroma2.it/) (Licata et al., 2020), a public repository that stores manually-annotated causal relationships between proteins and other biologically relevant entities, was performed to explore the signal transduction relationship of the former hub genes. Then, TF-miRNA coregulatory network analysis was implemented on the hub genes through curated regulatory interaction information collected from the RegNetwork repository (Liu et al., 2015) with Networkanalyst online tools (Zhou G. et al., 2019).
Analysis the heterogeneity of striata in the AD and PD single-cell transcriptome
The basal ganglia are frequently affected in AD and PD. However, there is limited understanding of the fundamental molecular processes involved. The striatum is crucial for motor learning and various cognitive functions. Therefore, the single-cell dataset GSE161045 was obtained from the GEO database (https://www.ncbi.nlm.nih.gov/), which contains 4 AD, 4 PD and 4 control striata samples. The data underwent thorough preprocessing, including the removal of cells with fewer than 200 genes, more than 5,000 genes, or more than 25% mitochondrial genes. After this process, a total of 31,030 filtered cells remained for analysis. To mitigate potential batch effects due to sample identity, we utilized the “Harmony” R package (version 0.1.0) for batch correction (Korsunsky et al., 2019). Subsequently, we used the FindNeighbors and FindClusters functions (with a resolution of 0.1) for cell cluster identification, following the uniform manifold approximation and projection (UMAP) analysis. Prior to this, we referenced preexisting markers from published literature to identify specific cell types and their respective markers through SingleCellBase (Meng et al., 2023). Finally, we present bubble plots depicting the expression of 26 AIDHS genes in different disease and cell groups.
Serum preparation containing YXKFY
The YXKFY samples used in this study were acquired from Hubei Fenghuang Baiyunshan Pharmaceutical Co., Ltd. Blood samples containing YXKFY were obtained from the YXKFY rat group. The original concentration of the YXKFY was 0.4 g/mL. The rats were administered a dosage of 2 mg/mL of the YXKFY. Each rat received 5 mL of the YXKFY orally, once in the morning and once in the afternoon, for a total of two administrations per day. This oral administration was continued for 5 consecutive days. Following the last administration, the abdominal aorta was punctured to collect blood 1 h later. The collected blood samples were centrifuged at 3000 r/min to separate the serum, which was then incubated at 56°C for 30 min for inactivation. The serum was subsequently filtered and sterilized using a 0.22 micron filter, and stored at −20°C for future use. All of the animal procedures, such as housing and care, and experimental protocols were approved by the Ethics Committee of Shanghai University of Traditional Chinese Medicine (No. PZSHUTCM210702015).
Cell culture and treatment
Human neuroblastoma SH-SY5Y cells were cultured in a medium composed of DMEM supplemented with 10% fetal bovine serum (FBS, Gibco, Australia) at a controlled environment of 37°C with 5% CO2. The cells were evenly distributed with a density of 10,000 cells per well in 96-well plates and allowed to grow for a period of 12 h. To determine the optimal concentration of the hypoxic model in SH-SY5Y cells, the viability of SH-SY5Y cells was assessed using the CCK8 assay after exposure to different concentrations of CoCl2 (0 μM, 25 μM, 50 μM, 100 μM, 200 μM, 400 μM, 500 μM, 600 μM, 800 μM, and 1000 μM) for 24 h. In order to evaluate the potential neuroprotective properties of YXKFY on SH-SY5Y cells, CoCl2 group was exposed to 800 μM for 24 h; and YXKFY group was subjected to 800 μM CoCl2 and serum containing 10% YXKFY for 24 h.
ROS measurement
ROS generation in SH-SY5Y cells was examined using the fluorescent dye 2′,7′-dichlorodihydrofluorescein diacetate (DCFH-DA; Beyotime, Shanghai, China). SH-SY5Y cells were firstly washed three times with phosphate-buffered saline (PBS). Subsequently, the cells were incubated in PBS supplemented with 10 μM of DCFH-DA at a temperature of 37°C for a duration of 20 min. Following the incubation period, cells were once again washed three times using PBS. To visualize the cells, the measurements and analysis were conducted using a FACS Calibur flow cytometer (Beckman, San Jose, United States).
Validation of YXKFY-associated AIDHS genes by RT-qPCR
The extraction of total RNA from cells was carried out using TRIzol (Invitrogen Corporation, CA, United States). Subsequently, cDNA was synthesized using the cDNA Synthesis SuperMix. In order to detect the mRNA expression levels of AIDHS, RT-qPCR was conducted using GAPDH as an internal reference. The amplification primers used in the experiment were obtained from Shanghai Sangon Biological Engineering Technology and are detailed in Table 1. To ensure consistency, the mRNA expression level of each target gene was normalized to that of GAPDH within the same sample. The relative expression of each target gene was determined using the 2−ΔΔCT method.
Analysis of YXKFY and drug serum via high performance liquid chromatography (HPLC)
To perform the analysis, an Acquity UPLC® BEH C18 column (100 mm × 2.1 mm, 1.7 µm) was utilized in conjunction with a Dionex Ultimate 3,000 high-pressure liquid system. The mobile phase consisted of a mixture of 0.1% formic acid in water and methanol. The elution procedure involved a stepwise gradient: from 0 to 4 min, 4% methanol; from 4 to 10 min, 4%–12% methanol; from 10 to 30 min, 12%–70% methanol; at 30 min, 70% methanol; from 35 to 38 min, 70%–95% methanol; at 38 min, 95% methanol; and from 42 to 45 min, 4% methanol. The flow velocity was maintained at 0.3 mL/min. The column temperature was set at 40°C, and the sample size used for each analysis was 5 μL. For the analysis and detection of the compounds, a Q Exactive quadrupole-electrostatic field orbitrap high-resolution mass spectrometer was employed. The ion source used was an electrospray ion source (H-ESI). The scanning mode utilized for this analysis was full MS/SIM mode (covered m/z 80–1,200), allowing for the detection of both positive and negative ions.
Molecular docking validation
The Protein Data Bank (PDB, https://www.rcsb.org/) was accessed to retrieve the crystal structures of the primary targets. The three-dimensional configurations of the active compounds were obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/). Utilizing AutoDock Vina (http://vina.scripps.edu/), molecular docking was carried out to compute the affinity between the active compounds and primary targets. The corresponding PDB identifiers were 5oog, 2 × 7 k and 3s3j for BLVRB, PPIL1, and TGM2, respectively. The AlphaFold algorithm was used to predict the structures of FKBPL and PDGFRB, which are proteins without a known 3D structure in the PDB database. Then, water molecules were eliminated from the proteins, and polar hydrogen atoms were added. Subsequently, active binding sites were constructed based on the ligand’s position in the PDB complex.
Statistical analysis
The data presented in this study were expressed as the mean ± standard deviation (SD), using the social science statistical software package for data analysis. An independent sample t-test was employed to compare the two groups of data, while a one-way analysis of variance was used to compare multiple sets of data. The column graphs were created using GraphPad Prism 6 software. A significance level of p < 0.05 was considered statistically significant.
Results
Identification of DEPs in AD and PD
Human serum protein expression dataset of GSE39087 (AD and PD) was applied for our differential protein expression analysis. After data preprocessing and differential expression analysis, we identified 1988 DEPs between AD and healthy controls (|Fold change| ≥ 1.2 and p-value < 0.05), including 723 upregulated (DNAJC8, IGLC1, LGALS1, ICAM4, and TRA2A) and 1,265 downregulated proteins (EDC3, HRH1, RP9, IL9, and HOXA1) (Figure 2A). The top 30 DEPs of AD (Figure 2C) were visualized in heatmap. In addition, we also identified 344 DEPs between PD and healthy controls (|Fold change| ≥ 1.2 and p-value < 0.05), including 90 upregulated (VEGF, RASL11B, CRB3, DBF4B, and PCBD2) and 254 downregulated proteins (ERO1LB, C1orf162, WDR22, MTMR2, and RABEP2) (Figure 2B). And the top 30 DEPs of PD were visualized in heatmap (Figure 2D).
Figure 2. Identification of differentially expressed proteins (DEPs) in AD and PD. Volcano plots of DEPs in AD (A) and PD (B). Heatmaps of the top 30 DEPs in AD (C) and PD (D). Red, upregulated DEPs; blue, downregulated DEPs.
WGCNA of the whole transcriptome expression matrix
We performed WGCNA analysis to find similar gene expression patterns tended to exert similar biological functions in GSE5281 (AD) and GSE26927 (PD). Figures 3A, B displayed the hierarchical clustering dendrograms of the samples. Subsequently, we examined the soft threshold powers of the network topology and determined β values of 5 and 12 to be the optimal soft-thresholding parameters for GSE5281 (AD) and GSE26927 (PD) datasets (Figures 3E, F). From the GSE5281 dataset, we identified 9 modules (Figure 3C), while from the GSE26927 dataset, we identified 15 modules (Figure 3D). To investigate the associations between gene modules and AD or PD, we constructed module-trait diagrams. In GSE5281, the four modules with the highest correlation to AD were MEyellow (r = −0.57, p = 5e−15), MEblack (r = 0.53, p = 3e−13), MEblue (r = 0.49, p = 6e−11), and MEred (r = 0.46, p = 8e−10). In GSE26927, the four modules with the highest correlation to PD were MEpink (r = 0.71, p = 5e−04), MEmagenta (r = 0.59, p = 0.006), MEpurple (r = 0.55, p = 0.01), and MEyellow (r = 0.55, p = 0.01).
Figure 3. Identification of modules correlated with AD and PD in transcriptome datasets. Hierarchical clustering dendrograms were constructed for the samples in the AD (A) and PD (B) datasets. Additionally, module-trait relationship diagrams were generated for the AD (C) and PD (D) datasets. In these diagrams, each row represents a color module, while each column represents a clinical trait. The correlation and p-value for each module-trait combination are displayed in the respective cells. Furthermore, an analysis was performed to assess the scale independence and mean connectivity for the optimal soft threshold powers in the AD (E) and PD (F) datasets.
Hypoxia-related genes involved in AD and PD shared by transcriptome and proteome
To discover functional genes concurrently involved in both AD and PD across the transcriptome and proteome, we employed two methods to identify key genes. Firstly, we used WGCNA to identify modules that are related to AD and PD. Secondly, we conducted differential expression analysis to identify the most dysregulated genes. In total, we identified 200 genes in AD and 17 genes in PD that were concurrently involved in the transcriptome and proteome, respectively (Figures 4A, B). Enrichment analysis of concurrent Genes Involved in AD showed that the terms were mainly enriched in cytokine signaling in immune system, PI3K-Akt signaling pathway, mitochondrial matrix, inflammatory response, and positive regulation of superoxide anion generation (Figure 4C), whereas PD was related to histone modifying activity, regulation of muscle system process, focal adhesion, regulation of metal ion transport and negative regulation of striated muscle cell differentiation (Figure 4D).
Figure 4. Hypoxia-related genes involved in AD and PD shared by transcriptome and proteome. Venn diagram of hypoxia-related genes shared by transcriptome and proteome in AD (A) and PD (B). Pathway enrichment analyses of the common hypoxia-related genes in AD (C) and PD (D).
In view of the important role of hypoxia-related genes in AD and PD, we extracted 106 and 9 hypoxia-related genes involved in AD and PD for further analysis by intersecting the DEPs and WGCNA related mRNAs with 8,371 hypoxia-related genes.
Integrated development of AD and PD consensus signature
A total of 106 and 9 hypoxia-related genes involved in AD and PD were further incorporated into our integration program to develop an AIDHS. In the protein training cohort, we applied 110 algorithm combinations via ten-fold cross-validation to construct prediction models and calculated the average AUC of each algorithm in the remaining testing cohorts. In AD, the combination of glmBoost and LDA with higher average AUC (0.764) and proper genes (N = 17) was selected as the final model, which contained BLVRB, C10orf54, CDC37, CHCHD6, FKBPL, FN1, GLUL, HOXC4, PFN1, PPIL1, PRDX4, PSMB1, PSMG1, SLC16A14, SLC17A7, SMYD3, and SYK. In PD, the combination of glmBoost and RF with the highest average AUC (0.722) was selected as the final model (N = 3), including FMO5, HDAC4, and TGM2 (Figure 5).
Figure 5. Generation of the artificial intelligence-derived hypoxia-related gene signature in AD and PD. The C-indices of 110 machine-learning algorithms in eleven cohorts of AD (A) and seven in PD (B).
Altered of immune cell infiltration in AD and PD
We observed 8 (GSE60630) and 4 (GSE60631) significant differences in the infiltration of cells between AD and healthy samples. In GSE60630, T cells CD4 naïve, NK cells resting, Macrophages M0, Mast cells activated and Neutrophils cell types were upregulated, whereas T cells CD4 memory resting, T cells gamma delta and Macrophages M2 cell types downregulated (Figure 6A). In GSE60631, NK cells resting cell type was upregulated, whereas B cells naive, T cells CD4 memory resting and Macrophages M2 cell types downregulated (Figure 6B).
Figure 6. The results of immune cell infiltration analysis. A box plot of AD infiltrated immune cells in GSE60630 (A) and GSE60631 (B). A box plot of PD infiltrated immune cells in GSE57475 (C) and GSE99039 (D). *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
We observed 4 (GSE57475) and 7 (GSE99039) significant differences in the infiltration of cells between PD and healthy samples. In GSE57475, T cells CD4 memory activated and NK cells resting cell types were upregulated, whereas B cells memory and T cells CD4 memory resting cell types downregulated (Figure 6C). In GSE99039, 4 cell types were upregulated, including T cells CD4 naive, Monocytes, Macrophages M0, and Neutrophils; the three downregulated cell types were: B cells memory, Macrophages M1, and Macrophages M2 (Figure 6D).
Correlation analysis between AIDHS markers and infiltrating cells
|r| > 0.2 was set as a correlation threshold for further analysis. In AD, T cells CD4 naïve cell was correlated with BLVRB, C10orf54, CDC37, CHCHD6, GLUL, HOXC4, and PPIL1 in GSE63060; and C10orf54, CHCHD6, GLUL, PFN1, PSMB1, and SYK in GSE63061. T cells CD4 memory resting cell was correlated with C10orf54, CDC37, GLUL, HOXC4, PPIL1, PRDX4, PSMB1, PSMG1, SLC16A14, and SYK in GSE63060; and BLVRB, C10orf54, CDC37, GLUL, HOXC4, PFN1, PPIL1, PRDX4, PSMB1, PSMG1, SLC17A7, and SYK in GSE63061. T cells gamma delta cell was correlated with C10orf54, PFN1, PPIL1, PRDX4, PSMB1, PSMG1, and SYK in GSE63060; and C10orf54, CHCHD6, GLUL, PFN1, PPIL1, PRDX4, PSMB1, PSMG1, and SYK in GSE63061. NK cells resting cell was correlated with BLVRB, CDC37, HOXC4, and PSMB1 in GSE63060; and C10orf54, CHCHD6, GLUL, HOXC4, PFN1, PRDX4, PSMB1, PSMG1, and SYK in GSE63061. Macrophages M0 cell was correlated with C10orf54, CHCHD6, GLUL, PPIL1, PSMG1, and SYK in GSE63060; and C10orf54, GLUL, PPIL1, PRDX4, PSMB1, PSMG1, and SYK in GSE63061.Macrophages M2 cell was correlated with C10orf54, PRDX4, PSMG1, and SYK in GSE63060; and GLUL in GSE63061. Mast cells activated cell was correlated with CDC37 and SMYD3 in GSE63060; and BLVRB, PPIL1, and SMYD3 in GSE63061. Neutrophils were correlated with C10orf54, CDC37, CHCHD6, GLUL, HOXC4, PPIL1, PRDX4, PSMB1, PSMG1, SMYD3, and SYK in GSE63060, and in GSE63061 (Figures 7A, B).
Figure 7. The correlation between AIDHS genes and infiltrated immune cells in AD ((A): GSE60630, (B) GSE60630) and PD ((C): GSE57475, (D) GSE99039). Nodes that appear redder indicate a stronger positive correlation, while nodes that appear bluer indicate a stronger negative correlation.
In PD, B cells memory cell was related to FMO5 and TGM2 in GSE57475. T cells CD4 naïve cell was related to HDAC4 in GSE99039. Monocytes was related to FMO5 in GSE57475. Macrophages M0 cell was related to TGM2 in GSE57475. Neutrophils was related to HDAC4 in GSE99039 (Figures 7C, D).
TF–miRNA coregulatory network
The NetworkAnalyst online platform was used to generate a TF-miRNA co-regulatory network. The analysis of this network provided insights on the interaction between TFs, miRNAs, and the selected hub genes, which could explain the underlying mechanism of how hub gene expression is regulated. Figure 8A shows the TF–miRNA co-regulatory network. The AD TF–miRNA co-regulatory network is consisted of 269 nodes and 343 edges. A total of 149 miRNAs and 61 TF genes interacted with the validated AIDHS genes, including MYC, MAX, TFAP2A, USF1, NFIC, hsa-miR-632, hsa-miR-1, hsa-miR-124, hsa-miR-127-5p and hsa-miR-140-5p. Functional enrichment analysis showed these nodes were related to Osteoclast differentiation, Th17 cell differentiation, TNF signaling pathway, IL-17 signaling pathway, and Longevity regulating pathway (Figure 8B). HOXC4 owned the largest number of neighbors, followed by FN1, PFN1, CDC37, SLC17A7, SLC16A14, CHCHD6, SMYD3, PPIL1, BLVRB, and PRDX4.
Figure 8. TF–miRNA coregulatory network and pathway enrichment analyses in AD (A,B) and PD (C,D). The round nodes are AIDHS genes, the diamond nodes are TF genes, and the square nodes are miRNA genes.
Figure 8C shows the TF–miRNA co-regulatory network. The PD TF–miRNA co-regulatory network is consisted of 108 nodes and 109 edges. A total of 46 miRNAs and 28 TF genes interacted with the validated AIDHS genes, including USF1, TFAP2A, TFAP2C, SP1, POU2F1, hsa-miR-19a, hsa-miR-19b, hsa-miR-1, hsa-miR-106b and hsa-miR-124. Functional enrichment analysis showed these nodes were related to Apelin signaling pathway, cGMP-PKG signaling pathway, Cellular senescence, Cell cycle and Osteoclast differentiation (Figure 8D). HDAC4 owned the largest number of neighbors.
Construct the signaling information network based on SIGNOR
We further built the signaling information network of AIDHS through SIGNOR2.0. The AD signaling network included six interaction mechanisms, such as ubiquitination, transcriptional regulation, polyubiquitination, phosphorylation, dephosphorylation and binding. SYK showed the highest number of interaction, followed by FN1, PFN1, CDC37, and PSMB1. Moreover, FN1 interacted with more complex than other hub genes, including SNAIL/RELA/PARP1, FN1/SDC4, Av/b6 integrin, A8/b1 integrin, A5/b1 integrin complex. In addition, the analysis also identified the effect of actin cytoskeleton reorganization as potential phenotype (Figure 9A).
The PD signaling network included five interaction mechanisms, such as transcriptional regulation, phosphorylation, dephosphorylation, deacetylation and binding. HDAC4 showed the highest number of interaction and was inhibited by AMPK complex. TGM2 can directly phosphorylate H3C1 and Histone H3 (Figure 9B).
Analysis the heterogeneity of striata in the AD and PD single-cell transcriptome
The quality control summary of our single-cell RNA sequencing (scRNA-seq) data is presented in Figure 10A. After preprocessing the data, we used the harmony algorithm to merge the samples and effectively eliminate any potential batch effects. Figure 10B displays the representation of the merged dataset, which includes AD, PD and control samples after the implementation of the harmony algorithm. By following the standard procedures of Seurat, a total of 12 clusters were successfully identified and visualized using UMAP, as shown in Figure 10C. Subsequently, we observed eight distinct cell clusters (Figure 10D), including astrocytes (AQP4, GFAP), CD8 T cells (CCL5, NKG7, GZMK, GZMA, GZMM, CCL4, CD69, LGALS1, LGALS3, PDCD1, PRDM1, and CXCR6), endothelial cells (FLT1, CLDN5), mesenchymal cells (CD44, CHI3L1, and HIF1A), microglia cells (CSF1R, CD74, C3), neural cells (GAD1, GAD2, SOX2, and SLC17A7), oligodendrocytes (MBP, PLP1), and Opcs (PDGFRA, VCAN, CSPG4). Figure 11A showed the distribution of 8 cell types in 12 samples. And the top 5 marker genes were identified within the cell-type populations (Figure 11B). Figure 11C showed that most of the 26 AIDHS genes were differentially expressed among the three groups. These genes were specifically expressed in different cells, such as TGM2 being highly expressed in endothelial cells, PDGFRB being highly expressed in endothelial and mesenchymal cells, and SYK being highly expressed in microglia cells (Figure 11D).
Figure 10. Preprocessing of scRNA-seq data was conducted to prepare for downstream analysis. (A) Following quality control steps, the counts, percentages, and features related to mitochondrial and ribosome genes were assessed for each sample under analysis. (B) Visualization of sample distribution was achieved using UMAP plots after integration of the datasets using the harmony algorithm. (C) A clustering algorithm with a resolution of 0.1 was applied, resulting in the identification of 12 clusters. (D) UMAP plot was used to show the distribution of the eight different cell types.
Figure 11. ScRNA-seq data was employed to annotate and scrutinize the distribution of cell proportions among the three groups. (A) The overall cell fraction comparison among the three groups. (B) A heatmap was produced to display the levels of expression of marker genes in the identified cells. Bubble plots depicting the expression of 26 AIDHS genes in different disease (C) and cell groups (D).
YXKFY restored the damaged SH-SY5Y cells caused by CoCl2 and reduce the ROS level
The viability of SH-SY5Y cells was assessed using the CCK8 assay after exposure to different concentrations of CoCl2 for 24 h. As depicted in Figure 12A, the viability of cells gradually declined as the concentration of CoCl2 increased from 400 to 1000 μM. Notably, when the CoCl2 concentration reached 800 μM, SH-SY5Y cell viability dropped to 57%. Subsequently, the concentration of 800 μM was selected for subsequent experiments as it induced significant damaging effects. In order to examine the potential of YXKFY in preventing CoCl2-induced cell death, the CCK8 assay was performed on SH-SY5Y cells treated with YXKFY-containing serum in the presence of CoCl2. As depicted in Figure 12B, the addition of YXKFY-containing serum led to a notable improvement in cell viability, with an increase to 67%. Moreover, YXKFY-containing serum can reduce the ROS level in the hypoxia model (Figure 12C).
Figure 12. YXKFY restored the damaged SH-SY5Y cells caused by CoCl2 and reduce the ROS level. (A) The viability of SH-SY5Y cells was assessed using the CCK8 assay after exposure to different concentrations of CoCl2. (B) YXKFY-containing serum restored the damaged SH-SY5Y cells caused by CoCl2. (C) YXKFY-containing serum reduced the ROS level in the hypoxia model. ***: p < 0.001.
Validation of hub genes
In order to validate the findings of the bioinformatics analysis, we identified the 17 AIDHS genes of AD and 9 AIDHS genes of PD. The mRNA expression levels of these genes were assessed using RT-qPCR. The hypoxia group exhibited decreased expression of BLVRB, C10orf54, CDC37, FKBPL, FN1, GLUL, HOXC4, PPIL1, PRDX4, PSMB1, SLC17A7, SMYD3, SYK, FMO5, HDAC4, TGM2, PDGFRB, MYBPH, and JUP and increased expression of FXYD5 and DNM2 compared to the control group. Additionally, the expression of FKBPL, TGM2, PPIL1, BLVRB, and PDGFRB showed substantial recovery after YXKFY treatment (Figure 13).
Figure 13. Validation of the mRNA expression levels of AIDHS genes. *: p < 0.05, **: p < 0.01, ***: p < 0.001.
Identification of the components in YXKFY and YXKFY-containing serum
Based on the information provided by high-resolution mass spectrometry, such as accurate molecular masses for quasi-molecular ions and loading ions, the precise relative molecular mass of the primary mass spectra is inferred. The preliminary estimation of the various components obtained under high-resolution mass spectrometry is then conducted using the Peakview 1.2 software for molecular formula fitting. The HPLC analysis identified a total of 107 chemical components in the YXKFY, and a further 9 prototype components were identified in the YXKFY-containing serum (Table 2).
Molecular docking validation
The present research investigated binding affinity between the above five genes and the 8 prototype ingredients in the YXKFY-containing serum. The molecular docking results revealed different levels of binding between natural products and hub genes, with the majority of natural products showing good binding capabilities to the hub genes. Particularly, riboflavin and lysicamine exhibited higher binding affinity than other compounds (Table 3).
Discussion
To this day, researchers continue to face the challenging task of finding an effective drug for treating AD and PD using modern medicine. However, there is a growing body of evidence suggesting that traditional Chinese medicine may offer promise in the prevention and treatment of AD and PD (Pei et al., 2020; Chen et al., 2022). Hypoxia, characterized by oxygen deprivation, is implicated in neuronal stress, initiating a cascade that culminates in neuronal dysfunction and demise, thereby exacerbating the pathogenesis and progression of AD and PD (March-Diaz et al., 2021; Guo et al., 2022). Our study, therefore, presents an integrative strategy, merging bidirectional transcriptome and proteome analyses, to delineate the hypoxia-related molecular signatures and pathways that underpin the neuroprotective efficacy of YXKFY in both AD and PD.
In our study, we first profiled hypoxia-related mRNA and protein expression in AD and PD. Our analysis revealed that the dysregulated mRNAs and proteins in both conditions predominantly converge on mitochondrial and inflammatory pathways in AD, and pathways related to the muscle system in PD. It is a recognized phenomenon that hypoxia triggers the activation of multiple signaling cascades, notably cytokine signaling within the immune system (Sturrock et al., 2018) and the PI3K-Akt pathway (Pan et al., 2021), which are integral to inflammatory responses and the modulation of mitochondrial function (Wang et al., 2021). The mitochondrial matrix is acutely vulnerable to hypoxia, which can augment the production of superoxide anions, thereby exacerbating oxidative stress and neuroinflammation. (Spooner et al., 2021). Furthermore, hypoxia has been implicated in the upregulation of pro-inflammatory cytokines and chemokines, perpetuating the chronic inflammatory milieu characteristic of AD and PD (Merelli et al., 2021). Additionally, hypoxia impacts muscle function, with disruptions in muscle system processes linked to the motor symptoms of PD (Miterko et al., 2021). Consequently, elucidating the molecular underpinnings of hypoxia’s impact on these pathways is essential for devising innovative therapeutic strategies for both AD and PD.
In this study, we conducted a comprehensive analysis using nine independent cohorts to identify stable AIDHS for the development of AD and PD. Recognizing the potential influence of researcher preferences and biases in algorithm selection, we employed a rigorous approach by integrating ten distinct machine learning algorithms to discern the optimal signature from a comprehensive set of 113 algorithmic combinations. This strategy culminated in the identification of AIDHS in AD (17 genes) and PD (3 genes). Our in vitro validation corroborated the machine learning predictions, with the majority of genes showing significant differential expression following hypoxia treatment. Notably, YXKFY was found to ameliorate the dysfunction of five key genes (FKBPL, TGM2, PPIL1, BLVRB, and PDGFRB). While some of the biomarkers associated with these genes have been previously implicated in AD and PD, others represent novel findings in the context of neurodegenerative diseases. PPIL1 plays a crucial role as an enzyme-substrate pair within the spliceosome, exerting its function in the facilitation of RNA splicing and ensuring the survival of neurons, and its mutations may cause neurodegenerative disease (Chai et al., 2021). The silencing of TGM2, implicated in the pathogenesis of AD and PD, inhibited mitochondrial calcium influx, accumulation of mtROS, phosphorylation of Tau, and ultimately protected neuronal cells from Aβ-induced cell death (Lee et al., 2021). The depletion of PDGFRB within the precuneus region in individuals with AD has been linked to fibrinogen leakage, a decline in oxygenation, and the accumulation of fibrillar Aβ (Miners et al., 2018). Therefore, PPIL1, TGM2 and PDGFRB may serve as therapeutic targets of YXKFY in AD and PD. However, FKBPL and BLVRB have never been investigated in AD and PD, suggesting a new therapeutic strategy through hypoxia-related mechanism in AD and PD.
In our study, HPLC were also utilized to pinpoint specific ingredients in YXKFY, which led to the identification of 107 chemical components. The 8 prototype ingredients in the YXKFY-containing serum were found to be quinic acid, lysicamine, riboflavin, cinnamic acid, p-hydroxybenzoic acid, quercetin-7-O-α-L-rhamnoside, tumerone and isoamyl senecioate. Several of these components, including quinic acid, riboflavin, and cinnamic acid, have demonstrated significant anti-inflammatory effects and upregulating antioxidant enzyme activities in neurodegenerative disease. For instance, in in vitro studies, quinic acid has demonstrated neuroprotective and neurotrophic effects on Aβ-induced toxicity, as well as enhancing the activity of neurite outgrowth in PC 12 cells (Hur et al., 2001). It is suggested that quinic acid exerts its neuroprotective effects through the PKA signaling pathway and has successfully restored catalase levels (Rebai et al., 2017). Patients with dementia and AD have shown low levels of riboflavin, and supplementing with riboflavin has been found beneficial in treating cognitive impairment, as it shields cells from oxidative stress by boosting antioxidant enzyme activities and the glutathione redox cycle, and diminishing pro-inflammatory responses in the brain (Zhang et al., 2023). Cinnamic acid, a naturally occurring antioxidant, triggers PPARα activation, promoting lysosomal biogenesis and reducing amyloid plaque pathology in an AD mouse model (Chandra et al., 2019), and protecting dopaminergic neurons in a PD mouse model (Prorok et al., 2019). Therefore, the combination of our discovered results that YXKFY can reduce the ROS levels in hypoxic SH-SY5Y cell injury, suggesting that these components may be effective in YXKFY’s prevention and treatment of AD and PD through antioxidant damage. Moreover, Riboflavin has the best binding affinity with PPIL1, followed by BLVRB, TGM2, FKBPL, and PDGFRB. Additionally, lysicamine and tumerone also exhibit good binding affinity with these proteins. Collectively, the potential therapeutic effects of these chemical monomers from YXKFY on AD and PD through these genes warrant further exploration in our future studies.
Conclusion
In conclusion, our study leveraged a synergistic approach combining transcriptome and proteome analyses with machine learning to innovatively identify 17 and 3 hypoxia-associated biomarkers for AD and PD, respectively. Furthermore, we elucidated the neuroprotective mechanisms of YXKFY, highlighting its antioxidant properties. This was achieved through the identification of eight bioactive compounds, engagement with five molecular targets, and modulation of several pathways pivotal in neuroprotection.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
All animal procedures, such as housing and care and experimental protocols were approved by the Ethics Committee of Shanghai University of Traditional Chinese Medicine (No. PZSHUTCM210702015).
Author contributions
X-YC: Writing–original draft. M-RC: Writing–review and editing. C-CT: Writing–review and editing. C-QX: Writing–review and editing. Y-LZ: Writing–review and editing. YG: Writing–review and editing. X-XC: Writing–review and editing. JC: Writing–review and editing, Writing–original draft.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (82104873, 82274408), the Shanghai Sailing Program (20YF1445900), Xinglin Young Talent training Program of Shanghai University of Traditional Chinese Medicine (RCPY0040), Preparatory discipline leader training program of Shanghai TCM-Integrated hospital (RCPY0071), school-based cooperation project (hxkt2023016).
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.
References
Aborode, A. T., Pustake, M., Awuah, W. A., Alwerdani, M., Shah, P., Yarlagadda, R., et al. (2022). Targeting oxidative stress mechanisms to treat Alzheimer's and Parkinson's disease: a critical review. Oxid. Med. Cell Longev. 2022, 7934442. doi:10.1155/2022/7934442
Al-Reza, S. M., Yoon, J. I., Kim, H. J., Kim, J. S., and Kang, S. C. (2010). Anti-inflammatory activity of seed essential oil from Zizyphus jujuba. Food Chem. Toxicol. 48 (2), 639–643. doi:10.1016/j.fct.2009.11.045
Baillieul, S., Chacaroun, S., Doutreleau, S., Detante, O., Pépin, J. L., and Verges, S. (2017). Hypoxic conditioning and the central nervous system: a new therapeutic opportunity for brain and spinal cord injuries? Exp. Biol. Med. (Maywood) 242 (11), 1198–1206. doi:10.1177/1535370217712691
Breijyeh, Z., and Karaman, R. (2020). Comprehensive review on Alzheimer's disease: causes and treatment. Molecules 25 (24), 5789. doi:10.3390/molecules25245789
Burtscher, J., Mallet, R. T., Burtscher, M., and Millet, G. P. (2021). Hypoxia and brain aging: neurodegeneration or neuroprotection? Ageing Res. Rev. 68, 101343. doi:10.1016/j.arr.2021.101343
Chai, G., Webb, A., Lee, C., Antaki, D., Lee, S., Breuss, M. W., et al. (2021). Mutations in spliceosomal genes PPIL1 and PRP17 cause neurodegenerative pontocerebellar hypoplasia with microcephaly. Neuron 109 (2), 241–256 e9. doi:10.1016/j.neuron.2020.10.035
Chandra, S., Roy, A., Jana, M., and Pahan, K. (2019). Cinnamic acid activates PPARα to stimulate Lysosomal biogenesis and lower Amyloid plaque pathology in an Alzheimer's disease mouse model. Neurobiol. Dis. 124, 379–395. doi:10.1016/j.nbd.2018.12.007
Chen, P., Zhang, J., Wang, C., Chai, Y. H., Wu, A. G., Huang, N. Y., et al. (2022). The pathogenesis and treatment mechanism of Parkinson's disease from the perspective of traditional Chinese medicine. Phytomedicine 100, 154044. doi:10.1016/j.phymed.2022.154044
Dinkar Gore, D., Ahmad, F., Tikoo, K., Kumar Bansal, A., Kumar, D., and Pal Singh, I. (2023). Comparative quantitative analysis of fruit oil from Hippophae rhamnoides (seabuckthorn) by qNMR, FTIR and GC-MS. Chin. Herb. Med. 15 (4), 607–613. doi:10.1016/j.chmed.2023.05.005
Grabska-Kobyłecka, I., Szpakowski, P., Król, A., Książek-Winiarek, D., Kobyłecki, A., Głąbiński, A., et al. (2023). Polyphenols and their impact on the prevention of neurodegenerative diseases and development. Nutrients 15 (15), 3454. doi:10.3390/nu15153454
Guo, M., Ji, X., and Liu, J. (2022). Hypoxia and alpha-synuclein: inextricable link underlying the pathologic progression of Parkinson's disease. Front. Aging Neurosci. 14, 919343. doi:10.3389/fnagi.2022.919343
Hayes, M. T. (2019). Parkinson's disease and parkinsonism. Am. J. Med. 132 (7), 802–807. doi:10.1016/j.amjmed.2019.03.001
Hur, J. Y., Soh, Y., Kim, B. H., Suk, K., Sohn, N. W., Kim, H. C., et al. (2001). Neuroprotective and neurotrophic effects of quinic acids from Aster scaber in PC12 cells. Biol. Pharm. Bull. 24 (8), 921–924. doi:10.1248/bpb.24.921
Kaeidi, A., Taati, M., Hajializadeh, Z., Jahandari, F., and Rashidipour, M. (2015). Aqueous extract of Zizyphus jujuba fruit attenuates glucose induced neurotoxicity in an in vitro model of diabetic neuropathy. Iran. J. Basic Med. Sci. 18 (3), 301–306.
Kim, M. J., Jung, J. E., Lee, S., Cho, E. J., and Kim, H. Y. (2021). Effects of the fermented Zizyphus jujuba in the amyloid β(25-35)-induced Alzheimer's disease mouse model. Nutr. Res. Pract. 15 (2), 173–186. doi:10.4162/nrp.2021.15.2.173
Korsunsky, I., Millard, N., Fan, J., Slowikowski, K., Zhang, F., Wei, K., et al. (2019). Fast, sensitive and accurate integration of single-cell data with Harmony. Nat. Methods 16 (12), 1289–1296. doi:10.1038/s41592-019-0619-0
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Lee, H. J., Jung, Y. H., Choi, G. E., Kim, J. S., Chae, C. W., Lim, J. R., et al. (2021). Urolithin A suppresses high glucose-induced neuronal amyloidogenesis by modulating TGM2-dependent ER-mitochondria contacts and calcium homeostasis. Cell Death Differ. 28 (1), 184–202. doi:10.1038/s41418-020-0593-1
Lee, J., Cho, E., Kwon, H., Jeon, J., Jung, C. J., Moon, M., et al. (2019). The fruit of Crataegus pinnatifida ameliorates memory deficits in β-amyloid protein-induced Alzheimer's disease mouse model. J. Ethnopharmacol. 243, 112107. doi:10.1016/j.jep.2019.112107
Li, L., Dong, L., Xiao, Z., He, W., Zhao, J., Pan, H., et al. (2020). Integrated analysis of the proteome and transcriptome in a MCAO mouse model revealed the molecular landscape during stroke progression. J. Adv. Res. 24, 13–27. doi:10.1016/j.jare.2020.01.005
Licata, L., Lo Surdo, P., Iannuccelli, M., Palma, A., Micarelli, E., Perfetto, L., et al. (2020). SIGNOR 2.0, the SIGnaling network open resource 2.0: 2019 update. Nucleic Acids Res. 48 (D1), D504–D510. doi:10.1093/nar/gkz949
Liu, Z. P., Wu, C., Miao, H., and Wu, H. (2015). RegNetwork: an integrated database of transcriptional and post-transcriptional regulatory networks in human and mouse. Database (Oxford) 2015, bav095. doi:10.1093/database/bav095
Macdonald, R., Barnes, K., Hastings, C., and Mortiboys, H. (2018). Mitochondrial abnormalities in Parkinson's disease and Alzheimer's disease: can mitochondria be targeted therapeutically? Biochem. Soc. Trans. 46 (4), 891–909. doi:10.1042/BST20170501
March-Diaz, R., Lara-Ureña, N., Romero-Molina, C., Heras-Garvin, A., Ortega-de San Luis, C., Alvarez-Vergara, M. I., et al. (2021). Hypoxia compromises the mitochondrial metabolism of Alzheimer's disease microglia via HIF1. Nat. Aging 1 (4), 385–399. doi:10.1038/s43587-021-00054-2
Meng, F. L., Huang, X. L., Qin, W. Y., Liu, K. B., Wang, Y., Li, M., et al. (2023). singleCellBase: a high-quality manually curated database of cell markers for single cell annotation across multiple species. Biomark. Res. 11 (1), 83. doi:10.1186/s40364-023-00523-3
Merelli, A., Repetto, M., Lazarowski, A., and Auzmendi, J. (2021). Hypoxia, oxidative stress, and inflammation: three faces of neurodegenerative diseases. J. Alzheimers Dis. 82 (1), S109–S126. doi:10.3233/JAD-201074
Miners, J. S., Schulz, I., and Love, S. (2018). Differing associations between Aβ accumulation, hypoperfusion, blood-brain barrier dysfunction and loss of PDGFRB pericyte marker in the precuneus and parietal white matter in Alzheimer's disease. J. Cereb. Blood Flow. Metab. 38 (1), 103–115. doi:10.1177/0271678X17690761
Miterko, L. N., Lin, T., Zhou, J., van der Heijden, M. E., Beckinghausen, J., White, J. J., et al. (2021). Neuromodulation of the cerebellum rescues movement in a mouse model of ataxia. Nat. Commun. 12 (1), 1295. doi:10.1038/s41467-021-21417-8
Moreira, J., Machado, M., Dias-Teixeira, M., Ferraz, R., Delerue-Matos, C., and Grosso, C. (2023). The neuroprotective effect of traditional Chinese medicinal plants-A critical review. Acta Pharm. Sin. B 13 (8), 3208–3237. doi:10.1016/j.apsb.2023.06.009
Mu, H. Q., Liang, Z. Q., Xie, Q. P., Han, W., Yang, S., Wang, S. B., et al. (2020). Identification of potential crucial genes associated with the pathogenesis and prognosis of prostate cancer. Biomark. Med. 14 (5), 353–369. doi:10.2217/bmm-2019-0318
Nowell, J., Blunt, E., and Edison, P. (2023). Incretin and insulin signaling as novel therapeutic targets for Alzheimer's and Parkinson's disease. Mol. Psychiatry 28 (1), 217–229. doi:10.1038/s41380-022-01792-4
Pan, F., Xu, X., Zhan, Z., and Xu, Q. (2021). 6-Gingerol protects cardiomyocytes against hypoxia-induced injury by regulating the KCNQ1OT1/miR-340-5p/PI3K/AKT pathway. Panminerva Med. 63 (4), 482–490. doi:10.23736/S0031-0808.20.03956-7
Pei, H., Ma, L., Cao, Y., Wang, F., Li, Z., Liu, N., et al. (2020). Traditional Chinese medicine for Alzheimer's disease and other cognitive impairment: a review. Am. J. Chin. Med. 48 (3), 487–511. doi:10.1142/S0192415X20500251
Prorok, T., Jana, M., Patel, D., and Pahan, K. (2019). Cinnamic acid protects the nigrostriatum in a mouse model of Parkinson's disease via peroxisome proliferator-activated receptorα. Neurochem. Res. 44 (4), 751–762. doi:10.1007/s11064-018-02705-0
Rebai, O., Belkhir, M., Sanchez-Gomez, M. V., Matute, C., Fattouch, S., and Amri, M. (2017). Differential molecular targets for neuroprotective effect of chlorogenic acid and its related compounds against glutamate induced excitotoxicity and oxidative stress in rat cortical neurons. Neurochem. Res. 42 (12), 3559–3572. doi:10.1007/s11064-017-2403-9
Spooner, R. K., Taylor, B. K., Moshfegh, C. M., Ahmad, I. M., Dyball, K. N., Emanuel, K., et al. (2021). Neuroinflammatory profiles regulated by the redox environment predicted cognitive dysfunction in people living with HIV: a cross-sectional study. EBioMedicine 70, 103487. doi:10.1016/j.ebiom.2021.103487
Sturrock, A., Woller, D., Freeman, A., Sanders, K., and Paine, R. (2018). Consequences of hypoxia for the pulmonary alveolar epithelial cell innate immune response. J. Immunol. 201 (11), 3411–3420. doi:10.4049/jimmunol.1701387
Su, W., Yang, Y., Zhao, X., Cheng, J., Li, Y., Wu, S., et al. (2024). Potential efficacy and mechanism of eight mild-natured and bitter-flavored TCMs based on gut microbiota: a review. Chin. Herb. Med. 16 (1), 42–55. doi:10.1016/j.chmed.2023.08.001
Su, Z., Zhang, G., Li, X., and Zhang, H. (2023). Inverse correlation between Alzheimer's disease and cancer from the perspective of hypoxia. Neurobiol. Aging 131, 59–73. doi:10.1016/j.neurobiolaging.2023.07.002
Wang, L., Liu, Z., Liang, R., Wang, W., Zhu, R., Li, J., et al. (2022). Comprehensive machine-learning survival framework develops a consensus model in large-scale multicenter cohorts for pancreatic cancer. Elife 11, e80150. doi:10.7554/eLife.80150
Wang, Q., Wang, P., Qin, Z., Yang, X., Pan, B., Nie, F., et al. (2021). Altered glucose metabolism and cell function in keloid fibroblasts under hypoxia. Redox Biol. 38, 101815. doi:10.1016/j.redox.2020.101815
Yang, H., Cheng, J., Zhuang, H., Xu, H., Wang, Y., Zhang, T., et al. (2024). Pharmacogenomic profiling of intra-tumor heterogeneity using a large organoid biobank of liver cancer. Cancer Cell 42 (4), 535–551 e8. doi:10.1016/j.ccell.2024.03.004
Yang, M. D., Zhou, W. J., Chen, X. L., Chen, J., Ji, Q., Li, Q., et al. (2021). Therapeutic effect and mechanism of bushen-jianpi-jiedu decoction combined with chemotherapeutic drugs on postoperative colorectal cancer. Front. Pharmacol. 12, 524663. doi:10.3389/fphar.2021.524663
Yang, W., Yin, H., Wang, Y., Wang, Y., Li, X., Wang, C., et al. (2023). New insights into effects of Kaixin Powder on depression via lipid metabolism related adiponectin signaling pathway. Chin. Herb. Med. 15 (2), 240–250. doi:10.1016/j.chmed.2022.06.012
Zhang, M., Chen, H., Zhang, W., Liu, Y., Ding, L., Gong, J., et al. (2023). Biomimetic remodeling of microglial riboflavin metabolism ameliorates cognitive impairment by modulating neuroinflammation. Adv. Sci. (Weinh) 10 (12), e2300180. doi:10.1002/advs.202300180
Zhang, S., Ta, N., Zhang, S., Li, S., Zhu, X., Kong, L., et al. (2024). Unraveling pancreatic ductal adenocarcinoma immune prognostic signature through a naive B cell gene set. Cancer Lett. 594, 216981. doi:10.1016/j.canlet.2024.216981
Zhao, H., Liu, J., Wang, Y., Shao, M., Wang, L., Tang, W., et al. (2023). Polysaccharides from sea buckthorn (Hippophae rhamnoides L.) berries ameliorate cognitive dysfunction in AD mice induced by a combination of d-gal and AlCl(3) by suppressing oxidative stress and inflammation reaction. J. Sci. Food Agric. 103 (12), 6005–6016. doi:10.1002/jsfa.12673
Zhou, G., Soufan, O., Ewald, J., Hancock, R. E. W., Basu, N., and Xia, J. (2019b). NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 47 (W1), W234–W241. doi:10.1093/nar/gkz240
Keywords: Alzheimer’s disease, Parkinson’s disease, transcriptome, proteome, signature
Citation: Chen X-Y, Cheng M-R, Tang C-C, Xu C-Q, Zhong Y-L, Gao Y, Cheng X-X and Chen J (2024) Integrative transcriptome-proteome approach reveals key hypoxia-related features involved in the neuroprotective effects of Yang Xue oral liquid on Alzheimer’s and Parkinson’s disease. Front. Pharmacol. 15:1411273. doi: 10.3389/fphar.2024.1411273
Received: 02 April 2024; Accepted: 13 June 2024;
Published: 09 July 2024.
Edited by:
Tianhao Liu, Affiliated Hospital of Jiangnan University, ChinaReviewed by:
Abdur Rauf, University of Swabi, PakistanKevin Spelman, Health, Education and Research, United States
Copyright © 2024 Chen, Cheng, Tang, Xu, Zhong, Gao, Cheng and Chen. 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: Yuan Gao, gaoacu@sina.com; Xue-Xiang Cheng, chengxuexiang@bjanting.com; Jian Chen, alexandercj@126.com
†These authors have contributed equally to this work