- Open Access
Understanding the molecular aspects of oriental obesity pattern differentiation using DNA microarray
Journal of Translational Medicine volume 13, Article number: 331 (2015)
Human constitution, the fundamental basis of oriental medicine, is categorized into different patterns for a particular disease according to the physical, physiological, and clinical characteristics of the individuals. Obesity, a condition of metabolic disorder, is classified according to six patterns in oriental medicine, as follows: spleen deficiency syndrome, phlegm fluid syndrome, yang deficiency syndrome (YDS), food accumulation syndrome (FAS), liver depression syndrome (LDS), and blood stasis syndrome. In oriental medicine, identification of the disease pattern for individual obese patients is performed on the basis of differentiation in obesity syndrome index and, accordingly, personalized treatment is provided to the patients. The aim of the current study was to understand the obesity patterns in oriental medicine from the genomic point of view via determining the gene expression signature of obese patients using peripheral blood mononuclear cells as the samples.
The study was conducted in 23 South Korean obese subjects (19 female and four male) with BMI ≥25 kg/m2. Identification of oriental obesity pattern was based on the software-guided evaluation of the responses of the subjects to a questionnaire developed by the Korean Institute of Oriental Medicine. The expression profiles of genes were determined using DNA microarray and the level of transcription of genes of interest was further evaluated using quantitative real-time PCR (qRT-PCR).
Results and conclusion
Gene clustering analysis of the microarray data from the FAS, LDS, and YDS subjects exhibited disease pattern-specific upregulation of expression of several genes in a particular cluster. Further analysis of transcription of selected genes using qRT-PCR led to identification of specific genes, including prostaglandin endoperoxide synthase 2, G0/G1 switch 2, carcinoembryonic antigen-related cell adhesion molecule 3, cystein-serine-rich nuclear protein 1, and interleukin 8 receptor, alpha which were highly expressed in LDS obesity constitution. Our current study can be considered as a valuable contribution to the understanding of possible explanation for obesity pattern differentiation in oriental medicine. Further studies can address a novel possibility that the genomic and oriental empirical approaches can be combined and implemented in systematic and synergistic development of personalized medicine.
This clinical trial was registered in Clinical Research Information Service of Korea National Institute of Health (https://cris.nih.go.kr/cris/index.jsp). Registration number: KCT0000387
In recent years, DNA microarray, one of the most advanced techniques for profiling global gene expression, has emerged as a popular tool for investigation of human diseases as it enables comparison of gene expression patterns between normal and pathological tissues and provides unique data about disease mechanisms, regulatory pathways, and gene functions [1–5].
Obesity is a kind of pathophysiological condition characterized by excess accumulation of adipose tissue in the body, which is determined by genetic factors, overeating, inadequate physical activity, and improper lifestyle [6–8]. Obesity is also known to be a prime factor for increasing the risk of other serious diseases, such as diabetes, hypercholesterolemia, fatty liver disease, hypertension, arteriosclerosis, angina, cardiac infarction, stroke, Pickwick syndrome, and gout . Because of these facts, the mortality of obese patients is much higher than that of normal individuals [10–13]. In parallel with the current global scenario, obesity has become one of the most life-threatening problems in South Korea. The Korea National Health and Nutrition Examination reported that the overall prevalence of Korean adult obesity was 30.7 % in 2008, compared with 21.8 % in 1998 .
Traditional medicine, which includes herbal medicine, acupuncture, and other non-medication therapies, represents a comprehensive system of medicine characterized by its inherent theoretical concept and practical experience . Because of its unique prototype and remarkable efficacy with fewer adverse side effects, this system of medicine has gained a considerable interest internationally. Oriental medicine, which is closely related to traditional Chinese medicine, focuses primarily on constitutional approaches . The term ‘body constitution’ (BC) is coined from the Chinese medicinal literature and culture to state one’s current health status which reflects the condition and function of his/her body . More specifically, each constitution is categorized by the characteristics of body shape, face, voice, and psychological and physiological aspects [17, 18], and those characteristics differ between and among the constitutional types . It has been stated that individuals representing different constitution types differ in their susceptibility to disease with a varied degree . It has been also described that in general, a normal or balanced constitution represents an overall healthy state of an individual, whereas individuals with unbalanced constitutions are more susceptible to certain diseases [20, 21]. People can be classified into different types of body characteristics such as strong, weak, ulcer, and bleeding types based on one’s risk or tendency of disease development ; and it has been opined that the relative risk of several chronic diseases differs across different constitution groups [18, 23]. In other words, human constitutions can be categorized into different patterns for a particular disease depending on the physical, physiological, and clinical characteristics of the individuals. Thus, in oriental medicine, subjects with disease are diagnosed based on the pattern identification , a comprehensive system characterized by its own theoretical basis and practical experience . More specifically, the extensive analyses of symptoms and signs performed through this unique system have implications in determining the cause, nature, and location of the illness and the patient’s physical condition, and their treatment approaches . In other words, pattern identification is used as a tool by the medical professionals of oriental medicine to diagnose patients, which ultimately allows personalized treatments according to the symptoms of individual patients.
In oriental medicine, obesity is considered a disease state of humans, who are categorized as ‘Obese People,’ Soft Extratic People,’ or ‘Flesh People.’  Six disease patterns in obesity among the Korean population have been identified and implemented by the Korea Institute of Oriental Medicine (KIOM, Daejeon, South Korea) . These are categorized as spleen deficiency syndrome (SDS), phlegm fluid syndrome (PFS); yang deficiency syndrome (YDS); food accumulation syndrome (FAS), also termed as indigestion or food stagnation syndrome; liver depression syndrome (LDS) or stagnation of liver qi syndrome; and blood stasis syndrome (BSS) . Doctors practicing oriental medicine make a diagnosis of obese patients on the basis of differentiation in obesity syndrome index and accordingly identify the disease pattern and prescribe medicines, providing personalized treatment to the patients. In the field of western medicine, in general, individual patterns are not considered in diagnosis of diseases and prescription of medicines to patients. In contrast, the oriental medicine field follows inherent or empirical traits of patients, sometimes resulting in better medicinal impact with fewer adverse side effects [26–28].
In the current study, we attempted to evaluate the oriental obesity patterns from the genomic point of view via determining the expression profiles of genes using DNA microarray and quantitative real-time PCR (qRT-PCR) techniques. Microarray-based techniques have been adopted previously for identification of the pattern of gene expression in traditional medicinal constitutions [29, 30]. In our study, gene clustering analysis of the microarray data of subjects representing FAS, LDS, and YDS patterns exhibited enriched expression profile of several genes. In addition, subsequent analysis of transcription of selected genes using qRT-PCR led to identification of marker genes that are highly expressed in LDS obesity constitution. Although our study did not categorize all oriental obesity patterns from the genomics point of view, it made a valuable contribution to the classification and acquisition of specific marker genes signifying an empirical science based oriental obesity pattern.
The scientific investigation reported here was an additional study of a clinical trial conducted in Dongguk University Oriental Hospital, South Korea [registered in Clinical Research Information Service of Korea National Institute of Health, (https://cris.nih.go.kr/cris/index.jsp), registration number: KCT0000387]. Briefly, 166 obese subjects were enrolled for the main trial who responded to an advertisement that explained the purpose of the study. They were provided with a brief outline of the trial along with the exclusion criteria. The recruitment procedure was continued for more than a year (October, 2009 to November, 2010). From the entire pool of the participants, 23 subjects (19 females and 4 males) were randomly chosen for the present microarray study. Additional obese subjects were also selectively or randomly chosen from the volunteer pool to provide samples for comparative gene expression analysis. To comply with this study, subjects should be ethnically Korean and within 18 to 65-years-of-age. They should be non-smoker and must possess in-general good physical condition with BMI ≥25 kg/m2. Their body weight should not have changed within ±3 kg during the last 3 months. Prior to the trial, all screened volunteers were informed in detail of the purpose and principles of the study, and the freedom to withdraw from the study. The subjects were then instructed to peruse the consent forms provided to them and sign the forms if they agreed to participate in the trial. Following this, the general health conditions of the volunteers were evaluated by a physician. The following abnormal health conditions were considered among the exclusion criteria: heart, renal, or hepatic disease, history of alanine aminotransferase (ALT) or aspartate aminotransferase (AST) enzyme activity >2.5 times the normal upper limit, or level of creatinine >2.0 mg/dL, as well as malignancy. In addition, subjects under medications that may interfere with the body composition and metabolism, and subjects who had exhibited a hypersensitive reaction to a drug were excluded. Pregnancy at any time point during the study period was also among the exclusion criteria. The trial was initiated following the approval of the medical, scientific, and ethical aspects of the study by the Institutional Review Board of Dongguk University Oriental Hospital (09-CR-12).
Identification of oriental obesity pattern
A computer-based self-reporting online questionnaire developed by Korean Institute of Oriental Medicine (KIOM, Daejeon, South Korea) was used for identification of oriental obesity pattern of the subjects (http://ammrc.kiom.re.kr/renewal/obesity/). During their first visit to the hospital, the subjects were asked to enter their responses to the queries depicted in the online form. The software program developed by KIOM analyzes the responses in the questionnaire and based on these designates every subject as one out of existing six obesity patterns in oriental medicine, as follows: SDS, PFS, YDS, FAS, LDS, and BSS. This questionnaire consists of 52 questions: 8 questions about general or systemic symptoms; 8 on personality and emotion, 18 on digestive function and 18 on circulatory function (see Additional file 1: Table S1 for further details). Each question which has different weight according to its significance for pattern diagnosis is valued on a 5-point scale (1, little; 2, mild; 3, moderate; 4, intensive; 5, very intensive). The software automatically calculates the points for each question after a logical evaluation of the relative significance of the response against the function under which the question belongs. Based on the points accumulated through the subject’s responses to the questionnaire, the KIOM software computes and assigns the score for each oriental obesity pattern. Finally, each subject is designated by a pattern which obtains the highest score among the six patterns in that particular subject. This questionnaire has been validated in several studies measuring the credibility and reliability [26, 28, 31, 32]. The agreement rate between the obesity patterns diagnosed by experts in the field of obesity in oriental medicine and the pattern diagnosis generated by the KIOM questionnaire was investigated [26, 28]. The reported Cronbach alpha coefficient ranged from 0.720 to 0.848.
Isolation of PBMCs
Shortly after their enrolment, blood was collected from the subjects following standard protocol. PBMCs were isolated from the whole heparinized blood by standard gradient centrifugation using Ficoll-Paque Plus reagent (GE Healthcare, Little Chalfont, Waukesha, WI, USA) according to the manufacturer’s instruction. Briefly, 8 mL of whole blood was diluted with phosphate-buffered saline (PBS, pH 7.4) and added carefully to the tube filled with 15 mL of Ficoll-Paque plus reagent. PBMCs were then separated into soluble interface between the upper plasma layer and lower Ficoll layer by centrifuging the samples at 400×g for 30 min at room temperature followed by collection of the resulting PMBC layer. The isolated PBMC fractions were washed twice by centrifugation at 100×g for 10 min at room temperature using PBS.
RNA preparation and DNA microarray
Total RNA of the PMBCs was extracted using TRI reagent (Ambion, Austin, TX, USA) and the RNeasy mini kit (Qiagen, Hilden, Germany) following the reagent and kit manufacturer’s instructions, respectively. The yield of RNA ranged from 5.02 to 15.37 μg with an average of 8.85 μg. The integrity of extracted RNA was verified by gel electrophoresis. DNA microarray of the samples and subsequent analysis of the data were performed as described previously . Briefly, 5 μg of total RNA was reverse transcribed for generation of double-stranded cDNA (dscDNA) using a SuperScript double-stranded cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA). Reactions were terminated by addition of EDTA followed by RNase A treatment. Samples were then subjected to ethanol-precipitation and finally rehydrated to make the stock solution of dscDNA at a concentration of 250 ng/μl. Next, 1 μg dscDNA was labeled with Cy3-conjugated random 9-mer (TriLink Biotechnologies, San Diego, CA, USA) using Klenow fragment (NEB, Beverly, MA, USA); the labeled samples were then subjected to isopropanol precipitation. Four micrograms of Cy3-labeled DNA (containing sample tracking control and alignment oligo) was then hybridized to NimbleGen, 12-plex, human microarray slides (Human Gene Expression 12 × 135 K Microarray, NimbleGene, Madison, WI, USA) for 18 h at 42 °C using the NimbleGen Hybridization system (NimbleGen). Subsequently, the array slides were washed by vigorous agitation in 1 × SSC + 0.1 % SDS for 5 min at 55 °C and in 0.1 × SSC + 0.1 % SDS for 5 min at room temperature. The slides were then rinsed with distilled water and dried by centrifugation. The array images were captured using an InnoScan 900 Series Microarray Scanner (Innopsys, Carbonne, France) and the signals extracted from the scanned images were then imported into NimbleScan software (version 2.5, Nimblegen) for grid alignment and analysis of gene expression data. Expression data were normalized using a quantile normalization method  and Robust Multichip Average (RMA) algorithms . Gene ontology (GO) analysis was performed using the software toolkit provided in the web-accessible Database for Annotation, Visualization, and Integrated Discovery (DAVID) programs (http://niaid.abcc.ncifcrf.gov/home.jsp) . For this analysis, the selected gene list was uploaded into the program and the gene list of Nimblegen microarray chip was used as the background. For clustering analysis, data from microarray was uploaded to the Cluster program and cluster analysis was performed by applying parameters of interest, as described previously . Clustered data were visualized using the TreeView program.
Quantitative real-time PCR (qRT-PCR)
Total RNA of the PMBCs was extracted using TRI reagent (Ambion) according to the reagent manufacturer’s instructions; 500 ng of RNA was processed for cDNA synthesis using the high capacity cDNA reverse transcription kit (Applied Biosystems, Foster City, CA, USA) following the instructions provided by the kit manufacturer. qRT-PCR of the samples was performed in a StepOne™ real-time PCR system (Applied Biosystems). The primer sequences used in our experiment were as follows: NM_000634 forward, 5′-ACC TGG CCG GTG CTT CAG TTA G-3′ and reverse, 5′-TCC CAG CAG GCT CAG CAG GAA-3′; NM_000963 forward, 5′-TTG CCC GAC TCC CTT GGG TGT-3′ and reverse, 5′-AAG TCC ACC CCA TGG CCC AGC-3′; NM_001815 forward, 5′-CAC CAA TGC ATC CCT GCT GAT CCA-3′ and reverse, 5′-GCT GGT TCT TCC AGT TTT GGC AAG G-3′; NM_015714 forward, 5′-GCG CCG TGC CAC TAA GGT CAT TC-3′ and reverse, 5′-AGA CGT CTG GCG GCC GTG AA-3′; NM_033027 forward, 5′-GCA GTG GAA GCT TTC GGC AGC-3′ and reverse, 5′-AGC TCC CGC TTC TCC TCC CG-3′.
The following conditions were applied for the PCR amplification reactions: an initial incubation step at 95 °C for 2 min, followed by 40 amplification cycles, each one consisting of a denaturation step at 95 °C for 10 s, an annealing step at 60 °C for 10 s, and an extension step at 72 °C for 30 s. Following this reaction, a melting curve analysis was performed in order to verify the specificity of the amplicon. The StepOne™ Real-Time PCR Systems software supplied by the instrument manufacturer (Applied Biosystems) was used for processing and analysis of the data. Diluted aliquots (1/10, 1/100) of the cDNA reaction mixture were analyzed for generation of a standard curve. Relative amount of mRNA was quantified using the freshly prepared standard curve.
Results are expressed as the mean ± standard deviation (SD). All statistical analyses were performed using SPSS software (Version 20, SPSS Inc, Chicago, IL, USA). One-way ANOVA was used for comparison of differences among the three oriental obesity patterns. Independent student’s t test was used for comparison of differences between two groups. Correlations between clinical parameters and the level of gene expression were assessed by Pearson’s correlation test. To determine correlation values, heat map analyses were performed using PermutMatrix software (Version 1.0.3 EN, PermutMatrix, Montpellier, France). P value <0.05 was considered statistically significant. As mentioned above, gene ontology (GO) term enrichment was analyzed using tools available on the DAVID website  and data from ‘GOTERM_BP_ALL’ database are presented in Table 1.
Anthropometric and body composition parameters, and oriental obesity patterns of the subjects
The identification numbers, age, sex, general characteristics, anthropometric and body composition parameters, and oriental obesity patterns of the 23 subjects (19 females and 4 males) involved in the present microarray study are shown in Table 2. Among the subjects, 9 showed symptoms of LDS (3 males and 6 females), 8 exhibited FAS pattern (1 male and 7 females) and 6 showed symptoms of YDS (all were females) (Table 3). Statistical analysis using one-way ANOVA revealed no significant differences in the obesity-related parameters among these three disease patterns (Table 3). However, since there are possibilities that age and sex of the subjects might influence the bias, we further extended our statistical evaluation using ANCOVA being adjusted by these two parameters as covariates. The results reconfirmed no significant differences in the obesity-related parameters (Table 3), consistent with our ANOVA analysis. This assures that in our study, the selection bias was minimized by appropriate randomization. Further, to judge whether the distributions of TG and HDL-C might be skewed, we also carried out the non-parametric method of Kruskal–Wallis test for these two parameters. The analysis (P values: 0.722 and 0.319 for TG and HDL-C, respectively) did not show any difference from our evaluation using parametric method.
Validation of specific syndrome differentiation by gene clustering of the obese patients
DNA microarray analysis of the PMBC of the subjects showed no marked difference in the average expression of the probes on the array among the three syndrome indexes mentioned above (Fig. 1), indicating that the quality of the performed microarray assay was satisfactory. When we performed the gene clustering analysis, those three oriental obesity patterns were categorized into three entities (Fig. 2). We explored various filtering conditions for clustering analysis. In most of these tested filtering conditions, clustering of genes in a syndrome index-dependent manner was found. Here we present the clustering results based on one of the filtering conditions as described above where we observed optimum syndrome index dependency. More specifically, three marked classifications were evident when whole genes were filtered according to a specific condition (at least one observation with value ≥2.0), similarity metric [correlation (centered)] and clustering method (centroid linkage). The categories were further divided into two parts. One part was comprised of seven subjects among whom five (71 %) represented LDS pattern (A) and two (29 %) belonged to the FAS pattern type (Fig. 2). The other part representing the remaining 16 subjects was further divided into two groups. One group consisted of 10 subjects among whom five represented the FAS pattern (50 %, B) and the other group consisted of six subjects among whom four (67 %) exhibited the YDS pattern (C).
Identification of differentially expressed genes (DEGs) in specific oriental obesity pattern
For identification of gene sets that are presumed to be differentially regulated among the classified oriental obesity patterns, we performed further gene clustering analysis using the microarray data. As shown in Fig. 2, significant gene sets for each pattern were not detected when all gene expression profiles of 23 subjects were processed for clustering analysis. This result suggests that each pattern may consist of subtype patterns which differ from each other in their gene expression profiles. In this regard, we narrow down numbers of expression profiles based upon clustering results in Fig. 2. More specifically, we performed gene clustering analysis using the microarray data of commonly upregulated genes from the categorized five subjects of LDS (A), five subjects of FAS (B), and four subjects of YDS patterns (C), as mentioned above (Fig. 2). A total of 382 genes were subjected to clustering analysis (Fig. 3a, left panel, the list of genes is depicted in Additional file 2: Table S2) which ultimately revealed that among these genes, 176 were specifically clustered (cluster 1).
Our analysis also demonstrated that cluster 1 was significantly upregulated (P < 0.0001) only in subjects representing LDS pattern compared to volunteers belonging to FAS and YDS pattern types (Fig. 3a, left panel, Fig. 3b). According to gene ontology analysis, genes that were highly enriched in this cluster were found to be related to defense response against pathogens or chemical insults (Table 1). Next, from all of the genes grouped in cluster 1, we sorted a total of 14 genes that were commonly up-regulated in LDS subjects (Fig. 3a, right panel). The identification number and name of these genes are shown in Table 4. Correlation between these 14 genes and the obesity parameters of the aforementioned 14 subjects as analyzed by heatmap plots is shown in Fig. 4. To determine whether the increase in expression of these 14 genes in LDS pattern was significant, we sorted them from the microarray data set of those 14 individuals and performed statistical analysis using student’s t-test. Accordingly, our data showed that among those 14 selected genes, the expression of 10 genes was significantly higher in LDS obesity pattern compared to FAS and YDS patterns (Fig. 5).
Confirmation of higher expression of particular genes in specific oriental obesity pattern
Next, we made an effort to further validate the higher expression of the identified 10 marker genes mentioned above in LDS obese patients (Fig. 5). For this, we selected five genes (NM_000963, prostaglandin endoperoxide synthase 2, PTGS2; NM_015714, G0/G1 switch 2, G0S2; NM_001815, carcinoembryonic antigen-related cell adhesion molecule 3, CEACAM3; NM_033027, cystein-serine-rich nuclear protein 1, CSRNP1; and NM_000634, interleukin 8 receptor, alpha, CXCR1) as representatives and measured their expression individually by quantitative RT-PCR using PBMC samples of the subjects. In this experiment, 15 LDS PBMC samples and 15 randomly selected PBMC samples (for PTGS2 and G0S2), and 10 LDS PBMC samples and 10 randomly selected PBMC samples (for CEACAM3, CSRNP1, and CXCR1) collected from the additional subjects were used for the gene expression profiling. The age, gender, and the clinical features of these subjects are depicted in Additional file 3: Table S3 and Additional file 4: Table S4, respectively. As shown in Fig. 6 and Additional file 5: Figure S1, significantly higher expression of all of the above mentioned five genes was observed in the LDS obese patients compared to the randomly selected obese subjects. Further to examine whether the distribution of the relative mRNA level of PTGS2 might be skewed, we carried out the statistical analysis of this parameter using non-parametric method of Mann–Whitney’s U test. The analysis further confirmed significantly higher PTGS2 expression in LDS compared to random samples (P = 0.03). The relative expression of PTGS2 and G0S2 genes in the 15 LDS and 15 randomly selected subjects are depicted individually in Additional file 6: Table S5 and Additional file 7: Table S6, respectively.
Obesity, which was officially considered a disease by the World Health Organization (WHO) in 1997 , is a serious and alarming global public health problem. Accumulating evidence suggests that obesity can increase the risk of many other life-threatening diseases, including certain types of cancers. WHO predicts that approximately 2.3 billion adults in the global population will be overweight by 2015 and 700 million of them will be obese . Notably, obesity is also becoming a major concern in South Korea as The Korea National Health and Nutrition Examination Survey reported an overall incidence of Korean adult obesity of 30.7 % in 2008, compared with 21.8 % in 1998 .
In the entire human population, each individual has his or her own entity, which differs from that of others. The fundamental concept in the practice of oriental medicine is based on the human constitution, which is classified according to the physical, physiological, and clinical characteristics of the individual. According to this concept, individuals from different constitutional patterns differ in their metabolic pattern and susceptibility to disease with a varying degree. As a norm, clinicians in oriental medicine focus on evaluation of the cause, nature, and location of the disease according to the signs and symptoms shown by the patients , and, based on this information prescribe individualized treatment to the patients . Interestingly, emerging scientific evidence indicates a genetic basis for oriental medicine, where it has been opined that further investigations using up to date genomic tools could help to introduce oriental medicine as an alternative to personalized treatment practiced in western medicine.
In general, oriental obesity patterns are clinically judged by the oriental medical doctors on the basis of the patient’s current symptoms, personality, digestive or circulatory conditions. According to the Korea Institute of Oriental Medicine, among the Korean population, there are six disease patterns in obese people , as follows: SDS, PFS, YDS, FAS, LDS, and BSS. Although the pattern identification is regarded as rather empirical from the scientific point of view, the questionnaire for obesity pattern identification has been validated in various studies [26, 28, 31, 32]. In general, SDS represents a condition for deficiency of the spleen, such as spleen qi deficiency, spleen yang deficiency, and spleen yin deficiency . PFS depicts (1) a condition for retention of phlegm and fluid in any part of the body; (2) a particular state for retention of fluid in the gastrointestinal tract, e.g., gastric retention in pyloric stenosis (WHO). YDS represents a pathological condition characterized by deficiency of the body’s yang qi, which leads to abatement of functions, decline in metabolic activity, and decreased body reactions as well as deficiency-cold manifestations . FAS denotes a pattern/syndrome marked by epigastric and abdominal distention, vomiting of sour matter, anorexia, offensive odor of stools, and curdy and slimy tongue coating . LDS represents a state of stagnation of qi in the liver caused by impairment of free coursing . BSS denotes a morbid condition of blood stagnancy in a particular area of the body resulting from sluggish flow of qi, deficiency of qi or blood, trauma, or yin-cold .
Inter-individual genetic differences in the human population, which are bestowed by a number of factors, such as single nucleotide polymorphisms , copy number variations , and epigenetic or gene expression modifications [41–43], play a vital role in the contribution to phenotypic variations. These variations influence an individual’s anthropometric parameters as well as susceptibility to disease and response to the environment . Therefore, in-depth study of genetic variation may facilitate understanding of the mechanisms underlying differential risk of disease and response to drugs in the human population. Indeed, rapid progress in ‘omics’ sciences has opened a new horizon in traditional medicine research. Accordingly, considerable interest has been paid worldwide, including Korea, India, China, and Japan to studies for understanding the genomic, proteomic, and metabolomic basis of disease patterns.
In recent years, DNA microarray has emerged as a popular technique for investigation of human diseases since it provides unique data about disease mechanisms, regulatory pathways, and gene functions . Microarrays have been employed for investigation of several pathological conditions, such as inflammation, atherosclerosis, breast cancer, colon cancer, and pulmonary fibrosis . In parallel, microarray-based techniques have also been applied in earlier studies to characterize the mode of gene expression in traditional medicinal constitutions and disease patterns. An earlier study attempted to address a scientific approach for oriental medicine treatment by combining oriental therapy for obesity, such as oral administration of Ephedra sinica extract or Sinomenium acutum extract, and western techniques including DNA microarray . However, the empirical and inherent properties of oriental medicine are not well understood, which prompted us to conduct the current study, where we performed DNA microarray analysis using PBMC samples from obese patients in order to categorize specific symptoms in accordance with oriental obesity patterns.
In this study, we acquired genome-wide transcription profiles of 23 subjects who represented the three most frequent oriental obesity pattern types: FAS, LDS, and YDS. According to our microarray data, no marked difference in the average expression of the probes on the array was evident among the three obesity patterns mentioned above (Fig. 1), indicating that the quality of the performed microarray assay was high. When the all genes and array data were filtered and clustered following a specific condition, the subjects were agglomerated into three main categories represented by LDS, FAS, and YDS patterns, respectively. Thus, our results may be considered as an evidence-based approach addressing a possible molecular explanation for obesity pattern differentiation which is normally practiced in oriental medicine.
A number of recent studies have been conducted to understand the possible genetic explanation for the differentiated constitution types or disease patterns in traditional medicine. These include a microarray survey which identified 785 up-regulated genes and 954 down-regulated genes in YDS in comparison with normal constitution . In another study, analysis of genome wide expression through cDNA microarrays revealed a number of differentially expressed genes in each category of individuals from the three most contrasting constitution types (Vata, Pitta, and Kapha) as defined by Ayurveda-based method in India . On the other hand, genes such as angiotensin I converting enzyme, parathyroid hormone, tumor necrosis factor, CD4, glucose phosphate isomerase, C-reactive protein, and vascular endothelial growth factor have been found to show strong relationships with Kidney-Yang Deficiency syndrome .
We further performed gene clustering analysis using the microarray data of commonly categorized five subjects representing LDS, five subjects exhibiting FAS, and four subjects of the YDS. As a result, 176 genes were clustered (cluster 1) and showed significant up-regulation (P < 0.0001) in subjects representing LDS pattern compared to subjects belonging to FAS and YDS pattern types, suggesting that stagnation of qi in the liver may contribute to the etiology of obesity. Indeed, liver stagnation, which is generally caused by prolonged strong emotions or depression, leads to disharmony between the spleen and liver, ultimately giving rise to fluid retention . Because liver becomes depressed in this condition, the gall bladder is also depressed and exhausted. As a result, the ebb and flow of these organs become unbalanced, and the qi mechanism does not flow freely. As a consequence, transformation of fat turbidity becomes difficult and over the course of time the condition leads to obesity .
In an earlier study, 12 proteins were differentially expressed in the serum of patients with liver stagnation syndrome, compared with the normal control group . Further analysis revealed that these differential proteins were mainly related to immunization, neuroendocrine and nutrient metabolism. While, in our study, gene ontology analysis of microarray data demonstrated that in cluster 1, where all 176 genes were significantly upregulated in LDS, the genes related to defense response against pathogens or chemical insults were highly enriched, suggesting a possible functional role of these genes in response to the onset and development of liver stagnation. As liver is involved in neutralization and detoxification of harmful chemicals or pathogens, the up-regulation of defense-system related genes in LDS is not an unexpected phenomenon.
For more extensive analysis of our microarray data, we sorted 14 genes from the aforementioned cluster 1 that were commonly up-regulated in LDS. Results showed that expression of 10 of those 14 genes was significantly elevated in LDS compared to FAS and YDS patterns. Correlation between these 14 genes and the obesity parameters of the aforementioned 14 subjects as analyzed by heatmap plots is shown in Fig. 4. Most notable finding in correlation analysis was that among 14 up-regulated genes in LDS subjects, 6 genes showed significant correlation with systolic blood pressure and most of the genes except NM024075 (TSEN34) displayed a tendency of negative correlation. It is widely known in oriental medicine that LDS, or liver qi stagnation in other word, is related to chronic restraint of stress [50, 51] which ultimately leads to liver disease and stress related symptoms including high blood pressure. Although these genes have not yet found to have a direct effect on blood pressure, this may provide a reasonable evidence for oriental medical theory on which further study will be needed.
Our further study revealed that the expression of five representative members of these 10 genes (NM_000963, PTGS2; NM_015714, G0S2; NM_001815, CEACAM3; NM_033027, CSRNP1; and NM_000634, CXCR1) as determined by quantitative RT-PCR was significantly higher in LDS obese patients compared to randomly selected samples. Notably, a number of studies have reported the impact of PTGS2 and G0S2 genes on hepatic function and pathophysiology. The cyclooxygenases 2-prostanoid pathway has been reported to play vital and complex roles in pathogenesis of various liver diseases . A number of studies have indicated that the COX2-prostanoid pathway may inhibit hepatic fibrogenesis by decreasing proliferation, migration, and contractility of hepatic stellate cells . On the other hand, recent studies have reported that targeted disruption of G0S2 alters hepatic energy balance  and overexpression of this gene induces triglyceride accumulation in the mouse liver leading to development of hepatic steatosis . While emerging evidence indicates possible linkage between CXCR1, CEACAM3, and CSRNP1 genes and liver function. CXCR1 represents the cellular targets of IL-8  which is known to be involved in the pathogenesis of chronic liver disease. Furthermore, it has been reported that CXCR1 expression is induced in hepatocytes after ischemia/reperfusion (I/R) injury . The results suggests that CXCR1 facilitates repair and regenerative responses of liver after I/R injury. It has been found that CXCR1/CXCR2 ligands promote neutrophil migration to the liver and regulate hepatocyte survival and proliferation . CEACAM3 belongs to CEACAMs which are a group of mammalian immunoglobulin related glycoproteins. They play an important role in cell–cell recognition and modulation of cellular processes such as the shaping of tissue architecture and neovascularization, as well as regulation of hepatic insulin sensitivity by facilitating insulin metabolism in liver [58, 59]. Accumulating evidence indicates that CEACAM3 is critically involved in the opsonin-independent recognition and elimination of several Human-specific bacterial pathogens [60, 61]. Considering the above findings on LDS together, our study can be considered as a probable molecular explanation of pattern-specific pathophysiological condition categorized in oriental obesity medication.
Our study has some limitations. First, although our study can be considered as exploratory, the number of subjects in the present microarray-based study is not large. Data based on larger samples are scientifically more valid compared to smaller samples (being provided that all other conditions are equal) because larger samples minimize the probability of errors and maximize the accuracy and precision of the measurements. Second, we cannot consider that the procedure adopted in this study for pattern identification is absolutely concrete as until now there is no gold standard for the pattern identification in oriental medicine. Third, the obesity patterns identified in this study cannot be explained from the western medicinal point of view as the principle of oriental medicine and conventional medicine are quite different. Fourth, there was a time lag between completion of the study and writing, submission and peer review of the manuscript, all of which can be considerable.
In conclusion, our study can be considered as an evidence-based approach providing a possible molecular explanation for obesity pattern differentiation practiced in oriental medicine. More specifically, our findings address a novel possibility that the genomic and oriental empirical approaches can be combined and implemented in systematic and synergistic development of personalized medicine. In our study, gene clustering analysis of the microarray data of obese subjects exhibited disease pattern-specific upregulation of expression of several genes in a particular cluster. Even though we were able to identify five marker genes that were specifically expressed in PBMC of patients from LDS obesity pattern, further studies will allow to find more clear evidences that can segregate other oriental obesity patterns. Thus, our study can contribute to the understanding and development of personalized medicine based on oriental obesity pattern.
Arcellana-Panlilio M, Robbins SM. I. Global gene expression profiling using DNA microarrays. Am J Physiol Gastrointest Liver Physiol. 2002;282:G397–402.
Asyali MH, Colak D, Demirkaya O, Inan MS. Gene expression profile classification: a review. Curr Bioinform. 2006;1:55–73.
Gresham D, Dunham MJ, Botstein D. Comparing whole genomes using DNA microarrays. Nat Rev Genet. 2008;9:291–302.
Schena M, Shalon D, Davis RW, Brown PO. Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science. 1995;270:467–70.
Tarca AL, Romero R, Draghici S. Analysis of microarray experiments of gene expression profiling. Am J Obstet Gynecol. 2006;195:373–88.
Parsons TJ, Manor O, Power C. Physical activity and change in body mass index from adolescence to mid-adulthood in the 1958 British cohort. Int J Epidemiol. 2006;35:197–204.
Gallagher D, Visser M, Sepulveda D, Pierson RN, Harris T, Heymsfield SB. How useful is body mass index for comparison of body fatness across age, sex, and ethnic groups? Am J Epidemiol. 1996;143:228–39.
Anuurad E, Shiwaku K, Nogi A, Kitajima K, Enkhmaa B, Shimono K, Yamane Y. The new BMI criteria for asians by the regional office for the western pacific region of WHO are suitable for screening of overweight to prevent metabolic syndrome in elder Japanese workers. J Occup Health. 2003;45:335–43.
Seifalian AM, Filippatos TD, Joshi J, Mikhailidis DP. Obesity and arterial compliance alterations. Curr Vasc Pharmacol. 2010;8:155–68.
Lee JAKK, Ko HY, Bae KH, Park SY. Recent topics of clinical trials in obesity and metabolic study. J Soci Kor Med Obesi Res. 2009;9:15–22.
Lenz TL, Hamilton WR. Supplemental products used for weight loss. J Am Pharm Assoc JAPhA. 2003;44:59–67 (quiz 67–58).
Park J, Lee M, Kim H, Hong S, Lee D, Yoo J, Choi S, Moon J, Lim C, Lee J. Efficacy and adverse events of bangpungtongseong-san (bofutsusho-san) and bangkihwangki-tang (boiogiot-tang) by oriental obesity pattern identification on obese subjects: randomized, double blind, placebo-controlled trial. J Oriental Rehab Med. 2011;21:265–78.
Park JHLM, Kim HJ. Review on the research related to Bangpungtongseong-san(Bofutshusho-san) in Korea. J Oriental REhab Med. 2011;21:189–202.
Korea Centers for Disease Control and Prevention (KCDC). National Health and Nutrition Survey Microdata 2008. Prevalence of Obestiy. Chungcheongbuk-do: Korea Centers for Disease Control and Prevention; 2009.
World Health Organization. WHO international standard terminologies on traditional medicine in the western pacific region. Western Pacific region: World Health Organization; 2007.
Sun LJ. Study and interpret of body constitution in Chinese medicine. Beijing: Zhongguo Zhongyiyao Chubanshe; 2008. p. 3–5.
Shim EB, Lee S, Kim JY, Earm YE. Physiome and sasang constitutional medicine. J Physiol Sci. 2008;58:433–40.
Kim JY, Pham DD. Sasang constitutional medicine as a holistic tailored medicine. Evid Based Complement Altern Med. 2009;6:11–9.
Lee T, Lee S, Choe B, Song I. A study on the prevalences of chronic diseases according to Sasang constitution at a health examination center. J Sasang Constitut Med. 2005;17:32–45.
Wu Y, Cun Y, Dong J, Shao J, Luo S, Nie S, Yu H, Zheng B, Wang Q, Xiao C. Polymorphisms in PPARD, PPARG and APM1 associated with four types of Traditional Chinese Medicine constitutions. J Genetics Genomics. 2010;37:371–9.
Sun Y, Liu P, Zhao Y, Jia L, He Y, Xue SA, Zheng X, Wang Z, Wang N, Chen J. Characteristics of TCM constitutions of adult Chinese women in Hong Kong and identification of related influencing factors: a cross-sectional survey. J Transl Med. 2014;12:140.
Nana W, Junmei Z. Constitution thoughts in treatise on febrile disease. Henan Tradit Chin Med. 2007;9:003.
Lee T-G, Koh B, Lee S. Sasang constitution as a risk factor for diabetes mellitus: a cross-sectional study. Evid Based Complement Altern Med. 2009;6:99–103.
Kim HJ, Bae HS, Park SU, Moon SK, Park JM, Jung WS. Clinical approach to the standardization of oriental medical diagnostic pattern identification in stroke patients. Evid Based Complement Altern Med. 2011. doi:10.1155/2011/768492.
Hd S. Basic concepts and treatment methods of obesity in oriental medicine. East Asian Soc Diet Life. 2005;29:1–15.
Moon JSKB, Rye EK, Choi SM. A study of syndrome index differentiation in obesity. J Soc Korean Med Obes Res. 2007;7:55–69.
Hong JA, Kim SS, Shin HD. A case report on types of identification of symptom in obese patients. J Oriental Rehab Med. 2011;11:63–71.
Kang KWMJ, Kang BG, Kim BY, Choi SM. The comparison of pattern identification diagnosis acoording to symptom scale based on obesity pattern identification questionnaire. Soc Kor Med Obes Res. 2009;37:37–44.
Wang Q, Yao S. Molecular basis for cold-intolerant yang-deficient constitution of traditional Chinese medicine. Am J Chin Med. 2008;36:827–34.
Prasher B, Negi S, Aggarwal S, Mandal AK, Sethi TP, Deshmukh SR, Purohit SG, Sengupta S, Khanna S, Mohammad F. Whole genome expression and biochemical correlates of extreme constitutional types defined in Ayurveda. J Transl Med. 2008;6:48.
Moon JSKB, Kang KW, Sim MS, Choi SM. Weighting method based on experts opinions for obesity syndrome differentiation questionnaire. J Soc Korean Med Obes Res. 2008;8:51–61.
Kang KWMJ, Kang BG, Kim BY, Kim NS, Yoo JH, Shin MS, Choi SM. The discrimination model for the pattern identification diagnosis of overweight patients. Korean J Oriental Med. 2008;14:41–6.
Hong SW, Kim S, Lee DK. The role of Bach2 in nucleic acid-triggered antiviral innate immune responses. Biochem Biophys Res Commun. 2008;365:426–32.
Bolstad BM, Irizarry RA, Åstrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–93.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4:P3.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998;95:14863–8.
Kopelman PG. Obesity as a medical problem. Nature. 2000;404:635–43.
World Health Organization. Preventing chronic diseases: a vital investment. Geneva: World Health Organization; 2005. http://www.who.int/chp/chronic_disease_report/en/index.html.
Suh Y, Vijg J. SNP discovery in associating genetic variation with human disease phenotypes. Mutation Res Fundam Mol Mech Mutagen. 2005;573:41–53.
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P. Origins and functional impact of copy number variation in the human genome. Nature. 2009;464:704–12.
Sharp AJ, Cheng Z, Eichler EE. Structural variation of the human genome. Annu Rev Genomics Hum Genet. 2006;7:407–42.
Hurles ME, Dermitzakis ET, Tyler-Smith C. The functional impact of structural variation in humans. Trends Genet. 2008;24:238–45.
Schones DE, Zhao K. Genome-wide approaches to studying chromatin modifications. Nat Rev Genet. 2008;9:179–91.
Sachidanandam R, Weissman D, Schmidt SC, Kakol JM, Stein LD, Marth G, Sherry S, Mullikin JC, Mortimore BJ, Willey DL. A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms. Nature. 2001;409:928–33.
Jayapal M, Melendez AJ. DNA microarray technology for target identification and validation. Clin Exp Pharmacol Physiol. 2006;33:496–503.
Yoo YG, Lee HJ, Lee D-K, Joh H-G, Kim D-I, Kim H. Effect of the oral administration of Ephedra Sinica extract on suppression of body weight gains. BioChip J. 2010;4:217–22.
Zhou X, Liu B, Wu Z, Feng Y. Integrative mining of traditional Chinese medicine literature and MEDLINE for functional gene networks. Artif Intell Med. 2007;41:87–104.
Wei S, Xie X. Acupuncture for the treatment of simple obesity: basic and clinical aspects. In: Agbo EC, editor. Innovations in biotechnology. INTECH Open Access Publisher; 2012. doi:10.5772/29783.
Lu CL, Qv XY, Jiang JG. Proteomics and syndrome of Chinese medicine. J Cell Mol Med. 2010;14:2721–8.
Luo H, Ding J, Yue G, Chen J. Metabonomic study of syndrome of liver qi stagnation and spleen deficiency in rats. Zhong xi yi jie he xue bao J Chin Integr Med. 2007;5:307–13.
Hwang MJ, Chung SH, Hwang DS, Song MY. Correlation between visceral adipose tissue and stagnation of liver qi in Korean perimenopausal women. J Oriental Rehab Med. 2008;18:143–56.
Hu K-Q. Cyclooxygenase 2 (COX2)-prostanoid pathway and liver diseases. Prostaglandins Leukot Essent Fatty Acids. 2003;69:329–37.
Wang Y, Zhang Y, Qian H, Lu J, Zhang Z, Min X, Lang M, Yang H, Wang N, Zhang P. The g0/g1 switch gene 2 is an important regulator of hepatic triglyceride metabolism. PLoS One. 2013;8:e72315.
Zhang X, Xie X, Heckmann BL, Saarinen AM, Czyzyk TA, Liu J. Targeted disruption of G0/G1 Switch Gene 2 enhances adipose lipolysis, alters hepatic energy balance, and alleviates high-fat diet-induced liver steatosis. Diabetes. 2014;63:934–46.
Murphy PM. Neutrophil receptors for interleukin-8 and related CXC chemokines. Semin Hematol. 1997;34(4):311–8.
Clarke C, Kuboki S, Sakai N, Kasten KR, Tevar AD, Schuster R, Blanchard J, Caldwell CC, Edwards MJ, Lentsch AB. CXC chemokine receptor-1 is expressed by hepatocytes and regulates liver recovery after hepatic ischemia/reperfusion injury. Hepatology. 2011;53:261–71.
Saiman Y, Friedman SL. The role of chemokines in acute liver injury. Front Physiol. 2012. doi:10.3389/fphys.2012.00213.
Kuespert K, Pils S, Hauck CR. CEACAMs: their role in physiology and pathophysiology. Curr Opin Cell Biol. 2006;18:565–71.
Poy MN, Yang Y, Rezaei K, Fernström MA, Lee AD, Kido Y, Erickson SK, Najjar SM. CEACAM1 regulates insulin clearance in liver. Nat Genet. 2002;30:270–6.
Schmitter T, Pils S, Sakk V, Frank R, Fischer K-D, Hauck CR. The granulocyte receptor carcinoembryonic antigen-related cell adhesion molecule 3 (CEACAM3) directly associates with Vav to promote phagocytosis of human pathogens. J Immunol. 2007;178:3797–805.
Schmitter T, Agerer F, Peterson L, Münzner P, Hauck CR. Granulocyte CEACAM3 is a phagocytic receptor of the innate immune system that mediates recognition and elimination of human-specific pathogens. J Exp Med. 2004;199:35–46.
Profs. DL, SK, and HK were responsible for the clinical study design and fund raising. JHP, KH conducted the experiment and collected data. SWH, JWY, and SB wrote the manuscript. All authors read and approved the final manuscript.
This study was supported by a Grant from the Traditional Korean Medicine R&D Project, Ministry for Health & Welfare & Family Affairs, Republic of Korea (B090014).
The authors declare that they have no competing interests.
Sun Woo Hong, Jae-Wook Yoo and Shambhunath Bose contributed equally to this work
Additional file 1: Table S1. Obesity pattern diagnosis questionnaire. The questionnaire consisting of 52 questions addressing systemic symptoms, personality and emotion, digestive function, and circulatory function of the subjects are enlisted.
Additional file 2: Table S2. List of commonly upregulated genes that were subjected to cluster analysis.
Additional file 3: Table S3. Clinical characteristics of 15 randomly selected subjects and 15 volunteers representing liver depression syndrome (LDS) pattern the samples of whom were used for the quantitative real-time PCR analyses of prostaglandin endoperoxide synthase 2 (PTGS2), G0/G1 switch 2 (G0S2) genes.
Additional file 4: Table S4. Clinical characteristics of 10 randomly selected subjects and 10 volunteers representing liver depression syndrome (LDS) pattern the samples of whom were used for the quantitative real-time PCR analyses of carcinoembryonic antigen-related cell adhesion molecule 3 (CEACAM3), cystein-serine-rich nuclear protein 1 (CSRNP1), and interleukin 8 receptor, alpha (CXCR1) genes.
Additional file 5: Figure S1. Expression levels of CEACAM3, CSRNP1, and CXCR1 genes in 10 randomly selected subjects (Random) or 10 volunteers representing LDS pattern (LDS) were measured using quantitative RT-PCR. P-value was determined using student’s t-test.
Additional file 6: Table S5. Expression levels of endogenous PTGS2 gene in 15 randomly selected subjects (Random) or 15 volunteers representing LDS pattern (LDS) as measured by quantitative RT-PCR.
Additional file 7: Table S6. Expression levels of endogenous G0S2 gene in 15 randomly selected subjects (Random) or 15 volunteers representing LDS pattern (LDS) as measured by quantitative RT-PCR.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Hong, S.W., Yoo, JW., Bose, S. et al. Understanding the molecular aspects of oriental obesity pattern differentiation using DNA microarray. J Transl Med 13, 331 (2015). https://doi.org/10.1186/s12967-015-0692-9
- Oriental medicine
- Pattern identification