- 1Medical College, Institute of Translational Medicine, Yangzhou University, Yangzhou, China
- 2Jiangsu Key Laboratory of Integrated Traditional Chinese and Western Medicine for Prevention and Treatment of Senile Diseases, Yangzhou University, Yangzhou, China
- 3Department of Neurology, Taizhou Second People's Hospital, Taizhou, China
- 4Department of Neurology, Yangzhou Hospital of Chinese Medicine, Yangzhou, China
Objective: We aimed to investigate the effect and mechanisms of action of two drug pairs [Huangqi-Chuanxiong and Sanleng-Ezhu Herb (HCSE)] on the treatment of ischemic stroke.
Materials and methods: We mined the current literature related to ischemic stroke and formulated a new formulation of Chinese herbs. Then, we identified the main candidate target genes of the new formulation by network pharmacology. Next, we performed enrichment analysis of the target genes to identify the potential mechanism of action of the new formulation in the treatment of ischemic stroke. Next, we experimentally validated the mechanism of action of the new formulation against ischemic stroke. Infarct volume and neurological deficits were evaluated by 2,3,5-triphenyltetrazolium (TTC) staining and Longa’s score, respectively. The predicted pathways of signal-related proteins were detected by western blotting.
Results: We mined the current literature and identified a new formulation of Chinese herbs for the treatment of ischemic stroke. The formulation included Huangqi, Chuanxiong, Sanleng and Ezhu. Next, we used network pharmacological analysis to identify 23 active compounds and 327 target genes for the new formulation. The key target genes were MAPK3, MAPK1, HSP90AA1, STAT3, PIK3R1, PIK3CA and AKT1. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed significant enrichment of the PI3K/AKT and MAPK/ERK signaling pathways. By performing experiments, we found that the new formulation reduced the infarct volume of middle cerebral artery occlusion (MCAO) induced mice and activated the PI3K/AKT and MAPK/ERK signaling pathways. These findings confirmed that the new formulation has a significant protective effect against ischemic stroke injury by activating the PI3K/AKT and MAPK/ERK signaling pathways.
Conclusion: We identified a new treatment formulation for ischemic stroke by data mining and network pharmacological target prediction. The beneficial effects of the new formulation act by regulating multiple target genes and pathways. The mechanism of action of the new formulation may be related to the AKT and ERK signaling pathways. Our findings provide a theoretical basis for the effects of the new formulation on ischemic stroke injury.
Introduction
Stroke has become the second leading cause of global death and is the leading cause of death and disability among adults in China. Ischemic stroke accounts for 69.6–77.8% of stroke cases and is associated with mortality and recurrence rates of 7 and 16%, respectively (Wang et al., 2017; Yiping Chen et al., 2020). Ischemic stroke causes irreversible damage to the brain tissue by influencing the distribution of affected blood vessels, thus leading to a variety of neurological symptoms and signs. Ischemic stroke adversely affects the health of patients and places significant socioeconomic burden on their families (Mukundan and Seidenwurm, 2018).
According to the theory of traditional Chinese medicine (TCM) (Liu et al., 2007), the main syndromes and causes of ischemic stroke are Qi deficiency and blood stasis syndrome (Li et al., 2014; Zhai et al., 2020). Therefore, the main principle of treatment for ischemic stroke is the invigoration of Qi and the promotion of blood circulation (Yu Wang et al., 2020). Based on this principle, practitioners of ancient Chinese medicine have suggested many treatments for stroke (Xu et al., 2019). The most commonly used and most effective treatment is Buyang Huanwu Decoction, invented by Wang Qingren, a famous doctor in the Qing Dynasty. Several studies have shown that the modified Buyang Huanwu Decoction has significant positive effects on the neurological deficits of patients with ischemic stroke (Xi Chen et al., 2020). When using Buyang Huanwu Decoction, the effect of invigoration of Qi and the promotion of blood circulation is mainly induced by two specific drugs: Huangqi and Chuanxiong. Some drugs mainly treated blood stasis, such as Taoren (Xi et al., 2013). However, research on the molecular effects associated with the invigoration of Qi and the promotion of blood circulation is scarce. Therefore, in the present study, we investigated the molecular mechanisms underlying the traditional methods used to prevent and treat ischemic stroke.
Data mining and network pharmacology can systematically evaluate the interactions between diseases and drugs and identify the specific mechanisms of action of drugs on their gene targets. Many previous studies have used data mining and network pharmacology to investigate the treatment of diseases by TCM (Li and Zhang, 2013; Sun et al., 2020). These strategies can provide new ideas for research and allow for the more accurate application of TCM.
In the present study, data mining and network pharmacology were used to evaluate the related mechanisms and effects of Huangqi-Chuanxiong and Sanleng-Ezhu Herb pairs (HCSE) on the treatment of patients suffering from ischemic stroke with Qi deficiency and blood stasis syndrome. Our findings provide scientific evidence to support the application of the new formulation for the treatment of ischemic stroke.
Materials and Methods
Data Mining
Literature Review
First, we searched the China National Knowledge Infrastructure, Wanfang Data, and PubMed using the keywords “stroke” and “ischemic stroke with Qi deficiency and blood stasis syndrome” to identify articles published the 1st September 2011 and the 1st September 2021. Then, we used the keywords “circulating blood and removing stasis” and “broking blood stasis” to search the clinical literature on the use of TCM to treat ischemic stroke with Qi deficit and blood stasis syndrome (Liu et al., 2012).
Data Screening
The prescriptions included in the retrieved literature were analyzed, and matrix distribution was performed.
Association Analysis
The drugs were analyzed using the a priori module included in SPSS Modeler software (version 18.0; IBM Corp., Armonk, NY, USA). The minimum number of conditional supports, indicating the number of drug combinations in the prescriptions, was set to 10. The minimum rule confidence, indicating the probabilities of occurrence for the first and second terms of the rule, was set to 80%. Drug pairs were then identified based on the association rules, and the new formulations were combined using TCM theory analysis.
Network Pharmacology
Prediction of Target Genes for the Drug Ingredients
The active compounds of the drugs were screened based on an oral bioavailability (OB) ≥ 30%, a drug-likeness (DL) ≥ 0.18, and a blood-brain barrier (BBB) ≥ −0.3 using the TCMSP database (https://tcmsp-e.com/) (Ru et al., 2014). Target proteins of the active ingredients were identified using SwissTargetPrediction (Swiss Institute of Bioinformatics, Basel, Switzerland; http://www.swisstargetprediction.ch/) (Daina et al., 2019).
Disease Target Identification
The keyword “ischemic stroke” was used to search disease targets in the GeneCards database (https://www.genecards.org/) (Safran et al., 2021). According to the principles of TCM, the main symptoms of ischemic stroke with Qi deficiency and blood stasis syndrome are a pale white complexion, a shortness of breath, palpitations, spontaneous sweating, loose stools, limb swelling, angular salivation, a purple tongue, a thin and whitish tongue coating, and a deep and faint pulse. We used these keywords as search terms in the GeneCards database. These symptoms were then selected as the targets for Qi deficiency and blood stasis syndrome. Overlapping targets were considered as potential targets of ischemic stroke with Qi deficiency and blood stasis syndrome.
Predicting Overlapping Drug and Disease Targets
Drug and disease targets were imported into Jvenn software (http://jvenn.toulouse.inra.fr/app/example.html) (Bardou et al., 2014) to construct a Venn diagram; overlapping targets were then considered to be potential therapeutic targets for the new treatment for ischemic stroke.
Construction of a Protein-Protein Interaction Network
A PPI network was constructed for the potential therapeutic target proteins of the new treatment formulation using the STRING database (https://string-db.org). The results were then visualized using Cytoscape software (version 3.8.2; https://cytoscape.org/). The degree of freedom was indicated by node size and color. The topological properties of the target genes were analyzed using the CytoNCA tool (Cytoscape software plugin) (Tang et al., 2015; Bei Yin et al., 2020). The PPI network was used to screen key targets of the new treatment based on a degree centrality (DC) > two-fold of the median and betweenness centrality (BC), closeness centrality (CC), and an eigenvector centrality (EC) > one-fold of the median. Higher quantitative values were correlated with a greater importance of the node.
Gene Ontology Enrichment and Kyoto Encyclopedia of Genes and Genomes Pathway Analyses
Metascape (https://metascape.org/) was used to perform GO and KEGG pathway enrichment analyses of 327 candidate target genes, according to the molecular function (MF), biological process (BP), and cellular component (CCT) categories (Zhou et al., 2019). p values were used to evaluate the proteins each GO annotation, thus reflecting the significance of the biological function. In addition, the FDR error control method (FDR <0.05) was used to test and correct the p value. The threshold value of p < 0.05 was finally used to screen the biological processes with significant differences. GO and KEGG enrichment analyses were visualized by the Bioinformatics platform (http://www.bioinformatics.com.cn/).
Experimental Verification
Experimental Animals
Healthy male Institute of Cancer Research mice (ICR; body weight: 32–35 g) were provided by the Comparative Medicine Center of Yangzhou University (Yangzhou, China) and used after 1 week of acclimation. The mice were housed in a controlled condition with a 12-h light/dark cycle at 23°C and 60% humidity with free access to food and water. All experimental investigation procedures for animals were permitted by the Yangzhou University Institutional Animal Care and Use Committee (Grant No. YIACUC-14-0015).
Middle Cerebral Artery Occlusion
The middle cerebral artery occlusion (MCAO) model was created using the following procedure. The mice were continuously anesthetized by the inhalation of a mixture of 30% oxygen, 70% nitrogen and 3–4% isoflurane (RWD Life Science, Guangdong Province, China). The mice were then fixed in the supine position. Next, we sterilized and dissected the skin of the middle of the neck. The left common carotid artery (CCA), the internal carotid artery (ICA), and the external carotid artery (ECA), were then separated with a blunt instrument. The ECA was ligated, the CCA and ICA were clamped with a bulldog clip, and a “V” incision was made in the ECA with ophthalmic scissors. A 0.23 mm monofilament nylon suture (Beijing Biotechnology Co., Ltd.) was then inserted through the “V” incision into the ICA until slight resistance was achieved. After 45 min of arterial occlusion, the monofilament nylon suture was removed, and the blood was reperfused for 24 h. In the Sham group, only the CCA, ICA and ECA were separated; no monofilament nylon suture was inserted. The body temperature of the mice was maintained at 37.0–37.5°C during surgery.
Experimental Groups and Drug Treatments
The mice were randomly divided into four groups as follows (n = 21 in each group): 1) a control group and 2) a pre-HCSE (13.65 g/kg) group. Each of these two groups was then further divided into a Sham group and a MCAO group. All animals were killed at 24 h after cerebral ischemia reperfusion. The method used to prepare HCSE involved mixing 225 g of Huangqi (100 g), Chuanxiong (50 g), Sanleng (45 g), and Ezhu (30 g) with eight times the volume of distilled water, followed by boiling for 2 h. Then, the filtrate was collected, and the filtrate was extracted again with three times the volume of distilled water. The filtrate was mixed twice and concentrated to 225 ml by a rotary evaporator. Therefore, the final crude drug concentration was 1 g/ml. All herbs were purchased from the Yangzhou Hospital of Chinese Medicine, Jiangsu Province, China. According to the dose used for ischemic stroke patients in the clinic (90 ml per day), we determined the dose for our experimental animals by extrapolation from the human dose in accordance with a previous study (Reagan-Shaw et al., 2008). Finally, a dose of 13.65 g/kg was obtained for administration. We chose 6.83 g/kg, 13.65 g/kg, and 27.3 g/kg as the HCSE doses in our preliminary study. We found that the neuroprotective effects (as determined by TTC staining) of the 13.65 g/kg and 27.3 g/kg doses were better than those with the 6.83 g/kg dose. Therefore, we selected a dose of 13.65 g/kg for use in the present study.
HCSE was administered intra-gastrically the same day, 24 h, and 48 h before ischemic surgery, and was administered twice a day. The Sham group and the ischemia group received the same amount of 0.9% saline intra-gastrically for the same durations. Following the last administration, the mice underwent ischemia/reperfusion or Sham operation.
Neurological Deficit Assessment
The degree of neurological deficit was assessed 24 h after reperfusion using the Longa score, as follows: 1) no neurological deficit: 0 points; 2) inability to fully extend the front paw on the paralyzed side: 1 point; 3) circling to the paralyzed side during walking: 2 points; 4) leaning towards the paralyzed side when walking: 3 points; and 5) inability to walk spontaneously, with loss of consciousness: 4 points. A score of more than 1 indicates that the MCAO model had been successfully established.
TTC Staining and the Quantification of Infarct Volume
Next, we performed TTC staining experiments with reference to relevant published literature (Liu et al., 2020; Zhang et al., 2021). After 24 h of ischemia/reperfusion, mice were anesthetized by isoflurane inhalation and killed by cervical dislocation. Brain tissue was then removed and cut into 2 mm thick coronal sections. The sections were then stained using 2% 2,3,5-triphenyltetrazolium chloride (Sigma-Aldrich, St. Louis, MO) in the dark at 37.0°C for 30 min. After staining, the brain sections were fixed in 4% paraformaldehyde buffer. The live portion of the brain section was red, and the infarcted portion was pale white. The infarct volume and whole volume of each brain slice were measured using Image Pro Plus 6.0 software (Media Cybernetics, Bethesda, MD, USA). The infarct volume of each section was multiplied by the layer thickness (2 mm) to calculate the total infarct volume. The ratio of total infarct volume/whole brain volume × 100% was used as the total infarct ratio.
Western Blot Analysis
Western blotting was carried out as described in our previous article (Liu et al., 2021). First, hippocampal and cortical tissues were lysed with a Whole Protein Extraction Kit (Solarbio, Beijing, China). The protein concentration of each sample was then measured using a BCA protein assay kit (Vazyme, Nanjing, China). Equal amounts of protein (30 μg) were then separated by 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to nitrocellulose membranes (Millipore, Bedford, USA). The membranes were cultured in 5% BSA in TBS containing 0.1% Tween 20 for 90 min to prevent non-specific protein binding sites from binding to antibodies. The membranes were cut according to molecular weight and then incubated overnight at 4°C with rabbit anti-AKT (1:1,000, Cell Signaling Technology), rabbit anti-p-AKT (1:2000, Cell Signaling Technology), rabbit anti-ERK 1/2 (1:1,000, Cell Signaling Technology), rabbit anti-p-ERK 1/2 (1:1,000, Cell Signaling Technology). Subsequently, the membrane was incubated with a corresponding secondary antibody at room temperature for 2 h. The protein membranes were then to Super Signal West Pico chemiluminescence substrates (CWBIO, Beijing, China) and observed by Fluor Chem M (Protein Simple, USA). Western blot results were analyzed by Image Pro Plus 6.0.
Data Analysis
Data are expressed as a mean ± standard deviation (‾X ± SD) and analyzed by GraphPad Prism 8.0. SPSS software (version 21.0; IBM, Armonk, NY, USA). The student’s two-tailed t-test was used for comparisons between two groups and one-way analysis of variance followed by Dunnett’s test was used for the comparison of three or more groups. p < 0.05 was considered statistically significant.
Results
Data Mining
Literature Analysis
We identified 177 articles related to TCM treatment of ischemic stroke with Qi deficiency and blood stasis syndrome, including 76 prescriptions based on 119 TCM drugs. In total, 77 articles were related to breaking blood and removing blood stasis for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome; these articles included 52 prescriptions based on 94 TCM drugs.
Types of TCM Drugs
The 10 most common TCM drugs used for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome were Huangqi (n = 69), Chuanxiong (n = 64), Honghua (n = 51), Danggui (n = 50), Dilong (n = 50), Chishao (n = 46), Taoren (n = 40), Danshen (n = 27), Shuizhi (n = 20), Jixueteng (n = 19), and Sanqi (n = 14). The 10 most common TCM drugs used for breaking blood and removing stasis were Chuanxiong (n = 46), Honghua (n = 39), Sanleng (n = 37), Ezhu (n = 37), Danshen (n = 35), Danggui (n = 35), Chishao (n = 29), Taoren (n = 29), Huangqi (n = 26), and Dilong (n = 19).
Association Analyses
The five most commonly used drug combinations for TCM treatment of ischemic stroke with Qi deficiency and blood stasis syndrome were Chuanxiong/Huangqi, Huangqi/Chuanxiong, Chishao/Honghua, Dilong/Honghua, and Danggui/Honghua (Table 1). The five most commonly used drug combinations for breaking blood and removing blood stasis for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome were Ezhu/Sanleng, Sanleng/Ezhu, Honghua/Danggui, Chuanxiong/Danggui, and Chuanxiong/Taoren (Table 2). Based on these results, we decided to combine Huangqi, Chuanxiong, Sanleng, and Ezhu to form a new formulation (HCSE).
Network Pharmacology
Chemical Component Analysis
After excluding components with undiscovered targets, 26 candidate components of Huangqi, Chuanxiong, Sanleng, and Ezhu were identified in the TCMSP database based on specific thresholds (OB ≥ 30%, DL ≥ 0.18, and BBB ≥ −0.3). In total, 12, 5, 5, and 3 compounds were identified from Huangqi, Chuanxiong, Sanleng and Ezhu, respectively (Table 3). Of these, Huangqi, Sanleng and Ezhu all contain Hederagenin as a key ingredient; Huangqi and Sanleng both contain Formononetin as a key ingredient.
Integration of Disease and Drug Target Genes
The GeneCards database was used to integrate 3,854 genes related to ischemic stroke and 5,103 genes related to the syndrome of Qi deficiency and blood stasis syndrome; these analyses identified a total of 2,696 genes. Then, the target genes of the drug components and disease were compared in a Venn diagram using Jvenn software. The Venn diagram showed that 327 genes were potential targets of the new formulation for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome (Figure 1).
FIGURE 1. Venn diagram of the overlapping drug and disease target genes. The red circle represents the targets for ischemic stroke with Qi deficiency and blood stasis syndrome. The blue circle represents the targets for the new formulation. The 327 overlapping genes are potential therapeutic targets for the new formulation against ischemic stroke with Qi deficiency and blood stasis syndrome.
PPI Network
We imported the 327 potential therapeutic target genes into the STRING database. Then, we imported target genes with a confidence >0.9 into Cytoscape 3.8.2 software to construct a PPI network. In the PPI network, the nodes represent the target genes, and the node size and color represent the degrees of freedom. As shown in Figure 2, the network consisted of 244 nodes and 1,081 edges. The central properties of the nodes were estimated by topological analysis. Sixteen key gene targets of the new formulation for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome were screened using the following criteria: DC > 11, BC > 728.62274, CC > 0.37587029, and EC > 0.09713785. In Figure 3, the size of the nodes is proportional to the DC. Notably, MAPK3, MAPK1, HSP90AA1, STAT3, PIK3R1, PIK3CA, and AKT1 were the main target proteins involved in the pathogenesis of ischemic stroke with Qi deficiency and blood stasis syndrome. Among these targets, MAPK3 (degree = 50) was shown to be the most important protein in the PPI network.
FIGURE 2. Protein-protein interaction network of the core proteins targeted by the new formulation for the treatment of ischemic stroke with Qi deficiency and blood stasis syndrome. The network included 244 nodes and 1,081 edges. Diamonds repr6+ the target proteins, with darker colors indicating increasing importance.
FIGURE 3. Key target genes of the new treatment formulation for ischemic stroke with Qi deficiency and blood stasis syndrome. Diamonds represent target genes, with darker colors indicating increasing importance.
GO Enrichment Analysis
To explore the functional distribution of the gene targets for the formulation, we imported the 327 predicted target genes into the Metascape database for GO enrichment analysis. These genes were found to be associated with multiple BP, CCT, and MF categories (p < 0.01; Figure 4). The main BP terms were involved in cellular response to nitrogen compound, cellular response to organonitrogen compound, circulatory system processes, the MAPK cascade, and ion homeostasis. The CCT terms included membrane rafts, membrane microdomains, post-synapse, the perinuclear region of the cytoplasm, and dendrites. The MF terms included protein kinase activity, phosphotransferase activity, alcohol group as acceptor kinase activity, transmembrane receptor protein tyrosine kinase activity, and transmembrane receptor protein kinase activity.
FIGURE 4. Gene Ontology (GO) enrichment analysis of the new treatment formulation for ischemic stroke with Qi deficiency and blood stasis syndrome. The top 10 biological processes, cellular components, and molecular functions are shown as a bar chart. The x-axis represents the GO terms and the y-axis represents the number of annotated genes.
KEGG Pathway Annotation
The 327 potential target genes were imported into the Metascape database for KEGG pathway enrichment analysis. Analysis showed that these targets were enriched in 191 pathways (p < 0.01). The top 20 pathways were analyzed using a bioinformatics database (Figure 5). Among these, the PI3K/AKT and MAPK signaling pathways were highly associated with the target genes, thus suggesting that the active ingredients of the new formulation exert beneficial effects on Qi and blood stasis by regulating the PI3K/AKT and MAPK signaling pathways.
FIGURE 5. The top 20 pathways identified by KEGG analysis of the new treatment formulation for ischemic stroke with Qi deficiency and blood stasis syndrome. The y-axis represents the name of each pathway, and the size of the dots corresponds to the number of genes annotated in the pathway. The x-axis represents the ratio corresponding to the significance of the differential protein enrichment, with higher ratios indicating higher reliability.
Experimental Verifications
The Neuroprotective Effects of HCSE on Ischemic Stroke
The volume of cerebral infarcts was determined by TTC staining. As shown in Figures 6A,B, no cerebral infarction was seen in the Sham and Sham + HCSE groups. However, the infarct volume in the MCAO group accounted for approximately 44% of the whole brain (p < 0.05). Compared with the MCAO group, the infarct volume in the MCAO + HCSE group was significantly smaller (approximately 23% of the brain volume; p < 0.05). In addition, the neurological function of the mice was assessed by applying the Longa (0–4) scale. As shown in Figure 6C, no neurological impairment was observed in the Sham and Sham + HCSE groups, although mice in the MCAO group had the highest score (approximately 3.64 points) (p < 0.05). Mice in the MCAO + HCSE group had significantly lower neurological function scores (approximately 2.35 points) compared to the MCAO group (p < 0.05).
FIGURE 6. The neuroprotective effects of HCSE. (A) TTC staining in the Sham, Sham + HCSE, MCAO and MCAO + HCSE groups. (B) The infarct volume was expressed as the ratio of (infarct volume/the whole brain volume) × 100%. (C) Longa score. (n = 7 per group; *p < 0.05, significantly different from the corresponding Sham group, #p < 0.05, significantly different from the MCAO group).
The Effect of HCSE on the Levels of AKT/P-AKT/ERK 1/2/p-ERK 1/2 Signaling-Related Proteins in the Hippocampus and Cortex With Ischemic Stroke in Mice
To investigate the possible mechanisms of HCSE against ischemic stroke, we investigated the expression levels of AKT/p-AKT/ERK 1/2/p-ERK 1/2 signaling related proteins by Western blotting. In hippocampal tissue, as shown in Figure 7, there was no significant difference in the p-AKT/AKT and p-ERK 1/2/ERK 1/2 protein ratios in the Sham + HCSE group when compared with the Sham group (p < 0.05). The p-AKT/AKT protein ratio in the MCAO group was lower than that in the Sham group (p < 0.05), while the p-ERK 1/2/ERK 1/2 protein ratio was higher than that in the Sham group (p < 0.05). However, the p-AKT/AKT protein ratio in the MCAO + HCSE group was significantly higher than that in the MCAO group (p < 0.05), and the p-ERK 1/2/ERK 1/2 protein ratio was significantly lower than that in the MCAO group (p < 0.05). Western blot analyses of cortical tissues revealed the same trend as the hippocampal tissue. As shown in Figure 8, the p-AKT/AKT protein ratio in the MCAO + HCSE group was significantly higher than that in the MCAO group (p < 0.05). In addition, the p-ERK 1/2/ERK 1/2 protein ratio was significantly lower than that in the MCAO group (p < 0.05).
FIGURE 7. Changes in the levels of proteins related to the AKT/p-AKT/ERK1/2/p-ERK 1/2 signaling pathway in the hippocampus. (A) Western blotting of related proteins. (B) The ratio of p-AKT/AKT. (C) The ratio of p-ERK 1/2/p-ERK 1/2. (n = 7 per group; *p < 0.05, significantly different from the corresponding Sham group, #p < 0.05, significantly different from the MCAO group).
FIGURE 8. Changes in the levels of proteins related to the AKT/p-AKT/ERK1/2/p-ERK1/2 signaling pathway in the cortex. (A) Western blotting of related proteins. (B) The ratio of p-AKT/AKT. (C) The ratio of p-ERK1/2/p-ERK 1/2. (n = 7 per group; *p < 0.05, significantly different from the corresponding Sham group, #p < 0.05, significantly different from the MCAO group).
Discussion
Although the treatment of stroke is rapidly evolving, research on the use of TCM to treat ischemic stroke is scarce. In the present study, we identified a new treatment formulation for ischemic stroke by the application of data mining and identified the mechanism of action of this treatment formulation using network pharmacology. The new formulation consists of Huangqi, Chuanxiong, SanLeng, and Ezhu. We concluded that the main active ingredients of the new formulation include trans-11-eicosenoic acid, ethyl linoleate, 7-O-methylisomucronulatol, myricanone, bisdemethoxycurcumin, wallichilide, astrapterocarpan, and kumatakenin. The active ingredients target a series of genes, including MAPK3, MAPK1, HSP90AA1, STAT3, PIK3R1, PIK3CA and AKT1, and can potentially regulate the PI3K/AKT, MAPK/ERK, Ras, cAMP and Rap1 signaling pathways to exert a neuroprotective role in patients with ischemic stroke, Qi deficiency and blood stasis syndrome. The cerebral protective effects of each of the drugs that constitutes HCSE have been extensively reported in previous studies; Huangqi promotes the proliferation of neural stem cells, Chuanxiong inhibits inflammation, Sanleng prevents the aggregation of platelets, and Ezhu inhibits autophagy (Huang et al., 2018; Fei Yin et al., 2020; Min Wang et al., 2020; Jia et al., 2021; Wang et al., 2021). In the present study, we found that HCSE effectively reduced cerebral infarcts and the neuronal dysfunction caused by ischemic stroke with Qi deficiency and blood stasis syndrome. KEGG pathway annotation suggested that the active ingredients of the new formulation exert beneficial effects by regulating the PI3K/AKT and MAPK signaling pathways. By analyzing the KEGG enrichment of these targets in ischemic stroke with Qi deficiency and blood stasis syndrome, we found that the PI3K/AKT and MAPK signaling pathway were the main pathways involved. Thus, by regulating the PI3K/AKT and MAPK signaling pathways, it may be possible to improve Qi and blood in patients with ischemic stroke.
The PI3K and MAPK signaling pathways are abnormally regulated in cerebral ischemia; by modulating these pathways, it is possible to alleviate the neuronal injury caused by ischemia (Tu et al., 2015). Research has also shown that HSP90 inhibitors protect against ischemia-induced neural progenitor cell death via the PI3K/AKT and MAPK/ERK pathways (Kwon et al., 2008; Wang et al., 2011; Bradley et al., 2014). On the one hand, the PI3K/AKT signaling pathway is known to play key roles in regulating cell proliferation, differentiation, apoptosis, and migration (Samakova et al., 2019) and is critical for neuronal growth and survival following cerebral ischemia (Zhao et al., 2016). Previous studies have reported that modulation of the PI3K/AKT/mTOR pathway upregulates the bcl-2 protein and increases ischemic tolerance in the semi-dark zone, thereby reducing apoptosis, increasing VEGF expression, and promoting cerebral angiogenesis (Hou et al., 2018; Liang et al., 2018). On the other hand, the extracellular signal-regulated kinases MAPK1/ERK2 and MAPK3/ERK1 are members of the MAP kinase family and are involved in cell proliferation, differentiation, transcriptional regulation, and apoptosis. A series of cascade reactions are known to be involved in ischemic stroke (Sun and Nan, 2016). Numerous previous studies have shown that the MAPK/ERK signaling pathway disrupts the blood-brain barrier, affects neurocyte apoptosis, and enhances the expression of neuronal inflammatory factors after ischemic stroke (Irving et al., 2000; Cao et al., 2016; Wang et al., 2019). The downregulation of MAPK1 has been shown to reduce the levels of TNF-α, IL-6, and reactive oxygen species, thereby reducing neuroinflammation, oxidative stress, and neuronal damage (Zhang et al., 2020). In a previous study, Mostajeran found that inhibition of the ERK signaling pathway by U0126 reduced neuronal death and significantly upregulated the expression of Tie-2, thereby promoting post-stroke vascular stabilization and angiogenesis (Mostajeran et al., 2017).
The cascade response of the PI3K and MAPK signaling pathways is regulated by complex feedback and crosstalk mechanisms. Zhou et al. reported negative crosstalk between the MAPK and PI3K/AKT signaling pathways, with AKT inhibiting the MAPK signaling pathway by phosphorylating and inhibiting the Raf1 node during cerebral ischemia (Zhou et al., 2015). This is consistent with the results of the present study, in which phosphorylated AKT levels were significantly lower, and phosphorylated ERK 1/2 levels were significantly higher, in mice suffering from ischemic stroke. Levels of phosphorylated AKT were increased while levels of phosphorylated ERK 1/2 were decreased in mice with ischemic stroke when treated with HCSE. We suggest that HCSE protects against post-ischemic injury not only by reducing ERK activity, but also by increasing crosstalk between AKT and ERK. The reduction in ERK activity exerts cerebral protective effects.
In summary, we used data mining and network pharmacological target prediction to identify a new treatment formulation (HCSE) for patients suffering from ischemic stroke with Qi deficiency and blood stasis syndrome. HCSE significantly reduced infarct volume and improved neurological function. We found that multiple target genes and pathways participated in the action of HCSE against cerebral ischemia stroke. The mechanisms underlying the neuroprotective effects of HCSE were closely related to activation of the PI3K/AKT and MAPK/ERK signaling pathways. Our results provided evidence for the positive effects of HCSE on ischemic stroke with Qi deficiency and blood stasis syndrome. However, further research now needs to verify the mechanisms underlying the action of HSCE on the improvement of cerebral ischemic injury via the regulation of the PI3K/AKT and MAPK/ERK pathways; it is possible that other neuroprotective mechanisms are involved. In addition, it is important to investigate the internal mechanisms that link Qi deficiency, blood stasis syndrome and cerebral ischemic injury, as this could provide a more enhanced theoretical basis for the prevention and treatment of stroke with TCM.
Data Availability Statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics Statement
The animal study was reviewed and approved by And all experimental investigation procedures were authorized by the Yangzhou University Institutional Animal Care and Use Committee (YIACUC-14–0015).
Author Contributions
LZ performed the experiments, LZ and LDD drafted the manuscript. ZHX, PS, MMS, LDD, and ZC analyzed the data; BCY designed the experiments. All authors have discussed the results, reviewed the manuscript and approved the final version for publication.
Funding
This work was partially supported by the National Natural Science Foundation of China (81401005).
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
Bardou, P., Mariette, J., Escudié, F., Djemiel, C., and Klopp, C. (2014). Jvenn: an Interactive Venn Diagram Viewer. BMC Bioinformatics 15, 293. doi:10.1186/1471-2105-15-293
Bradley, E., Zhao, X., Wang, R., Brann, D., Bieberich, E., and Wang, G. (2014). Low Dose Hsp90 Inhibitor 17AAG Protects Neural Progenitor Cells from Ischemia Induced Death. J. Cel. Commun Signal 8 (4), 353–362. doi:10.1007/s12079-014-0247-5
Cao, G., Jiang, N., Hu, Y., Zhang, Y., Wang, G., Yin, M., et al. (2016). Ruscogenin Attenuates Cerebral Ischemia-Induced Blood-Brain Barrier Dysfunction by Suppressing TXNIP/NLRP3 Inflammasome Activation and the MAPK Pathway. Int. J. Mol. Sci. 17 (9), 1418. doi:10.3390/ijms17091418
Daina, A., Michielin, O., and Zoete, V. (2019). SwissTargetPrediction: Updated Data and New Features for Efficient Prediction of Protein Targets of Small Molecules. Nucleic Acids Res. 47 (W1), W357–W364. doi:10.1093/nar/gkz382
Fei Yin, F., Zhou, H., Fang, Y., Li, C., He, Y., Yu, L., et al. (2020). Astragaloside IV Alleviates Ischemia Reperfusion-Induced Apoptosis by Inhibiting the Activation of Key Factors in Death Receptor Pathway and Mitochondrial Pathway. J. Ethnopharmacol. 248, 112319. doi:10.1016/j.jep.2019.112319
Hou, Y., Wang, K., Wan, W., Cheng, Y., Pu, X., and Ye, X. (2018). Resveratrol Provides Neuroprotection by Regulating the JAK2/STAT3/PI3K/AKT/mTOR Pathway after Stroke in Rats. Genes Dis. 5 (3), 245–255. doi:10.1016/j.gendis.2018.06.001
Huang, L., Chen, C., Zhang, X., Li, X., Chen, Z., Yang, C., et al. (2018). Neuroprotective Effect of Curcumin against Cerebral Ischemia-Reperfusion via Mediating Autophagy and Inflammation. J. Mol. Neurosci. 64 (1), 129–139. doi:10.1007/s12031-017-1006-x
Irving, E. A., Barone, F. C., Reith, A. D., Hadingham, S. J., and Parsons, A. A. (2000). Differential Activation of MAPK/ERK and P38/SAPK in Neurones and Glia Following Focal Cerebral Ischaemia in the Rat. Brain Res. Mol. Brain Res. 77 (1), 65–75. doi:10.1016/s0169-328x(00)00043-7
Jia, J., Li, X., Ren, X., Liu, X., Wang, Y., Dong, Y., et al. (2021). Sparganii Rhizoma: A Review of Traditional Clinical Application, Processing, Phytochemistry, Pharmacology, and Toxicity. J. Ethnopharmacol. 268, 113571. doi:10.1016/j.jep.2020.113571
Kwon, H. M., Kim, Y., Yang, S. I., Kim, Y. J., Lee, S. H., and Yoon, B. W. (2008). Geldanamycin Protects Rat Brain through Overexpression of HSP70 and Reducing Brain Edema after Cerebral Focal Ischemia. Neurol. Res. 30 (7), 740–745. doi:10.1179/174313208X289615
Li, S., and Zhang, B. (2013). Traditional Chinese Medicine Network Pharmacology: Theory, Methodology and Application. Chin. J. Nat. Med. 11 (2), 110–120. doi:10.1016/S1875-5364(13)60037-0
Li, S. M., Xu, H., and Chen, K. J. (2014). The Diagnostic Criteria of Blood-Stasis Syndrome: Considerations for Standardization of Pattern Identification. Chin. J. Integr. Med. 20 (7), 483–489. doi:10.1007/s11655-014-1803-9
Liang, Z., Chi, Y. J., Lin, G. Q., Luo, S. H., Jiang, Q. Y., and Chen, Y. K. (2018). MiRNA-26a Promotes Angiogenesis in a Rat Model of Cerebral Infarction via PI3K/AKT and MAPK/ERK Pathway. Eur. Rev. Med. Pharmacol. Sci. 22 (11), 3485–3492. doi:10.26355/eurrev_201806_15175
Liu, S. Z., Bao, Z. X., and Zhang, R. L. (2007). Theory Discuss of the Pathogenesis Theory of Deficiency of Qi and Blood Stasis in Lschemic Stroke. Chin. Arch. Tradit. Chin. 25, 97–98. doi:10.13193/j.archtcm.2007.01.99.liuszh.040
Liu, Y., Yin, H. J., Shi, D. Z., and Chen, K. J. (2012). Chinese Herb and Formulas for Promoting Blood Circulation and Removing Blood Stasis and Antiplatelet Therapies. Evid. Based Complement. Alternat Med. 2012, 184503. doi:10.1155/2012/184503
Liu, D., Bai, X., Ma, W., Xin, D., Chu, X., Yuan, H., et al. (2020). Purmorphamine Attenuates Neuro-Inflammation and Synaptic Impairments after Hypoxic-Ischemic Injury in Neonatal Mice via Shh Signaling. Front. Pharmacol. 11, 204. doi:10.3389/fphar.2020.00204
Liu, J., Wang, F., Sheng, P., Xia, Z., Jiang, Y., and Yan, B. C. (2021). A Network-Based Method for Mechanistic Investigation and Neuroprotective Effect on Treatment of Tanshinone Ⅰ against Ischemic Stroke in Mouse. J. Ethnopharmacol. 272, 113923. doi:10.1016/j.jep.2021.113923
Wang, M., Yao, M., Liu, J., Takagi, N., Yang, B., Zhang, M., et al. (2020). Ligusticum Chuanxiong Exerts Neuroprotection by Promoting Adult Neurogenesis and Inhibiting Inflammation in the hippocampus of ME Cerebral Ischemia Rats. J. Ethnopharmacol. 249, 112385. doi:10.1016/j.jep.2019.112385
Mostajeran, M., Edvinsson, L., Warfvinge, K., Singh, R., and Ansar, S. (2017). Inhibition of Mitogen-Activated Protein Kinase 1/2 in the Acute Phase of Stroke Improves Long-Term Neurological Outcome and Promotes Recovery Processes in Rats. Acta Physiol. (Oxf) 219 (4), 814–824. doi:10.1111/apha.12632
Mukundan, G., and Seidenwurm, D. J. (2018). Economic and Societal Aspects of Stroke Management. Neuroimaging Clin. N. Am. 28 (4), 683–689. doi:10.1016/j.nic.2018.06.009
Reagan-Shaw, S., Nihal, M., and Ahmad, N. (2008). Dose Translation from Animal to Human Studies Revisited. FASEB J. 22 (3), 659–661. doi:10.1096/fj.07-9574LSF
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a Database of Systems Pharmacology for Drug Discovery from Herbal Medicines. J. Cheminform 6, 13. doi:10.1186/1758-2946-6-13
Safran, M., Rosen, N., Twik, M., BarShir, R., Stein, T. I., Dahary, D., et al. (2021). “The GeneCards Suite,” in Practical Guide to Life Science Databases. Editors I. Abugessaisa, and T. Kasukawa. (Singapore: Springer), 27–56. doi:10.1007/978-981-16-5812-9_2
Samakova, A., Gazova, A., Sabova, N., Valaskova, S., Jurikova, M., and Kyselovic, J. (2019). The PI3k/Akt Pathway Is Associated with Angiogenesis, Oxidative Stress and Survival of Mesenchymal Stem Cells in Pathophysiologic Condition in Ischemia. Physiol. Res.68, S131–S138. doi:10.33549/physiolres.934345
Sun, J., and Nan, G. (2016). The Mitogen-Activated Protein Kinase (MAPK) Signaling Pathway as a Discovery Target in Stroke. J. Mol. Neurosci. 59 (1), 90–98. doi:10.1007/s12031-016-0717-8
Sun, J. H., Sun, F., Yan, B., Li, J. Y., and Xin, L. (2020). Data Mining and Systematic Pharmacology to Reveal the Mechanisms of Traditional Chinese Medicine in Mycoplasma Pneumoniae Pneumonia Treatment. Biomed. Pharmacother. 125, 109900. doi:10.1016/j.biopha.2020.109900
Tang, Y., Li, M., Wang, J., Pan, Y., and Wu, F. X. (2015). CytoNCA: a Cytoscape Plugin for Centrality Analysis and Evaluation of Protein Interaction Networks. Biosystems 127, 67–72. doi:10.1016/j.biosystems.2014.11.005
Tu, X. K., Yang, W. Z., Chen, J. P., Chen, Y., Chen, Q., Chen, P. P., et al. (2015). Repetitive Ischemic Preconditioning Attenuates Inflammatory Reaction and Brain Damage after Focal Cerebral Ischemia in Rats: Involvement of PI3K/Akt and ERK1/2 Signaling Pathway. J. Mol. Neurosci. 55 (4), 912–922. doi:10.1007/s12031-014-0446-9
Wang, G., Krishnamurthy, K., and Tangpisuthipongsa, D. (2011). Protection of Murine Neural Progenitor Cells by the Hsp90 Inhibitor 17-Allylamino-17-Demethoxygeldanamycin in the Low Nanomolar Concentration Range. J. Neurochem. 117 (4), 703–711. doi:10.1111/j.1471-4159.2011.07239.x
Wang, W., Jiang, B., Sun, H., Ru, X., Sun, D., Wang, L., et al. (2017). Prevalence, Incidence, and Mortality of Stroke in China: Results from a Nationwide Population-Based Survey of 480 687 Adults. Circulation 135 (8), 759–771. doi:10.1161/CIRCULATIONAHA.116.025250
Wang, Q. C., Lu, L., and Zhou, H. J. (2019). Relationship between the MAPK/ERK Pathway and Neurocyte Apoptosis after Cerebral Infarction in Rats. Eur. Rev. Med. Pharmacol. Sci. 23 (12), 5374–5381. doi:10.26355/eurrev_201906_18206
Wang, Y., Liu, X., Hu, T., Li, X., Chen, Y., Xiao, G., et al. (2021). Astragalus Saponins Improves Stroke by Promoting the Proliferation of Neural Stem Cells through Phosphorylation of Akt. J. Ethnopharmacol 277, 114224. doi:10.1016/j.jep.2021.114224
Chen, X., Chen, H., He, Y., Fu, S., Liu, H., Wang, Q., et al. (2020). Proteomics-Guided Study on Buyang Huanwu Decoction for its Neuroprotective and Neurogenic Mechanisms for Transient Ischemic Stroke: Involvements of EGFR/PI3K/Akt/Bad/14-3-3 and Jak2/Stat3/Cyclin D1 Signaling Cascades. Mol. Neurobiol. 57 (10), 4305–4321. doi:10.1007/s12035-020-02016-y
Xi, S., Qian, L., Tong, H., Yue, L., Zhao, H., Wang, D., et al. (2013). Toxicity and Clinical Reasonable Application of Taoren (Semen Persicae) Based on Ancient and Modern Literature Research. J. Tradit Chin. Med. 33 (2), 272–279. doi:10.1016/s0254-6272(13)60139-9
Xu, H. Y., Zhang, Y. Q., Liu, Z. M., Chen, T., Lv, C. Y., Tang, S. H., et al. (2019). ETCM: an Encyclopaedia of Traditional Chinese Medicine. Nucleic Acids Res. 47 (D1), D976–D982. doi:10.1093/nar/gky987
Chen, Y., Wright, N., Guo, Y., Turnbull, I., Kartsonaki, C., Yang, L., et al. (2020). Mortality and Recurrent Vascular Events after First Incident Stroke: a 9-year Community-Based Study of 0·5 Million Chinese Adults. Lancet Glob. Health 8 (4), e580–e590. doi:10.1016/S2214-109X(20)30069-3
Wang, Y., Zhang, L., Pan, Y. J., Fu, W., Huang, S. W., Xu, B., et al. (2020). Investigation of Invigorating Qi and Activating Blood Circulation Prescriptions in Treating Qi Deficiency and Blood Stasis Syndrome of Ischemic Stroke Patients: Study Protocol for a Randomized Controlled Trial. Front. Pharmacol. 11, 892. doi:10.3389/fphar.2020.00892
Yin, B., Bi, Y. M., Fan, G. J., and Xia, Y. Q. (2020). Molecular Mechanism of the Effect of Huanglian Jiedu Decoction on Type 2 Diabetes Mellitus Based on Network Pharmacology and Molecular Docking. J. Diabetes Res. 2020, 5273914. doi:10.1155/2020/5273914
Zhai, X., Wang, X., Wang, L., Xiu, L., Wang, W., and Pang, X. (2020). Treating Different Diseases with the Same Method-A Traditional Chinese Medicine Concept Analyzed for its Biological Basis. Front. Pharmacol. 11, 946. doi:10.3389/fphar.2020.00946
Zhang, L., Wang, J., Liu, Q., Xiao, Z., and Dai, Q. (2020). Knockdown of Long Non-coding RNA AL049437 Mitigates MPP+ -induced Neuronal Injury in SH-Sy5y Cells via the microRNA-205-5p/MAPK1 axis. Neurotoxicology 78, 29–35. doi:10.1016/j.neuro.2020.02.004
Zhang, J., Rui, Y., Gao, M., Wang, L., and Yan, B. C. (2021). Expression of Long Non-coding RNA RGD1566344 in the Brain Cortex of Male Mice after Focal Cerebral Ischemia-Reperfusion and the Neuroprotective Effect of a Non-coding RNA RGD1566344 Inhibitor. Cell Mol Neurobiol 41 (4), 705–716. doi:10.1007/s10571-020-00877-4
Zhao, E. Y., Efendizade, A., Cai, L., and Ding, Y. (2016). The Role of Akt (Protein Kinase B) and Protein Kinase C in Ischemia-Reperfusion Injury. Neurol. Res. 38 (4), 301–308. doi:10.1080/01616412.2015.1133024
Zhou, J., Du, T., Li, B., Rong, Y., Verkhratsky, A., and Peng, L. (2015). Crosstalk between MAPK/ERK and PI3K/AKT Signal Pathways during Brain Ischemia/Reperfusion. ASN Neuro 7 (5), 1759091415602463. doi:10.1177/1759091415602463
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun. 10 (1), 1523. doi:10.1038/s41467-019-09234-6
Glossary
AKT Protein Kinase B
BBB Blood-Brain Barrier
BC Betweenness Centrality
BCA Bicinchonininc Acid
bcl-2 B-Cell Lymphoma-2
BP Biological Process
BSA Bovine Serum Albumin
cAMP Cyclic Adenosine Monophosphate
CC Cellular Component Closeness Centrality
CC Cellular Component Closeness Centrality
CCA Common Carotid artery
DC Degree Centrality
DL Drug-Likeness
EC Eigenvector Centrality
ECA External Carotid artery
ERK Extracellular Regulated Protein Kinases 1 and 2
FDR False Discovery Rate
GO Gene Ontology
HCSE Huangqi-Chuanxiong and Sanleng-Ezhu
HSP90 Heat Shock Protein 90
HSP90AA1 Heat Shock Protein 90 Alpha Family Class A Member 1
ICA Internal Carotid artery
ICR Institute of Cancer Research
IL6 Interleukin 6
KEGG Kyoto Encyclopedia of Genes and Genomes
MAPK Mitogen-Activated Protein Kinases
MAPK1 Mitogen-Activated Protein Kinase 1
MAPK3 Mitogen-Activated Protein Kinase 3
MCAO Middle Cerebral artery Occlusion
MF Molecular Functions
mTOR Mechanistic Target of Rapamycin Kinase
OB Oral Bioavailability
p-AKT Phospho-Protein kinase B
p-ERK1/2 Phospho-Extracellular Regulated Protein Kinases 1 and 2
PI3K Phosphoinositide 3-kinase
PIK3CA Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha
PIK3R1 Phosphoinositide-3-Kinase Regulatory Subunit 1
PPI Protein-Protein Interaction
Raf1 Raf-1 Proto-Oncogene, Serine/Threonine Kinase
Ras Rat Sarcoma
SDS-PAGE Sodium Dodecyl Sulfate Polyacrylamide Gel Electrophoresis
STAT3 Signal Transducer and Activator of Transcription 3
TBS Tris Buffered Saline
TCM Traditional Chinese Medicine
Tie-2 Tyrosine-Protein Kinase
TNF-α Tumor Necrosis Factor-α
TTC 2, 3, 5-Triphenyl Tetrazolium Chloride
VEGF Vascular Endothelial Growth Factor
Keywords: huangqi-chuanxiong, sanleng-ezhu herb, ischemic stroke, clinical data mining, network pharmacology
Citation: Zhao L, Ding LD, Xia ZH, Sheng P, Shen MM, Cai ZM and Yan BC (2022) A Network-Based Approach to Investigate the Neuroprotective Effects and Mechanisms of Action of Huangqi-Chuanxiong and Sanleng-Ezhu Herb Pairs in the Treatment of Cerebral Ischemic Stroke. Front. Pharmacol. 13:844186. doi: 10.3389/fphar.2022.844186
Received: 27 December 2021; Accepted: 18 February 2022;
Published: 23 March 2022.
Edited by:
Li-Nan Zhang, Hebei Medical University, ChinaReviewed by:
Yucong Peng, Zhejiang University, ChinaFeng Zhang, Third Hospital of Hebei Medical University, China
Copyright © 2022 Zhao, Ding, Xia, Sheng, Shen, Cai and Yan. 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: Bing Chun Yan, YmN5YW5AeXp1LmVkdS5jbg==
†These authors have contributed equally to this work