Identification of differentially expressed non‑coding RNAs and mRNAs involved in Qi stagnation and blood stasis syndrome
- Authors:
- Published online on: December 7, 2018 https://doi.org/10.3892/etm.2018.7068
- Pages: 1206-1223
-
Copyright: © Chen et al. This is an open access article distributed under the terms of Creative Commons Attribution License.
Abstract
Introduction
It is a common phenomenon that the symptoms and signs are slightly different among patients with the same disease (1). For instance, coronary heart disease (CHD) manifests as chest pain (2). However, certain patients suffer from a stabbing pain, while others have complaints of a burning pain (3). Within the framework of Traditional Chinese Medicine (TCM) theory, practitioners always prescribe a medical formula based on symptoms and signs (4). In fact, TCM practitioners may divide patients with the same disease into several subtypes according to these slight differences and prescribe the corresponding formulae. In TCM theory, these subtypes are known as Zheng (5). It is renowned that the advancement of Zheng depends on the genes inherited from the parents, the lifestyle and the environmental conditions, which may differ between individual patients, and therefore, different types of Zheng exist (6). The first rule of Zheng is that one type of disease may be divided into several different Zheng (7).
A disequilibrium of the body's original condition is accountable for the development of Zheng, and certain deviations from the equilibrium may lead to certain types of disease and manifestations thereof (8). For instance, Qi stagnation and blood stasis syndrome (QSBSS), a type of Zheng, describes the condition of unsmooth flow of Qi and blood, which may occur in diseases including CHD, chronic gastritis (CG) and rheumatoid arthritis (RA) (9). This exemplifies the second rule of Zheng, namely that one type of Zheng may be the basis of several different types of disease (7). However, the western scientific basis for this theory, including the substantial essence of this Zheng, remains to be fully elucidated. Previous studies in the field mainly focus on the subtypes of one disease based on Zheng, while those on the biological basis of the same Zheng in different diseases are scarce (10,11).
The genetic profile of an individual provides information on what could theoretically happen in the body, while protein expression profiles are representative of the processes that are active in the body at the time-point of determination (12). Of note, RNAs are the connection between DNAs and proteins, which reflect the interaction of the genotype and the environment (8). Among the different types of RNAs, messenger RNAs (mRNAs) are transcripts of DNA and the information they carry is transferred by their translation into proteins in accordance with the central dogma (13). Long non-coding RNAs (lncRNAs) are a class of non-protein-coding RNAs of >200 nucleotides in length, and they are known to have roles in controlling chromatin structure, transcriptional regulation and post-transcriptional processing (14,15). Furthermore, circular RNAs (circRNAs), which have drawn an increasing amount of attention recently, are a type of non-canonical form of alternative splicing and are more stable than linear RNAs (16,17). Network pharmacology studies based on the theory of systems biology may be used for the preliminary validation of the interaction between TCM herbs and RNA targets (18). Due to its core methods based on the holistic view and the dynamic balance of the human body, network pharmacology is useful for the study of TCM herbs with complex mechanisms of action (19). In the present study, a network pharmacology analysis was employed to preliminarily validate the differentially expressed (DE) RNAs in subjects with QSBSS. The purpose of the present study was to explore the biological basis of one common Zheng, QSBSS, from the perspective of ncRNAs and mRNAs in a diversity of associated diseases.
Materials and methods
Patients and samples
The protocol of the present study was approved by the Ethics Committee of Guang'anmen Hospital (Beijing, China). The study included 5 patients with CHD and QSBSS (group a), 5 patients with CG and QSBSS (group b), 5 patients with RA and QSBSS (group c) and 5 completely healthy subjects (control) who presented at Guang'anmen Hospital affiliated to the China Academy of Chinese Medical Sciences (Beijing, China) between October 2016 and December 2016. Informed consent was provided by all of the participants prior to the study.
The diagnosis of QSBSS was based on the judgement of an expert clinician according to the standard criteria approved by the China Food and Drug Administration (20) and the scaled diagnosis criteria, which is used to diagnose Qi Stagnation and Blood Stasis Syndrome (21). The criteria for QSBSS included primary symptoms and signs [i) distending pain; ii) tingling sensation anywhere on the body; iii) pain in a fixed place; iv) pain aggravated by touch; v) lumps in the body; vi) self-reported irritability; vii) clinical depression diagnosed by a psychiatrist], secondary symptoms and signs [i) loss of appetite; ii) paleness of the face; iii) hemorrhagic spots on any part of the body], TCM tongue presentations [i) purple, dusky tongue; ii) stasis macules or stasis spots on the tongue] and pulse presentation [i) choppy pulse; ii) wiry pulse] (22). A patient was diagnosed with QSBSS if they had >2 of the primary symptoms and signs, >1 of the secondary symptoms and signs, and 1 of the tongue or pulse presentations. The diagnosis of CHD, RA and CG were based on guidelines for the diagnosis and management of patients with stable ischaemic heart disease (23), the diagnosis and treatment for rheumatoid arthritis (24) and the consensus of chronic gastritis (25), respectively. From each of the patients, blood samples (3–4 ml) were collected in EDTA-containing tubes.
RNA isolation and library preparation
RNA degradation and contamination was monitored on 1% agarose gels. RNA purity was confirmed using the Agencourt AMPure XP (cat. no. A63881; Beckman Coulter, Brea, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit for the Bioanalyzer 2100 system (cat. no. 5067-1511; Agilent Technologies, Inc., Santa Clara, CA, USA). The methods of lncRNA quantification were identical to the conventional methods for mRNAs. Quantification of circRNAs was performed with exonuclease to exclude non-circRNAs. After RNA extraction, the total RNA was digested with TruSeq Stranded Total RNA with Ribo-Zero Gold (cat. no. 15021048; Illumina, Inc., San Diego, CA, USA). The sample was then subjected to reverse transcription using SuperScript II Reverse Transcriptase (cat. no. 18064014; Invitrogen; Thermo Fisher Scientific, Inc., Waltham, MA, USA).
A total of 3 µg RNA per sample was used as input material for the RNA sample preparations. First, ribosomal (r)RNA was removed using Agencourt RNAClean XP (cat. no. A63987; Beckman Coulter), and rRNA-free material was cleaned up by ethanol precipitation. Subsequently, sequencing libraries were generated using the rRNA-depleted RNA. Library preparation was performed according to manufacturer's protocol (Illumina). In brief, fragmentation was performed using divalent cations at 94°C in 5 µl rRNA Binding Buffer and 5 µl rRNA Removal Mix-Gold. First-strand complementary (c)DNA was synthesized using 8 µl First Strand Synthesis Act D Mix and SuperScript II Reverse Transcriptase. Second-strand cDNA synthesis was performed using 5 µl End Repair Control (2 µl End Repair Control + 98 µl Resuspension Buffer). The reaction buffer contained deoxynucleoside triphosphate where thymine was replaced by uracil. After adenylation of the 3′ends of the DNA fragments, an adaptor with a hairpin loop structure was ligated to prepare for hybridization. To select cDNA fragments of preferentially 150–200 bp in length, the library fragments were purified using the AMPure XP system (cat. no. A63881; Beckman Coulter). The DNA fragment was enriched by using 5 µl polymerase chain reaction (PCR) Primer Cocktail and 25 µl PCR Master Mix prior to performing PCR. The mixture was incubated in the following conditions: 98°C for 30 sec, 15 cycles 98°C for 10 sec, 60°C for 30 sec, 72°C for 30 sec, then 72°C for 5 min and held at 10°C. Finally, products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system.
Sequencing of RNAs and DE RNAs analysis
The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3cBot-HS (Illumina) according to the manufacturer's protocols. After cluster generation, the libraries were sequenced on an Illumina HiSeq 2,500 platform, and 125 bp paired-end and 50 bp single-end reads were generated. The number of reads for each sample was 8×107. The number of clean bases of lncRNAs, circRNAs and mRNAs for each sample in filtered data was 12–13 Gb. The error rate of each base in lncRNAs, circRNAs and mRNAs for each sample was <0.001. The spliced transcription products of each sample were combined and screened as lncRNAs with Cuffmerge program in the Cufflinks 2.2.1 software package. All transcripts that overlapped with known mRNAs, other non-coding RNA and non-lncRNA were discarded. Next, the transcripts that longer than 200 bp and >2 exons were obtained, and CPC (version 0.9-r2) (26), PLEK (version 1.2) (27), CNCI (version 1.0) (28) and Pfam (version 30) (29) software packages were used to predict transcripts with coding potential. The characteristics (including length, type and number of exons) of lncRNA were analysed after screening by the Cuffcompare program in the Cufflinks 2.2.1 software package (30). CircRNAs were identified using the CircRNA Identifier tool (31). The expression of lncRNAs were calculated using fragments per kb per million reads (FPKM), while the expression of circRNAs were calculated by spliced reads per million reads (RPM). Then, DESeq2 with negative binomial distribution was used to analyse DE of RNAs (32). All sequencing processes and analyses were performed by Shanghai OE Biotech Co., Ltd. (Shanghai, China).
Gene ontology (GO) annotations and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis
GO annotations and KEGG pathway analysis were performed to investigate the possible roles of the DE ncRNAs. GO annotations were performed to identify regulatory networks of DE genes (http://geneontology.org). KEGG analysis was also performed to explore the enriched pathways of the DE RNAs based on the KEGG database (http://www.genome.jp/kegg/). To reveal the role and interactions among the DE ncRNAs and mRNAs in QSBSS patients, an ncRNAs regulatory network was constructed using Cytoscape software version 3.2.1 (http://www.cytoscape.org/).
Network pharmacology validation
The process of network pharmacology validation was performed in three steps. First, according to the sequencing results of DE RNAs, the RNA targets associated with QSBSS were prepared by ID normalization at http://www.uniprot.org/ (33). Next, the active components of Xuefu Zhuyu decoction, a classic and CFDA-approved TCM formula for treating QSBSS, were retrieved from the TCM Systems Pharmacology Database (TCMSP) and the TCM Integrated Database (34,35). The compounds were then screened according to parameters including drug-likeness (DL) and oral bioavailability (OB), and the specific ingredients were selected if their DL ≥0.18 and OB ≥30%, suggested criteria given by the TCMSP database (34). Finally, the targets of the compounds were predicted using the TCMSP analysis platform and the merged compound-RNA target network was constructed using Cytoscape (version 3.2.1), in order to determine whether any of the DE RNAs associated with QSBSS were targets of the components of Xuefu Zhuyu decoction.
Statistical analysis
The data regarding the clinical characteristics of the participants are expressed as the mean ± standard error of the mean and the corresponding results were statistically analyzed using one-way analysis of variance. The results of DE RNAs were analyzed by the algorithm of DESeq2 in R software (32). The regulatory network for ncRNAs was analyzed by the algorithm of ClueGo in Cytoscape software. P<0.05 was considered to indicate a statistically significant difference.
Results
Clinical characteristics of participants
A total of 20 participants were recruited for the present study and were divided into 4 groups. The clinical characteristics of these 4 groups were listed in Table I. The 4 groups were matched in terms of sex, as well as body height and weight, and there were no significant differences in terms of the frequency of hypertension, diabetes mellitus and cancer between the groups (P>0.05).
DE ncRNAs and mRNAs, and their overlap among patient groups with different types of disease
The results of the DE ncRNAs and mRNAs in the comparisons were as follows: A total of 678 DE lncRNAs (448 upregulated and 230 downregulated), 166 DE circRNAs (59 upregulated and 107 downregulated), and 2,676 DE mRNAs (1,287 upregulated and 1,389 downregulated) were identified in group a (CHD with QSBSS vs. controls), respectively. In group b (CG with QSBSS vs. control), 1,158 DE lncRNAs (901 upregulated and 257 downregulated), 116 DE circRNAs (38 upregulated and 78 downregulated) and 2,880 DE mRNAs (1,336 upregulated and 1,544 downregulated) were detected. A total of 697 DE lncRNAs (470 upregulated and 227 downregulated), 120 DE circRNAs (49 upregulated and 71 downregulated) and 2,776 DE mRNAs (1,420 upregulated and 1,356 downregulated) were identified in group c (RA with QSBSS vs. control). DE ncRNAs (Figs. 1 and 2) and mRNAs (Fig. 3) in the samples were displayed using Volcano plot and clustering maps for groups a, b and c individually. The top 10 upregulated and top 10 downregulated lncRNAs, circRNAs and mRNAs in groups a, b and c are listed in Tables II, III and IV, respectively. Overlapped DE lncRNAs, circRNAs and mRNAs among groups a, b and c were illustrated in a Venn diagram (Fig. 4). A total of 374 up-regulated DE mRNAs, 323 down-regulated DE mRNAs, 81 up-regulated DE lncRNAs, 23 down-regulated DE lncRNAs and 2 down-regulated DE circRNAs overlapped. No up-regulated DE circRNAs overlapped.
Functional prediction of ncRNAs
To further determine the possible functions of the overlapped DE RNAs and investigate the association between the functional groups and their underlying mechanisms in the biological networks, the Cytoscape plugin ClueGo was utilized. Significant GO terms in the categories MF and BP were specifically associated with nucleic acid metabolic process, regulation of macromolecule metabolic process, peptidyl-lysine modification and mitotic cell cycle (Figs. 5 and 6). Of note, these DE RNAs coordinate with each other to perform their functions and the relevant major biochemical pathways and signal transduction pathways were predicted by enrichment analysis based on the KEGG database. A scatterplot of the DE RNAs was employed to demonstrate the KEGG enrichment analysis. In this plot, the degree of KEGG enrichment is assessed by the enrichment factor, P-value and number of genes. A greater enrichment factor is associated with a P-value closer to zero and a greater number of genes, and therefore a more significantly enriched pathway. When the data were analyzed regarding the co-expression of genes of DE lncRNAs, DE mRNAs and target genes of DE circRNAs, the most significantly enriched KEGG pathways were sphingolipid, neurotrophin, AMP-activated protein kinase (AMPK), mRNA surveillance pathway and endocytosis (Fig. 7). These major BPs and signal transduction pathways provide insight regarding further directions of research on the biological basis of QSBSS.
Merged network of DE RNAs and target fishing
As presented in the network for compounds and targets (Fig. 8) and the merged network (Fig. 9), the targets of active components of Xuefu Zhuyu decoction partly covered the targets in the DE RNAs of QSBSS discovered by sequencing. The shared targets of QSBSS and Xuefu Zhuyu Decoction were interleukin 1β, DNA topoisomerase I, potassium voltage-gated channel subfamily H member 2, nuclear receptor coactivator 1, mitogen-activated protein kinase 8, insulin-like growth factor 1 receptor and calcium/calmodulin-dependent protein kinase kinase 2.
Discussion
In the present study, expression profiles of lncRNAs, circRNAs and mRNA were obtained by next-generation sequencing and compared in patients with QSBSS and three different types of disease (CHD, CG or RA vs. normal controls). A total of 104 lncRNAs, 2 circRNAs and 697 mRNAs were finally identified to be significantly DE among groups a, b and c in a Venn diagram. These DE RNAs were subsequently displayed in a Volcano Plot and integrated into hierarchical categories according to heat maps and hierarchical clustering. GO and KEGG analyses were used to predict their relevant functions by ClueGo. A network pharmacology study was also performed to preliminarily confirm the DE RNAs with the hypothesis that the DE RNAs of QSBSS may be regulated by the TCM herbs with efficacy in treating QSBSS. These DE RNAs and their corresponding pathways may be considered highly relevant to QSBSS, and may contribute to the explanation of Zheng and prediction of therapeutic targets of the TCM formula.
QSBSS is a common Zheng that is thought to be associated with the manifestation of numerous types of disease, including CHD, CG and RA, and it is always used to define the condition of unsmooth flow of Qi and blood (36). In addition to distending pain in a fixed place, one symptom of QSBSS is emotional problems, including irritability and depression (21). The development of QSBSS is complex (37) and the mechanism remains largely elusive. The competing endogenous (ce)RNA hypothesis posits that miRNA may bind to an argonaute protein and form the miRNA-guided RNA-induced silencing complex (RISC) to target the relevant mRNA, leading to accelerated degradation, blocked translation and decreased expression, while ceRNA with shared target sites competes for miRNA targeting and binding in the RISC, thus sequestering the miRNA activity, resulting in repression of the expression the target gene of the respective miRNA (14). Furthermore, lncRNA, circRNA and pseudogene transcripts have gained considerable attention as they may work together to regulate ceRNA mechanisms. Despite the growing interest in exploring the biological basis of Zheng from DNA and RNA (38–41), the molecular mechanisms underlying the interaction of ncRNAs and mRNAs in QSBSS have remained largely elusive. Therefore, for the first time, the present study provided and compared expression profiles of lncRNAs, circRNAs and mRNAs in QSBSS obtained by next-generation sequencing. These pioneering discoveries may deepen the current understanding of the mechanisms of Zheng as well as provide a foundation for further scientific research into the molecular basis of TCM theory.
Previous studies have indicated the links between QSBSS and the processes of inflammation, immune response and sympathetic regulation (42–44). Of note, glycoprotein containing disulfide bonds, G protein-coupled receptor signaling pathways and catecholamine neurotransmitter activity were identified to constitute the biological basis of QSBSS through methods of data mining, which integrated the results of nearly all of the previous research on QSBSS (45). Another previous study suggested that QSBSS patients suffered from an imbalance of lower pneumogastric nerve activity and higher sympathetic nerve activity, as indicated by a heart rate variability analysis (43). Although evidence for the biological basis of QSBSS has accumulated, no comprehensive analysis of the ncRNA and mRNA profile in QSBSS has been performed to date, to the best of our knowledge. Thus, next-generation sequencing was applied in the present study to analyze DE ncRNAs and mRNAs in patients with QSBSS.
Bioinformatics analysis of the RNA sequencing results of the present study indicated that DE RNAs in QSBSS were most significantly enriched in pathways including the sphingolipid signaling pathway, the neurotrophin signaling pathway, AMPK and endocytosis. These results imply that QSBSS is not only associated with an imbalance of neural function, but also with a deregulation of energy metabolism. Specifically, sphingomyelin and its metabolic products are known to have second messenger functions in a variety of cellular signaling pathways (46,47). Among its metabolic products, ceramide and S1P regulate cellular responses to stress, with generally opposing effects (48,49). Neurotrophins are a family of trophic factors involved in the differentiation and survival of neural cells, which are central for the development of diseases including hereditary sensory and autonomic neuropathy, as well as congenital pain insensitivity with anhidrosis (50). Endocytosis is a mechanism for cells to remove ligands, nutrients, plasma membrane proteins and lipids from the cell surface, which is essential for the advance of hereditary sensory and autonomic neuropathy (51). AMPK is a serine threonine kinase that is highly conserved through evolution and the AMPK system acts as a sensor of the cellular energy status (52,53). These functions of enriched pathways may be indicative of why QSBSS frequently manifests as pain accompanied with obscure sensations, including a distending and tingling sensation, as well as emotional problems, including irritability and depression (17). Of note, QSBSS is a stagnation of Qi and the flow of blood, which is known to carry the substance, including oxygen and energy in the body (54). Hence, it is likely that the biological basis of QSBSS may be directly associated with a dysfunction of energy metabolism and the autonomic nervous system, which requires more in-depth study in the future.
To the best of our knowledge, the present study was the first to explore the biological basis of Zheng in subjects with different diseases rather than different Zheng in one disease and it was attempted to interpret QSBSS, one of the common Zheng, from the perspective of ncRNAs and mRNAs, which were determined by next-generation sequencing technology. It is worth mentioning that the grouping design in the present study may be implemented to assess any type of Zheng in further studies. However, one of the limitations of the present study is the lack of validation of DE RNAs by quantitative (q)PCR as the cohort size was very small and qPCR analysis was not performed, which is likely to be resolved in a future study. Besides, the inter-group and inter-patient variability should also be analysed in the future study. Furthermore, one intended utility of these DE RNAs and enriched pathways is to predict the potential targets of TCM formulae that treating QSBSS, which mainly include interleukin 1β, DNA topoisomerase I and potassium voltage-gated channel subfamily H member 2. These putative targets may be candidates for experimental validation of the mechanisms of action of the TCM formulae, which may be essential for the optimization of ancient formulae and the discovery of novel drugs.
Acknowledgements
Not applicable.
Funding
This work was supported by the China Food and Drug Administration Project (grant no. 201207009).
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Authors' contributions
GC designed the present study, and analyzed and interpreted the data. JLG, HQH, CL, YML, JL and JW enrolled the participants and examination of the blood samples. GC was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
Ethical approval and consent to participate
The protocol of the present study was approved by the Ethics Committee of Guang'anmen Hospital (Beijing, China). Informed consent was provided by all of the participants.
Patient consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
References
Yuan R and Lin Y: Traditional Chinese medicine: An approach to scientific proof and clinical validation. Pharmacol Ther. 86:191–198. 2000. View Article : Google Scholar : PubMed/NCBI | |
Balta S, Cakar M, Demırkol S, Kucuk U, Ay SA and Unlu M: May chest pain describe coronary heart disease? Croat Med J. 54:4112013. View Article : Google Scholar : PubMed/NCBI | |
Wang C, Niimi M, Watanabe T, Wang Y, Liang J and Fan J: Treatment of atherosclerosis by traditional Chinese medicine: Questions and quandaries. Atherosclerosis. 277:136–144. 2018. View Article : Google Scholar : PubMed/NCBI | |
Zha LH, He LS, Lian FM, Zhen Z, Ji HY, Xu LP and Tong XL: Clinical strategy for optimal traditional chinese medicine (TCM) herbal dose selection in disease therapeutics: Expert consensus on classic TCM herbal formula dose conversion. Am J Chin Med. 43:1515–1524. 2015. View Article : Google Scholar : PubMed/NCBI | |
Wang YH and Xu AL: Zheng: A systems biology approach to diagnosis and treatments. Sci Art Sci Traditional Med (Suppl). S13–S15. 2014. | |
Blázovics A: The interpretation and integration of traditional Chinese phytotherapy into Western-type medicine with the possession of knowledge of the human genome. Orv Hetil. 159:696–702. 2018.(In Hungarian). View Article : Google Scholar : PubMed/NCBI | |
Josephine PB: A global scientific challenge: Learning the right lessons from ancient healing practices. Sci Art Scie Traditional Med (Suppl). S7–S9. 2014. | |
Chen G, Xue WD and Zhu J: Full genetic analysis for genome-wide association study of Fangji: A powerful approach for effectively dissecting the molecular architecture of personalized traditional Chinese medicine. Acta Pharmacol Sin. 39:906–911. 2018. View Article : Google Scholar : PubMed/NCBI | |
Gao JL, Chen G, He QY, Li J and Wang J: Analysis of chinese patent medicine prescriptions for qi stagnation and blood stasis syndrome. Zhongguo Zhong Yao Za Zhi. 42:187–191. 2017.(In Chinese). PubMed/NCBI | |
Li Y, Li R, Ouyang Z and Li S: Herb Network Analysis for a Famous TCM Doctor's prescriptions on treatment of rheumatoid arthritis. Evid Based Complement Alternat Med. 2015:4513192015.PubMed/NCBI | |
Liu C, Chen G, Liu D, Nie S, Wang X, Huang X, Li H, Li Y, Luo S, Zhao G, et al: Observational study of the association between TCM zheng and types of coronary artery stenosis: Protocol of a multicenter case series study. Evid Based Complement Alternat Med. 2018:25649142018. View Article : Google Scholar : PubMed/NCBI | |
Guttman M and Rinn JL: Modular regulatory principles of large non-coding RNAs. Nature. 482:339–346. 2012. View Article : Google Scholar : PubMed/NCBI | |
Gilbert WV, Bell TA and Schaening C: Messenger RNA modifications: Form, distribution, and function. Science. 352:1408–1412. 2016. View Article : Google Scholar : PubMed/NCBI | |
Salmena L, Poliseno L, Tay Y, Kats L and Pandolfi PP: A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell. 146:353–358. 2011. View Article : Google Scholar : PubMed/NCBI | |
Thomson DW and Dinger ME: Endogenous microRNA sponges: Evidence and controversy. Nat Rev Genet. 17:272–283. 2016. View Article : Google Scholar : PubMed/NCBI | |
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al: Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 495:333–338. 2013. View Article : Google Scholar : PubMed/NCBI | |
Liu J, Liu T, Wang X and He A: Circles reshaping the RNA world: From waste to treasure. Mol Cancer. 16:582017. View Article : Google Scholar : PubMed/NCBI | |
Hopkins AL: Network pharmacology: The next paradigm in drug discovery. Nat Chem Biol. 4:682–690. 2008. View Article : Google Scholar : PubMed/NCBI | |
Li S and Zheng B: Traditional Chinese medicine network pharmacology: Theory, methodology and application. Chin J Nat Med. 11:110–120. 2013. View Article : Google Scholar : PubMed/NCBI | |
Zheng XY: Guiding Principles for Clinical Research of New Drugs in Traditional Chinese Medicine (Trial) 2002. 2. 1. China Med Sci Technol Press; BeiJing, BJ: pp. 69–126. 2002 | |
Wang J, Gao JL, Chen G and Haoqiang H: Diagnosis criteria on Qi Stagnation and Blood Stasis Syndrome. Chin J Experimental Traditional Med Formulae. 24:16–20. 2018.(In Chinese). | |
Shen CT, Zhang L, Wang Z, Chen QG, Chen BW, Yuan Y, Zhong Y, Zhu J, Chen Y and Wang YY: Standard for symptoms and signs of Traditional Chinese Medicine in new drug clinical trials. J Chin Med. 54:1265–1267. 2013. | |
American College of Cardiology Foundation/American Heart Association: 2012 ACCF/AHA/ACP/AATS/PCNA/SCAI/STS guideline for the diagnosis and management of patients with stable ischemic heart disease. J Am College Cardiololy. 60:2564–603. 2012. View Article : Google Scholar | |
Rheumatology Branch of Chinese Medical Association: The diagnosis and treatment for rheumatoid arthritis. Chin J Rheumatol. 14:265–269. 2010.(In Chinese). | |
Fang JY, Liu WZ, Li ZS, Du YQ, Ji XL, Ge ZZ, Li YQ, Ni JM, Lv NH, Wu KC, et al: Digestive Disease Branch of Chinese Medical Association: Consensus on chronic gastritis in China (Shanghai 2012). Zhong Guo Yi Xue Qian Yan Za Zhi. 5:44–55. 2013. | |
Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L and Gao G: CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 35:W345–W349. 2007. View Article : Google Scholar : PubMed/NCBI | |
Li A, Zhang J and Zhou Z: PLEK: A tool for predicting long non-coding RNAs and messenger RNAs based on an improved k-mer scheme. BMC Bioinformatics. 15:3112014. View Article : Google Scholar : PubMed/NCBI | |
Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C and Zhao Y: Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 41:e1662013. View Article : Google Scholar : PubMed/NCBI | |
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al: Pfam: The protein families database. Nucleic Acids Res. 42:(Database Issue). D222–D230. 2014. View Article : Google Scholar : PubMed/NCBI | |
Ghosh S and Chan CK: Analysis of RNA-Seq data using TopHat and Cufflinks. Methods Mol Biol. 1374:339–361. 2016. View Article : Google Scholar : PubMed/NCBI | |
Gao Y, Wang J and Zhao F: CIRI: An efficient and unbiased algorithm for de novo circular RNA identification. Genome Biol. 16:42015. View Article : Google Scholar : PubMed/NCBI | |
Love MI, Huber W and Anders S: Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:5502014. View Article : Google Scholar : PubMed/NCBI | |
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, et al: UniProt: The universal protein knowledgebase. Nucleic Acids Res. 32:(Database Issue). D115–D119. 2004. View Article : Google Scholar : PubMed/NCBI | |
Ru JL, Li P, Wang J, Zhou W, Li B, Huang C, Li P, Guo ZH, Tao WY, Yang YF, et al: TCMSP: A database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 6:132014. View Article : Google Scholar : PubMed/NCBI | |
Xue R, Fang Z, Zhang M, Yi Z, Wen C and Shi T: TCMID: Traditional Chinese medicine inegrative database for herb molecular mechanism analysis. Nucleic Acids Res. 41:(Database Issue). D1089–D1095. 2013. View Article : Google Scholar : PubMed/NCBI | |
Liu L, Liu JX, Guo H and Ren JX: Recent advances on pericytes in microvascular dysfunction and traditional Chinese medicine prevention. Zhongguo Zhong Yao Za Zhi. 42:3072–3077. 2017.(In Chinese). PubMed/NCBI | |
Wang T, Jia C, Chen Y, Li X and Cheng J: Analysis on establishment and affecting factors of qi stagnation and blood stasis rat model. Zhongguo Zhong Yao Za Zhi. 37:1629–1633. 2012.(In Chinese). PubMed/NCBI | |
Chen CS, Lin LW, Hsieh CC, Chen GW, Peng WH and Hsieh MT: Differential gene expression in hemodialysis patients with ‘cold’ zheng. Am J Chin Med. 34:377–385. 2006. View Article : Google Scholar : PubMed/NCBI | |
Lu YY, Zhao Y, Song YN, Dong S, Wei B, Chen QL, Hu YY and Su SB: Serum cytokine profiling analysis for zheng differentiation in chronic hepatitis B. Chin Med. 10:242015. View Article : Google Scholar : PubMed/NCBI | |
Kiyama R: DNA microarray-based screening and characterization of traditional chinese medicine. Microarrays (Basel). 6:42017. View Article : Google Scholar | |
Hu XQ and Su SB: An overview of epigenetics in Chinese medicine researches. Chin J Integr Med. 23:714–720. 2017. View Article : Google Scholar : PubMed/NCBI | |
Li W, Liu P, Duan J, Ma CH, Shi XQ and Guo JM: Effects of Xiangfu Siwu decoction on nerve endocrine immune system in female rats with Qi Stagnation and Blood Stasis Syndrome. Chin J Exp Trad Med Formulae. 20:99–104. 2014. | |
Yin CE, Zhang NN, Lv FF, Guan L and Wei XJ: Study on the correlation of Qi Stagnation and Blood Stasis Syndrome and heart rate variability in patients with Coronary Heart Disease. J Sichuan Traditional Chin Med. 21:15–16. 2003.(In Chinese). | |
Ren JX, Liu JX and Lin CR: Comparative analysis on the biological basis of blood stasis syndrome induced by Qi Stagnation and Qi deficiency in patients with unstable angina pectoris. Zhongguo Zhong Xi Yi Jie He Za Zhi. 30:352–356. 2010.(In Chinese). PubMed/NCBI | |
Liu JW: A preliminary bionetwork research of Qi deficiency and Qi stagnation in the context of coronary heart disease angina pectoris. Beijing Univ Chin Med,. 2:29–34. 2015. | |
Kolesnick RN: Sphingomyelin and derivatives as cellular signals. Prog Lipid Res. 30:1–38. 1991. View Article : Google Scholar : PubMed/NCBI | |
Hannun YA: The sphingomyelin cycle and the second messenger function of ceramide. J Biol Chem. 269:3125–3128. 1994.PubMed/NCBI | |
Liu J, Beckman BS and Foroozesh M: A review of ceramide analogs as potential anticancer agents. Future Med Chem. 5:1405–1421. 2013. View Article : Google Scholar : PubMed/NCBI | |
Brunkhorst R, Vutukuri R and Pfeilschifter W: Fingolimod for the treatment of neurological diseases-state of play and future perspectives. Front Cell Neurosci. 8:2832014. View Article : Google Scholar : PubMed/NCBI | |
Arevalo JC and Wu SH: Neurotrophin signaling: Many exciting surprises. Cell Mol Life Sci. 63:1523–1537. 2006. View Article : Google Scholar : PubMed/NCBI | |
Grant BD and Donaldson JG: Pathways and mechanisms of endocytic recycling. Nat Rev Mol Cell Biol. 10:597–608. 2009. View Article : Google Scholar : PubMed/NCBI | |
Steinberg GR and Kemp BE: AMPK in health and disease. Physiol Rev. 89:1025–1078. 2009. View Article : Google Scholar : PubMed/NCBI | |
Sanchez AM, Candau RB, Csibi A, Pagano AF, Raibon A and Bernardi H: The role of AMP-activated protein kinase in the coordination of skeletal muscle turnover and energy homeostasis. Am J Physiol Cell Physiol. 303:C475–C485. 2012. View Article : Google Scholar : PubMed/NCBI | |
Wang YY and Zhang HM: The scientific explanation of Qi. J Chin Med. 58:811–813. 2017.(In Chinese). |