Elucidating the anti-atherosclerosis mechanism of Si Jun Zi Decoction by integrating network pharmacology and transcriptomic with experimental validation in vivo

Elucidating the anti-atherosclerosis mechanism of Si Jun Zi Decoction by integrating network pharmacology and transcriptomic with experimental validation in vivo | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Elucidating the anti-atherosclerosis mechanism of Si Jun Zi Decoction by integrating network pharmacology and transcriptomic with experimental validation in vivo Weiyan Chen, Tianmin Ji, Ying Yang, Zhuo Zhao, Hao Gao, Lijiang Zhou, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4872904/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Objective: Sijunzi Decoction (SJZD) is a kind of traditional Chinese medicine (TCM) formula,which has the contribution to anti-atherosclerosis.This study aims to explore the potential mechanism of the treatment of atherosclerosis (AS) with SJZD. Methods: Comprehensive analysis with transcriptomics and network pharmacology,combining with in vivo experiment and therapeutic targets. Results: Through the in vivo experiment,it was found that SJZD could significantly improve the blood lipid level and lipid deposition in ApoE−/− mice.Furthermore,4 biological processes(lipid localization,specific granule,positive regulation of lipid localization,positive regulation of lipid transport),7 targets (SPP1,EGF,OLR1,LDLR,PON1,SLP1,PLAU), and 6 active ingredients(including beta-sitosterol,kaempferol,quercetin,stigmasterol, naringenin and isorhamnetin) play an extremely important role in the treatment of AS with SJZD.We also found that the genes of OLR,SPP1 and EGF were involved in regulating oxidative stress during the progress of atherosclerosis and the improvement of endothelial dysfunction.Furthermore,SJZD could interfere the expression of mRNA and protein of OLR1,SPP1 and EGF,and respectively reduce and increase the density level of MDA and SOD in serum significantly. Conclusion: Inhibiting excessively oxidative stress and improving endothelial dysfunction by regulating OLR1/SPP1/EGF pathway could be the mechanism,by which Sijunzi decoction resists AS. network pharmacology transcriptomic techniques TCM Sijunzi Decoction atherosclerosis OLR1/SPP1/EGF pathway Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 1. Introduction Atherosclerosis (AS) is one of the most common vascular disfunction all over the world and could lead to stroke and coronary heart disease.It is associated with health risks that are major causes of death [ 1 ] .Previous studies have demonstrated that AS has a connection with endothelial cell damage,inflammation,lipid metabolism disorders,and immune dysfunction [ 2 ] .Endothelial cell injury,which results in plaque, leads to further instability and rupture of it, and secondary thrombosis will accelerate the progression of the disease [ 3 ] .Vascular endothelial cell injury caused by oxidative stress is an significant factor and an early pathological change in the development of AS [ 4 – 5 ] .In China,Chinese herbal medicine (CHM) has been used to treat atherosclerotic complications for a long time,such as angina pectoris,myocardial infarction and stroke [ 6 ] .However,the molecular mechanism is still unclear.Sijunzi Decoction (SJZD),consisting of Atractylodes macrocephala Panax ginseng,Glycyrrhiza uralensis and Poria cocos,has been applied for thousands of years [ 7 ] . SJZD as a classic prescription for treating spleen deficiency syndrome,has obvious anti-oxidative stress effect [ 8 ] .However,the composition of traditional Chinese medicine (TCM) is complex,with many active components and function through multiple pathways. To further explore its mechanism, network pharmacology and transcriptomic techniques are essential to be introduced. Network pharmacology,which is similar to the overall concept of TCM,constructs a drug-component-target-disease regulatory network through joint analysis of disease-related genes [ 9 ] .It is helpful to elucidate the mechanism of drug comprehensively as much as possible and interpret the synergistic effect of multiple components and targets of TCM [ 10 ] .It is also used to evaluate multi-target drugs and complex diseases [ 11 ] .Transcriptomics is a subject that elaborates gene transcription and regulation in tissues and cells with a holistic view,and could comprehensively explain the regulation mechanism of gene expression in complex diseases by integrating a variety of calculation ways.For the study of gene transcription and regulation in cells or tissues [ 12 ] , the combined application of network pharmacology and transcriptomics could preliminarily screen the active ingredients and molecular targets of drugs,and provide convenient conditions for subsequent research [ 10 ] .Based on the above, the aim of our research is to find the possible mechanism and molecular targets of SJZD,which could improve vascular endothelial injury and prevent AS,mainly through network pharmacology and transcriptomic technology. 2. Materials and Methods 2.1 Experimental animal 16 SPF grade male ApoE − / − mice (body weight, 18–20 g) and 8 male C57BL/6 mice were provided b y Beijing Vital River Laboratory Animal Technology Co.Ltd. The Committee on the Ethics of Animal Experiments at Liaoning University of Traditional Chinese Medicine approved this study(permitted number SCXK(Jing) 2012-0001).All mice were raised in the Animal Experimental Center of Liaoning University of Traditional Chinese Medicine and allowed free movement with food and water ad libitum. 2.2 Establishment of AS model mice and treatments After 1 week of adaptive feeding in SPF grade mice,8 C57BL/6 mice were considered as a normal group(NG),and 16 ApoE −/− mice were randomly assigned to the model group(MG) and SJZ Decoction group(SJZD).The normal group was fed with basal diet daily,while the other groups were fed with high-fat diet(containing 21% fat, 10% lard, 1% cholesterol, etc.)with free access to food and water for 12 weeks.The drug intervention was given to the group,which the ApoE −/− mice treated with 20 g/kg/d SJZ Decoction for the last 4 weeks.The normal group and model group were given an equal volume of normal saline daily by gavage once a day. 2.3 Serum and tissue samples After the treatment period,all mice were fasted overnight and anesthetized with chloral hydrate.Eyeball blood was collected from sacrificed mice after resting for 30 min and subjected to centrifugation at 3000× g and 4°C for 15 min in order to obtain serum, and stored in aliquots at -80°C for further analysis.Quickly separate the aortic arch,thoracic-abdominal aorta,and rinse 3 times with ice-cold saline.The aorta was fixed in 4% paraformaldehyde for atherosclerotic lesion analysis.The remaining part is frozen at -80°C until it is used for biochemical determination,iron content analysis and RNA extraction. 2.4 Serum lipids The levels of serum TC,TG,HDL-C,and LDL-C were measured using an automated biochemical analyzer(SIEMENS,Munich,DE,Germany). 2.5 Histopathology test The arterial tissue was fixed in 4% paraformaldehyde solution overnight,and stained using the conventional method of hematoxylin-eosin(HE)staining.The specific steps are as follows:ⅰ.The tissue is dehydrated,transparent,embedded in paraffin, and then sectioned(thickness:~5 microns). Ⅱ.Paraffin sections are dewaxed with xylene, rehydrated with gradient alcohol,stained with hematoxylin,differentiated with hydrochloric acid alcohol,counterstained with eosin, then dehydrated,transparent with xylene,and mounted with neutral gum. Ⅲ.Observe the morphology of rat aorta under a microscope. Oil Red O staining of aorta: complete aorta was isolated and fixed in 10% neutral formalin for 24 h.The complete aorta was in 60% isopropanol for 10min,stained in oil red O (Solarbio,China) working fluid (3:2) for 3h,washed in 60% isopropanol for 6 times till the background color became white,and then photographed. 2.6 Transcriptomic sequencing and data analysis 2.6.1 RNA extraction and library preparation The aorta tissues samples were randomly pooled from the mice in the NG,MG,and SJZD (n = 3).The total RNA was extracted using 1mL of Trizol Reagent (Invitrogen,Carlsbad,CA,USA) according to the manufacturer’s instructions.The total RNA was reverse transcribed into cDNA.NEBNext® UltraTM RNA Library Preparation Kit for Illumina(NEB,USA) was used to generate sequencing libraries.HiSeq X instrument was used for whole transcriptome sequencing.Cut adapt was used to filter the raw data to get clean data for further analysis. 2.6.2 Differential expression genes (DEGs) analysis Two comparisons-MG versus NG, and SJZD vs MG were executed application of the DESeq2 R package (1.16.1) to determine differential expression genes.DESeq2 provided statistical routines for defining various expression in digital gene expression data application model according to the negative binomial distribution.The significance of DESeq2 was determined by p-value 1,which was implemented to define DEGs. 2.6.3 PPI network construction and analysis We used the Venn map to scalp the DEGs,which were down-regulated by MG and up-regulated by SJZD, or DEGs up-regulated by MG and down-regulated by SJZD. The screened DEGs might be the latent targets of SJZD to treat AS. 2.6.4 Functional enrichment analysis The GO and KEGG pathways analysis were enriched by the R software package to realize the functions of candidate DEGs.In our research,we used the cluster Profiler to analyze the GO and KEGG pathways which was focused on were presented. 2.6 Transcriptomic sequencing and data analysis 2.6.1 RNA extraction and library preparation The aorta tissues samples were randomly pooled from the mice in the NG,MG,and SJZD (n = 3).The total RNA was extracted using 1mL of Trizol Reagent (Invitrogen,Carlsbad,CA,USA) according to the manufacturer’s instructions.The total RNA was reverse transcribed into cDNA.NEBNext® UltraTM RNA Library Preparation Kit for Illumina(NEB,USA) was used to generate sequencing libraries.HiSeq X instrument was used for whole transcriptome sequencing.Cut adapt was used to filter the raw data to get clean data for further analysis. 2.6.2 Differential expression genes (DEGs) analysis Two comparisons-MG versus NG, and SJZD vs MG were executed application of the DESeq2 R package (1.16.1) to determine differential expression genes.DESeq2 provided statistical routines for defining various expression in digital gene expression data application model according to the negative binomial distribution.The significance of DESeq2 was determined by p-value 1,which was implemented to define DEGs. 2.6.3 PPI network construction and analysis We used the Venn map to scalp the DEGs,which were down-regulated by MG and up-regulated by SJZD, or DEGs up-regulated by MG and down-regulated by SJZD. The screened DEGs might be the latent targets of SJZD to treat AS. 2.6.4 Functional enrichment analysis The GO and KEGG pathways analysis were enriched by the R software package to realize the functions of candidate DEGs.In our research,we used the cluster Profiler to analyze the GO and KEGG pathways which was focused on were presented. 2.7 Network pharmacology 2.7.1 Gatheration of active compounds and acquisition of corresponding targets of SJZ decoction The active ingredients of four herbal medicines in SJZ decoction were obtained from TCMSP and the BATMAN-TCM database.According to the recommended drug screening criteria of the TCMSP database,chemical compounds with OB ≥ 30% and DL ≥ 0.18 may considered active,and they were choosed as effective constituents for further analysis through the TCMSP database.Targets of active ingredients were acquired from the TCMSP database.The target protein names corresponding to active compounds were normalized in UniProt. 2.7.2 Collect the common targets of SJZ Decoction and AS Using “atherosclerosis” as a keyword,the predicted genes of AS were collected from Gene Cards database, CTD database and TTD database.The search results found in three database were merged and deleted duplicate targets to obtain all target genes of AS.The common targets of SJZ Decoction and AS were identified by the bioinformatics online tool ( http://www.bioinformatics.com.cn/).I n short,the Venn map was used to draw the targets of SJZ Decoction and the disease targets of AS.Next,the core targets were considered as the potential targets for SJZ Decoction in treatment of AS. 2.7.3 Network construction and Functional enrichment The core targets were entered into String online tool to obtain the interrelated information of protein-protein interactions (PPI).Next, Cytoscape 3.7.2 was used to establish the visualized PPI network and calculate the degree centrality. Moreover, the key proteins of the PPI network were acquired. The GO and KEGG pathways of the potential proteins were enriched and analyzed by R 3.6.3. The results were analized by using p-value and count values. 2.8 Structure ingredient-gene-pathway network (IGP) We assessed pathways in the GO-biological process,which appears together in network pharmacology and transcriptome research.The genes related to this category were found in the transcriptome results,and the gene-pathway network was established.According to the network pharmacology results,the related active molecules of genes in this part were analyzed.Ingredients,genes,and GO-biological process were also put into Cytoscape tool to demonstrate a visualized Ingredient-gene-pathway (IGP) network. 2.9 Validation of compound-target interaction To validate the compound-target association,the molecular docking program was performed with Auto Dock Tools software(version 4.2).All the 3D structures of core target proteins were obtained from the PDB database ( http://www.rcsb.org/),s o that proteins and ligands could be prepared in the Auto Dock before performing the docking progress.The whole molecular docking process included setup of proteins,determination of docking sites and docking of proteins to small molecules.Firstly, during the process of protein preparation,water molecules were removed from the protein structures as well as the region arround the amide portion and atomic groups.Next,the protein binding sites were defined and edited.Finally,insert the compounds into the protein-binding site.If the small molecule docked with the protein successfully,it showed a labelled interaction site and a corresponding docking score.Hence,target proteins and compounds with high docking scores were choosed for further analysis. 2.10 Quantitative real-time PCR to examine mRNA expression TRIZOL reagent (Solarbio,China) was used to extract total RNA from the sample.The cDNA template was synthesized using a commercially available reverse transcription reaction (Bio-Rad,USA),and the quantitative RT-PCR experiment was performed using an ABI7500 rapid quantitative PCR instrument (Applied Biosystems, USA).The reaction procedure is:pre-denaturation at 95°C for 5 min,denaturation at 95°C for 20 s,annealing at 58°C for 30 s,extension at 72°C for 20 s,40 cycles.The relative quantification (2 −ΔΔCt ) method was used to calculate the relative amount of mRNA. The PCR primers were as follows:EGF,5’-ACAGAAGGAGTAGATACGCTTG-3’-(forward),and 5’-GATTATTCGATGATGCTTCCCG-3’-(reverse).OLR1,5’-GAAGCCTGCGAATGACGAGC-3’-(forward),5’-ACACCAGGCAGAGGATGACC-3’-(reverse).COL1,5’-TGAACGTGGTGTACAAGGTC-3’-(forward),and 5’-CCATCTTTACCAGGAGAACCAT-3’-(reverse).PON1,5’-AAGAGGAAAGATCACTCTTGCA-3’-(forward),and 5’-GGTCCAATAGCAGCTATATCGT-3’-(reverse).SLP1,5’-ATGTATGATGCTTAACCCTCCC-3’-(forward),and 5’-AGGCAGACTTTCCCACATATAC-3’-(reverse).SPP1,5’-AAACACACAGACTTGAGCATTC-3’-(forward),and 5’-TTAGGGTCTAGGACTAGCTTGT-3’-(reverse).PLAU,5’-GCTTGTTTCTCATGAACAGTGT-3’-(forward),and 5’-TTCGATGTTACAGATAAGCGGT-3’-(reverse).MMP2,5’-ACTTTGAGAAGGATGGCAAGTA-3’-(forward),and 5’-CTTCTTATCCCGGTCATAGTCC-3’-(reverse).LDLR1,5’-GAAGCCTGCGAATGACGAGC-3’-(forward),and 5’-ACACCAGGCAGAGGATGACC-3’-(reverse).And GAPDH,5’-GACATGCCGCCTGGAGAAAC-3’-(forward),and 5’-AGCCCAGGATGCCCTTTAGT-3’-(reverse). 2.11 Western blot analysis Total proteins were extracted from tissues using RIPA lysis buffer containing 1% phenylmethylsulphonyl fuoride (PMSF) for 30 minute on ice.The supernatants were collected after being clarified by centrifuging at 12000 rpm for 15 minutes at 4°C. Protein concentrations were measured by the BCA protein assay.The same amounts of protein (60 µg) were separated by electrophoresis on 10%-12% sodium dodecyl sulphate polyacrylamide gels and transferred onto PVDF membranes.Next,PVDF membranes were incubated overnight at 4°C with the antibody(proteintech,China). After washing with TBST 3 times,the membranes were incubated for 1 h with fluorescent secondary antibody at room temperature.Then be washed with TBST again, the membranes were scanned using a fluorescent scanner. 2.12 Immunohistochemistry(IHC)analysis IHC was performed on aorta tissue sections.Each sample was fixed in formalin and embedded in paraffin.The blocks were sliced into 5 µm-thick sections,After the tissue sections were dewaxed,antigen retrieval was performed,hydrated in a graded series of alcohols,and subjected to heat-activated antigen retrieval.After blocking endogenous peroxidase activity,the tissue was incubated with OLR1,SPP1 and EGF antibodies at room temperature for 4 h.Subsequently,the sections were washed and incubated with biotinylated secondary antibody at room temperature for 30 min.The reaction complexes were visualized with diaminobenzidine and counterstained with hematoxylin.Observe the morphology under the microscope. 2.13 Enzyme-linked immunosorbent assay (ELISA) analysis The serum SOD and MDA were determined using an ELISA kit according to the manufacturer's instructions. Specifically, after separating and collecting mice blood sample, the serum was obtained by centrifugation. Use the following commercial kits to evaluate various substances in serum by ELISA analysis: MDA (Mouse MDA ELISA KIT; Shanghai Enzyme Linked Biotechnology Co., Ltd.), SOD (Mouse SOD ELISA KIT; Shanghai Enzyme Linked Biotechnology Co., Ltd.). 2.14 Statistical analysis Unless otherwise specified, all experiments are conducted in no less than three independent parallel studies, and each parallel study consists of three repeated measurements. All data were analyzed using IBM's SPSS Statistics 21.0 software and expressed as mean ± standard error. One-way analysis of variance (ANOVA) is used for statistical analysis of comparisons between different groups. For all tests, a P-value < 0.05 was considered statistically significant. Unless marked with an asterisk, differences between groups can be considered statistically insignificant. 3. Results 3.1 SJZD regulates Serum lipid profile and Aortic lesion areas in ApoE −/− mice. In order to observe the anti-atherosclerotic effect of SIZD, we first measured the level of blood lipids.In the MG, the serum levels of TC and LDL-C increased significantly than those in the NG( P < 0.01, Fig. 1 A), while those of HDL-C decreased significantly( P < 0.01).The ApoE −/− mice in SJZD exhibited significantly lower serum levels of TC and LDL-C than the ApoE −/− mice in the model group ( P < 0.01). Then we carried out HE and Oil Red O staining.HE staining results (Fig. 1 B) showed that the endothelium of NG arranged regularly. No obvious atherosclerotic lesion was observed. In the MG, a greater amount of the number of aortic atherosclerotic plaque was observed. However, SJZD treatment significantly reduced the lipid deposition in the aorta of ApoE −/− mice and the area of atherosclerotic plaque was smaller. Furthermore, Oil Red O staining (Fig. 1 C). In the MG, atherosclerotic plaque development was observed, with the intimal thickening of the aorta, and an increased amount of massive Oil Red O-stained lipid. The amount of Oil Red O-stained lipids in the SJZD decreased significantly. 3.2Analysis of DEGs altered after treatment by Transcriptomic sequencing A total of 631.24 million raw reads were obtained in the sequencing. After filtering, 619.38 (98.12%) clean reads were obtained for following analysis. 600.13 million (95.07%) clean reads mapped to reference genome and 19.25 million clean reads were uniquely mapped. In order to further explore the mechanism of SJZD anti-AS, we performed transcriptome sequencing of mouse aorta samples from NG, MG and SJZD. Consequently,in order to identify the potential genes associated with AS, we analyzed DEGS. As shown in Fig. 2 A, a heat map distribution of those DEGs. The DEGs shown in red indicate a high expression, whereas those in blue indicate a low expression. The heat map demonstrated that several regions with highly altered (up-regulated or down-regulated) expression were shown in SJZD, NG, and MG groups. DEGs were further filtered by applying volcano plots, where the x-axis shows the Log2FC value and the y-axis corresponds to the value of log10 (p-value). The red dots correspond to significantly up-regulated DEGs ( P 1); the green dots correspond to significantly downregulated DEGs ( P 1); and the black dots indicate DEGs that were no statistical significance ( P > 0.05; |Log2FC|<1). There were 2754 DEGs in the MG/NG comparison (total gene counts: DEG-up = 2162; DEG-down = 592; Fig. 2 B). Moreover, in the SJZD/AS comparison, there were 1680 DEGs (total gene counts: DEG-up = 540, DEG-down = 1140; Fig. 2 C). 3.3 PPI network analysis of common SJZT targets to determine the core targets and bioinformatics analysis of predictive targets. As indicated in the Venn diagrams, in the mice treatment by SJZD, 33 genes downregulated by AS and upregulated by SJZD, and 306 genes upregulated by AS and downregulated by SJZD. Fortunately, 339 genes were rescued, indicating that SJZD reversed some AS-induced genes. Those 339 overlapping genes that indicate the opposite trends were marked as latent genes of interest (Fig. 3 A-B). The PPI network of the 339 latent targets was established by Cytoscape 3.7.2. As shown in Fig. 3 C, some of the remarkable genes in the PPI network included PON1, EGF, SLPI were initially upregulated in AS, but later afterwards reversed by SJZD. The GO enrichment results showed that 339 latent targets of SJZD were involved in biological processes, including lipid localization, eukocyte cell-cell adhesion, and leukocyte proliferation. The cellular component involved, including the membrane raft, nuclear chromatin, and nuclear transcription factor complex. It was closely related to molecular function, such as receptor ligand activity and cytokine receptor binding. As shown in Fig. 3 D-F, the GO enrichment results were scalped. The results of KEGG pathways enrichment showed that 339 DEGs were enriched into nine pathways, in which the pathways with P < 0.05 were defined as the key pathways (Fig. 3 G). Among them, Including AGE-RAGE signaling pathway in diabetic complications, Th17 cell differentiation, and Non-alcoholic fatty liver disease are AS-related signal pathways. 3.4 Determination of effective components and common anti-AS indexes of SJZD. A total of 136 active ingredients of four herbal medicines in SJZ decoction were identified based on threshold values of OB > 30% and DL > 0.18, including 22 constituents in ginseng, 15 constituents in poria cocos, 7 constituents in atractylodes, 92 constituents in licorice after deleting the duplicate data. Details of the active ingredients are listed in (supplement Table S1). In addition, 258 targets of active constituents of SJZ decoction were obtained from the TCMSP database and the gene names of these targets were collected via the Uniprot database (Supplementary Table S2). Furthermore, we had collected AS-related genes/targets from Gene Cards database, CTD database, TTD database. And we used TCMSP to supplement them and then added targets. Finally, we totally collected 3729 AS-related targets. Information on these targets is provided in Supplementary Table S3. After mapping SJZD targets and AS targets, a total of 171 common targets were obtained (Fig. 4 A). To better understand the interactions between common targets, PPI network was established to better interpret the mechanisms of SJZD in AS treatment by using STRING software (combined score ≥ 0.9). As shown in Fig. 4 B, the network consisted of 165 nodes and 1283 edges with an average node degree of 7.776, network centralization of 0.342, and an average number of neighbors of 15.552. According to the degree principle of each target, AKT1, IL6, MAPK1, JUN, and MAPK8 were determined as the hub targets (Fig. 4 C). To reflect the interactions intuitively between the active compounds of SJZD and their potential targets, the C–T network was constructed by mapping 114 active compounds to their 208 corresponding potential targets. As shown in Fig. 4 D, the network consisted of 322 nodes (114 active compound nodes, and 208 compound associated target nodes) and 1400 interaction edges.In the network, quercetin (degree 136), kaempferol (degree 85), and 7-Methoxy-2-methyl isoflavone (degree 31) presented the maximum interactions with potential targets, indicating that these active compounds with high degree values could play an significant role in the latent pharmacological effects of SJZD. The C-T network revealed intimate communications between active constituents and related targets, which provided a reference to further investigate the pharmacological mechanisms of SJZD. 3.5 Gene ontology enrichment analysis and targetpathway network construction The GO functional analyses is carried out in the David database to reveal the function of 171overlapping targets. GO enrichment analysis included 197biological processes (BP), 59 cellular components (CC) and 126 molecular functions (MF) with a threshold value of P < 0.05. The GO analysis results are shown in Fig. 5 A–C. The BP results mainly comprised of positive regulation of transcription from RNA polymerase II promoter, positive regulation of transcription, DNA-templated, response to drug. The CC analysis indicated that the overlapping targets were mainly related to cytosol, nucleoplasm, extracellular space. The MF results included protein binding, enzyme binding, identical protein binding. The results of the KEGG enrichment analysis indicated that 124 pathways meeting the threshold value of P < 0.05 were significantly enriched. As depicted in Fig. 5 D, the KEGG pathways of SJZD against atherosclerosis were mainly related to P13K-Akt signaling pathway, TNF signaling pathway, MAPK signaling pathway, HIF-1 signaling pathway, and FpxO signaling pathway. 3.6 Structuration of ingredient-gene-pathway (IGP) network and verification of core genes by qRT-PCR. In order to find a reliable signal pathway for SZJD to improve AS, we chose DEGs and GO-biological process pathway, in which the pathways with P < 0.05, by intersecting results of network pharmacology with transcriptome (Fig. 6 A). The results of GO-biological processes pathway revealed that 339 DEGs were enriched including lipid localization. Ingredients, genes, and pathways were also introduced into Cytoscape to establish a visualized Ingredient-gene-pathway (IGP) network for further analysis (Fig. 6 B). We performed qRT-PCR analyses to investigate whether SJZD affected the transcriptional regulation of the genes screened above. As shown in Fig. 2 , EGF, OLR1, COL1, SLP, SPP1, PLAU and MMP2 mRNA were significantly increased in ApoE −/− mice, while PON-1mRNA significantly reduced. compared with those in normal mice ( P 0.05).The expressions of PON-1and LDLR mRNA was significantly increased and the expressions of EGF, OLR1, COL1, SPP1, PLAU and MMP2 mRNA was significantly reduced under SJZD treatment ( P < 0.05 or ( P 0.05, Fig. 6 C). 3.7SJZD regulates OLR1/SPP1/EGF pathway to affect endothelial injury and prevent atherosclerosis Among the genes screened above, we found that there is a link between OLR1,SPP1 and EGF, which is related to endothelial injury and oxidative stress. To examine whether antioxidative effect of SJZD is mediated by OLR1,SPP1 and EGF, we detected their protein expression in mice aorta. The protein levels of OLR1,SPP1 and EGF were significantly increased in the model group, compared to the normal mice ( P < 0.01). SJZD intervention decreased the expression levels of these proteins( P < 0.01, Figs. 7 A). The above results suggest that SJZD can improve vascular endothelial injury in atherosclerotic mice by down-regulating the expression of OLR1,SPP1 and EGF protein. In order to further explore whether SJZD can improve endothelial function and reduce oxidative stress in ApoE −/− mice, the expression of OLR1, SPP1 and EGF were detected by immunehistochemistry. As shown in Figs. 7 B, immunohistochemical(IHC)analysis of OLR1,SPP1 and EGF showed that the area occupied by these proteins in the lesions of ApoE −/− mice was significantly higher than that of normal mice ( OLR1 was30.21vs. 18.39%, SPP1 was 25.40 vs.10.16%, EGF was 29.13 vs.11.80%, P < 0.01) However, after SJZD administration, this situation was obviously reversed, and plaque deposition and these proteins occupied area of aorta decreased significantly ( P < 0.01 or P < 0.05,Fig. 7 C). These results were consistent with the changes in protein levels in mouse aortas, as revealed by Western Blot. All in all, these results indicate that SJZD can reduce the level of OLR1, SPP1 and EGF in atherosclerotic lesions and improve oxidative stress. According to the above results, we speculate that the therapeutic effect of Sijunzi decoction on ApoE −/− mice may be related to oxidative stress. In order to further verify the effect of SJZD treatment on oxidative stress in ApoE −/− mice, we detected the levels of SOD and MDA in serum of mice. compared to that from normal mice, the serum levels of anti-oxidative SOD was significantly reduced and pro-oxidative MDA was increased in ApoE −/− mice ( P < 0.01). However, the serum level of SOD increased and the level of MDA decreased in ApoE −/− mice treated with SJZD ( P < 0.01, Figs. 7 D). These results indicated that SJZD treatment results in a marked suppression in oxidative stress in ApoE −/− mice through OLR1/SPP1/EGF pathway. 3.8 Validation of compound-target interaction To improve the accuracy of the connection between small molecular compounds and the core target proteins, we used molecular docking to evaluate the interactions between compounds and target proteins. It is generally believed that when the docking binding energy less than − 4.0 kcal/mol, it could be considered to have potential binding activity between the ligand and receptor [ 13 ] . A binding energy score less than − 5 kcal/mol suggests strong binding between the ligand and receptor [ 14 , 15 ] . Molecular docking indicated that kaempferol (A), quercetin (B) and naringenin (C) had good affinity for the three important compounds, SPP1, EGF, OLR1, in SJZT. 4. Discussion Many studies have shown that SJZD is a promising formula for treatment of cardiovascular related diseases. Preliminary basic research shows that SJZD has the effect of anti-atherosclerosis, however, the pharmacological mechanism of SJZD against AS is not yet fully understood, and the effects of its bioactive compounds and targets remain elusive. Therefore, by applying network pharmacology and transcriptomic methods, this study, which is combined with in vivo experiment, was carried out to identify the underlying mechanism and therapeutic targets of SJZD against AS effect. Firstly,the results about serum lipids and atherosclerotic plaques of aorta by HE and Oil Red O suggested that ApoE −/− AS mice model is conducted successfully and SJZT has a regulatory effect on them,which is consistent with previous studies.Then, transcriptomics sequencing was used to obtain active targets of AS.Through transcriptomics results of aorta tissue,339 genes were screened out,which indicated that SJZD may reverse some AS-induced genes.The PPI network analysis showed that some of the notable genes were reversed by SJZD,which means these genes may serve as potential targets and play an important role in the development of AS.The data of GO and KEGG showed that the importance of AS related terms was extremely high, which supports our previous hypothesis that the reliability of SJZD in the treatment of AS.Furthermore,a network pharmacology approach was used to identify 114 active compounds and 208 hub targets.The component-target network indicated that SJZD compounds might affect multiple targets and possess overlapping some,which results in synergistic effects.GO enrichment analysis showed that 208 targets were enriched in some biological processes,such as inflammatory response,angiogenesis and anti-apoptosis.KEGG analysis demonstrated that through the inhibition of some signal pathways,such as PI3K-Akt,TNF, MAPK and HIF-1, SJZD could resist AS. Finally, to make further analysis, we coupled network pharmacology with known transcriptomic results, and a visualized IGP network was established. According to the results,4 biological process (lipid localization, specific granule, positive regulation of lipid transport, positive regulation of lipid localization), 7 targets (SPP1, EGF, OLR1, LDLR, PON1, SLP1, PLAU), and 6 active ingredients (including beta-sitosterol, kaempferol, quercetin, stigmasterol, naringenin and isorhamnetin) play an extremely important role in the treatment of AS with SJZD. From above, we find that four biological processes involved in AS are mainly related to lipids and specific granules. As it is known, AS is a process caused by multiple factors, and lipid accumulation and subsequent inflammation have been shown to play a critical role. During further progression of atherosclerotic lesion, many other specific granules, such as low-density lipoprotein (LDL), come into play, which results in, locally or systemically, chronic low-grade inflammatory state. The infiltration and retention of atheroprone lipoproteins, predominately LDL, in the vessel wall, is considered as an initial step in the development of AS [ 16 ] . The transportation of atheroprone particles, like LDL, across the endothelium, is critical to the initiation, progression, and regression of AS [ 17 ] . In this study, we attempted to explore a prospective novel therapeutic for AS based on the bioinformatic methods. In vivo and some animal models have shown that quercetin has a wide range of biological actions, including anti-carcinogenic, antiviral,anti-inflammatory,antioxidant,and psychostimulant activities,as well as the ability to inhibit lipid peroxidation,capillary permeability and platelet aggregation, and to stimulate mitochondrial biogenesis [ 21 ] .Recently,the mechanism of the anti-atherosclerotic effects of quercetin have been investigated extensively [ 22 ] .In addition,quercetin suppressed vascular ROS formation endothelial dysfunction in HFD-fed ApoE-/- mice,with beneficial effects on atherosclerotic plaque formation. Naringenin is a major flavanone in citrus fruits that has multiple pharmacological attributes, such as anticancer and antiatherogenic antioxidant activity and free radical scavenging effects.Naringenin plays an important role in resisting cardiovascular diseases [ 23 ] .Naringenin attenuates endothelial dysfunction by reducing ROS accumulation and increasing nitric oxide production in endothelial cells [ 24 ] .Finally, we performed molecular docking to identify the reliability of the interactions of 3 active compounds and 3 core targets based on the previous study.We found that quercetin, kaempferol and naringenin bind most closely to the EGF, SPP1 and OLR1 proteins. In terms of secondary structure,quercetin and kaempferol are similar in structure.This shows that the binding pockets of kaempferol and quercetin are almost the same on protein,but the binding free energy and binding site are not necessarily the same. These results, which are consistent with our previous findings, show that the strategy of applying network pharmacology to find potential active compounds is reliable and feasible. 5. Conclusions In conclusion,benefiting by the integration of network pharmacology, transcriptomics,further experimental verification and validation of compound-target interaction,our study clarified that LOX-1/SPP1/EGF pathway may be one of the possible targets of SJZD to antiatherogenic.Kaempferol,quercetin and naringenin are the effective constitutes of SJZD to antiatherogenic.These findings are expected to guide the therapeutic research of AS in the future,and to provide a methodological reference for the pharmacodynamic material basis and targets of TCM in treatment of diseases by the integrated analysis combined with network pharmacology and transcriptomics.Although we have examined the possible role of SJZD in treating AS,there are some limitations in the present study.Further in vivo experimental validation is needed to support our research. Declarations Compliance with ethical standards All animal experiments and methods were performed in accordance with the Health Guide for the Care and Use of Laboratory Animals published by the National Institute of Health. Financial support and sponsorship This study was supported by National Natural Science Foundation of China grant (No. 81874372 to Wenna Chen);the Natural Science Foundation of Liaoning Province (20180550767). Conflicts of interest There are no conflicts of interest. Author Contribution Weiyan Chen were responsible for conception and design. Tianmin Ji and Ying Yang was responsible for manuscript writing. Zhuo Zhao, Hao Gao, were responsible for collection and assembly of data. Lijiang Zhou and Ying Wang were responsible for data analysis and interpretation. All authors were responsible for manuscript writing. All authors were responsible for the final approval of the manuscript. References Zhang Y, Tang Y, Yan J (2022) LncRNA-XIST Promotes Proliferation and Migration in ox-LDL Stimulated Vascular Smooth Muscle Cells through miR-539-5p/SPP1 Axis. Oxid Med Cell Longev. Jan 4; 2022:9911982 Wang R, Wang M, Ye J, Sun G, Sun X (2021) Mechanism overview and target mining of atherosclerosis: Endothelial cell injury in atherosclerosis is regulated by glycolysis (Review). Int J Mol Med 47(1):65–76 Wang R, Wang M, Ye J, Sun G, Sun X (2021) Mechanism overview and target mining of atherosclerosis: Endothelial cell injury in atherosclerosis is regulated by glycolysis (Review). Int J Mol Med 47(1):65–76 Hou X, Yang S, Yin J (2019) Blocking the REDD1/TXNIP axis ameliorates LPS-induced vascular endothelial cell injury through repressing oxidative stress and apoptosis. Am J Physiol Cell Physiol 316(1):C104–C110 Meng B, Li Y, Ding Y, Xu X, Wang L, Guo B, Zhu B, Zhang J, Xiang L, Dong J, Liu M, Xiang L, Xiang G (2021) Myeloid-derived growth factor inhibits inflammation and alleviates endothelial injury and atherosclerosis in mice. Sci Adv 7(21):eabe6903. 10.1126/sciadv.abe6903 Wang C, Niimi M, Watanabe T, Wang Y, Liang J, Fan J (2018) Treatment of atherosclerosis by traditional Chinese medicine: Questions and quandaries. Atherosclerosis 277:136–144 Huang C, Zhu Z, Cao X, Chen X, Fu Y, Chen Z, Li L, Song X, Jia R, Yin Z, Ye G, Feng B, Zou Y (2017) A Pectic Polysaccharide from Sijunzi Decoction Promotes the Antioxidant Defenses of SW480 Cells. Molecules 22(8):1341. 10.3390/molecules22081341 Ren Xiuyu K, Ziwei Z, Yanchun Q, Chaoyi L, Yafu, Li Jingke (2021) The effect of Sijunzi decoction on the apoptosis of photoreceptor cells after retinal detachment in rats with spleen deficiency syndrome. China J Integr Traditional Chin Western Med 41(02):204–211 Zhang R, Zhu X, Bai H, Ning K (2019) Network Pharmacology Databases for Traditional Chinese Medicine: Review and Assessment. Front Pharmacol 10:123 Li AP, He SS, Zhang WN, Zhang LC, Liu YT, Li K, Qin XM (2020) Exploration the active compounds of Astragali Radix in treatment of adriamycin nephropathy by network pharmacology combined with transcriptomic approach. J Ethnopharmacol 258:112537. 10.1016/j Dai S, Wang H, Wang M, Zhang Y, Zhang Z, Lin Z (2021) Comparative transcriptomics and network pharmacology analysis to identify the potential mechanism of celastrol against osteoarthritis. Clin Rheumatol 40(10):4259–4268 He S, Li A, Zhang W et al (2020) An integrated transcriptomics and network pharmacology approach to exploring the mechanism of adriamycin-induced kidney injury. Chem Biol Interact 325:109096 Li J, Luo H, Liu X, Zhang J, Zhou W, Guo S, Chen X, Liu Y, Jia S, Wang H, Li B, Cheng G, Wu J (2020) Dissecting the mechanism of Yuzhi Zhixue granule on ovulatory dysfunctional uterine bleeding by network pharmacology and molecular docking. Chin Med 15:113. 10.1186/s13020-020-00392-0 Guo C, Kang X, Cao F, Yang J, Xu Y, Liu X, Li Y, Ma X, Fu X (2021) Network Pharmacology and Molecular Docking on the Molecular Mechanism of Luo-hua-zi-zhu (LHZZ) Granule in the Prevention and Treatment of Bowel Precancerous Lesions. Front Pharmacol 12:629021. 10.3389/fphar Liu J, Liu M, Yao Y, Wang J, Li Y, Li G, Wang Y (2012) Identification of novel potential β-N-acetyl-D-hexosaminidase inhibitors by virtual screening, molecular dynamics simulation and MM-PBSA calculations. Int J Mol Sci 13(4):4545–4563. 10.3390/ijms13044545 Tabas I, Williams KJ, Boren J (2007) Subendothelial lipoprotein retention as the initiating process in atherosclerosis: update and therapeutic implications. Circulation 116:1832–1844 Zhang X, Fernández-Hernando C (2020) Transport of LDLs into the arterial wall: impact in atherosclerosis. Curr Opin Lipidol 31(5):279–285 Li Y, Zhao L, Qi W (2022) Uric acid, as a double-edged sword, affects the activity of epidermal growth factor (EGF) on human umbilical vein endothelial cells by regulating aging process. Bioengineered 13(2):3877–3895 Barreto J, Karathanasis SK, Remaley A, Sposito AC (2021) Role of LOX-1 (Lectin-Like Oxidized Low-Density Lipoprotein Receptor 1) as a Cardiovascular Risk Predictor: Mechanistic Insight and Potential Clinical Use. Arterioscler Thromb Vasc Biol 41(1):153–166 Tian K, Ogura S, Little PJ, Xu SW, Sawamura T (2019) Targeting LOX-1 in atherosclerosis and vasculopathy: current knowledge and future perspectives. Ann N Y Acad Sci 1443(1):34–53 Li Y, Yao J, Han C, Yang J, Chaudhry MT, Wang S, Liu H, Yin Y, Quercetin (2016) Inflamm Immun Nutrients 8(3):167. 10.3390/nu8030167 Luo M, Tian R, Lu N (2020) Quercetin Inhibited Endothelial Dysfunction and Atherosclerosis in Apolipoprotein E-Deficient Mice: Critical Roles for NADPH Oxidase and Heme Oxygenase-1. J Agric Food Chem 68(39):10875–10883 Orhan IE, Nabavi SF, Daglia M, Tenore GC, Mansouri K, Nabavi SM (2015) Naringenin and atherosclerosis: a review of literature. Curr Pharm Biotechnol 16(3):245–251 Heidary Moghaddam R, Samimi Z, Moradi SZ, Little PJ, Xu S, Farzaei MH (2020) Naringenin and naringin in cardiovascular disease prevention: A preclinical review. Eur J Pharmacol 887:173535 Supplementary Tables Supplementary Tables S1 to S3 are not available with this version. Additional Declarations No competing interests reported. Supplementary Files GA.png Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4872904","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":347971524,"identity":"2cf026d3-1ca1-413f-b3ba-f6bd9bf26d06","order_by":0,"name":"Weiyan Chen","email":"","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine，Liaoning","correspondingAuthor":false,"prefix":"","firstName":"Weiyan","middleName":"","lastName":"Chen","suffix":""},{"id":347971525,"identity":"d6b358f3-2340-48f7-8b59-3e435412b6e8","order_by":1,"name":"Tianmin Ji","email":"","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine，Liaoning","correspondingAuthor":false,"prefix":"","firstName":"Tianmin","middleName":"","lastName":"Ji","suffix":""},{"id":347971526,"identity":"3c917777-8ba1-4432-ac3f-92ec7d168e20","order_by":2,"name":"Ying Yang","email":"","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine，Liaoning","correspondingAuthor":false,"prefix":"","firstName":"Ying","middleName":"","lastName":"Yang","suffix":""},{"id":347971527,"identity":"25b83c9f-27de-4f54-b3b7-c1791b6847a5","order_by":3,"name":"Zhuo Zhao","email":"","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine，Liaoning","correspondingAuthor":false,"prefix":"","firstName":"Zhuo","middleName":"","lastName":"Zhao","suffix":""},{"id":347971528,"identity":"893d0981-99d7-452c-b833-53fc490a457f","order_by":4,"name":"Hao Gao","email":"","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine，Liaoning","correspondingAuthor":false,"prefix":"","firstName":"Hao","middleName":"","lastName":"Gao","suffix":""},{"id":347971529,"identity":"77b26bb7-d2f9-43de-86fc-d22e54c8049d","order_by":5,"name":"Lijiang Zhou","email":"","orcid":"","institution":"Affiliated Hospital of Liaoning University of Traditional Chinese Medicine","correspondingAuthor":false,"prefix":"","firstName":"Lijiang","middleName":"","lastName":"Zhou","suffix":""},{"id":347971530,"identity":"b403620a-8cfb-4df2-97ba-3719742ef245","order_by":6,"name":"Ying Wang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA+ElEQVRIiWNgGAWjYBACPhiDX/7xwQcfKiTk5AlpYYMxJBvSkg1nnLEwNmwgVovBgRwzad62ikSGA4S0sPcefs1TcSex4cCxZAPeeRIJjA3MDx/dwKeF51yaNc+ZZ4mNjc0HH0huk8hjZ2AzNs7Bp0Uix8yYt+1wYjMzW7KB4TaJYsYGHjZpvFrk3wC1/Duc2MbGYyaROEcC6EJCWiR4jB/zNhxO7OEBajnYQIwWnhwzxjnHnhnPkGBLNmw4JmFs2EzAL/zsZ4w/vKm5I7v/BvPBx39q6uTk2ZsfPsanBWSRFA9KXDDjVw5W8vEHwegbBaNgFIyCEQ0Ao55NwOTg7foAAAAASUVORK5CYII=","orcid":"","institution":"Liaoning University of Traditional Chinese Medicine","correspondingAuthor":true,"prefix":"","firstName":"Ying","middleName":"","lastName":"Wang","suffix":""}],"badges":[],"createdAt":"2024-08-07 07:56:24","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4872904/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4872904/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":63905008,"identity":"ad83d425-10b4-433e-bfa2-3b0f4ef4ee46","added_by":"auto","created_at":"2024-09-03 15:03:34","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":340578,"visible":true,"origin":"","legend":"\u003cp\u003eSJZD regulates Serum lipid profile and Aortic lesion areas in ApoE-/-mice.\u0026nbsp; (A) Detection of TC, TG and LDL-C, HDL-C in serum of ApoE-/-mice. The results showed that SJZD could significantly down-regulate the level of serum TC,TG,LDL-C and up-regulate the level of HDL-C in ApoE\u003csup\u003e-/-\u003c/sup\u003e mice, (\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01, n=6);(B,C) Representative histology of H\u0026amp;E and Oil Red O staining in aortic Aorta. Original magnification,200um.\u0026nbsp; In MG group, the inner wall of plaque thickened and the infiltration of plaque and inflammatory cells increased. The plaque area of SJZD is lower than that in MG. SJZD treatment significantly reduced the lipid deposition in the aorta of ApoE\u003csup\u003e−/− \u003c/sup\u003emice and the area of atherosclerotic plaque was smaller. Data are presented as the means ± standard deviation. NG：normal group; MG; model group ;SJZD: SJZ Decoction group；MG was compared to NG ,\u003csup\u003e*\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e**\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01;SJZD was compared with MG， \u003csup\u003e▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e▲▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01. TC, total cholesterol; TG, triglyceride; LDL-C, low-density-lipoprotein-cholesterol; HCL-C, high-density- lipoprotein-cholesterol. The levels of TC TG, LDL-C and HDL-C are shown in mmol/L.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/7d39b73110e3413254066ae7.png"},{"id":63905582,"identity":"4f7edfc4-f968-43e3-9f5e-068a1aa9f4b2","added_by":"auto","created_at":"2024-09-03 15:11:34","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":242727,"visible":true,"origin":"","legend":"\u003cp\u003eDEGs were analyzed to determine latent genes associated with AS. (A) heat map distribution of those DEGs. The heat map demonstrated that several regions with highly altered (up-regulated or down-regulated) expression were shown in SJZD, NG, and MG groups. Red represents high expression and blue represents low expression. (B) DEGs in the MG/NG comparison. There were 2754 DEGs in the MG/NG comparison, total gene counts: DEG-up = 2162; DEG-down=592; (C) DEGs in the SJZD/AS comparison. there were 1680 DEGs (total gene counts: DEG-up = 540, DEG-down = 1140\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/ba18e287e0edfbacbf64cb33.png"},{"id":63906119,"identity":"aada49f5-3ff0-4c15-a544-cb086f80a8a2","added_by":"auto","created_at":"2024-09-03 15:19:34","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":531161,"visible":true,"origin":"","legend":"\u003cp\u003ebioinformatics analysis of the core targets. (A,B) Venn diagram of potential genes of interest. Venn diagram showed that 339 genes affected by AS could be saved by SJZD treatment, of which 33 genes were up-regulated and 306 genes were down-regulated. (C) PPI network analysis of SJZD anti-AS related protein from STRING database. The results showed that PON1, EGF and SLPI genes were initially up-regulated in AS, but were subsequently reversed by SJZD. (D-G) Enrichment analysis of core targets. 339 DEGs were enriched into nine pathways, in which the key pathway is \u003cem\u003eP\u003c/em\u003e\u0026lt;0.05. Including AGE-RAGE signaling pathway in diabetic complications, Th17 cell differentiation, and Non-alcoholic fatty liver disease are AS-related signal pathways.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/6fd1fbcc6db4ff10bb6eaf2b.png"},{"id":63905009,"identity":"277c8ed2-e1b0-4e66-b517-cf0e477fbb55","added_by":"auto","created_at":"2024-09-03 15:03:34","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":525732,"visible":true,"origin":"","legend":"\u003cp\u003eDetermination of effective components and common anti-AS indexes of SJZD. (A) mapping SJZD targets and AS targets. According to the threshold of OB \u0026gt; 30% and DL \u0026gt; 0.18, 136 active components of 4 herbs in Sijunzi decoction were identified, 258 active compounds of SJZD were obtained, and the gene names of these targets were collected to predict the related targets of AS. The result was a total of 171 common goals of SJZD and AS. (B) PPI network of SJZD anti-AS related protein from STRING database. The network consisted of 165 nodes and 1283 edges with an average node degree of 7.776, network centralization of 0.342, and an average number of neighbors of 15.552. (C) A shared target network between SJZD potential targets and AS targets. AKT1, IL6, MAPK1, JUN, and MAPK8 were determined as the hub targets. (D) The C-T network was constructed to reflect the interaction between the active constituents of SJZD and its potential targets. the network consisted of 322 nodes (114 active compound nodes, and 208 compound associated target nodes) and 1400 interaction edges. In the network, quercetin (degree 136), kaempferol (degree 85), and 7-Methoxy-2-methyl isoflavone (degree 31) presented the maximum interactions with potential targets, indicating that these active compounds with high degree values could play an significant role in the latent pharmacological effects of SJZD.\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/d9f47a27554bd60b9fe39078.png"},{"id":63905012,"identity":"7eb78d94-b250-4a53-9b8c-50da86b1d958","added_by":"auto","created_at":"2024-09-03 15:03:34","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":334863,"visible":true,"origin":"","legend":"\u003cp\u003ePathway enrichment analysis. (A-C) The distribution of GO entries in biological process, molecular function and cell composition are shown. The BP results mainly comprised of positive regulation of transcription from RNA polymerase II promoter, positive regulation of transcription, DNA-templated, response to drug. The CC analysis indicated that the overlapping targets were mainly related to cytosol, nucleoplasm, extracellular space. The MF results mainly included protein binding, enzyme binding, identical protein binding. (D) KEGG enrichment analysis for key pathways. The color scale indicates the different thresholds for the p-values, and the size of the dots represents the number of genes corresponding to each term. The results of the KEGG enrichment analysis showed that 124 pathways meeting the threshold value of \u003cem\u003eP\u003c/em\u003e\u0026lt;0.05 were significantly enriched. The KEGG pathways of SJZD against atherosclerosis were mainly related to P13K-Akt signaling pathway, TNF signaling pathway, MAPK signaling pathway, HIF-1 signaling pathway, and FoxO signaling pathway.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/2deca9f70b2849d858d6761c.png"},{"id":63905016,"identity":"0af40972-846a-49a4-8220-a077e04231af","added_by":"auto","created_at":"2024-09-03 15:03:35","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":356549,"visible":true,"origin":"","legend":"\u003cp\u003eStructuration of ingredient-gene-pathway (IGP) network\u0026nbsp;and verification of core genes by qRT-PCR.(A) Screening of anti-AS signal pathways by crossing DEGS and GO biological process pathways. We found a reliable anti-AS pathway of SJZD. (B) Ingredient-gene-pathway (IGP) network for analysis. We screened for genes such as EGF, OLR1, COL1, SLP, SPP1, PLAU and MMP2. (C)Aorta gene expression of EGF, OLR1, COL1, PON-1, SLP1, SPP1, PLAU, MMP2 and LDLR mRNA were detected by RT-qPCR method. In addition to SLP1( \u003cem\u003eP\u003c/em\u003e＞0.05,n=3), SJZD can significantly affect the transcriptional regulation of the above screened genes(\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01 or \u003cem\u003eP\u003c/em\u003e\u0026lt;0.05,n=3). data are presented as the means±standard deviation. NG：normal group; MG; model group; SJZD: SJZ Decoction group; MG was compared to NG, \u003csup\u003e*\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e**\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01; SJZD was compared with MG，\u003csup\u003e▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e▲▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/0df7a8c4a7beb0967140e10c.png"},{"id":63906120,"identity":"18b16a87-154b-459b-aab7-a6ac67036905","added_by":"auto","created_at":"2024-09-03 15:19:34","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":278242,"visible":true,"origin":"","legend":"\u003cp\u003eSJZD alleviates oxidative stress and the expression of AS-related genes and proteins. (A) Protein was extracted from mice aorta. Then the levels of OLR1, SPP1 and EGF were determined by western blotting assay. (B) Quantification and statistical analysis of western blotting results. SJZD intervention decreased the expression levels of these proteins (\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01, n =3). (C) Immunohistochemical staining on sections of the thoracic artery with antibodies against OLR1, SPP1 and EGF (100um). (D) Quantitation of OLR1, SPP1 and EGF positive areas is expressed as a percentage of the total lesion area. SJZD could significantly reverse the deposition of aortic plaques and the aortic area occupied by these proteins in ApoE\u003csup\u003e\u003cstrong\u003e−/−\u003c/strong\u003e\u003c/sup\u003e\u003cstrong\u003e \u003c/strong\u003emice (\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01 or \u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, n=3). (E) Production of MDA and SOD in serum were measured by ELISA kit. After treatment with SJZD, the level of serum SOD increased and the level of MDA decreased in ApoE-/-mice (\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01, n=8). These results suggest that SJZD treatment ameliorates endothelial damage and inhibits oxidative stress in ApoE\u003csup\u003e-/-\u003c/sup\u003e mice through OLR1/SPP1/EGF pathway. The quantitative results were depicted. data are presented as the means±standard deviation.NG: normal group; MG: model group; SJZD:SJZ Decoction group; MG was compared to NG ,\u003csup\u003e*\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e**\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01;SJZD was compared with MG，\u003csup\u003e▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.05, \u003csup\u003e▲▲\u003c/sup\u003e\u003cem\u003eP\u003c/em\u003e\u0026lt;0.01.The levels of MDA and SOD are shown in nmol/mL\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/ea44eadb85f11b6b87929dc2.png"},{"id":63905581,"identity":"2d58e397-e64f-4b21-a30c-1deb3188302d","added_by":"auto","created_at":"2024-09-03 15:11:34","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":262427,"visible":true,"origin":"","legend":"\u003cp\u003eMolecular docking between small molecule ligands and core targets protein. (A) Quercetin with EGF, (B) Quercetin with OLR1, (C) Quercetin with SPP1, (D) Naringenin with OLR1, (E) Naringenin with SPP1. The results of molecular docking showed that kaempferol, quercetin and Naringenin had potential binding activity to three important compounds SPP1, EGF and OLR1 in SJZD.\u003c/p\u003e","description":"","filename":"8.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/14d058096bbda56558b617fe.png"},{"id":80051410,"identity":"83667d70-35b3-4057-b036-c93a290d9d77","added_by":"auto","created_at":"2025-04-07 10:23:51","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4220239,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/786a4122-49da-42b0-911b-eb7744c9ac58.pdf"},{"id":63905579,"identity":"7424e21f-2211-4670-bb3d-7909f22462f4","added_by":"auto","created_at":"2024-09-03 15:11:34","extension":"png","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":561385,"visible":true,"origin":"","legend":"","description":"","filename":"GA.png","url":"https://assets-eu.researchsquare.com/files/rs-4872904/v1/623cf87ab55450024728c5a3.png"}],"financialInterests":"No competing interests reported.","formattedTitle":"Elucidating the anti-atherosclerosis mechanism of Si Jun Zi Decoction by integrating network pharmacology and transcriptomic with experimental validation in vivo","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eAtherosclerosis (AS) is one of the most common vascular disfunction all over the world and could lead to stroke and coronary heart disease.It is associated with health risks that are major causes of death \u003csup\u003e[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]\u003c/sup\u003e.Previous studies have demonstrated that AS has a connection with endothelial cell damage,inflammation,lipid metabolism disorders,and immune dysfunction\u003csup\u003e[\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]\u003c/sup\u003e.Endothelial cell injury,which results in plaque, leads to further instability and rupture of it, and secondary thrombosis will accelerate the progression of the disease\u003csup\u003e[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]\u003c/sup\u003e.Vascular endothelial cell injury caused by oxidative stress is an significant factor and an early pathological change in the development of AS\u003csup\u003e[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e–\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]\u003c/sup\u003e.In China,Chinese herbal medicine (CHM) has been used to treat atherosclerotic complications for a long time,such as angina pectoris,myocardial infarction and stroke\u003csup\u003e[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]\u003c/sup\u003e.However,the molecular mechanism is still unclear.Sijunzi Decoction (SJZD),consisting of Atractylodes macrocephala Panax ginseng,Glycyrrhiza uralensis and Poria cocos,has been applied for thousands of years\u003csup\u003e[\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]\u003c/sup\u003e. SJZD as a classic prescription for treating spleen deficiency syndrome,has obvious anti-oxidative stress effect\u003csup\u003e[\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]\u003c/sup\u003e.However,the composition of traditional Chinese medicine (TCM) is complex,with many active components and function through multiple pathways.\u003c/p\u003e \u003cp\u003eTo further explore its mechanism, network pharmacology and transcriptomic techniques are essential to be introduced. Network pharmacology,which is similar to the overall concept of TCM,constructs a drug-component-target-disease regulatory network through joint analysis of disease-related genes\u003csup\u003e[\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e]\u003c/sup\u003e.It is helpful to elucidate the mechanism of drug comprehensively as much as possible and interpret the synergistic effect of multiple components and targets of TCM\u003csup\u003e[\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]\u003c/sup\u003e.It is also used to evaluate multi-target drugs and complex diseases\u003csup\u003e[\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]\u003c/sup\u003e.Transcriptomics is a subject that elaborates gene transcription and regulation in tissues and cells with a holistic view,and could comprehensively explain the regulation mechanism of gene expression in complex diseases by integrating a variety of calculation ways.For the study of gene transcription and regulation in cells or tissues\u003csup\u003e[\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]\u003c/sup\u003e, the combined application of network pharmacology and transcriptomics could preliminarily screen the active ingredients and molecular targets of drugs,and provide convenient conditions for subsequent research\u003csup\u003e[\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]\u003c/sup\u003e.Based on the above, the aim of our research is to find the possible mechanism and molecular targets of SJZD,which could improve vascular endothelial injury and prevent AS,mainly through network pharmacology and transcriptomic technology.\u003c/p\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003cdiv id=\"Sec13\" class=\"Section3\"\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section3\"\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section3\"\u003e \u003c/div\u003e \u003c/div\u003e "},{"header":"2. Materials and Methods","content":"\u003ch2\u003e2.1 Experimental animal\u003c/h2\u003e\u003cp\u003e16 SPF grade male ApoE \u003csup\u003e− / −\u003c/sup\u003e mice (body weight, 18–20 g) and 8 male C57BL/6 mice were provided b\u003cem\u003ey Beijing\u003c/em\u003e Vital River Laboratory \u003cem\u003eAnimal\u003c/em\u003e Technology Co.Ltd. The Committee on the Ethics of Animal Experiments at Liaoning University of Traditional Chinese Medicine approved this study(permitted number SCXK(Jing) 2012-0001).All mice were raised in the Animal Experimental Center of Liaoning University of Traditional Chinese Medicine and allowed free movement with food and water ad libitum.\u003c/p\u003e\u003ch2\u003e2.2 Establishment of AS model mice and treatments\u003c/h2\u003e\u003cp\u003eAfter 1 week of adaptive feeding in SPF grade mice,8 C57BL/6 mice were considered as a normal group(NG),and 16 ApoE\u003csup\u003e−/−\u003c/sup\u003emice were randomly assigned to the model group(MG) and SJZ Decoction group(SJZD).The normal group was fed with basal diet daily,while the other groups were fed with high-fat diet(containing 21% fat, 10% lard, 1% cholesterol, etc.)with free access to food and water for 12 weeks.The drug intervention was given to the group,which the ApoE\u003csup\u003e−/−\u003c/sup\u003e mice treated with 20 g/kg/d SJZ Decoction for the last 4 weeks.The normal group and model group were given an equal volume of normal saline daily by gavage once a day.\u003c/p\u003e\u003ch2\u003e2.3 Serum and tissue samples\u003c/h2\u003e\u003cp\u003eAfter the treatment period,all mice were fasted overnight and anesthetized with chloral hydrate.Eyeball blood was collected from sacrificed mice after resting for 30 min and subjected to centrifugation at 3000× g and 4°C for 15 min in order to obtain serum, and stored in aliquots at -80°C for further analysis.Quickly separate the aortic arch,thoracic-abdominal aorta,and rinse 3 times with ice-cold saline.The aorta was fixed in 4% paraformaldehyde for atherosclerotic lesion analysis.The remaining part is frozen at -80°C until it is used for biochemical determination,iron content analysis and RNA extraction.\u003c/p\u003e\u003ch2\u003e2.4 Serum lipids\u003c/h2\u003e\u003cp\u003eThe levels of serum TC,TG,HDL-C,and LDL-C were measured using an automated biochemical analyzer(SIEMENS,Munich,DE,Germany).\u003c/p\u003e\u003ch2\u003e2.5 Histopathology test\u003c/h2\u003e\u003cp\u003eThe arterial tissue was fixed in 4% paraformaldehyde solution overnight,and stained using the conventional method of hematoxylin-eosin(HE)staining.The specific steps are as follows:ⅰ.The tissue is dehydrated,transparent,embedded in paraffin, and then sectioned(thickness:~5 microns). Ⅱ.Paraffin sections are dewaxed with xylene, rehydrated with gradient alcohol,stained with hematoxylin,differentiated with hydrochloric acid alcohol,counterstained with eosin, then dehydrated,transparent with xylene,and mounted with neutral gum. Ⅲ.Observe the morphology of rat aorta under a microscope.\u003c/p\u003e\u003cp\u003eOil Red O staining of aorta: complete aorta was isolated and fixed in 10% neutral formalin for 24 h.The complete aorta was in 60% isopropanol for 10min,stained in oil red O (Solarbio,China) working fluid (3:2) for 3h,washed in 60% isopropanol for 6 times till the background color became white,and then photographed.\u003c/p\u003e\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e2.6 Transcriptomic sequencing and data analysis\u003c/h2\u003e \u003cdiv id=\"Sec8\" class=\"Section3\"\u003e \u003ch2\u003e2.6.1 RNA extraction and library preparation\u003c/h2\u003e \u003cp\u003e The aorta tissues samples were randomly pooled from the mice in the NG,MG,and SJZD (n = 3).The total RNA was extracted using 1mL of Trizol Reagent (Invitrogen,Carlsbad,CA,USA) according to the manufacturer’s instructions.The total RNA was reverse transcribed into cDNA.NEBNext® UltraTM RNA Library Preparation Kit for Illumina(NEB,USA) was used to generate sequencing libraries.HiSeq X instrument was used for whole transcriptome sequencing.Cut adapt was used to filter the raw data to get clean data for further analysis.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003ch2\u003e2.6.2 Differential expression genes (DEGs) analysis\u003c/h2\u003e \u003cp\u003eTwo comparisons-MG versus NG, and SJZD vs MG were executed application of the DESeq2 R package (1.16.1) to determine differential expression genes.DESeq2 provided statistical routines for defining various expression in digital gene expression data application model according to the negative binomial distribution.The significance of DESeq2 was determined by p-value \u0026lt; 0.05 and |log2(fold-change)|\u0026gt;1,which was implemented to define DEGs.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section3\"\u003e \u003ch2\u003e2.6.3 PPI network construction and analysis\u003c/h2\u003e \u003cp\u003eWe used the Venn map to scalp the DEGs,which were down-regulated by MG and up-regulated by SJZD, or DEGs up-regulated by MG and down-regulated by SJZD. The screened DEGs might be the latent targets of SJZD to treat AS.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section3\"\u003e \u003ch2\u003e2.6.4 Functional enrichment analysis\u003c/h2\u003e \u003cp\u003eThe GO and KEGG pathways analysis were enriched by the R software package to realize the functions of candidate DEGs.In our research,we used the cluster Profiler to analyze the GO and KEGG pathways which was focused on were presented.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e\u003ch2\u003e2.6 Transcriptomic sequencing and data analysis\u003c/h2\u003e\u003ch2\u003e2.6.1 RNA extraction and library preparation\u003c/h2\u003e\u003cp\u003e The aorta tissues samples were randomly pooled from the mice in the NG,MG,and SJZD (n = 3).The total RNA was extracted using 1mL of Trizol Reagent (Invitrogen,Carlsbad,CA,USA) according to the manufacturer’s instructions.The total RNA was reverse transcribed into cDNA.NEBNext® UltraTM RNA Library Preparation Kit for Illumina(NEB,USA) was used to generate sequencing libraries.HiSeq X instrument was used for whole transcriptome sequencing.Cut adapt was used to filter the raw data to get clean data for further analysis.\u003c/p\u003e\u003ch2\u003e2.6.2 Differential expression genes (DEGs) analysis\u003c/h2\u003e\u003cp\u003eTwo comparisons-MG versus NG, and SJZD vs MG were executed application of the DESeq2 R package (1.16.1) to determine differential expression genes.DESeq2 provided statistical routines for defining various expression in digital gene expression data application model according to the negative binomial distribution.The significance of DESeq2 was determined by p-value \u0026lt; 0.05 and |log2(fold-change)|\u0026gt;1,which was implemented to define DEGs.\u003c/p\u003e\u003ch2\u003e2.6.3 PPI network construction and analysis\u003c/h2\u003e\u003cp\u003eWe used the Venn map to scalp the DEGs,which were down-regulated by MG and up-regulated by SJZD, or DEGs up-regulated by MG and down-regulated by SJZD. The screened DEGs might be the latent targets of SJZD to treat AS.\u003c/p\u003e\u003ch2\u003e2.6.4 Functional enrichment analysis\u003c/h2\u003e\u003cp\u003eThe GO and KEGG pathways analysis were enriched by the R software package to realize the functions of candidate DEGs.In our research,we used the cluster Profiler to analyze the GO and KEGG pathways which was focused on were presented.\u003c/p\u003e\u003ch2\u003e2.7 Network pharmacology\u003c/h2\u003e\u003ch2\u003e2.7.1 Gatheration of active compounds and acquisition of corresponding targets of SJZ decoction\u003c/h2\u003e\u003cp\u003eThe active ingredients of four herbal medicines in SJZ decoction were obtained from TCMSP and the BATMAN-TCM database.According to the recommended drug screening criteria of the TCMSP database,chemical compounds with OB ≥ 30% and DL ≥ 0.18 may considered active,and they were choosed as effective constituents for further analysis through the TCMSP database.Targets of active ingredients were acquired from the TCMSP database.The target protein names corresponding to active compounds were normalized in UniProt.\u003c/p\u003e\u003ch2\u003e2.7.2 Collect the common targets of SJZ Decoction and AS\u003c/h2\u003e\u003cp\u003eUsing “atherosclerosis” as a keyword,the predicted genes of AS were collected from Gene Cards database, CTD database and TTD database.The search results found in three database were merged and deleted duplicate targets to obtain all target genes of AS.The common targets of SJZ Decoction and AS were identified by the bioinformatics online tool (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.bioinformatics.com.cn/).I\u003c/span\u003e\u003cspan address=\"http://www.bioinformatics.com.cn/).I\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003en short,the Venn map was used to draw the targets of SJZ Decoction and the disease targets of AS.Next,the core targets were considered as the potential targets for SJZ Decoction in treatment of AS.\u003c/p\u003e\u003ch2\u003e2.7.3 Network construction and Functional enrichment\u003c/h2\u003e\u003cp\u003eThe core targets were entered into String online tool to obtain the interrelated information of protein-protein interactions (PPI).Next, Cytoscape 3.7.2 was used to establish the visualized PPI network and calculate the degree centrality. Moreover, the key proteins of the PPI network were acquired. The GO and KEGG pathways of the potential proteins were enriched and analyzed by R 3.6.3. The results were analized by using p-value and count values.\u003c/p\u003e\u003ch2\u003e2.8 Structure ingredient-gene-pathway network (IGP)\u003c/h2\u003e\u003cp\u003eWe assessed pathways in the GO-biological process,which appears together in network pharmacology and transcriptome research.The genes related to this category were found in the transcriptome results,and the gene-pathway network was established.According to the network pharmacology results,the related active molecules of genes in this part were analyzed.Ingredients,genes,and GO-biological process were also put into Cytoscape tool to demonstrate a visualized Ingredient-gene-pathway (IGP) network.\u003c/p\u003e\u003ch2\u003e2.9 Validation of compound-target interaction\u003c/h2\u003e\u003cp\u003eTo validate the compound-target association,the molecular docking program was performed with Auto Dock Tools software(version 4.2).All the 3D structures of core target proteins were obtained from the PDB database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.rcsb.org/),s\u003c/span\u003e\u003cspan address=\"http://www.rcsb.org/),s\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003eo that proteins and ligands could be prepared in the Auto Dock before performing the docking progress.The whole molecular docking process included setup of proteins,determination of docking sites and docking of proteins to small molecules.Firstly, during the process of protein preparation,water molecules were removed from the protein structures as well as the region arround the amide portion and atomic groups.Next,the protein binding sites were defined and edited.Finally,insert the compounds into the protein-binding site.If the small molecule docked with the protein successfully,it showed a labelled interaction site and a corresponding docking score.Hence,target proteins and compounds with high docking scores were choosed for further analysis.\u003c/p\u003e\u003ch2\u003e2.10 Quantitative real-time PCR to examine mRNA expression\u003c/h2\u003e\u003cp\u003eTRIZOL reagent (Solarbio,China) was used to extract total RNA from the sample.The cDNA template was synthesized using a commercially available reverse transcription reaction (Bio-Rad,USA),and the quantitative RT-PCR experiment was performed using an ABI7500 rapid quantitative PCR instrument (Applied Biosystems, USA).The reaction procedure is:pre-denaturation at 95°C for 5 min,denaturation at 95°C for 20 s,annealing at 58°C for 30 s,extension at 72°C for 20 s,40 cycles.The relative quantification (2\u003csup\u003e−ΔΔCt\u003c/sup\u003e) method was used to calculate the relative amount of mRNA. The PCR primers were as follows:EGF,5’-ACAGAAGGAGTAGATACGCTTG-3’-(forward),and 5’-GATTATTCGATGATGCTTCCCG-3’-(reverse).OLR1,5’-GAAGCCTGCGAATGACGAGC-3’-(forward),5’-ACACCAGGCAGAGGATGACC-3’-(reverse).COL1,5’-TGAACGTGGTGTACAAGGTC-3’-(forward),and 5’-CCATCTTTACCAGGAGAACCAT-3’-(reverse).PON1,5’-AAGAGGAAAGATCACTCTTGCA-3’-(forward),and 5’-GGTCCAATAGCAGCTATATCGT-3’-(reverse).SLP1,5’-ATGTATGATGCTTAACCCTCCC-3’-(forward),and 5’-AGGCAGACTTTCCCACATATAC-3’-(reverse).SPP1,5’-AAACACACAGACTTGAGCATTC-3’-(forward),and 5’-TTAGGGTCTAGGACTAGCTTGT-3’-(reverse).PLAU,5’-GCTTGTTTCTCATGAACAGTGT-3’-(forward),and 5’-TTCGATGTTACAGATAAGCGGT-3’-(reverse).MMP2,5’-ACTTTGAGAAGGATGGCAAGTA-3’-(forward),and 5’-CTTCTTATCCCGGTCATAGTCC-3’-(reverse).LDLR1,5’-GAAGCCTGCGAATGACGAGC-3’-(forward),and 5’-ACACCAGGCAGAGGATGACC-3’-(reverse).And GAPDH,5’-GACATGCCGCCTGGAGAAAC-3’-(forward),and 5’-AGCCCAGGATGCCCTTTAGT-3’-(reverse).\u003c/p\u003e\u003ch2\u003e2.11 Western blot analysis\u003c/h2\u003e\u003cp\u003eTotal proteins were extracted from tissues using RIPA lysis buffer containing 1% phenylmethylsulphonyl fuoride (PMSF) for 30 minute on ice.The supernatants were collected after being clarified by centrifuging at 12000 rpm for 15 minutes at 4°C. Protein concentrations were measured by the BCA protein assay.The same amounts of protein (60 µg) were separated by electrophoresis on 10%-12% sodium dodecyl sulphate polyacrylamide gels and transferred onto PVDF membranes.Next,PVDF membranes were incubated overnight at 4°C with the antibody(proteintech,China). After washing with TBST 3 times,the membranes were incubated for 1 h with fluorescent secondary antibody at room temperature.Then be washed with TBST again, the membranes were scanned using a fluorescent scanner.\u003c/p\u003e\u003ch2\u003e2.12 Immunohistochemistry(IHC)analysis\u003c/h2\u003e\u003cp\u003eIHC was performed on aorta tissue sections.Each sample was fixed in formalin and embedded in paraffin.The blocks were sliced into 5 µm-thick sections,After the tissue sections were dewaxed,antigen retrieval was performed,hydrated in a graded series of alcohols,and subjected to heat-activated antigen retrieval.After blocking endogenous peroxidase activity,the tissue was incubated with OLR1,SPP1 and EGF antibodies at room temperature for 4 h.Subsequently,the sections were washed and incubated with biotinylated secondary antibody at room temperature for 30 min.The reaction complexes were visualized with diaminobenzidine and counterstained with hematoxylin.Observe the morphology under the microscope.\u003c/p\u003e\u003ch2\u003e2.13 Enzyme-linked immunosorbent assay (ELISA) analysis\u003c/h2\u003e\u003cp\u003eThe serum SOD and MDA were determined using an ELISA kit according to the manufacturer's instructions. Specifically, after separating and collecting mice blood sample, the serum was obtained by centrifugation. Use the following commercial kits to evaluate various substances in serum by ELISA analysis: MDA (Mouse MDA ELISA KIT; Shanghai Enzyme Linked Biotechnology Co., Ltd.), SOD (Mouse SOD ELISA KIT; Shanghai Enzyme Linked Biotechnology Co., Ltd.).\u003c/p\u003e\u003ch2\u003e2.14 Statistical analysis\u003c/h2\u003e\u003cp\u003eUnless otherwise specified, all experiments are conducted in no less than three independent parallel studies, and each parallel study consists of three repeated measurements. All data were analyzed using IBM's SPSS Statistics 21.0 software and expressed as mean ± standard error. One-way analysis of variance (ANOVA) is used for statistical analysis of comparisons between different groups. For all tests, a P-value \u0026lt; 0.05 was considered statistically significant. Unless marked with an asterisk, differences between groups can be considered statistically insignificant.\u003c/p\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec24\" class=\"Section2\"\u003e \u003ch2\u003e3.1 SJZD regulates Serum lipid profile and Aortic lesion areas in ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003emice.\u003c/h2\u003e \u003cp\u003eIn order to observe the anti-atherosclerotic effect of SIZD, we first measured the level of blood lipids.In the MG, the serum levels of TC and LDL-C increased significantly than those in the NG(\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01, Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA), while those of HDL-C decreased significantly(\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01).The ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice in SJZD exhibited significantly lower serum levels of TC and LDL-C than the ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice in the model group (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01).\u003c/p\u003e \u003cp\u003eThen we carried out HE and Oil Red O staining.HE staining results (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB) showed that the endothelium of NG arranged regularly. No obvious atherosclerotic lesion was observed. In the MG, a greater amount of the number of aortic atherosclerotic plaque was observed. However, SJZD treatment significantly reduced the lipid deposition in the aorta of ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice and the area of atherosclerotic plaque was smaller.\u003c/p\u003e \u003cp\u003eFurthermore, Oil Red O staining (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). In the MG, atherosclerotic plaque development was observed, with the intimal thickening of the aorta, and an increased amount of massive Oil Red O-stained lipid. The amount of Oil Red O-stained lipids in the SJZD decreased significantly.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec25\" class=\"Section2\"\u003e \u003ch2\u003e3.2Analysis of DEGs altered after treatment by Transcriptomic sequencing\u003c/h2\u003e \u003cp\u003eA total of 631.24\u0026nbsp;million raw reads were obtained in the sequencing. After filtering, 619.38 (98.12%) clean reads were obtained for following analysis. 600.13\u0026nbsp;million (95.07%) clean reads mapped to reference genome and 19.25\u0026nbsp;million clean reads were uniquely mapped. In order to further explore the mechanism of SJZD anti-AS, we performed transcriptome sequencing of mouse aorta samples from NG, MG and SJZD.\u003c/p\u003e \u003cp\u003eConsequently,in order to identify the potential genes associated with AS, we analyzed DEGS. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA, a heat map distribution of those DEGs. The DEGs shown in red indicate a high expression, whereas those in blue indicate a low expression. The heat map demonstrated that several regions with highly altered (up-regulated or down-regulated) expression were shown in SJZD, NG, and MG groups.\u003c/p\u003e \u003cp\u003eDEGs were further filtered by applying volcano plots, where the x-axis shows the Log2FC value and the y-axis corresponds to the value of log10 (p-value). The red dots correspond to significantly up-regulated DEGs (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05; |Log2FC|\u0026gt;1); the green dots correspond to significantly downregulated DEGs (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05; | Log2FC|\u0026gt;1); and the black dots indicate DEGs that were no statistical significance (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026gt;\u0026thinsp;0.05; |Log2FC|\u0026lt;1). There were 2754 DEGs in the MG/NG comparison (total gene counts: DEG-up =\u0026thinsp;2162; DEG-down =\u0026thinsp;592; Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB). Moreover, in the SJZD/AS comparison, there were 1680 DEGs (total gene counts: DEG-up =\u0026thinsp;540, DEG-down =\u0026thinsp;1140; Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003e3.3 PPI network analysis of common SJZT targets to determine the core targets and bioinformatics analysis of predictive targets.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eAs indicated in the Venn diagrams, in the mice treatment by SJZD, 33 genes downregulated by AS and upregulated by SJZD, and 306 genes upregulated by AS and downregulated by SJZD. Fortunately, 339 genes were rescued, indicating that SJZD reversed some AS-induced genes. Those 339 overlapping genes that indicate the opposite trends were marked as latent genes of interest (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA-B). The PPI network of the 339 latent targets was established by Cytoscape 3.7.2. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC, some of the remarkable genes in the PPI network included PON1, EGF, SLPI were initially upregulated in AS, but later afterwards reversed by SJZD.\u003c/p\u003e \u003cp\u003eThe GO enrichment results showed that 339 latent targets of SJZD were involved in biological processes, including lipid localization, eukocyte cell-cell adhesion, and leukocyte proliferation. The cellular component involved, including the membrane raft, nuclear chromatin, and nuclear transcription factor complex. It was closely related to molecular function, such as receptor ligand activity and cytokine receptor binding. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eD-F, the GO enrichment results were scalped. The results of KEGG pathways enrichment showed that 339 DEGs were enriched into nine pathways, in which the pathways with \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05 were defined as the key pathways (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eG). Among them, Including AGE-RAGE signaling pathway in diabetic complications, Th17 cell differentiation, and Non-alcoholic fatty liver disease are AS-related signal pathways.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec26\" class=\"Section2\"\u003e \u003ch2\u003e3.4 Determination of effective components and common anti-AS indexes of SJZD.\u003c/h2\u003e \u003cp\u003eA total of 136 active ingredients of four herbal medicines in SJZ decoction were identified based on threshold values of OB\u0026thinsp;\u0026gt;\u0026thinsp;30% and DL\u0026thinsp;\u0026gt;\u0026thinsp;0.18, including 22 constituents in ginseng, 15 constituents in poria cocos, 7 constituents in atractylodes, 92 constituents in licorice after deleting the duplicate data. Details of the active ingredients are listed in (supplement Table S1). In addition, 258 targets of active constituents of SJZ decoction were obtained from the TCMSP database and the gene names of these targets were collected via the Uniprot database (Supplementary Table S2).\u003c/p\u003e \u003cp\u003eFurthermore, we had collected AS-related genes/targets from Gene Cards database, CTD database, TTD database. And we used TCMSP to supplement them and then added targets. Finally, we totally collected 3729 AS-related targets. Information on these targets is provided in Supplementary Table S3.\u003c/p\u003e \u003cp\u003eAfter mapping SJZD targets and AS targets, a total of 171 common targets were obtained (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). To better understand the interactions between common targets, PPI network was established to better interpret the mechanisms of SJZD in AS treatment by using STRING software (combined score\u0026thinsp;\u0026ge;\u0026thinsp;0.9). As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB, the network consisted of 165 nodes and 1283 edges with an average node degree of 7.776, network centralization of 0.342, and an average number of neighbors of 15.552. According to the degree principle of each target, AKT1, IL6, MAPK1, JUN, and MAPK8 were determined as the hub targets (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003eTo reflect the interactions intuitively between the active compounds of SJZD and their potential targets, the C\u0026ndash;T network was constructed by mapping 114 active compounds to their 208 corresponding potential targets. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eD, the network consisted of 322 nodes (114 active compound nodes, and 208 compound associated target nodes) and 1400 interaction edges.In the network, quercetin (degree 136), kaempferol (degree 85), and 7-Methoxy-2-methyl isoflavone (degree 31) presented the maximum interactions with potential targets, indicating that these active compounds with high degree values could play an significant role in the latent pharmacological effects of SJZD. The C-T network revealed intimate communications between active constituents and related targets, which provided a reference to further investigate the pharmacological mechanisms of SJZD.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec27\" class=\"Section2\"\u003e \u003ch2\u003e3.5 Gene ontology enrichment analysis and targetpathway network construction\u003c/h2\u003e \u003cp\u003eThe GO functional analyses is carried out in the David database to reveal the function of 171overlapping targets. GO enrichment analysis included 197biological processes (BP), 59 cellular components (CC) and 126 molecular functions (MF) with a threshold value of \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05. The GO analysis results are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA\u0026ndash;C. The BP results mainly comprised of positive regulation of transcription from RNA polymerase II promoter, positive regulation of transcription, DNA-templated, response to drug. The CC analysis indicated that the overlapping targets were mainly related to cytosol, nucleoplasm, extracellular space. The MF results included protein binding, enzyme binding, identical protein binding.\u003c/p\u003e \u003cp\u003eThe results of the KEGG enrichment analysis indicated that 124 pathways meeting the threshold value of \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05 were significantly enriched. As depicted in Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eD, the KEGG pathways of SJZD against atherosclerosis were mainly related to P13K-Akt signaling pathway, TNF signaling pathway, MAPK signaling pathway, HIF-1 signaling pathway, and FpxO signaling pathway.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec28\" class=\"Section2\"\u003e \u003ch2\u003e3.6 Structuration of ingredient-gene-pathway (IGP) network and verification of core genes by qRT-PCR.\u003c/h2\u003e \u003cp\u003eIn order to find a reliable signal pathway for SZJD to improve AS, we chose DEGs and GO-biological process pathway, in which the pathways with \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05, by intersecting results of network pharmacology with transcriptome (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eA). The results of GO-biological processes pathway revealed that 339 DEGs were enriched including lipid localization. Ingredients, genes, and pathways were also introduced into Cytoscape to establish a visualized Ingredient-gene-pathway (IGP) network for further analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eB). We performed qRT-PCR analyses to investigate whether SJZD affected the transcriptional regulation of the genes screened above. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, EGF, OLR1, COL1, SLP, SPP1, PLAU and MMP2 mRNA were significantly increased in ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003emice, while PON-1mRNA significantly reduced. compared with those in normal mice (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0,01). The expression of LDLR mRNA was not statistically significant. The expression of LDLR mRNA was no significant difference (\u003cem\u003eP\u003c/em\u003e\u0026gt;0.05).The expressions of PON-1and LDLR mRNA was significantly increased and the expressions of EGF, OLR1, COL1, SPP1, PLAU and MMP2 mRNA was significantly reduced under SJZD treatment (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05 or (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01). The expression of SLP1 mRNA was no significant difference (\u003cem\u003eP\u003c/em\u003e\u0026gt;0.05, Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eC).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec29\" class=\"Section2\"\u003e \u003ch2\u003e3.7SJZD regulates OLR1/SPP1/EGF pathway to affect endothelial injury and prevent atherosclerosis\u003c/h2\u003e \u003cp\u003eAmong the genes screened above, we found that there is a link between OLR1,SPP1 and EGF, which is related to endothelial injury and oxidative stress. To examine whether antioxidative effect of SJZD is mediated by OLR1,SPP1 and EGF, we detected their protein expression in mice aorta. The protein levels of OLR1,SPP1 and EGF were significantly increased in the model group, compared to the normal mice (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01). SJZD intervention decreased the expression levels of these proteins(\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01, Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eA). The above results suggest that SJZD can improve vascular endothelial injury in atherosclerotic mice by down-regulating the expression of OLR1,SPP1 and EGF protein.\u003c/p\u003e \u003cp\u003eIn order to further explore whether SJZD can improve endothelial function and reduce oxidative stress in ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice, the expression of OLR1, SPP1 and EGF were detected by immunehistochemistry. As shown in Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eB, immunohistochemical(IHC)analysis of OLR1,SPP1 and EGF showed that the area occupied by these proteins in the lesions of ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice was significantly higher than that of normal mice ( OLR1 was30.21vs. 18.39%, SPP1 was 25.40 vs.10.16%, EGF was 29.13 vs.11.80%, \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01) However, after SJZD administration, this situation was obviously reversed, and plaque deposition and these proteins occupied area of aorta decreased significantly (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01 or \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.05,Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eC). These results were consistent with the changes in protein levels in mouse aortas, as revealed by Western Blot. All in all, these results indicate that SJZD can reduce the level of OLR1, SPP1 and EGF in atherosclerotic lesions and improve oxidative stress.\u003c/p\u003e \u003cp\u003eAccording to the above results, we speculate that the therapeutic effect of Sijunzi decoction on ApoE\u003csup\u003e\u003cb\u003e\u0026minus;/\u0026minus;\u003c/b\u003e\u003c/sup\u003emice may be related to oxidative stress. In order to further verify the effect of SJZD treatment on oxidative stress in ApoE\u003csup\u003e\u003cb\u003e\u0026minus;/\u0026minus;\u003c/b\u003e\u003c/sup\u003e mice, we detected the levels of SOD and MDA in serum of mice. compared to that from normal mice, the serum levels of anti-oxidative SOD was significantly reduced and pro-oxidative MDA was increased in ApoE\u003csup\u003e\u003cb\u003e\u0026minus;/\u0026minus;\u003c/b\u003e\u003c/sup\u003e mice (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01). However, the serum level of SOD increased and the level of MDA decreased in ApoE\u003csup\u003e\u003cb\u003e\u0026minus;/\u0026minus;\u003c/b\u003e\u003c/sup\u003e mice treated with SJZD (\u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u0026lt;\u0026thinsp;0.01, Figs.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eD). These results indicated that SJZD treatment results in a marked suppression in oxidative stress in ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e mice through OLR1/SPP1/EGF pathway.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec30\" class=\"Section2\"\u003e \u003ch2\u003e3.8 Validation of compound-target interaction\u003c/h2\u003e \u003cp\u003eTo improve the accuracy of the connection between small molecular compounds and the core target proteins, we used molecular docking to evaluate the interactions between compounds and target proteins. It is generally believed that when the docking binding energy less than \u0026minus;\u0026thinsp;4.0 kcal/mol, it could be considered to have potential binding activity between the ligand and receptor \u003csup\u003e[\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]\u003c/sup\u003e. A binding energy score less than \u0026minus;\u0026thinsp;5 kcal/mol suggests strong binding between the ligand and receptor \u003csup\u003e[\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e]\u003c/sup\u003e. Molecular docking indicated that kaempferol (A), quercetin (B) and naringenin (C) had good affinity for the three important compounds, SPP1, EGF, OLR1, in SJZT.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eMany studies have shown that SJZD is a promising formula for treatment of cardiovascular related diseases. Preliminary basic research shows that SJZD has the effect of anti-atherosclerosis, however, the pharmacological mechanism of SJZD against AS is not yet fully understood, and the effects of its bioactive compounds and targets remain elusive. Therefore, by applying network pharmacology and transcriptomic methods, this study, which is combined with in vivo experiment, was carried out to identify the underlying mechanism and therapeutic targets of SJZD against AS effect.\u003c/p\u003e \u003cp\u003eFirstly,the results about serum lipids and atherosclerotic plaques of aorta by HE and Oil Red O suggested that ApoE\u003csup\u003e\u0026minus;/\u0026minus;\u003c/sup\u003e AS mice model is conducted successfully and SJZT has a regulatory effect on them,which is consistent with previous studies.Then, transcriptomics sequencing was used to obtain active targets of AS.Through transcriptomics results of aorta tissue,339 genes were screened out,which indicated that SJZD may reverse some AS-induced genes.The PPI network analysis showed that some of the notable genes were reversed by SJZD,which means these genes may serve as potential targets and play an important role in the development of AS.The data of GO and KEGG showed that the importance of AS related terms was extremely high, which supports our previous hypothesis that the reliability of SJZD in the treatment of AS.Furthermore,a network pharmacology approach was used to identify 114 active compounds and 208 hub targets.The component-target network indicated that SJZD compounds might affect multiple targets and possess overlapping some,which results in synergistic effects.GO enrichment analysis showed that 208 targets were enriched in some biological processes,such as inflammatory response,angiogenesis and anti-apoptosis.KEGG analysis demonstrated that through the inhibition of some signal pathways,such as PI3K-Akt,TNF, MAPK and HIF-1, SJZD could resist AS. Finally, to make further analysis, we coupled network pharmacology with known transcriptomic results, and a visualized IGP network was established. According to the results,4 biological process (lipid localization, specific granule, positive regulation of lipid transport, positive regulation of lipid localization), 7 targets (SPP1, EGF, OLR1, LDLR, PON1, SLP1, PLAU), and 6 active ingredients (including beta-sitosterol, kaempferol, quercetin, stigmasterol, naringenin and isorhamnetin) play an extremely important role in the treatment of AS with SJZD.\u003c/p\u003e \u003cp\u003eFrom above, we find that four biological processes involved in AS are mainly related to lipids and specific granules. As it is known, AS is a process caused by multiple factors, and lipid accumulation and subsequent inflammation have been shown to play a critical role. During further progression of atherosclerotic lesion, many other specific granules, such as low-density lipoprotein (LDL), come into play, which results in, locally or systemically, chronic low-grade inflammatory state. The infiltration and retention of atheroprone lipoproteins, predominately LDL, in the vessel wall, is considered as an initial step in the development of AS \u003csup\u003e[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]\u003c/sup\u003e. The transportation of atheroprone particles, like LDL, across the endothelium, is critical to the initiation, progression, and regression of AS \u003csup\u003e[\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e]\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn this study, we attempted to explore a prospective novel therapeutic for AS based on the bioinformatic methods. In vivo and some animal models have shown that quercetin has a wide range of biological actions, including anti-carcinogenic, antiviral,anti-inflammatory,antioxidant,and psychostimulant activities,as well as the ability to inhibit lipid peroxidation,capillary permeability and platelet aggregation, and to stimulate mitochondrial biogenesis\u003csup\u003e[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]\u003c/sup\u003e.Recently,the mechanism of the anti-atherosclerotic effects of quercetin have been investigated extensively\u003csup\u003e[\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]\u003c/sup\u003e.In addition,quercetin suppressed vascular ROS formation endothelial dysfunction in HFD-fed ApoE-/- mice,with beneficial effects on atherosclerotic plaque formation. Naringenin is a major flavanone in citrus fruits that has multiple pharmacological attributes, such as anticancer and antiatherogenic antioxidant activity and free radical scavenging effects.Naringenin plays an important role in resisting cardiovascular diseases\u003csup\u003e[\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]\u003c/sup\u003e.Naringenin attenuates endothelial dysfunction by reducing ROS accumulation and increasing nitric oxide production in endothelial cells \u003csup\u003e[\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]\u003c/sup\u003e.Finally, we performed molecular docking to identify the reliability of the interactions of 3 active compounds and 3 core targets based on the previous study.We found that quercetin, kaempferol and naringenin bind most closely to the EGF, SPP1 and OLR1 proteins. In terms of secondary structure,quercetin and kaempferol are similar in structure.This shows that the binding pockets of kaempferol and quercetin are almost the same on protein,but the binding free energy and binding site are not necessarily the same. These results, which are consistent with our previous findings, show that the strategy of applying network pharmacology to find potential active compounds is reliable and feasible.\u003c/p\u003e"},{"header":"5. Conclusions","content":"\u003cp\u003eIn conclusion,benefiting by the integration of network pharmacology, transcriptomics,further experimental verification and validation of compound-target interaction,our study clarified that LOX-1/SPP1/EGF pathway may be one of the possible targets of SJZD to antiatherogenic.Kaempferol,quercetin and naringenin are the effective constitutes of SJZD to antiatherogenic.These findings are expected to guide the therapeutic research of AS in the future,and to provide a methodological reference for the pharmacodynamic material basis and targets of TCM in treatment of diseases by the integrated analysis combined with network pharmacology and transcriptomics.Although we have examined the possible role of SJZD in treating AS,there are some limitations in the present study.Further in vivo experimental validation is needed to support our research.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eCompliance with ethical standards\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAll animal experiments and methods were performed in accordance with the Health Guide for the Care and Use of Laboratory Animals published by the National Institute of Health.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFinancial support and sponsorship\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis study was supported by National Natural Science Foundation of China grant (No. 81874372 to Wenna Chen);the Natural Science Foundation of Liaoning Province (20180550767).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConflicts of interest\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThere are no conflicts of interest.\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eWeiyan Chen were responsible for conception and design. Tianmin Ji and Ying Yang was responsible for manuscript writing. Zhuo Zhao, Hao Gao, were responsible for collection and assembly of data. Lijiang Zhou and Ying Wang were responsible for data analysis and interpretation. All authors were responsible for manuscript writing. All authors were responsible for the final approval of the manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eZhang Y, Tang Y, Yan J (2022) LncRNA-XIST Promotes Proliferation and Migration in ox-LDL Stimulated Vascular Smooth Muscle Cells through miR-539-5p/SPP1 Axis. Oxid Med Cell Longev. Jan 4; 2022:9911982\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang R, Wang M, Ye J, Sun G, Sun X (2021) Mechanism overview and target mining of atherosclerosis: Endothelial cell injury in atherosclerosis is regulated by glycolysis (Review). Int J Mol Med 47(1):65\u0026ndash;76\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang R, Wang M, Ye J, Sun G, Sun X (2021) Mechanism overview and target mining of atherosclerosis: Endothelial cell injury in atherosclerosis is regulated by glycolysis (Review). Int J Mol Med 47(1):65\u0026ndash;76\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHou X, Yang S, Yin J (2019) Blocking the REDD1/TXNIP axis ameliorates LPS-induced vascular endothelial cell injury through repressing oxidative stress and apoptosis. Am J Physiol Cell Physiol 316(1):C104\u0026ndash;C110\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMeng B, Li Y, Ding Y, Xu X, Wang L, Guo B, Zhu B, Zhang J, Xiang L, Dong J, Liu M, Xiang L, Xiang G (2021) Myeloid-derived growth factor inhibits inflammation and alleviates endothelial injury and atherosclerosis in mice. Sci Adv 7(21):eabe6903. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1126/sciadv.abe6903\u003c/span\u003e\u003cspan address=\"10.1126/sciadv.abe6903\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang C, Niimi M, Watanabe T, Wang Y, Liang J, Fan J (2018) Treatment of atherosclerosis by traditional Chinese medicine: Questions and quandaries. Atherosclerosis 277:136\u0026ndash;144\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHuang C, Zhu Z, Cao X, Chen X, Fu Y, Chen Z, Li L, Song X, Jia R, Yin Z, Ye G, Feng B, Zou Y (2017) A Pectic Polysaccharide from Sijunzi Decoction Promotes the Antioxidant Defenses of SW480 Cells. Molecules 22(8):1341. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3390/molecules22081341\u003c/span\u003e\u003cspan address=\"10.3390/molecules22081341\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRen Xiuyu K, Ziwei Z, Yanchun Q, Chaoyi L, Yafu, Li Jingke (2021) The effect of Sijunzi decoction on the apoptosis of photoreceptor cells after retinal detachment in rats with spleen deficiency syndrome. China J Integr Traditional Chin Western Med 41(02):204\u0026ndash;211\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang R, Zhu X, Bai H, Ning K (2019) Network Pharmacology Databases for Traditional Chinese Medicine: Review and Assessment. Front Pharmacol 10:123\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi AP, He SS, Zhang WN, Zhang LC, Liu YT, Li K, Qin XM (2020) Exploration the active compounds of Astragali Radix in treatment of adriamycin nephropathy by network pharmacology combined with transcriptomic approach. J Ethnopharmacol 258:112537. \u003cdiv class=\"ExternalRefDOI\"\u003e10.1016/j\u003c/div\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDai S, Wang H, Wang M, Zhang Y, Zhang Z, Lin Z (2021) Comparative transcriptomics and network pharmacology analysis to identify the potential mechanism of celastrol against osteoarthritis. Clin Rheumatol 40(10):4259\u0026ndash;4268\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHe S, Li A, Zhang W et al (2020) An integrated transcriptomics and network pharmacology approach to exploring the mechanism of adriamycin-induced kidney injury. Chem Biol Interact 325:109096\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi J, Luo H, Liu X, Zhang J, Zhou W, Guo S, Chen X, Liu Y, Jia S, Wang H, Li B, Cheng G, Wu J (2020) Dissecting the mechanism of Yuzhi Zhixue granule on ovulatory dysfunctional uterine bleeding by network pharmacology and molecular docking. Chin Med 15:113. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1186/s13020-020-00392-0\u003c/span\u003e\u003cspan address=\"10.1186/s13020-020-00392-0\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGuo C, Kang X, Cao F, Yang J, Xu Y, Liu X, Li Y, Ma X, Fu X (2021) Network Pharmacology and Molecular Docking on the Molecular Mechanism of Luo-hua-zi-zhu (LHZZ) Granule in the Prevention and Treatment of Bowel Precancerous Lesions. Front Pharmacol 12:629021. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3389/fphar\u003c/span\u003e\u003cspan address=\"10.3389/fphar\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu J, Liu M, Yao Y, Wang J, Li Y, Li G, Wang Y (2012) Identification of novel potential β-N-acetyl-D-hexosaminidase inhibitors by virtual screening, molecular dynamics simulation and MM-PBSA calculations. Int J Mol Sci 13(4):4545\u0026ndash;4563. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3390/ijms13044545\u003c/span\u003e\u003cspan address=\"10.3390/ijms13044545\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTabas I, Williams KJ, Boren J (2007) Subendothelial lipoprotein retention as the initiating process in atherosclerosis: update and therapeutic implications. Circulation 116:1832\u0026ndash;1844\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang X, Fern\u0026aacute;ndez-Hernando C (2020) Transport of LDLs into the arterial wall: impact in atherosclerosis. Curr Opin Lipidol 31(5):279\u0026ndash;285\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Y, Zhao L, Qi W (2022) Uric acid, as a double-edged sword, affects the activity of epidermal growth factor (EGF) on human umbilical vein endothelial cells by regulating aging process. Bioengineered 13(2):3877\u0026ndash;3895\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBarreto J, Karathanasis SK, Remaley A, Sposito AC (2021) Role of LOX-1 (Lectin-Like Oxidized Low-Density Lipoprotein Receptor 1) as a Cardiovascular Risk Predictor: Mechanistic Insight and Potential Clinical Use. Arterioscler Thromb Vasc Biol 41(1):153\u0026ndash;166\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTian K, Ogura S, Little PJ, Xu SW, Sawamura T (2019) Targeting LOX-1 in atherosclerosis and vasculopathy: current knowledge and future perspectives. Ann N Y Acad Sci 1443(1):34\u0026ndash;53\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Y, Yao J, Han C, Yang J, Chaudhry MT, Wang S, Liu H, Yin Y, Quercetin (2016) Inflamm Immun Nutrients 8(3):167. \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3390/nu8030167\u003c/span\u003e\u003cspan address=\"10.3390/nu8030167\" targettype=\"DOI\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLuo M, Tian R, Lu N (2020) Quercetin Inhibited Endothelial Dysfunction and Atherosclerosis in Apolipoprotein E-Deficient Mice: Critical Roles for NADPH Oxidase and Heme Oxygenase-1. J Agric Food Chem 68(39):10875\u0026ndash;10883\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eOrhan IE, Nabavi SF, Daglia M, Tenore GC, Mansouri K, Nabavi SM (2015) Naringenin and atherosclerosis: a review of literature. Curr Pharm Biotechnol 16(3):245\u0026ndash;251\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeidary Moghaddam R, Samimi Z, Moradi SZ, Little PJ, Xu S, Farzaei MH (2020) Naringenin and naringin in cardiovascular disease prevention: A preclinical review. Eur J Pharmacol 887:173535\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"},{"header":"Supplementary Tables","content":"\u003cp\u003eSupplementary Tables S1 to S3 are not available with this version.\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"network pharmacology, transcriptomic techniques, TCM, Sijunzi Decoction, atherosclerosis, OLR1/SPP1/EGF pathway","lastPublishedDoi":"10.21203/rs.3.rs-4872904/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4872904/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eObjective:\u003c/strong\u003e Sijunzi Decoction (SJZD) is a kind of traditional Chinese medicine (TCM) formula,which has the contribution to anti-atherosclerosis.This study aims to explore the potential mechanism of the treatment of atherosclerosis (AS) with SJZD.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMethods:\u003c/strong\u003e Comprehensive analysis with transcriptomics and network pharmacology,combining with in vivo experiment and therapeutic targets.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResults: \u003c/strong\u003eThrough the in vivo experiment,it was found that SJZD could significantly improve the blood lipid level and lipid deposition in ApoE−/− mice.Furthermore,4 biological processes(lipid localization,specific granule,positive regulation of lipid localization,positive regulation of lipid transport),7 targets (SPP1,EGF,OLR1,LDLR,PON1,SLP1,PLAU), and 6 active ingredients(including beta-sitosterol,kaempferol,quercetin,stigmasterol, naringenin and isorhamnetin) play an extremely important role in the treatment of AS with SJZD.We also found that the genes of OLR,SPP1 and EGF were involved in regulating oxidative stress during the progress of atherosclerosis and the improvement of endothelial dysfunction.Furthermore,SJZD could interfere the expression of mRNA and protein of OLR1,SPP1 and EGF,and respectively reduce and increase the density level of MDA and SOD in serum significantly.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConclusion:\u003c/strong\u003eInhibiting excessively oxidative stress and improving endothelial dysfunction by regulating OLR1/SPP1/EGF pathway could be the mechanism,by which Sijunzi decoction resists AS.\u003c/p\u003e","manuscriptTitle":"Elucidating the anti-atherosclerosis mechanism of Si Jun Zi Decoction by integrating network pharmacology and transcriptomic with experimental validation in vivo","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-03 15:03:29","doi":"10.21203/rs.3.rs-4872904/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"a4f5e732-c922-4416-ac95-fb9b5c9f54bd","owner":[],"postedDate":"September 3rd, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2025-04-07T10:23:24+00:00","versionOfRecord":[],"versionCreatedAt":"2024-09-03 15:03:29","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4872904","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4872904","identity":"rs-4872904","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}