Comparative transcriptomic analysis revealed novel potential therapeutic targets of traditional Chinese medicine (Pinggan-Qianyang decoction) on vascular remodeling in spontaneously hypertensive rats
Chinese Medicine volume 16, Article number: 21 (2021)
Both experimental and clinical studies have revealed satisfactory effects with the traditional Chinese formula Pinggan Qianyang decoction (PGQYD) for improving vascular remodeling caused by essential hypertension. The present study explored various therapeutic targets of PGQYD using mRNA transcriptomics.
In this study, rats were randomly divided into three groups: Wistar-Kyoto (WKY; normal control), spontaneously hypertensive (SHR), and PGQYD-treated rat groups. After 12 weeks of PGQYD treatment, behavioral tests were employed and the morphology of thoracic aortas were examined with hematoxylin–eosin (HE) and Masson staining and electron microscopy. The mRNA expression profiles were identified with RNA-Seq and quantitative real-time PCR to validate changes in gene expression observed with microarray analysis. The gene ontology and pathway enrichment analyses were carried out to predict gene function and gene co-expressions. Pathway networks were constructed to identify the hub biomarkers, which were further validated by western blotting and immunofluorescence analysis.
After PGQYD treatment, the behavioral tests and histological and morphological findings of vascular remodeling were obviously meliorated compared with the SHR group. In the rat thoracic aorta tissues, 626 mRNAs with an exact match were identified. A total of 129 of mRNAs (fold change > 1.3 and P-value < 0.05) were significantly changed in the SHR group compared to the WKY group. Among them, 16 mRNAs were markedly regulated by PGQYD treatment and validated with quantitative real-time PCR. Additionally, target prediction and bioinformatics analyses revealed that these mRNAs could play therapeutic roles through biological processes for regulating cell metabolic processes (such as glycation biology), biological adhesions, rhythmic processes, and cell autophagy. The cellular signaling pathways involved in autophagy may be AGE–RAGE/PI3K/Akt/mTOR signaling pathway.
The present study provides novel insights for future investigations to explore the mechanisms by which PGQYD may effectively inhibit vascular remodeling by activating the AGE–RAGE/PI3K/Akt/mTOR signal pathway in cell autophagy biology.
Essential hypertension (EH) is a major public health problem both in middle-aged and elderly people. It is both a complex disease and an important risk factor for other cardiovascular outcomes, such as sudden death, stroke, myocardial infarction, heart failure, and renal diseases . Unfortunately, the control of arterial hypertension is far from optimal and has improved only minimally over the last few decades . Side effects of anti- hypertensive drugs, complaints due to their blood pressure lowering effect, and inadequate compliance are the key factors in the background of inadequate control of hypertension [3, 4]. Therefore, it is urgent to find new antihypertensive therapy in the future. The spontaneously hypertensive (SHR) rat model is a suitable model for studying the development and consequences of hypertension. The development of vascular remodeling is an early and important consequence of hypertension. Vascular remodeling is mainly characterized by vascular smooth muscle cell hypertrophy and increased production of the extracellular matrix . Vascular remodeling is initially an adaptive process that evolves in response to long-term pressure overload, but later it can contribute to the development of hypertensive target organ damage .
The development of new effective therapies for EH is urgently required. Fortunately, traditional Chinese medicine (TCM), with more than a 2000-year history, includes some of the oldest herbal medicines in the world . Based on a multi-component and multi-target approach, TCM has become an important source for new drug development to treat essential hypertension [8, 9]. The World Health Organization encourages the incorporation of herbal remedies into main stream medical systems . Pinggan Qianyang decoction (PGQYD), a famous TCM formula, is by far the most frequently used for the treatment of EH (87.9%) . Accumulating evidence has suggested that PGQYD improves target organ damage of EH through vascular protection, vascular anti-aging effects, and improvement in vascular remodeling . Moreover, our previous study confirmed that PGQYD tends to improve vascular remodeling by enhancing HSP27 expression through the p38MAPK signaling pathway . However, the multifaceted regulatory mechanisms of PGQYD acting on EH have not been fully elucidated due to the lack of appropriate methods.
Advanced omics technology is considered to be a holistic and efficient tool to study the use of TCM, and it can be used as a bridge between TCM and western medicine . Transcriptomics is a powerful tool for detecting global alterations in RNA expression and, consequently, changes in the corresponding proteins [15, 16]. A transcriptomics study of Buyang-Huanwu Decoction (BYHWD) for treating intracerebral hemorrhage showed that three tRNAs (rno-tRFi-Ser-25a, rno-tRF5-Ala-16a, and rno-tRF5-Glu-29a) are specifically regulated by these drugs . Another transcriptomics study of a TCM formula (Baoyuan decoction, BYD) for treating myocardial ischemia (MI) disclosed that the cardioprotection of BYD is mainly involved in the regulation of energy homeostasis, oxidative stress, apoptosis, inflammation, cardiac contractile dysfunction, and extracellular matrix remodeling . Therefore, it is rational to believe that transcriptomics analysis of mRNAs may help shed light on multifaceted mechanisms of TCM formulas by identifying exact therapeutic targets and their interactions . For that purpose, we aimed to utilize RNAsequencing technologies to evaluate the mRNA expression levels in Wistar-Kyoto (WKY; normal control) rats and the spontaneously hypertensive (SHR) and PGQYD-treated rats groups. Then hypertension-induced changed mRNAs (SHR vs WKY) and PGQYD-regulated altered mRNAs (PGQYD vs SHR) were obtained, and the key genes and pathways were identified with bioinformatics. Finally, the hub biomarkers were further validated by investigating protein expression levels. We evaluated their biological functions to reveal the therapeutic mechanisms of PGQYD (Fig. 1). These new findings may provide a novel perspective to illustrate the molecular mechanisms of PGQYD treatment for vascular remodeling after hypertension.
Materials and methods
Preparation of PGQYD
The ingredients of PGQYD included five medicinal herbs: Tian Ma (Gastrodiae Rhizoma), Gou Teng (Uncaria Rhynchophylla), Shi Jue Ming (Haliotis Discus Hannai), Mao Li (Ostrea Gigas Thunberg), and Niu Xi (Cyathulae Radix) at a ratio of 10:20:30:30:20, respectively (dry weight). All herbs were obtained from Xiangya Hospital, Central South University (CSU; Changsha, China). Furthermore, each herb was authenticated by the herbal medicinal botanist, Professor Shao Liu, in the Department of Pharmacy at Xiangya Hospital, CSU. The decoction was prepared and subjected to quality control as previously described . Finally, the powder was dissolved in distilled water at a concentration of 0.2 g/ml for intragastric administration.
Animals and experimental groups
Twenty healthy male spontaneous hypertensive rats [14 weeks old, SCXK (Hu) 2018-009] and 10 healthy male WKY rats [14 weeks old, SCXK (Hu) 2018-009] weighing 220 − 240 g were purchased from Shanghai Slake Experimental Animal Co., Ltd. All the animals were reared under identical conditions (room temperature at 25 °C and 12-h light–dark cycle) and had free access to standard rodent water and food. All animal protocols were approved by the Committee on the Use and Care of Animals of CSU and conformed to the Guidelines for the Care and Use of Laboratory Animals. After 1 week of adaptive feeding, the SHR rats with systolic blood pressure (SBP) more than 165 mmHg and the WKY rats with SBP less than 120 mmHg were used in this investigation. The first group was an age-matched normotensive control group (WKY, n = 10). The SHR rats were randomly divided into two groups. One group received no treatment (SHR, n = 10), whereas the other group (PGQYD, n = 10) received 0.7 g/kg/day as the dosage of PGQYD for 4 weeks.
Assessment of behaviors and measurement of blood pressure
The degree of irritability was observed and graded : 1 point: screamed and jumped when holding the neck; 2 points: bit while holding the neck; 3 points: when the tail was raised, the rat screamed, jumped, and even bit the rats in the same cage, frequently fought, and bit the iron cage. If the above situation was not obvious, it was scored as level 0. The rotation tolerance time was measured. The rats were placed on a home-made rotating platform with a rotating speed of 45 r·min-1. If the rats did not fall after rotating for 2 min, the experiment was stopped. SBP was determined in conscious rats using the tail-cuff method . The rat tails were occluded with the appropriate size tube-shaped tail cuff connected to the tail cuff device. The SBP and HR of the rats were measured via an XH200 thermostatic non-invasive blood pressure meter and an MD3000 biological signal collection and processing system (Huaibei Zhenghua Biological Instrument Equipment Co., Ltd., Anhui, China). Each result was averaged from three repeats at the beginning of the study, and at 1 week, 2 weeks, 3 weeks, and 4 weeks subsequently until euthanasia. On the 28th day, rats were deeply anesthetized by intraperitoneal injection of pentobarbital and perfused with ice-cold saline. Then, thoracic aortic tissues were harvested for subsequent analysis.
Histological and morphological assay
Rats were anesthetized with an intraperitoneal injection of 5 mL/kg 20% urethane at the end of each week of whole-day drug administration. The thoracic aorta below the aortic arch of each rat was stripped and clipped. A portion was fixed in 8% neutral formaldehyde, embedded in paraffin, sectioned at 5 µm, and stained with hematoxylin–eosin (HE) and Masson, and then observed with electron microscopy . Light microscopy was used to image each cross-sectional slice, of which there were five per rat. Each vascular ring in the perpendicular position and the vessel media wall was observed. The images were observed under a Leica imaging system (Leica Microsystems GmbH, Wetzlar, Germany). The media thickness (MT) and inner diameter (LD) were measured, and the ratio of media thickness to inner diameter (MT/LD) was calculated. Other parts of the thoracic aorta were removed from the adventitia and were promptly refrigerated at − 80 °C for processing in an mRNA assay.
Transcript profile analysis
The total RNA of thoracic aorta tissues from nine rats (WKY group, n = 3, SHR group, n = 3, and PGQYD group, n = 3) were extracted using TRIzol (Invitrogen, Carlsbad, CA, USA) reagent for RNA sequencing and were purified according to the manufacturer’s instructions. Total RNA was extracted from the right globus pallidus and purified with an RNeasy Mini Kit (Qiagen, Redwood City, CA, USA). Arraystar Rat mRNA microarrays (v2.0, containing 24,626 coding transcripts) were used to detect the expression of mRNAs in a total of nine rats (three groups with three replicates). Tissue preparation and microarray hybridization were performed with an Agilent Gene Expression Hybridization Kit (Agilent Technology, CA, USA) . After washing, the arrays were scanned with an Agilent Microarray Scanner and finally analyzed with Agilent Feature Extraction software (version 10.5.1.1). Differentially expressed transcripts (thresholds of ≥ 1.3-fold and P values of < 0.05) were identified by comparing the SHR model and normal control groups (SHR vs WKY) and the SHR + PGQYD-treated and SHR groups (PGQYD vs SHR). Then, we determined the intersection between the upregulated transcripts in the SHR vs WKY groups and the downregulated transcripts in the PGQYD vs SHR groups as well as the intersection between the downregulated transcripts in the SHR vs WKY groups and the upregulated transcripts in the PGQYD vs SHR groups to identify the mechanism by which PGQYD reverses pathophysiological changes in SHR.
Analysis of differentially expressed genes (DEGs)
The expression level for each gene was determined by the number of reads uniquely mapped to the specific gene and the total number of uniquely mapped reads in the sample. The gene expression level was calculated with the reads per kb per million (RPKM) method . The formula to calculate FPKM was as follows: FPKM = (number of mapping fragments) × 103 × 106/[(length of transcript) × (number of total fragments)]. Then, the NOI seq method was applied to screen DEGs between two groups, with the threshold of significance as fold change of RPKM ≥ 3 and probability ≥ 0.8 .
Pathway enrichment analysis and GO analysis
GO (http://www.geneontology.org) analysis was performed to determine biological roles based on the molecular functions, biological processes, and cellular components of the aberrantly expressed mRNAs; P < 0.05 and false discovery rate (FDR) < 0.05 were used as thresholds to define markedly enriched GO terms/pathways . Pathway analysis (based on KEGG, http://www.genome.jp/kegg/) was performed to explore the pathways significantly enriched in DEGs .
Real-time quantitative PCR validation
The procedure of real-time quantitative PCR was based on a previous study . Total RNA was extracted from tissues in each group and 5 µg of total RNA was reversely transcribed to cDNA (Tiangen Biotech, Shanghai, China) based on the manufacturer’s manuals. RT-qPCR was performed in a reaction system (25 µL) that contained SYBR Green/Fluorescein qPCR Master Mix, forward primers, reverse primers, and cDNA. The ABI7500 realtime PCR system (provided by Applied Biosystems, CA, USA) was employed for PCR. Alterations in mRNA expression in thoracic aortas tissues were assessed with the 2−△△Ct method normalized with endogenous control GAPDH. The PCR primers that were used are listed in Additional file 1: Table S1.
Immunofluorescence (IFC) analysis
The staining procedure for IFC analysis was based on a previous study . The 5 µm acetone-fixed thoracic aorta frozen sections were incubated for 1 h in 10% bovine serum to block nonspecific protein–protein interactions. They were cultured in PBS (overnight at 4 °C) and incubated with primary antibodies (rabbit anti-rat RAGE antibodies, 1:200, Abcam, Cambridge, MA, UK; rabbit anti-rat EGR-1 antibodies, 1:200, Abcam, Cambridge, UK). After washing three times for 10 min, the sections were incubated in PBS for 1 h at room temperature with secondary antibodies [goat anti-rabbit Alexa Fluor 594 (1:200)]. The nuclei were counterstained with 4-0-6-diamidino-2-phenylindole (DAPI; Invitrogen, Carlsbad, CA, USA). The immunofluorescence staining was analyzed using a laser-scanning confocal microscope (SLM 510, Carl Zeiss Meditec, Inc., Jena, Germany).
Western blotting analysis
From the frozen thoracic aorta tissues, proteins were extracted and the concentrations were measured with a bicinchonininc acid (BCA) protein assay. Samples were subjected to incubation with primary antibodies antiCaSR (provided by Abcam; 1:800), renin (provided by Bioss, Beijing, China; 1:250), AT2R (provided by Abcam; 1:800), angiotensin II type 1 receptor (AT1R) (provided by Abcam; 1:800), MasR (Abcam; 1:800), B cell leukemia-2 (Bcl-2) (provided by Cell Signaling Technology, Danvers, USA; 1:1000), Bcl-2- associated X protein (Bax; provided by Cell Signaling Technology; 1:800), and caspase-3 (provided by Abcam; 1:500), and normalized to GAPDH (provided by Sugisuke Bridge, China; 1:1000). Bio-Rad Quantity One software (provided by Bio-Rad) was used to quantify the intensity of the protein bands. The primary antibodies used are listed in Additional file 1: Table S2.
Statistical analyses were carried out with SPSS 19.0 (SPSS Inc., Chicago, IL, USA). Assays were conducted at least three times unless otherwise stated. All values were expressed as either the mean ± standard deviation (SD) or the mean ± the standard error (SE), except for the neurobehavioral scores. The comparison of data between groups was assessed using a one-way analysis of variance (ANOVA) followed by Tukey’s multiple-comparison test. Neurobehavioral scores were expressed as the median (interquartile range, IQR) and were analyzed using the Mann–Whitney U test. P < 0.05 was considered as statistical significance.
Effects of PGQYD on behaviors and blood pressure in SHRs
As shown in Fig. 2a, b, compared with the WKY group, the degree of irritability in the SHR group was significantly increased, while the rotation tolerance time was dramatically decreased (P < 0.01). Compared with the SHR group, the degree of irritability in the PGQYD group was significantly decreased (P < 0.05), and the rotation tolerance time was obviously prolonged (P < 0.05). The results showed that PGQYD could improve the behavior of SHR rats. As shown in Fig. 2c, before treatment, the blood pressure of the two groups of SHRs was significantly increased compared with that of the WKY group (P < 0.01). Compared with the SHR group, systolic blood pressure was reduced in the PGQYD group after 2 weeks (P < 0.05), and the efficacy was stable with the extension of the course of treatment (5th week). It was suggested that PGQYD could reduce the systolic blood pressure of SHR rats.
Effect of PGQYD on vascular remodeling in the morphology of SHRs
Masson and HE staining showed that the aortic tunica media of the SHR group was thicker than that of WKY group, and the aortic tunica media of PGQYD-treated rats was thinner than that of the SHR rats in the SHR group (Fig. 3a, b). As shown in Fig. 3d, e, both MT and MT/LD were higher in the SHR group than in the WKY group (MT: 118.6 ± 10.3 µm vs 85.3 ± 9.6 µm, respectively, P = 0.009; MT/LD: 1.83 ± 0.18 vs 1.25 ± 0.19, respectively, P = 0.02). However, both MT and MT/LD were significantly lower in the PGQYD group than in the SHR group (MT: 98.3 ± 11.2 µm vs 118.6 ± 10.3 µm, respectively, P = 0.02; MT/LD: 1.47 ± 0.21 vs 1.83 ± 0.18, respectively, P = 0.04).
Comparison of rat aortic electron microscopy results in each group are shown in Fig. 3c. In the WKY group, organelles in the cytoplasm of vascular smooth muscle cells were not well developed, mitochondria were occasionally seen, and the endoplasmic reticulum was visible. In the SHR group, organelles in the cytoplasm of vascular smooth muscle cells were more developed, more mitochondria were seen, and the endoplasmic reticulum was visible. In the PGQYD group, organelles in the plasma of vascular smooth muscle cells were more developed, and the mitochondria and endoplasmic reticulum were visible. It was suggested that PGQYD could improve morphological changes in the aortic vascular tissue and reverse vascular remodeling in SHR rats.
Classification of the transcriptomic phenotypes and identification of the altered transcripts with PGQYD treatment
A non-target global gene expression analysis was employed to explore how vascular remodeling altered aortic tissue mRNA expression, as well as evaluate the intervention with the PGQYD treatment on SHR rats. According to the cut-off criteria of |log2-fold change (FC)|> 1 and FDR < 0.05, a total of 465 DEGs were detected in the SHR group when compared to the WKY group (Fig. 4a, c), and a total of 511 DEGs were detected in the PGQYD group when compared to the SHR group (Fig. 4b, d). There were 129 and 70 DEGs (upregulated and downregulated, respectively) in the SHR group (relative to the WKY group) and PGQYD-treat (relative to the SHR group) group (Fig. 4e, f, respectively). Interestingly, PGQYD treatment could reverse these changes in the SHR group and created expression profiles close to WKY levels (Fig. 4f). Additionally, we identified 22 significantly dysregulated mRNAs in the SHR group: 16 were upregulated, while six were downregulated (SHR vs WKY). After PGQYD treatment, 22 mRNAs were obviously changed: six were upregulated, while 16 were downregulated (PGQYD vs SHR).
All expression levels of 22 candidate transcripts (Table 1) were validated at the transcript level by qRT-PCR (Fig. 5). The qRT-PCR results indicated statistically that the SHR group had increases in 11 gene expression levels (SMAD2, GSK-3, METTL3, CASP3, SIRT1, RAGE, HMGN4, PTGS1, JAK2, EIF4B, and DNMT1; Figs. 5a, b) and decreases in five gene expression levels (Egr-1, MFN2, HSP27, PTEN, and SEPN1; Fig. 5c).
Enrichment analysis of GO and KEGG pathways
GO and KEGG assignments were used to classify the genes associated with the response of vascular remodeling to PGQYD treatment. GO classification analysis provided the percentage and distribution of top-level GO terms that were portrayed in three categories: (A) cellular component; (B) molecular function, and (C) biological process. As shown in Figs. 6a and 7a, we found that transcripts regulated in the thoracic aorta tissues of these model rats, as compared with the WKY group and PGQYD treatment group, could be mapped to biological processes (BP) for regulating cell metabolic processes, biological adhesions, rhythmic processes, cell autophagy, cell components (CC) for cell junctions, membrane parts and extracellular regions, molecular functions (MF) for receptor binding, signal transducer activity, antioxidant activity, transcription factor activity, and catalytic activity. For further characterization of the DEGs, we performed pathway enrichment analysis with GSEA. As shown in Fig. 6b–d, we found that transcripts regulated in the thoracic aorta tissues of these model rats, compared with the WKY group, could be mapped to signaling pathways, such as the AGE–RAGE signaling pathway, oxidative phosphorylation, MAPK signaling pathway, and PI3K–AKT signaling pathway. As shown in Fig. 7b–d, we found that transcripts regulated in the thoracic aorta tissues of these model rats, compared with the PGQYD treatment group, could be mapped to signaling pathways, such as the AGE–RAGE signaling pathway, phagosomes, ECM-receptor interactions, and the PI3K–AKT signaling pathway.
Gene co-expression network analysis of the differences in the WKY, SHR, and PGQYD groups
According to transcriptome analysis of the three groups, 16 gene expressions were dysregulated (including 11 upregulated genes and six down regulated genes) and affirmed by qRT-PCR detection. These dysregulated genes were closely related to the vascular remodeling of hypertension treated by PGQYD. ESR1 was identified as the hub gene (Fig. 8a). Combined with the prediction of the PI3K–Akt signal pathway, AGE–RAGE signal pathway, and mTOR signal pathway were closely related to the vascular remodeling of hypertension treated by PGQYD (Fig. 8b, Additional file 1: Table S3). It was suggested that GSK-3, eIF4B, PTEN, CCND1, SMAD2, JAK2, CASP3, EGR-1, RAGE, and FOXO1 were closely related to the above-mentioned signaling pathways.
Validation of targeted genes enriched in the AGE–RAGE signaling pathway after PGQYD treatment
Six genes were enriched in the AGE–RAGE signaling pathway (SMAD2, JAK2, RAGE, CASP3, and EGR-1). These results indicated that the role of the AGE–RAGE signaling pathway was important in the SHR model rats after PGQYD treatment. Thus, the western blot method and immunofluorescence analysis were utilized to determine the expression levels of SMAD2, JAK2, RAGE, CASP3, and EGR-1 in thoracic aorta tissues. The results showed that SMAD2, CASP3, JAK2, and RAGE expression levels were higher in the SHR group compared with the WKY group. However, after vascular remodeling injury, the expression levels of EGR-1 were remarkably decreased (Fig. 9). To different extents, the PGQYD treatment improved the expression levels of these five proteins. Among them, EGR-1 had remarkably increased protein expression levels (Fig. 10a). SMAD2, CASP3, JAK2, and RAGE had significantly decreased protein expression levels compared with the SHR group (Figs. 10b, c). The above results indicated that PGQYD effectively regulated the AGE–RAGE signaling pathway in vascular remodeling. Regulation of the AGE–RAGE signaling pathway was involved in the PGQYD pharmaceutic effect on hypertension, which agreed with the predicted mechanism mentioned previously (Fig. 11).
In this study, the SHR rats displayed remarkably enhanced vascular remodeling. Moreover, rat aortic electron microscopy images indicated that organelles in the cytoplasm of vascular smooth muscle cells were more developed, more mitochondria were seen, and the endoplasmic reticulum was visible compared to the WKY rats. Gastrodin can inhibit the proliferation of vascular smooth muscle cells, and delay aging . In addition, Gastrodin can promote the endogenous vasodilator and inhibit the release of vasoconstrictor, thus showing a anti-hypertensive effect by increasing the synthesis of NO and decreasing the contents of ET and Ang-II . The alkaloids in U. rhynchophylla can inhibit the proliferation and senescence of endothelial cells, thus protecting vascular endothelial cells and delaying vascular remodeling . it can also increase the apoptosis rate of vascular smooth muscle cells, but inhibit its proliferation, showing an indirect inhibitory effect on vascular remodeling induced by Ang-II . Studies had shown that ethanol extract of Achyranthes bidentata and Cassia obtusifolia can inhibit the activity of angiotensin-converting enzyme involved in RAAS system [35, 36]. It can affect the concentration of serum calcium ion, thus affecting the calcium channel, and has a strong and lasting anti-hypertensive effect, especially for those caused by long-term tension . In the present study, our findings showed that PGQYD, a specific TCM for treating essential hypertension, could substantially improve blood pressure, behavioral, pathology (vascular remodeling) and endothelial dysfunction in essential hypertension, indicating there were unique advantages in the herb combinations and synergistic actions of multiple ingredients in the TCM.
Whether it is combined with western medicine to reduce blood pressure or alone, PGQYD had the advantages of multi-target, multi-channel, safe and effective. At present, the research on the mechanism of anti-hypertension has not been elucidated, mainly focusing on improving endothelial function, inhibiting the excessive activation of RAAS system, and improving the function of autonomic nervous system. We can further study its anti-hypertensive mechanism, anti-hypertensive pathway and the protective effect of other target organs by genomics and metabonomics. In the past, we had found that PGQYD could regulate the expression of some specific proteins(TpxII, HSP27 and ANXA1) and miRNAs(miR-20a, miR-145, miR-30, and miR-98) by used proteomics and miRNA gene chip technology [38, 39]. In this study, we first described the transcriptomic features of PGQYD-treated SHR rats, and KEGG analysis suggested that the regulated transcripts were attributed to multiple functions, such as cell metabolic processes, biological adhesions, rhythmic processes, and cell autophagy. Most genes involved in the three cellular functions were dysregulated. According to the pathway and gene co-expression network of PGQYD-specific genes, the AGE–RAGE signaling pathway, PI3K–Akt signal pathway, and SMAD2 were at the core of the networks.
Some evidence had supported the view that the glycation of proteins was one of the main factors contributing to vascular remodling and cardiac target organ damage in mild hypertensives [40, 41]. Treatment with sRAGE had alleviated vascular adverse remodeling in SHR, possibly via suppression of oxidative stress and inflammation, improvement in RAS balance, and activation of PPAR-γ pathway . The present study had showed that the AGE–RAGE may have close relationship with the pathogenesis of Qi and Yin deficiency and blood stasis in diabetic nephropathy patients . The natural product gastrodin had potent activities in lowering blood glucose, improving insulin resistance and ameliorating diabetic nephropathy by reducing the contents of advanced glycation end product (AGE) in renal tissues . In this study, the results suggested that the process of cell metabolism in glycation may play an important role in the vascular remodeling with PGQYD treatment.
Blood pressure exhibits a robust circadian rhythm in health. This balanced rhythm in hypertension is perturbed via elevations in night-time blood pressure, inflicting silent damage to the vasculature and body organs. These data suggest that circadian clock genes(Per2, bmal1) may act to inhibit or control renin/angiotensin signaling . Chinese herbal medicine for PGQYD may lower the blood pressure smoothly and recover the circadian rhythm of blood pressure in essential hypertension patients related to the regulation of vascular endothelium function . Tianma Gouteng Yin can regulate the circadian rhythm of clock gene in vascular smooth muscle cells induced by AngII through AT1R signaling pathway . Isorhynchophylline was the main alkaloids in U. rhynchophylla had certain effect on regulating the disordered circadian rhythm of clock genes in vascular smooth muscle cells (A7r5 cell line) induced by Ang-II through AT1R-Src signaling pathway . In this study, the results suggested that the process of rhythmic processes may play an important role in the vascular remodeling with PGQYD treatment.
In our study, these findings showed that the expressions of RAGE, SMAD2, JAK2, and CASP3 were upregulated and the expression of EGR-1 was downregulated in the SHR model rats. Among them, five mRNAs were markedly regulated with PGQYD treatment. The results indicated that the mechanism for PGQYD reversing vascular remodeling may be involved in the process of glycation regulation, which was carried out by RAGE, SMAD2, JAK2, CASP3, and EGR-1. Through bioinformatics analysis, we found that five genes are involved in regulation through the AGE–RAGE signaling pathway, PI3K–Akt signaling pathway, and mTOR signaling pathway. The study showed that paeoniflorin inhibited autophagy at least partly by inhibiting RAGE and upregulating the level of p-mTOR to act against AGE-induced mesangial cell dysfunction.
Autophagy plays an important role in cellular waste removal and structural reconstruction . Tanshinone IIA can inhibit Angiotensin II-induced proliferation and autophagy of vascular smooth muscle cells via regulating the MAPK signaling pathway . The total alkaloids of U. rhynchophylla can reduce the apoptosis of endothelial cells by enhancing autophagy in SHR . Gastrodin suppresses apoptosis possibly mediated by inhibiting autophagy in serum deprivation-stimulated H9C2 cardiomyocytes .The mTOR signaling pathway is involved in vascular remodeling against stress, such as pressure overload in VSMCs. Jin X et al. have shown that activation of the PI3K/Akt signaling pathway can reduce levels of reactive oxygen species and lipid deposition, thereby inhibiting endothelial dysfunction and reversing the progression of vascular remodeling . Nε-Carboxymethyl-Lysine (the key active component of AGEs) significantly decreased the phosphorylation of PI3K/AKT signaling and restoration of PI3K/AKT signaling in cell apoptosis. The results showed RAGE induced foam cell apoptosis in diabetic atherosclerosis by inhibiting the PI3K/AKT pathway . Rhynchophylline and Isorhynchophylline may promote the apoptosis and inhibit rhe proliferation and migration of VSMC by inducing Ang-II . Gastrodin inhibits PDGF-BB induced VSMC migration, its mechanisms may be associated with the inhibition of the JNK signaling pathway activation . The rhynchophylline can inhibite MPP+-triggered neurotoxicity by stimulating MEF2D via activating PI3-K/Akt/GSK3β signaling pathway .
The PI3K/Akt/mTOR signaling pathway is thought to be a crucial regulator of autophagy and apoptosis . It has been confirmed that the mechanisms of apoptosis and autophagy in vascular smooth muscle cells were involved in vascular remodeling of hypertension . It has also been confirmed that activation of the PI3K/Akt/mTOR pathway was detected in the proliferation, migration, and autophagy of VSMCs, which enhanced phenotype transformation . Tianma Gouteng Decoction can significantly regulate the autophagy-related pathway Ca2 + /AMPK/mTOR of vascular endothelial cells in spontaneous hypertensive rats, which is probably one of the important mechanisms of its antihypertensive effect . In our study, the results indicated that the PI3K/Akt/mTOR signal pathway was involved in the process of vascular remodeling, and this regulated the mechanism of PGQYD in reversing vascular remodeling. These findings suggested that the mechanism of PGQYD in reversing vascular remodeling was involved in the AGE–RAGE/PI3K/Akt/mTOR signaling pathway, which regulated the target gene expression and regulated the process of cell autophagy biology.
In conclusion, we performed a transcriptomics analysis of gene expression changes in SHR rats under treatment with PGQYD using a RNA-seq technique. Overall, a total of 465 predicted protein-encoding genes were identified. Further gene expression analysis revealed a total of 16 DEGs, including 11 upregulated genes and five downregulated genes. Intensive bioinformatics analysis identified cell regulation, cell metabolic processes, biological adhesions, rhythmic processes, and cell autophagy as significant biological processes. The AGE–RAGE signaling pathway, PI3K–Akt signal pathway, and mTOR signal pathway were mainly involved. Overall, the results from this study might provide insights into the understanding of the mechanisms for the response to PGQYD. Furthermore, this work demonstrated the potential utility of the RNA-seq technique in ant- hypertension studies.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request.
Ioannidis JPA. Diagnosis and treatment of hypertension in the 2017 ACC/AHA guidelines and in the real world. JAMA. 2018;319(2):115–6.
Xie X, Atkins E, Lv J, Bennett A, Neal B, Ninomiya T, et al. Effects of intensive blood pressure lowering on cardiovascular and renal outcomes: updated systematic review and meta-analysis. Lancet. 2016;387(10017):435–43.
Zhao W, Katzmarzyk PT, Horswell R, Wang Y, Li W, Johnson J, et al. Aggressive blood pressure control increases coronary heart disease risk among diabetic patients. Diabetes Care. 2013;36(10):3287–96.
Basu S, Sussman JB, Rigdon J, Steimle L, Denton BT, Hayward RA. Benefit and harm of intensive blood pressure treatment: Derivation and validation of risk models using data from the SPRINT and ACCORD trials. PLoS Med. 2017;14(10):e1002410.
Zhou TL, Henry RMA, Stehouwer CDA, van Sloten TT, Reesink KD, Kroon AA. Blood pressure variability, arterial stiffness, and arterial remodeling. Hypertension. 2018;72(4):1002–10.
Pires PW, Rogers CT, McClain JL, Garver HS, Fink GD, Dorrance AM. Doxycycline, a matrix metalloprotease inhibitor, reduces vascular remodeling and damage after cerebral ischemia in stroke-prone spontaneously hypertensive rats. Am J Physiol Heart Circ Physiol. 2011;301(1):H87-H97.0.
Chen B, Wang Y, He Z, Wang D, Yan X, Xie P. Tianma Gouteng decoction for essential hypertension: Protocol for a systematic review and meta-analysis. Medicine (Baltimore). 2018;97(8):e9972.
Xiong X, Yang X, Liu Y, Zhang Y, Wang P, Wang J. Chinese herbal formulas for treating hypertension in traditional Chinese medicine: perspective of modern science. Hypertens Res. 2013;36(7):570–9.
Wang J, Xiong X. Control strategy on hypertension in chinese medicine. Evid Based Complement Alternat Med. 2012;2012:284847.
He SM, Li CG, Liu JP, Chan E, Duan W, Zhou SF. Disposition pathways and pharmacokinetics of herbal medicines in humans. Curr Med Chem. 2010;17(33):4072–113.
Zhong GW, Luo Y, Xiang L, Xie Y, Xie QY, Li YH. Clinical efficacy study on calming liver and restraining Yang formula in treating patients with mild or moderate degree of essential hypertension. Zhongguo Zhong Yao Za Zhi. 2010;35(6):776–81 (Chinese).
Zhong GW, Chen MJ, Luo YH, Xiang LL, Xie QY, Li YH, et al. Effect of Chinese herbal medicine for calming Gan (肝) and suppressing hyperactive yang on arterial elasticity function and circadian rhythm of blood pressure in patients with essential hypertension. Chin J Integr Med. 2011;17(6):414–20.
Zhong GW, Li W, Luo YH, Chen MJ, Yi ZJ, Chen GL, et al. Effects of the calming liver and suppressing yang method on proliferation and the expression of heat shock protein 27 in vascular smooth muscle cells of spontaneously hypertensive rats. Zhong Guo Lao Nian Yi Xue Za Zhi. 2009;29(2):385–8 (Chinese).
Eichner J, Rosenbaum L, Wrzodek C, Häring HU, Zell A, Lehmann R. Integrated enrichment analysis and pathway-centered visualization of metabolomics, proteomics, transcriptomics, and genomics data by using the InCroMAP software. J Chromatogr B Analyt Technol Biomed Life Sci. 2014;966:77–82.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.
Lappalainen T, Sammeth M, Friedländer MR, et al. Transcriptome and genome sequencing uncovers functional variation in humans. Nature. 2013;501(7468):506–11.
Li P, Tang T, Liu T, et al. Systematic analysis of tRNA-derived small RNAs reveals novel potential therapeutic targets of traditional chinese medicine (Buyang-Huanwu-Decoction) on intracerebral hemorrhage. Int J Biol Sci. 2019;15(4):895–908.
Du Z, Shu Z, Lei W, Li C, Zeng K, Guo X, et al. Integration of metabonomics and transcriptomics reveals the therapeutic effects and mechanisms of Baoyuan decoction for myocardial ischemia. Front Pharmacol. 2018;9:514.
Zhou W, Wang Y. A network-based analysis of the types of coronary artery disease from traditional Chinese medicine perspective: potential for therapeutics and drug discovery. J Ethnopharmacol. 2014;151(1):66–77.
Guang-wei Z, Wei L, Min-jing C, Zhen-jia Y, Guo-lin C, Ling-li X. Effects on the vascular remodeling and adiponectin expression in aorta in the spontaneously hypertensive rats by chinese herb mixture method. Zhong Hua Gao Xue Ya Za Zhi. 2008;16(9):812–6 (Chinese).
Zhong GW, Xiang L, Hu J, Yin Y, Chen Q, Fang X. Effect of Pinggan Qianyang recipe on the expression of Tpx II HSP27 and ANXA1 in the hypothalamus of spontaneously hypertensive rats with hyperactivity of liver-YANG syndrome. Zhong Nan Da Xue Xue Bao Yi Xue Ban. 2015;40(2):136–43 (Chinese).
Lvhua Y, Shazhou Z, Yunxia L. New equipments and methods in measuring systoIic blood pressure and diastolic blood pressure of rats without wounds. Zhong Guo Ying Yong Sheng Li Za Zhi. 1991;7(1):62–4 (Chinese).
Zhu X, Zhou Z, Zhang Q, Cai W, Zhou Y, Sun H, et al. Vaccarin administration ameliorates hypertension and cardiovascular remodeling in renovascular hypertensive rats. J Cell Biochem. 2018;119(1):926–37.
Ren S, Peng Z, Mao JH, Yu Y, Yin C, Gao X, et al. RNA-seq analysis of prostate cancer in the Chinese population identifies recurrent gene fusions, cancer-associated long noncoding RNAs and aberrant alternative splicings. Cell Res. 2012;22(5):806–21.
Wagner GP, Kin K, Lynch VJ. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012;131(4):281–5.
Tarazona S, García-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21(12):2213–23.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Zhang T, Tang N, Xi D, Zhao Y, Liu Y, Wang L, et al. Calcimimetic R568 improved cardiac remodeling by classic and novel renin-angiotensin system in spontaneously hypertensive rats. Exp Biol Med (Maywood). 2019;244(10):789–801.
Duan J, Han X, Ling S, Gan W, Sun L, Ni RZ, et al. Aortic Remodelling Is Improved by 2,3,5,4’- tetrahydroxystilbene-2-O-β-D-glucoside involving the Smad3 pathway in spontaneously hypertensive rats. Evid Based Complement Alternat Med. 2015;2015:789027.
Gen-meng Y, Juan L, Shi-jun H, et al. Effects of gastrodin against SH-SY5Y cell autophagy induced by methamphetamine and underlying mechanisms. Zhong Guo Yao Li Xue Tong Bao. 2018;34(4):577–83 (Chinese).
Wu SY. Effects of gastrodin injection combined with RAAS inhibitor on arterial elasticity, vascular overload index and oxidation antioxidant system in elderly patients with essential hypertension. J Mod Chin West Med. 2017;26:3242–5.
Chao L, Wen-qing Y, Yu W, et al. Total alkaloids from Uncaria rhynchophylla enhance autophagy and reduce oxidative stress-induced inflammatory response of endothelial cells. Zhong Hua Zhong Yi Yao Za Zhi. 2017;32(5):2137–41 (Chinese).
Tingting S. Effect and mechanism of isorhynchophylline on angiotensin II induced cardiomyocyte hypertrophy in rats. Shan Dong Yi Yao. 2017;57(38):44–7 (Chinese).
Ai-cui M, Xue Y, Ying J, et al. Inhibitory effect of Traditional Chinese medicine for Pinggan Qianyang method ofhaliotidis concha on angiotensin converting enzyme (ACE) in vitro. Zhong Guo Hai Yang Da Xue Xue Bao Nat Sci. 2012;42(z2):135–7 (Chinese).
Ming Q, Guo X, Guo-qi Z. Effect of ethanol extract of Achyranthes bidentata on blood pressure, myocardial ACE activity and myocardial cell diameter in spontaneously hypertensive rats. Zhong Guo Xian Dai Zhong Yao. 2010;12(6):34–7 (Chinese).
Wei J, Jing-feng L, Jiu-tang G, et al. Chemical constituents and pharmacological effects of haliotidis concha. Ji Lin Zhong Yi Yao. 2015;35(3):218–21 (Chinese).
Fan R, He F, Wang Y, Zhong GW, Li YH. Changes of protein expression profile in vascular tissues of spontaneously hypertensive rats treated by a compound Chinese herbal medicine. Zhong Xi Yi Jie He Xue Bao. 2011;9(6):643–50 (Chinese).
Zhong G, Fang X, Wang D, Chen Q, Tang T. Effects of the Pinggan Qianyang recipe on MicroRNA gene expression in the aortic tissue of spontaneously hypertensive rats. Evid Based Complement Alternat Med. 2015;2015:154691.
Bai Y, Shi X, Ke Y, Lin X, Hong H. Hypertension accelerates age-related intrarenal small artery (IRSA) remodelling and stiffness in rats with possible involvement of AGEs and RAGE. Histol Histopathol. 2020;35(1):97–109.
Maresca AM, Guasti L, Bozzini S, Mongiardi C, Tandurella N, Corso R, Zerba FG, Squizzato A, Campiotti L, Dentali F, Klersy C, Grandi AM, Falcone C. sRAGE and early signs of cardiac target organ damage in mild hypertensives. Cardiovasc Diabetol. 2019;18(1):17.
Liu Y, Yu M, Zhang L, Cao Q, Song Y, Liu Y, Gong J. Soluble receptor for advanced glycation end products mitigates vascular dysfunction in spontaneously hypertensive rats. Mol Cell Biochem. 2016;419(1–2):165–76.
Yueguang D, Zonglei G, Kefu C. The study of relationship between material bases of early diabetic nephropathy as deficiency of both Qi and Yin and obstruction of collaterals and AGE-RAGE. Zhejiang Chinese Medical University Xue Bao. 2016;40(2):75–80 (Chinese).
Yong Z, Yan-bin W, Wei-jia K. Experimental study of anti-diabetic effects of gastrodin in KK-Ay mice. Zhong Guo Yao Li Xue Tong Bao. 2018;34(7):917–24.
Pati P, Fulton DJ, Bagi Z, Chen F, Wang Y, Kitchens J, et al. Low-salt diet and circadian dysfunction synergize to induce angiotensin II-dependent hypertension in mice. Hypertension. 2016;67(3):661–8.
Liu Y, Wu Z, Zhenfu Z, Zhuli Y, Zhangyu J, Ran Y, et al. The Effect and mechanism of Tianma Gouteng decoction on regulating circadian rhythm of clock gene in A7r5 cells. Trad Chin Drug Res Clin Pharmacol. 2018;29(3):307–12 (Chinese).
Zhu-Li Y, Wu Z, Zhen-Fu Z, Liu Y, Zhang-Yu J, Rong Z, et al. Regulatory effect of isorhynchophylline on circadian rhythm of clock genes in vascular smooth muscle cells through AngiotensinIIType 1 receptor. J Guangzhou Univ Trad Chin Med. 2018;35(4):699–705 (Chinese).
Du J, Zhang C, Zhao W. Autophagy and hypertension. Adv Exp Med Biol. 2020;1207:213–6.
Cha L, Wen-qin Y, Yu W, Yan-peng T, Fen J, Yun-lun L. Rhynchophylla total alkaloids reduce endothelial cell inflammatory response induced by oxidative stress through enhancing autophagy. Zhong Hua Zhong Yi Yao Za Zhi. 2017;32(5):2137–41 (Chinese).
Lu J, Shan J, Liu N, Ding Y, Wang P. Tanshinone IIA can inhibit angiotensin II-induced proliferation and autophagy of vascular smooth muscle cells via regulating the MAPK signaling pathway. Biol Pharm Bull. 2019;42(11):1783–8.
Yao L, Ping Y, Hua H, Fei W, Jun-quan Y, Lin S, et al. Anti-apoptotic effect of gastrodin through inhibiting serum deprivation- induced autophagy in rat H9C2 cardiomyocytes. Zhong Guo Xin Yao Za Zhi. 2016;25(3):328–33 (Chinese).
Jin X, Xu Z, Fan R, et al. HO-1 alleviates cholesterol-induced oxidative stress through activation of Nrf2/ERK and inhibition of PI3K/AKT pathways in endothelial cells. Mol Med Rep. 2017;16(3):3519–27.
Wang Z, Bao Z, Ding Y, et al. Nε-carboxymethyl-lysine-induced PI3K/Akt signaling inhibition promotes foam cell apoptosis and atherosclerosis progression. Biomed Pharmacother. 2019;115:108880.
Jin-guo Y, Yun-lun L, Hong-lei Z. Protective effect of total alkaloids from Uncariae Ramulus cure Uncis combined with Raphani Semen on vascular endothelial cell injury by hypertension. Zhong Cheng Yao. 2013;35(5):889–93 (Chinese).
Li-hua Z, Hong-jing G, Lang W, Song T, Da Y, Ming-yue F, et al. Effect of gastrodin on rat vascular smooth muscle cell migration induced by platelet-derived growth factor-BB. Int J Cerebrovasc Dis. 2012;20(3):189–92.
Hu S, Mak S, Zuo X, Li H, Wang Y, Han Y. Neuroprotection against MPP+-induced cytotoxicity through the activation of PI3-K/Akt/GSK3β/MEF2D signaling pathway by rhynchophylline, the major tetracyclic oxindole alkaloid isolated from Uncaria rhynchophylla. Front Pharmacol. 2018;9:768.
Zhang M, Li F, Wang X, et al. MiR-145 alleviates Hcy-induced VSMC proliferation, migration, and phenotypic switch through repression of the PI3K/Akt/mTOR pathway. Histochem Cell Biol. 2020;153(5):357–66. https://doi.org/10.1007/s00418-020-01847-.
Pi S, Mao L, Chen J, Shi H, Liu Y, Guo X, et al. The P2RY12 receptor promotes VSMC-derived foam cell formation by inhibiting autophagy in advanced atherosclerosis. Autophagy. 2020. https://doi.org/10.1080/15548627.2020.1741202.
Ciceri P, Falleni M, Tosi D, et al. Therapeutic effect of iron citrate in blocking calcium deposition in high Pi-calcified VSMC: role of autophagy and apoptosis. Int J Mol Sci. 2019;20(23):5925.
Yun-pan L, Mei-yun L. Effect of Tianma Gouteng Decoction on Ca2+/AMPK/mTOR signaling pathway of vascular endothelial cells in spontaneous hypertensive rats. Zhong Guo Xian Dai Yi Xue Za Zhi. 2018;28(22):38–43 (Chinese).
This work was financially supported by Natural Science Foundation of China (Grant Nos. 81673858, 81704062 and 30500644) and the project of Natural Science Foundation of Hunan Province (Grant No. 2020JJ8073).
This work was supported by the Natural Science Foundation of Hunan (Grant Number 2020JJ8073) and National Natural Science Foundation of China (Grant Numbers 81673858, 81704062 and 30500644).
Ethics approval and consent to participate
The study was review and approved by the Ethics Committee of Xiangya Hospital, Central South University.
Consent for publication
The authors have declared that no competing interest exists.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sequences of the Primers and Product size in this Study. Table S2. Primary antibody information. Table S3. Pathways enrichment analysis of differently expressed genes in SHR of PGQYD-treated rats.
About this article
Cite this article
Yao, J., Zhang, C., Yang, Y. et al. Comparative transcriptomic analysis revealed novel potential therapeutic targets of traditional Chinese medicine (Pinggan-Qianyang decoction) on vascular remodeling in spontaneously hypertensive rats. Chin Med 16, 21 (2021). https://doi.org/10.1186/s13020-021-00431-4
- Traditional Chinese medicine
- Comparative transcriptomic analysis
- Vascular remodeling
- Spontaneously hypertensive rats
- AGE–RAGE signal pathway
- Cell autophagy