Discrete derivative: a data slicing algorithm for exploration of sharing biological networks between rheumatoid arthritis and coronary heart disease
- Guang Zheng†^{1, 2},
- Miao Jiang^{1},
- Xiaojuan He†^{1},
- Jing Zhao^{1},
- Hongtao Guo^{1, 3},
- Gao Chen^{1},
- Qinglin Zha^{1} and
- Aiping Lu†^{1}Email author
DOI: 10.1186/1756-0381-4-18
© Zheng et al; licensee BioMed Central Ltd. 2011
Received: 25 December 2010
Accepted: 23 June 2011
Published: 23 June 2011
Abstract
Background
One important concept in traditional Chinese medicine (TCM) is "treating different diseases with the same therapy". In TCM practice, some patients with Rheumatoid Arthritis (RA) and some other patients with Coronary Heart Disease (CHD) can be treated with similar therapies. This suggests that there might be something commonly existed between RA and CHD, for example, biological networks or biological basis. As the amount of biomedical data in leading databases (i.e., PubMed, SinoMed, etc.) is growing at an exponential rate, it might be possible to get something interesting and meaningful through the techniques developed in data mining.
Results
Based on the large data sets of Western medicine literature (PubMed) and traditional Chinese medicine literature (SinoMed), by applying data slicing algorithm in text mining, we retrieved some simple and meaningful networks. The Chinese herbs used in treatment of both RA and CHD, might affect the commonly existed networks between RA and CHD. This might support the TCM concept of treating different diseases with the same therapy.
Conclusions
First, the data mining results might show the positive answer that there are biological basis/networks commonly existed in both RA and CHD. Second, there are basic Chinese herbs used in the treatment of both RA and CHD. Third, these commonly existed networks might be affected by the basic Chinese herbs. Forth, discrete derivative, the data slicing algorithm is feasible in mining out useful data from literature of PubMed and SinoMed.
Background
Traditional Chinese Medicine (TCM), is one of China's splendid cultural heritages [1, 2] with various intelligent theoretical thinking. One important concept in TCM is called "Treating Different Diseases with the Same Therapy" (TDDST), which can be explained as that different diseases might be shown with similar TCM patterns based on TCM diagnostic information (such as symptoms, pulse feelings and tongue appearance). Therefore, they could be treated with similar therapies in TCM. For example, in Western medicine, Rheumatoid Arthritis (RA) and Coronary Heart Disease (CHD) are recognized as different diseases because they are different in etiology and pathology. However, in TCM pattern classification, they share similar TCM patterns during their development. RA and CHD, named as Impediment Pattern (Bi Zheng) and Palpitations/Angina Pectoris (Xin Ji/Xin Tong) respectively in TCM, could show similar TCM patterns in Qi deficiency and blood stasis based TCM pattern classification [2], thus the two diseases can be treated with similar therapy in this case (reinforcing the deficient Qi and dissolving the blood stasis), which has been practiced during the long-term TCM activities [2].
The successful clinical practise of TDDST in TCM may suggest the fact that there may be similar (or same) biological networks/basis in RA and CHD. Due to the complexity of biological networks/basis, the common features are impossible to be revealed by single experiment or research. Thus, they keep concealed hitherto. It can be supposed that some potential regularity might be discovered by integrated analysis of the global literatures on the two diseases, in both English and Chinese.
Although forming different theoretical systems, both TCM and Western Medicine are aiming at human health services. Thus the interdisciplinary research might lighten the cognition of health and diseases. Based on this, we retrieved data from both PubMed and SinoMed. In analyzing data from PubMed, we can calculate the biological networks/basis on Western Medicine commonly existed in both RA and CHD. From SinoMed, we can filter out the basic Chinese herbs used for the treatment of both RA and CHD. Therefore, we try to find out some common regularities between RA and CHD through the overall literatures in both PubMed (English literature) and SinoMed (TCM literature) by data mining technique. Progress in digital data acquisition and storage technology has resulted in the growth of huge databases. Thus, it is impossible for anyone to read them line by line, or record by record. Based on this, we turn to the technique of data mining. Data mining is the analysis of (often large) observational data sets to find unsuspected relationships and to summarize the data in novel ways that are both understandable and meaningful/useful to the data owner and users [3–12]. During the process of data mining, we explored databases of PubMed for English literature and SinoMed for TCM literature in Chinese.
Results and Discussion
After data mining, we got some well-structured networks, i.e., wheel-shaped, wheel-wheel shaped, and wheels-connected shaped, etc. What's more, these networks are not only in good shapes, they might also demonstrate something meaningful in researches.
Wheel-shaped Networks
To be more precise, take "Interleukin" confined within rectangles for example. If only for RA, besides interleukin-1, 6, and 10, we can retrieve interleukin-2, 5, 7, 8, 11, 12, 13, 15, 16, 17, 18, and 23. It is similar with CHD, apart from interleukin-1, 6, and 10, there are interleukins of interleukin-2, 8, and 18 are concerned.
At the first glance, one may wonder that why interleukin-2, 8, and 18 are not covered by the intersection of CHD and RA. The reason is: for DescriptorNames in CHD_RA, things under considerations are the co-occurrent both in CHD and RA. Interleukin-2, 8, and 18 do exist in data sets of CHD and RA. However, there is no co-occurrent DescriptorName pairs link them into each other on frequency "19". What's more, under the framework of TDDST in TCM theory, networks of this kinds or alike can suggest something new for the research both on treatment and pathogenesis of CHD and RA.
That is to say, through the calculation of data mining, we can get something which might be useful for further research on the treatment and pathogenesis of CHD and RA.
In Figure 2., among all the nodes connected with center node "Apoptosis", they can be grouped into two classes. One group is those nodes confined by rectangles. These nodes are more specific, usually indicating single substance which might be more meaningful in practice. The other group is those common nodes, tags of them are much more abstract. Each of them indicates a concept of a biological process, or a set of chemical material. Nodes within this group indicate knowledge in a much more abstract way, usually, they just indicate common sense.
Wheel-wheel Shaped Networks
In the left wheel, nodes around the center node tagged "Leukocytes" are all tagged with abstract concepts, i.e., "Blood Platelets", "Cell Migration Inhibition", "Cell Adhesion", "Cell Movement", and so on. As they cannot indicate much meaningful information, they are taken as noises and ignored.
In the right wheel, apart from abstract concepts (i.e., "Biological Markers", "Mast Cells", and "T-Lymphodytes"), we have "Interleukin"-1, 6, and 8, "Tumor Necrosis Factor - alpha", and "C-Reactive Protein" connected with center node with tag "Inflammation Mediators".
There is one interesting thing in this figure: an intersection composed of five nodes between two wheel-shaped networks centered by "Leukocytes" and "Inflammation Mediators". The intersection nodes are "Cell Adhesion Molecules", "Inflammation", "Macrophages", "Cytokines", and "Monocytes". What's more, the intersection nodes and the two centered nodes within wheel-shaped networks can also form a network. This network can be dug out from the inter-connected wheel-shaped network, which can is shown in lower part in Figure 3.
There are two interesting thing in this sub-network in Figure 3. One is the meaning of this sub-network which has been described above. The other is that we can get a third sub-network out from a network of two interconnected wheel-shaped networks.
Thus, we know that wheel-wheel shaped networks contain more information than single-wheel shaped network. They can also suggest much more in medical research.
Wheels-connected Shaped Networks
This figure can tell us one important principle: all complex networks can be taken as inter-connection of single wheel-shaped networks. Of course, there might be some sub-networks which are isolated, and have no connections with others.
There are two types of inter-connections among three wheel shaped networks who are centered by "Phenotype", "Membrane Glycoproteins", and "Genetic Predisposition to Disease" respectively. One type is those nodes connects two wheel-shaped networks, i.e., "Biological Markers" connects wheels centered by "Phenotype" and "Membrane Glycoproteins", "Tumor Necrosis Factor-alpha" and "Carrier Proteins" connect wheels centered by "Membrane Glycoproteins" and "Genetic Predisposition to Disease". "Mutation", "Base Sequence", "Linkage (Genetics)", "Haplotypes", and "Genetic Markers" connect wheels centered by "Phenotype" and "Genetic Predisposition to Disease". The other type is those who connect all the three wheel-shaped networks. For example, nodes with tag "Gene Frequency" and "Genotype" confined within a red rectangle, each of them has three edges out to the three wheel-shaped networks. Apart from the two types of inter-connecting nodes, there are also meaningful nodes tagged with single chemical materials. Each of them is confined within a rectangle, i.e., "Toll-Like Receptor 4", "CD40 Ligand", "Fas Ligand Protein", and "Tumor Necrosis Factor-alpha".
Overview of These Networks in Systems Biology
From Figure 1, 2, 3, and 4, we can find that the pathogenesis of RA and CHD are commonly associated with "inflammation", "apoptosis", "cytokines" and "macrophage" et., which also have been verified by some experimental research [13–15].
Take "inflammation" in Figure 1 for example. In the past decade, many studies reported that inflammation was the key pathogenesis in both RA and CHD [16, 17]. Some immune cells and molecules such as macrophages, neutrophils, monocytes, cytokines and so on were all involved in the occurrence and development of these two diseases. The abundance and activation of macrophages in the inflamed synovial membrane/pannus significantly correlated with the severity of RA [14]. Reactive oxygen, free radicals and lipase produced by macrophges played an important role in the development of CHD [18–20]. High levels of Interleukin-1 (IL-1), Interleukin-6 (IL-6), Interleukin-10 (IL-10) and Tumor necrosis factor-α (TNF-α) in the serum were observed in RA and CHD patients [21–23]. In Figure 1, we can also find these useful information, such as "macrophage", "neutrophils", "monocytes", "IL-1", "IL-6", "IL-10" and "TNF-α" etc. Therefore, these data demonstrate that our method to search the commonly existed biological basis/pathogenesis between RA and CHD is feasible.
Cross Query between SinoMed and PubMed
Results of data mining on RA and CHD in PubMed and SinoMed with Chinese herbs Angelica and Salvia
SinoMed | PubMed | ||
---|---|---|---|
Chinese Herb | Biology Basis | Detailed Information | |
Angelica | Inflammation | Apoptosis | Caspase 3, Caspases, bcl-2-Associated × Protein, Proto-Oncogene Proteins c-bcl-2, Transforming Growth Factor beta1, CD40 Ligand, Tumor Necrosis Factor-alpha, Tumor Suppressor Protein p53, Vascular Endothelial Growth Factor A |
Leukocytes | Cyclooxygenase 1, Phospholipases A, Prostaglandin-Endoperoxide Synthases | ||
Monocytes | Interleukin-8 | ||
Cytokines | Fibroblast Growth Factor 2, Intercellular Adhesion Molecule-1, Interleukin-2, Interleukin-6, NF-kappa B, Toll-Like Receptor 4, Transforming Growth Factor beta, Transforming Growth Factor beta1 | ||
Macrophages | Nitric Oxide Synthase Type II, Cyclooxygenase 2, Interleukin-2, NF-kappa B, Nitric Oxide Syn-thase, Tumor Necrosis Factor-alpha, Interleukin- 6, Cyclooxygenase 1, Interleukin-1, Prostaglandin- Endoperoxide Synthases | ||
Salvia | Inflammation | Apoptosis | Caspase 3, Caspases, Proto-Oncogene Proteins c-bcl-2, bcl-2-Associated × Protein, Tumor Suppressor Protein p53, L-Lactate Dehydrogenase, Proto-Oncogene Proteins c-akt, Proto-Oncogene Proteins, Tumor Necrosis Factor-alpha, Cyclooxygenase 2, Intercellular Adhesion Molecule-1, Interleukin-2, Interleukin-6, Interleukin-8, Platelet-Derived Growth Factor, Vascular Endothelial Growth Factor A |
Leukocytes | Intercellular Adhesion Molecule-1, Tumor Necrosis Factor-alpha | ||
Monocytes | Proto-Oncogene Proteins, Transforming Growth Factor beta1, Tumor Necrosis Factor-alpha | ||
Cytokines | Alanine Transaminase, Tumor Necrosis Factor- alpha, Interleukin-6, Interleukin-8, Nitric Oxide Synthase Type II, Platelet-Derived Growth Factor, Prostaglandin- Endoperoxide Synthases, Transforming Growth Factor beta1 | ||
Macrophages | Nitric Oxide Synthase Type II, Tumor Necro-sis Factor-alpha, Cyclooxygenase 2, NF-kappa B, Interleukin-1, Interleukin-6, Nitric Oxide Synthase, Transforming Growth Factor beta1, Caspases, Chemokine CCL2, E-Selectin, Intercellular Adhesion Molecule-1, Matrix Metalloproteinase 9 |
Results of data mining on RA and CHD in PubMed and SinoMed with Chinese herbs Safflower and Astragalus
SinoMed | PubMed | ||
---|---|---|---|
Chinese Herb | Biology Basis | Detailed Information | |
Safflower | Inflammation | Apoptosis | Proto-Oncogene Proteins c-bcl-2, Caspase 3, bcl-2- Associated × Protein, Caspases, Nitric Oxide Synthase, Prostaglandin-Endoperoxide Synthases, Vascular Endothelial Growth Factors |
Leukocytes | Interleukin-2, Leukotriene B4, Linoleic Acid, Throm-boxane B2, Triglycerides | ||
Monocytes | Interleukin-1, Tumor Necrosis Factor-alpha, Interleukin-6, NF-kappa B, Intercellular Adhesion Molecule-1 | ||
Cytokines | Interleukin-6, Tumor Necrosis Factor-alpha, Interleukin-1, Interleukin-10, Thromboxane B2, 6-Ketoprostaglandin F1 alpha, Cyclooxygenase 1, Cyclooxygenase 2, Leukotriene B4, NF-kappa B, Prostaglandin-Endoperoxide Synthases, Toll-Like Receptor 4 | ||
Macrophages | Eicosapentaenoic Acid, Interleukin-2, Prostaglandins E, Tumor Necrosis Factor-alpha, Interleukin-1, Leukotriene B4, Endothelial Growth Factors, Prostaglandin-Endoperoxide Synthases, Thromboxane B2, Vascular Endothelial Growth Factor A, Vascular Endothelial Growth Factors | ||
Astragalus | Inflammation | Apoptosis | Caspase 3, Tumor Necrosis Factor-alpha, Caspases, NF-kappa B, Proto-Oncogene Proteins c-bcl-2, bcl-2-Associated × Protein, Interleukin-10, Interleukin-6, Nitric Oxide Synthase Type II, Transforming Growth Factor beta1, Tumor Suppressor Protein p53 |
Leukocytes | Cyclooxygenase 1, E-Selectin, Intercellular Adhesion Molecule-1, NF-kappa B, Phospholipases A, Tumor Necrosis Factor-alpha | ||
Monocytes | Tumor Necrosis Factor-alpha | ||
Cytokines | Interleukin-2, Interleukin-10, Interleukin-6, Interleukin-1, Interleukin-8, Matrix Metalloproteinase 9, Toll-Like Receptor 4, Transforming Growth Factor beta1, Tumor Necrosis Factor-alpha, Vascular Endothelial Growth Factor A | ||
Macrophages | Interleukin-1, Interleukin-2, NF-kappa B, Toll-Like Receptor 4 |
Take "apoptosis" for example. As we know, apoptosis is the process of programmed cell death, which plays a pivotal role in tissue homoeostasis. Apoptosis disorders can lead to some serious diseases such as RA and CHD [26, 27], which is also demonstrated in our study, as shown in Figure 2. We queried "Angelica", "Salvia", "Safflower" and "Astragalus" in PubMed and find "apoptosis" again and some specific nodes such as "Caspase 3", "Caspases", "Proto-Oncogene Proteins c-bcl-2", "bcl-2-Associated × Protein", "Tumor Suppressor Protein p53", et al, which are related with apoptosis and previously emerged in Figure 2. These data indicate that the mechanism of these four Chinese herbs treating RA and CHD is partly through affecting cell apoptosis by regulating apoptosis-related signal pathway, which are consistent with some experimental results [28, 29].
- 1.
For the treatment of RA and CHD, some basic Chinese herbal medicines are used. They are Angelica, Salvia, Safflower and Astragalus. In TCM theory, Chinese herbs of Angelica, Salvia and Safflower can be grouped into the class of reinforcing blood, while Astragalus can be grouped into the class of reinforcing Qi;
- 2.
The mechanisms of these commonly used traditional Chinese herbal medicines for the treatment of both RA and CHD are mostly associated with "inflammation", "apoptosis", "cytokines", "monocytes", "macrophages", etc., which indicate that these herbs can affect the same biological networks that are commonly existed in RA and CHD;
- 3.
Data mining results demonstrate that our method is feasible, and they may support the concept of TDDST in TCM;
- 4.
These results can also give us some useful tips for the future research of these two diseases.
Conclusions
First, common biological networks/basis may exist
Through the calculation, we get partial results shown in Figure 1, 2, 3, and 4. These figures indicate the existence of common biological networks and biological basis in both RA and CHD;
Second, basic Chinese herbal medicines may affect these common biological networks/basis
Through the calculation of derivatives on different orders (i.e., primary, secondary, etc.) between PubMed and SinoMed, back and forth, we get the data that support the TCM concept of TDDST through data mining. In brief, the biological basis/networks commonly existed in RA and CHD, and they might be affected by the Chinese herbal medicines which are used in TCM therapies for both RA and CHD;
Third, our data slicing algorithm works
Through the above two items, we have reason to believe that our data slicing algorithm, together with other skills adopted in data mining, can dig out simple and meaningful information from large data sets from both PubMed and SinoMed. What's more, these results may support the concept of TDDST in TCM.
Methods
For data mining, data preparation, pretreatment, and treatment is fundamental and necessary [30–33] for the final result.
The whole process of data mining is scheduled as described below. First, mine the data retrieved from Pubmed. Through this process, we get the common biological networks/basis existing in both RA and CHD. Then, mine the data retrieved from SinoMed. By doing this, we get the basic Chinese herbal medicines (with the highest frequencies occurring in literature) which are used for the treatment of both RA and CHD. At last, we turn to PubMed again with the basic Chinese herbal medicines for verification. After analyzing the data retrieved from PubMed with these basic Chinese herbs, we will verify our hypothesis: these basic Chinese herbs might affect those biological networks/basis commonly existed in both RA and CHD. In brief, we can test whether or not these Chinese herbs used against RA and CHD in TCM can affect the biological networks/basis existed in both RA and CHD.
Data from PubMed and SinoMed
PubMed is a free database accessing the MEDLINE database of citations, abstracts and some full text articles on life sciences and biomedical topics. The United States National Library of Medicine (NLM) at the National Institutes of Health (NIH) maintains PubMed as part of the Entrez information retrieval system [34]. We retrieved data of RA and CHD from PubMed as basic material for data mining. SinoMed http://sinomed.cintcm.ac.cn/index.jsp is also a database like PubMed. The characteristic of SinoMed is that this database focus on Chinese literature in the fields of TCM, biological and medicine. We only retrieve TCM data in SinoMed for basic Chinese herbs used for the treatment of both RA and CHD.
Query and Download Data
As mentioned above, we queried the keywords of "Rheumatoid Arthritis" and "Coronary Heart Disease" in the database of PubMed http://www.ncbi.nlm.nih.gov/pubmed/ on May 9, 2010. To be more specific, we searched MeSH for "Rheumatoid Arthritis" [MeSH Major Topic] with Restrict Search to Major Topic headings only retrieved a record set of 67,049 papers. Search MeSH for "Coronary Heart Disease" [MeSH Major Topic] with Restrict Search to Major Topic headings only retrieved records of 115,757. When this query was done, we downloaded these retrieved data sets in the type of XML into local computer system for the pretreatment.
In SinoMed http://sinomed.imicams.ac.cn/index.jsp, we queried the MeSH terms of "Rheumatoid" and "Arthritis" for RA, "Coronary Heart Disease" for CHD in Chinese on Aug. 24, 2010. By querying terms of "Default: Coronary Heart Disease or Rheumatoid or Arthritis" in Chinese, we retrieved 60,967 records for CHD, and 13,686 records for RA. The difference between PubMed and SinoMed is that PubMed can download all data in one time, while in SinoMed, we can only download 500 record one time.
Pretreat Data
After retrieving data from PubMed (the same with SinoMed), we listed the data order by PMID and found that, for each PMID, there were several DescriptorNames associated with it. What's more, these DescriptorNames are also the keywords of this paper. Observed this, it is natural and intuitive to construct pairs of co-occurrent DescriptorNames from the retrieved data [35]. When the tables of co-occurrent DescriptorName pairs are constructed, it is natural to calculate the frequencies of these pairs [8, 36, 37]. When the tables of frequency pairs are built, our data slicing algorithm get its input raw material.
In order to analyze these data retrieved from PubMed and SinoMed by these given keywords, the best way is to get interesting data into the framework of structured file system, i.e., structured databases [5, 7, 9, 38]. Because only in the structured database, large amount of optimization techniques can be applied and taking effects. Then, we can get efficient processing of our algorithm on large data sets [5].
For PubMed data sets, we transferred the XML data to the structured database of Microsoft^{®} SQL^{®} 2000. In order to find new pattern and rules in data mining and text mining, we should understand our data first. The more we understand our data (especially with specified knowledge), the more accurately and efficiently we can do in further analysis [38]. For SinoMed data sets, we developed a tool to transfer its plain TXT data into Microsoft^{®} SQL^{®} 2000.
A Paper in PubMed with Its PMID and DescriptorNames
PMID | DescriptorName |
---|---|
20464912 | physiopathology |
20464912 | rehabilitation |
20464912 | Evidence-Based Medicine |
20464912 | Humans |
20464912 | Muscle Stretching Exercises |
20464912 | Physical Fitness |
20464912 | Resistance Training |
20464912 | Treatment Outcome |
Algorithm of Calculating Co-occurrent DescriptorNames
USE Table I |
---|
FOR each P MID |
k = Number_of_DescriptorName(PMID) |
j = 1 |
FOR DescriptorNames(i) (i =, 1, 2,..., k) |
DO while j ≤ k |
DescriptorNames Pair = DescriptorNames(i)+ |
DescriptorNames (j) |
j = j + 1 |
OUTPUT DescriptorName_Pair INTO |
table DM_pairs |
ENDDO |
j = 1 |
ENDFOR |
ENDFOR |
Results of Co-occurrent DescriptorName Pairs Calculated by Algorithm in Table 3 with Raw Data Listed in Table 2
DescriptorName_1 | DescriptorName_2 |
---|---|
physiopathology | rehabilitation |
physiopathology | Evidence-Based Medicine |
physiopathology | Humans |
physiopathology | Muscle Stretching Exercises |
physiopathology | Physical Fitness |
physiopathology | Resistance Training |
physiopathology | Treatment Outcome |
rehabilitation | Evidence-Based Medicine |
rehabilitation | Humans |
rehabilitation | Muscle Stretching Exercises |
rehabilitation | Physical Fitness |
rehabilitation | Resistance Training |
rehabilitation | Treatment Outcome |
Evidence-Based Medicine | Humans |
Evidence-Based Medicine | Muscle Stretching Exercises |
Evidence-Based Medicine | Physical Fitness |
Evidence-Based Medicine | Resistance Training |
Evidence-Based Medicine | Treatment Outcome |
Humans | Muscle Stretching Exercises |
Humans | Physical Fitness |
Humans | Resistance Training |
Humans | Treatment Outcome |
Muscle Stretching Exercises | Physical Fitness |
Muscle Stretching Exercises | Resistance Training |
Muscle Stretching Exercises | Treatment Outcome |
Physical Fitness | Resistance Training |
Physical Fitness | Treatment Outcome |
Resistance Training | Treatment Outcome |
By the algorithm in Table 4 we can build the table of co-occurrent DescriptorName pairs. By applying this algorithm with input of Table 3 we can get the result demonstrated in Table 5.
Based on Table 5 our algorithm of discrete derivatives on the distribution frequency can get its turn for calculations.
Design Algorithm to Treat Data Sets
In the process of data mining, we construct an data slicing algorithm called discrete derivatives. This algorithm is based on the calculation of frequency. First, it splits data into different "layers", i.e., slicing actions are executed based on the frequency distribution of DescriptorName pairs. When data slices are calculated, we construct the function of discrete derivatives f _{ dd } . This is based on the principle of derivative in advanced calculus [39] yet in discrete space. Function f _{ dd } can calculate the difference between two adjacent data slices. When discrete derivatives f _{ dd } are calculated, we will discuss one property of function f _{ dd } .
PubMed Data
Among the tables of co-occurrent DescriptorName pairs in CHD, RA and CHD_RA (intersection of CHD and RA), there might be something interesting existed. In the table of CHD, there are 7,277 lines. In RA, there are 11,952 lines, and in CHD_RA, there are 1,474 lines. That is, focused on the DescriptorName pairs within both RA and CHD, there are only 1,474 pairs left. Table CHD_RA forms the initial data of un-directed graph on the DescriptorName.
Algorithm of Calculating Frequency of Co-occurrent DescriptorName Pairs
USE table CHD_RA |
---|
k = max_line_number |
DO while k ≥ 1 |
GO top |
FOR DescriptorName_Pair(1)//The 1st pairs in CHD_RA |
COUNT its Frequency |
EndFor |
OUTPUT DescriptorName Pairs, Frequency INTO table |
CHD_RA_Frqncy |
DELETE all DescriptorName_Pair(1) from table |
CHD_RA |
k = max_line_number |
ENDDO |
When table CHD_RA_Frqncy is constructed, we can classify the data of co-occurrent DescriptorName pairs into different slices according to their frequency distributions. Thus, we get the table of data slices as CHD_RA_Frqncy _i where i = 1, 2, 3,···, 30 with 30 is the greatest frequency existed in table CHD_RA_Frqncy.
From Figure 6, we can see that there is a trend of simplification as frequency goes from 1 to 30. For example, numbers within parentheses keep descending as frequency increases.
What's more, as frequency goes up from 1 to 18, the main part of these networks are too complicate to check for naked eye. Then, it is hard for one to get useful information on the first view. However, there are some notable "wheel" shape sub-networks of existed in graphs on frequencies of 1, 2, 6, 8, 9, and 10. These wheel-shaped sub-networks might have some meaningful cue. As frequency goes from 19 to 24, some well structured sub-networks can be observed in these graphs. For example, on frequencies of 19, 20 and 21, there are distinct wheel-shape sub-networks in the left-top position of the graphs. As frequency goes beyond 25, it is clear that we can get wonderful relationships among DescriptorNames.
Check graphs on different frequencies, we can get that, apart from the fact that most of them are distinct, there are also some similar even equivalent graphs, i.e., 25 and 26.
Now, we can calculate the distribution of co-occurrent DescriptorNames pairs on each distinct frequency. We call this distribution as level distribution, because it slices original data set into different subsets according to the frequency distribution, and these distributions are calculated from co-occurrent DescriptorNames pairs. Thus, level distribution is the difference between two adjacent slices of initial data set. Formally, we have the following formal definition.
Definition 1 (Level Distribution) For given data set , let k be the hierarchical parameter, where k is the maximum frequency of occurrence and k = max(frequency(e _{ i })). For slices of data sets slice _{ i }(i = 1, 2,···, k) in , level distribution ld can be expressed by formula ld _{ i }= slice _{ i }- slice _{ i+1}where i = 1, 2,···, k - 1.
We calculate level distributions between different initial frequencies, i.e., top-left cell in Figure 7 with subscription of "1(179)" means the visualized data of level distribution "1". "1" means the data slice ID, and it represent the data between initial data slices of frequencies "1" and "2". "179" in parentheses means the number of nodes(DescriptorName) represented in level distribution "1". Figure 7 shows all level distributions of frequencies, i.e., from "1" to "30".
- 1.
Common knowledge already existed, useful but meaningless for research;
- 2.
Rarely involved knowledge, can be taken as noises or false positive findings;
- 3.
Emerging knowledge, useful and meaningful for research.
Graphs on frequencies "1" and "2" are both of this kind. As for graphs on frequencies "3"-"8" and "15", there is no specific structures exist. To get better understanding of these data, we have to check these graphs manually. Graphs on frequencies "10"-"14", "16" and "23", are clearly structured on topology, nodes in these graphs indicate explicit relationships among DescriptorNames. As for graphs on frequencies "18", "19", "21", "22", "24", and "26", nodes inside them are wheel-shaped with a centered node, all other nodes are around the centered one. These graphs indicate some knowledge to a large probability of popular science in professional eyes. Graphs on frequencies "17", "20", "25", and "27"-"29" are blank, this means that there is no content in the level of these frequencies. On "30", as there is nothing exists higher than it but empty set ∂, so, it is meaningless to calculate its content.
As we have calculated level distribution as shown in Figure 7, it is natural for us to think of calculating the discrete derivative of these level distributions for even more simple results. The discrete derivative came from the principle of derivative in advanced calculus [39], which means calculating the changing-rate/velocity of the data in a continuous space. Formally, we have .
Based on this, we give out the formal definition of discrete derivative.
Definition 2 (Discrete Derivative) For a given serious of data slices slice _{ i } with i = 1, 2,···, k where k is a natural number, discrete derivative is .
Definition 2 gives out the definition of the first order of discrete derivative on the slices of initial data set. Enlightened by the calculation of higher order derivatives in advanced calculus [39], we give out the definitions of higher order discrete derivatives.
Definition 3 (Higher Order of Discrete Derivative)
As ld' = slice _{ i }- slice _{ i+1}is the first order discrete derivative, and ld" = (ld") is the second discrete derivative. Then, higher order of discrete derivative is ld ^{(j)}= (ld ^{(j- 1)})' where j ∈ INT and j > 2.
In order to get the full view of the discrete derivatives, we calculate several orders of derivative (i.e., the first order, the second order, etc.) based on Definition 2. Most of them are similar and we show the first derivative in Figure 8.
Checking graphs in Figure 8, we can see that most of the graphs in Figure 8 are exactly the same as Figure 7, except on frequency "29". It is blanked by the definition of the algorithm on discrete derivative. It is meaningless to calculate the derivative of the last point, which is the same as frequency "30" in Figure 7. What's more, we continue our test of calculating the second, the third, and even the forth discrete derivatives, all these results demonstrate one interesting phenomenon: for a given frequency, discrete derivatives all orders are constant. That is to say, discrete derivatives of given frequency keep unchanged and independent of the orders of discrete derivative.
Now, take Figure 7 and 8 for example, it is very interesting that there are some slight differences between them. For example, pairs on frequencies of "2", "4", "9", "11", "12", "13",···. However, check them carefully, we can find that their nodes and edges are all the same, this means their topology of nodes and relationships are all the same.
Then, naturally, comes another question. Why graphs within these two figures under these frequencies are all the same? We all know that most derivative of functions change in the continuous space [39], i.e., , and so on. Of course, there are exceptions, i.e., c' = 0 where c is constant, (e ^{ x } )' = e ^{ x } and so on. These phenomena are all have their own explanations. As to the reason of our question, we believe that it is based on the essence of our data type: discrete characteristic data in PubMed. If we have all the DescriptorNames assigned with specific numbers in the continuous space, then, our discrete derivatives will have their numeric results. What's more, these results will change as those examples mentioned before.
Thus, we have our theorem of this interesting phenomenon.
Theorem 4 (Discrete Derivatives Keep Constant)
Given a data set with elements e _{ i } (i = 1, 2,···, n) where is the size of the data set. can be made hierarchical by k where k = max (frequency(e _{ i } )), frequency(e _{ i } ) is the occur times of e _{ i } within . Pairwise elements {e _{ i } , e _{ j } } stands for the co-occurrence within where i ≠ j. Then, discrete derivatives ld ^{(m)}= constant where m ≥ 1 is the order of discrete derivative.
where 1 ≤ i ≤ k - 1, and let .
By Definition 1, we have . What's more, by algorithm in Table 4 we know that e _{ i } is distinct co-occurrent DescriptorName pair. By algorithm in Table 6 we know that e _{ i } i = 1, 2,···, n are tagged with frequency frequency(e _{ i } ). Then, slice _{ i } is the collection of co-occurrent DescriptorName pairs with frequency(e _{ i } ) = 1.
In this way, we can prove that . Together with , we have hold, and the proof is done.
SinoMed Data
In order to mine out the basic Chinese herbs used in the treatment of RA and CHD in TCM, we turned to SinoMed for data sets. First, we calculated their commonly existed DescriptorNames in the same way as described in PubMed. Then, we built a table of Chinese herbs composed of 539 records as filter to filtrate the data retrieved from SinoMed. At last, we calculated the number of occurrences of Chinese herbs in the data of SinoMed. When this was done, we checked the list with TCM professionals. They found that the first four Chinese herbs listed at the top of number of occurrences are reasonable. These four Chinese herbs are "Angelica", "Salvia", "Safflower" and "Astragalus". These four Chinese herbs can be found in Figure 5. By now, we ascertain that these four Chinese herbs can be used to treat RA and CHD in TCM. However, we still do not know whether or not they have association with the knowledge hidden in PubMed which are studied by modern medical researchers. In order to verify our idea that the Chinese herbs used in both RA and CHD can affect the biological networks/basic in both RA and CHD, we explore PubMed with these four Chinese herbs.
When these data set typed XML downloaded from PubMed, we transfer them into Microsoft^{®} SQL^{®} 2000 as we did before.
Through the cross queries from SinoMed to PubMed with these four Chinese herbs mentioned above, we get meaningful data in Table 1 and 2 which can support our assumption that these basic Chinese herbs used to treat RA and CHD in TCM can actually affect the biological networks/basis that exist in RA and CHD.
Notes
Declarations
Acknowledgements
This work was partially supported by the National Eleventh Five Year Support Project of China (2006BAI04A10), the Innovative Methodology Project supported by MOST of China(2008IM020900), National Science Foundation of China (No. 30902003, 30973975, and 90709007), and Fundamental Research Funds for the Central Universities (No. lzujbky-2009-127).
Authors’ Affiliations
References
- Li Zhaoguo, Zhang Qingrong: English Textbook for Traditional Chinese Medicine. Shanghai Scientific and Technical Publishers. 2009Google Scholar
- Zhou Zhongying : Internal Medicine of Traditional Chinese Medicine. China Press of Traditional Chinese Medicine. 2007Google Scholar
- Li Jinyan, Wong Limsoon, Yang Qiang: Guest Editors' Introduction: Data Mining in Bioinformatics. IEEE Intelligent Systems. 2005, 20 (6): 16-18. 10.1109/MIS.2005.108.View ArticleGoogle Scholar
- Tzanis George, Berberidis Christos, Vlahavas Ioannis: Biological Data Mining. Encyclopedia of Database Technologies and Applications. 2005, 35-41.View ArticleGoogle Scholar
- Seifert Jeffrey: Data Mining: An Overview. CRS Report RL31798. 2004Google Scholar
- Nassif Houssam, Woodsz Ryan, Burnsideyz Elizabeth, Ayvacix Mehmet, Shavlik Jude, Page David: Information Extraction for Clinical Data Mining: A Mammography Case Study. Proceedings of the 2009 IEEE International Conference on Data Mining Workshops table of contents. 2009, 37-42.View ArticleGoogle Scholar
- Mathiak Brigitte, Eckstein Silke: Five Steps to Text Mining in Biomedical Literature. Proceedings of the Second European Workshop on Data Mining and Text Mining for Bioinformatics, held in Conjunction with ECML/PKDD in Pisa, Italy. 2004, 24: 47-50.Google Scholar
- Geng Liqiang, Hamilton Howard: Choosing the Right Lens: Finding What is Interesting in Data Mining. Studies in Computational Intelligence. 2007, Springer-Verlag Berlin Heidelberg, 43: 3-24. 10.1007/978-3-540-44918-8_1.Google Scholar
- Schmidt Sam, Vuillermin Peter, Jenner Bernard, Ren Yongli, Li Gang, Chen Phoebe Yi-Ping: Mining Medical Data: Bridging the Knowledge Divide. Proceedings of eResearch Australasia. 2008Google Scholar
- Hotho Andreas, Nürnberger Andreas, Paaß Gerhard: A Brief Survey of Text Mining. LDV Forum - GLDV Journal for Computational Linguistics and Language Technology. 2005, 20 (1): 19-62.Google Scholar
- Wu Xindong, Kumar Vipin, Ross Quinlan J: Top 10 algorithms in data mining. Knowl Inf Syst. 2007, Springer-Verlag London Limited, 14: 1C37-Google Scholar
- Heinen Stephanie, Thielen Bernhard, Schomburg Dietmar: KID- an algorithm for fast and efficient text mining used to automatically generate a database containing kinetic information of enzymes. BMC Bioinformatics. 2010, 11: 375-10.1186/1471-2105-11-375.View ArticlePubMedPubMed CentralGoogle Scholar
- Okamoto H, Hoshi D, Kiire A, Yamanaka H, Kamatani N: Molecular targets of rheumatoid arthritis. Inflamm Allergy Drug Targets. 2008, 7 (1): 53-66. 10.2174/187152808784165199.View ArticlePubMedGoogle Scholar
- Kinne RW, Brüer R, Stühlmuller B, Palombo-Kinne E, Burmester GR: Macrophages in rheumatoid arthritis. Arthritis Res. 2000, 2 (3): 189-202. 10.1186/ar86.View ArticlePubMedPubMed CentralGoogle Scholar
- Bucova M, Bernadic M, Buckingham T: C-reactive protein, cytokines and inflammation in cardiovascular diseases. Bratisl Lek Listy. 2008, 109 (8): 333-40.PubMedGoogle Scholar
- Candore G, Balistreri CR, Caruso M, Grimaldi MP, Incalcaterra E, Listi F, Vasto S, Caruso C: Pharmacogenomics: a tool to prevent and cure coronary heart disease. Curr Pharm Des. 2007, 13 (36): 3726-34. 10.2174/138161207783018617.View ArticlePubMedGoogle Scholar
- Smolen JS, Aletaha D: The assessment of disease activity in rheumatoid arthritis. Clin Exp Rheumatol. 2010, 28 (3 Suppl 59): S18-27.PubMedGoogle Scholar
- Amirbekian V, Lipinski MJ, Briley-Seabo KC: Detecting and assessing macrophages in vivo to evaluate atherosclerosis noninvasively using molecular MRI. Proc Natl Acad Sci USA. 2007, 104 (3): 961-966. 10.1073/pnas.0606281104.View ArticlePubMedPubMed CentralGoogle Scholar
- Tacke F, Alvarez D, Kaplan TJ: Monocyte subsets differentially employ CCR2, CCR5, and CX3CR1 to accumulate within atherosclerotic plaques. J Clin Invest. 2007, 117 (1):Google Scholar
- Moore KJ, Freeman MW: Scavenger receptors in atherosclerosisbeyond lipid uptake. Arterioscler Thromb Vasc Biol. 2006, 26 (8): 1702-1711. 10.1161/01.ATV.0000229218.97976.43.View ArticlePubMedGoogle Scholar
- Hussein MR, Fathi NA, El-Din AM, Hassan HI, Abdullah F, Al-Hakeem E, Backer EA: Alterations of the CD4(+), CD8 (+) T cell subsets, interleukins-1beta, IL-10, IL-17, tumor necrosis factor-alpha and soluble intercellular adhesion molecule-1 in rheumatoid arthritis and osteoarthritis: preliminary observations. Pathol Oncol Res. 2008, 14 (3): 321-8. 10.1007/s12253-008-9016-1.View ArticlePubMedGoogle Scholar
- Gotsman I, Stabholz A, Planer D, Pugatsch T, Lapidus L, Novikov Y, Masrawa S, Soskolne A, Lotan C: Serum cytokine tumor necrosis factor-alpha and interleukin-6 associated with the severity of coronary artery disease: indicators of an active inflammatory burden?. Isr Med Assoc J. 2007, 10 (7): 494-8.Google Scholar
- Danesh J, Kaptoge S, Mann AG, Sarwar N, Wood A, Angleman SB, Wensley F, Higgins JP, Len-non L, Eiriksdottir G, Rumley A, Whincup PH, Lowe GD, Gudnason V: Long-term interleukin-6 levels and subsequent risk of coronary heart disease: two new prospective studies and a systematic review. PLoS Med. 2007, 5 (4): e78-View ArticleGoogle Scholar
- WHO: WHO International Standard Terminologies on Traditional Medicine in the Western Pacific Region. WHO Library Cataloguing in Publication Data. 2007, (ISBN 978 92 9061 248 7)Google Scholar
- Xu Qinghua: Zhong Yi Yao Dui Da Quan. China Press of Traditional Chinese Medicine. 1996Google Scholar
- Liu H, Pope RM: Apoptosis in rheumatoid arthritis: friend or foe. Rheum Dis Clin North Am. 2004, 30 (3): 603-25. 10.1016/j.rdc.2004.04.010.View ArticlePubMedGoogle Scholar
- Xu F, Sun Y, Chen Y, Sun Y, Li R, Liu C, Zhang C, Wang R, Zhang Y: Endothelial cell apoptosis is responsible for the formation of coronary thrombotic atherosclerotic plaques. Tohoku J Exp Med. 2009, 218 (1): 25-33. 10.1620/tjem.218.25.View ArticlePubMedGoogle Scholar
- Han SY, Li HX, Ma X, Zhang K, Ma ZZ, Tu PF: Protective effects of purified Safflower extract on myocardial ischemia in vivo and in vitro. Phytomedicine. 2009, 16 (8): 694-702. 10.1016/j.phymed.2009.02.019.View ArticlePubMedGoogle Scholar
- Zhang JG, Yang N, He H: Effect of astraglus injection on serum apoptosis relevant factors in patients with chronic heart failure. Zhongguo Zhong Xi Yi Jie He Za Zhi. 2005, 25 (5): 400-3.PubMedGoogle Scholar
- Razali Mahir Ahmad, Ali Shahriyah: Generating Treatment Plan in Medicine: A Data Mining Approach. Shahriyah Ali Journal: American Journal of Applied Sciences. 2009, 6 (2): 345-351.Google Scholar
- Dubey , Machchhar Avinava, Bhattacharyya Jinesh, Chakrabarti Chiranjib, Soumen : Conditional Models for Non-smooth Ranking Loss Functions. ICDM '09: Proceedings of the 2009 Ninth IEEE International Conference on Data Mining. 2009, 129-138. isbn: 978-0-7695-3895-2View ArticleGoogle Scholar
- Achuthan NR, Gopalan Raj, Rudra Amit: Mining Value-Based Item Packages C An Integer Programming Approach. Data Mining. Edited by: Williams GJ, Simoff SJ. 2006, LNAI 3755, Springer-Verlag Berlin Heidelberg, 78-89.View ArticleGoogle Scholar
- Wang Xingwen, Huang Zhexue Joshua: A Case-Based Data Mining Platform. Data Mining. Edited by: Williams GJ, Simoff SJ. 2006, Springer-Verlag Berlin Heidelberg, LNAI 3755, 28 C 38Google Scholar
- Wikipedia: PubMed from Wikipedia. 2010, [http://en.wikipedia.org/wiki/PubMed]Google Scholar
- Campagna Andrea, Pagh Rasmus: Finding Associations and Computing Similarity via Biased Pair Sampling. 2009 Ninth IEEE International Conference on Data Mining. 2009, 61-70.View ArticleGoogle Scholar
- Abdel Badeeh M, Salem MR, Mahmoud SA: Mining Pattern Data Based on Rough Set Theory to Determine Thrombosis Disease. The International Journal of Artificial Intelligence and Machine Learning. 2005, 05 (I): 27-31.Google Scholar
- Boyen Peter, Neven Frank, Van Dyck Dries: SLIDER: Mining Correlated Motifs in Protein-Protein Interaction Networks. ICDM '09: Proceedings of the 2009 Ninth IEEE International Conference on Data Mining. 2009, 716-721. isbn: 978-0-7695-3895-2View ArticleGoogle Scholar
- Hand DJ, Mannila Heikki, Smyth Padhraic: Principles of data mining, ISBN 026208290X, 9780262082907. 2001, MIT PressGoogle Scholar
- Wrede Robert, Spiegel Murray: Schaum's outline of theory and problems of advanced calculus. The McGraw-Hill Companies. 2002Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.