- Research article
- Open Access
Transcriptomic profiles of peripheral white blood cells in type II diabetes and racial differences in expression profiles
BMC Genomics volume 12, Article number: S12 (2011)
Along with obesity, physical inactivity, and family history of metabolic disorders, African American ethnicity is a risk factor for type 2 diabetes (T2D) in the United States. However, little is known about the differences in gene expression and transcriptomic profiles of blood in T2D between African Americans (AA) and Caucasians (CAU), and microarray analysis of peripheral white blood cells (WBCs) from these two ethnic groups will facilitate our understanding of the underlying molecular mechanism in T2D and identify genetic biomarkers responsible for the disparities.
A whole human genome oligomicroarray of peripheral WBCs was performed on 144 samples obtained from 84 patients with T2D (44 AA and 40 CAU) and 60 healthy controls (28 AA and 32 CAU). The results showed that 30 genes had significant difference in expression between patients and controls (a fold change of <-1.4 or >1.4 with a P value <0.05). These known genes were mainly clustered in three functional categories: immune responses, lipid metabolism, and organismal injury/abnormaly. Transcriptomic analysis also showed that 574 genes were differentially expressed in AA diseased versus AA control, compared to 200 genes in CAU subjects. Pathway study revealed that "Communication between innate and adaptive immune cells"/"Primary immunodeficiency signaling" are significantly down-regulated in AA patients and "Interferon signaling"/"Complement System" are significantly down-regulated in CAU patients.
These newly identified genetic markers in WBCs provide valuable information about the pathophysiology of T2D and can be used for diagnosis and pharmaceutical drug design. Our results also found that AA and CAU patients with T2D express genes and pathways differently.
Type 2 diabetes (T2D) is the most common form diabetes, and it accounts for more than 90% of all diagnosed diabetic cases . It is a heterogeneous and multi-factorial disease. Although the exact pathogenesis has not been fully understood, older age, family history, physical inactivity and race/ethnicity are important risk factors[1, 2]. African American (AA) is the ethnicity with the highest incidence and mortality rate of T2D. According to 2010 statistics of adults aged 20 years or older, 15.7 million non-Hispanic whites have diabetes (10.2% of the ethnic group), compared to 4.9 million non-Hispanic blacks (18.7%). After adjusting for population age differences, AAs are 1.8 times as likely to have diabetes as CAUs at similar ages [3–5]. In addition, obesity and immune/inflammatory processes contribute to the onset of T2D [6–9]. For example, elevated levels of circulating IL-1β, IL-6, TNFα, and acute phase proteins in T2D have been extensively described in several studies. That may reflect the activation of innate immune cells by excessive levels of nutrient concentration, including glucose and free fatty acids. In turn, both IL-6 and IL-1β act on the liver to produce a characteristic dyslipidemia in metabolic syndrome, one which manifests as increased very low density lipoprotein (VLDL) and decreased high density lipoprotein (HDL) [10, 11]. IL-6 and tumor necrosis factor TNFα also participate in the impairment of insulin signaling in adipocytes .
Many clinical and epidemiological studies have been done to elucidate environmental contribution to T2D [13–15]. Microarray analysis will allow us to study genetic contribution and patterns of altered gene expression related to T2D. Thus far, human skeletal muscles, liver, adipose tissues, and pancreatic β-cells [16–22] have been chosen for gene expression assay in T2D. These studies have yielded important insights into the pathogenesis of T2D and its relationship to obesity. For example, adipose tissue gene expression has been found to show more differences between T2D and healthy control than muscle tissue gene expression. Many genes are inflammatory and metabolically important . Because it is difficult to obtain permission from a clinic to use tissue samples such as liver, adipose, and muscle for genetic studies, the peripheral blood could be a convenient source of cells for this study. To date, only one microarray study of human peripheral blood has been reported by Grayson et al . Their finding indicated that genes differentially expressed in T2D have key roles in T-cell activation and signaling, but this study used a limited sample size (n = 6) and lacked ethnicity information. In this study, we measured the up/down regulation of some T2D-related genes using white blood cells from T2D patients, and we compared the gene expression profiles between two ethnic groups. The information obtained from our study would provide further understanding of the pathogenesis of T2D and the disparities in expression profiles between AA and CAU.
Characteristics of studied subjects
This study comprised a balanced distribution of the studied subjects in gender and ethnicity: among 60 controls, 28 were AA including 14 females and 14 males, and 32 were CAU including 14 females and 20 males. Among 84 patients with T2D, 44 were AA including 22 females and 22 males, and 40 were CAU including 23 females and 17 males. As compared to AA, CAU had a significantly higher level of blood triglycerides (TG) in both the controls (106 ± 54.3 mg/dl in AA vs. 153 ± 77.8 mg/dl in CAU, p = 0.0009), and the patients (157 ± 128 mg/dl in AA vs. 207 ± 98.3 mg/dl in CAU, p = 0.037). There were no significant differences in other studied clinical parameters between two races (data for ethnicity differences were not shown). As compared to all controls (mixed), the patients group was 4.5 years older, had significantly higher body mass index (BMI), blood TG, and fasting glucose, and had lower high density apolipoprotein (HDL). There were no differences in low density apolipoproteins (LDL) and total cholesterols (Table 1) between controls and T2D patients.
Identification of differentially expressed genes between all T2D patients and healthy controls
Whole human genome (4 × 44K) oligo-microarrays of peripheral WBCs from 60 healthy subjects and 84 T2D subjects were performed. Of 41,000 probes available on the Agilent arrays, we identified 30 genes that were differentially expressed (a fold change of <-1.4 or >1.4 with a P value <0.05). Six known genes were up-regulated, including erythroid associated factor (ERAF), delta-aminolevulinate synthase 2 (ALAS2), oxysterol binding protein 2 (OSBP2), carbonic anhydrase I (CA1), serine/threonine/tyrosine kinase 1 (STYK1), and Zinc finger protein 2 (ZIC2). Sixteen known genes were down-regulated, including G0/G1 switch 2 (GOS2); telomerase-associated protein 1 (TEP1); prostaglandin-endoperoxide synthase 2 (PTGS2); interleukin 4 (IL4); interleukin 8 (IL8); interferon alpha-inducible protein 27 (IFI27); interferon-induced protein with tetratricopeptide repeats 3 (IFIT3); interferon-induced protein with tetratricopeptide repeats 2 (IFIT2); tumor necrosis factor alpha-induced protein 6T (NFAIP6); radical S-adenosyl methionine domain containing 2 (RSAD2); apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3A (APOBEC3A); ATP/GTP binding protein-like 4 (ABGL4); ATP-binding cassette sub-family A transporter (ABCA1); epithelial stromal interaction 1 isoform 2 (EPSTI1); Ephrin type-A receptor 6 (EPHA6); and leucine rich repeat neuronal protein 3 (LRRN3). Most of these genes belonged to three functional categories: immune and inflammatory responses, such as IL-4, IL-8, TNFAIP6, and PTGS2; lipid and glucose metabolism, such as ABCA1, OSBP2, and GOS2; and organismal injury and abnormalities, such as APOBEC3A. Thus, these genes might be the genetic signatures for T2D.
Ethnic differences in gene expression profiles
In order to study ethnic disparities, gene expression profiles between AA and CAU participants were analyzed and compared. Based on the same criteria described before, many more genes were found to be differentially expressed. 574 genes were differentially expressed in AA (28 healthy controls and 44 T2D patients) compared to 200 genes in CAU (32 healthy controls and 40 T2D patients). In the AA group, there were 294 down-regulated and 269 up-regulated genes in T2D patients; however, we only saw 134 down-regulated and 47 up-regulated genes in CAU T2D patients (Figure 1a and 1b). It was obvious that AA T2D patients had 2-6 times more altered genes than CAU T2D patents. Table 2 summarizes the top 10 up-regulated and down-regulated genes in different ethnic groups (a fold change of <-1.4 or >1.4 with a P value <0.05). Overlap analysis showed that AA had an entirely different spectrum of genetic signatures for T2D, because only one down-regulated gene and one up-regulated gene were common between the two ethnicities (Figure 1a and 1b).
Ethnic differences in pathway enrichment analysis
Gene expression profiling study showed ethnic disparities in this study cohort as described above. Because one pathway usually involves many different genes, it is not clear whether the differentially expressed genes in the two ethnic groups involved the same pathways. Using the tool of Ingenuity Canonical Pathways Analysis, the differentially expressed genes in two ethnic groups were clustered into various pathways. The top 10 down-regulated and up-regulated pathways were identified for T2D patients in AA and in CAU as shown in Figure 2. Similar to the results of differentially expressed single genes, the two T2D patient groups had different spectra of down- and up-regulated pathways. For example, significant enrichments of down-regulated pathways in AA patients were "communication between innate and adaptive immune cells" and "primary immunodeficiency signaling", whereas those in CAU patients were "interferon signaling" and "complement system" pathways. The main up-regulated pathways in AA included "inositol metabolism", "Wnt/β-catenin signaling", and "LXR/RXR activation", but for CAU they were "tight junction signaling", "cleavage and polyadenyation of pre-mRNA", and "methane metabolism".
Microarray validation by real-time PCR
To confirm the accuracy of the results in microarray study, real-time PCR were performed on 5 randomly selected genes (ABCA1, IL-4, IL-8 and PTGS2 except ALAS2) that were found to be differentially expressed in microarray analysis. All of these genes have a known function in glucose and lipid metabolism and in immune responses. The fold changes for the transcripts of these genes were normalized to the transcripts of glyceraldehyde-3-phosphate dehydrogenase (GAPDH); information of primers for the 5 selected genes and GAPDH was shown in table 3. The fold changes of those five gene expressions, as determined by microarray and real-time PCR, are shown in Figure 3. The results validated our microarray data.
Previous studies of gene expression in T2D have been restricted to adipose tissue and muscle or to animal models with limited sample size [16–22]. Before our study, the comparison of blood gene expression in AA and CAU with T2D patients was lacking. Our analysis of gene profiling in peripheral WBCs from a larger number of sample and ethnic groups begins to fill the gaps. We showed that 30 genes are differentially expressed in T2D patients when compared with healthy controls. Furthermore, we demonstrated that different gene expression signatures are present in AA and CAU patients.
In T2D, we see increased expression of six genes in which five of those have not been previously reported in gene expression studies. ERAF is associated with haemoglobin stabilization , ALAS2 is involved in immune responses , OSBP2 in lipid metabolism , STYK1 in signal transduction , and ZIC2 in transcriptional regulation . Carbonic anhydrase I, encoded by CA1, has been reported to be progressively increased in the urine of T2D patients . Interestingly, in our study, we found no evidence for up-regulation of classical inflammatory markers including the interleukin family, the TNF family, or interferon family which have been shown in the literature [5–9]; however, a panel of genes belonging to those families is down-regulated in our study cohort, including PTGS2, IL4, IL8, IFI27, IFIT3, IFIT2, TNFAIP6, RSAD2, and AGBL4. All of those are associated with the interferon signaling pathway, immune signaling pathways involved in inflammatory processes . Our results suggested that the innate immunity was actually compromised, and significantly lower inflammation was observed in T2D patients. These novel findings are inconsistent with the results from other studies, which may indicate an important difference in the pathophysiology of T2D between the tissues. A recent study  also showed that elevated levels of circulating inflammatory markers in T2D may not necessarily reflect the degree of inflammation in individual tissues because of the mass difference. In summary, peripheral WBC gene expression profiling may not mirror the changes of gene expression within a specific tissue such as adipose or muscle.
In addition to inflammatory markers, we found down-regulation of several genes (for example ABCA1 and GOS2) which are reported to be associated with lipid metabolism and T2D in gene knockout cell lines and/or animal models [33–35]. For example, recent studies have shown that ATP-binding cassette transporters ABCA1 and ABCG1 play important roles in macrophage cholesterol efflux to serum or HDL in macrophage foam cells . Mice lacking ABCA1 and ABCG1 accumulate inflammatory macrophage foam cells in various tissues such as lung, liver, spleen, or myocardium . Therefore, ABCA1 may have a unique role in dampening inflammation. Consistent with their finding, our data showed that ABCA1 gene was down-regulated in T2D patients. In addition, we also see that patients had significant lower HDL when compared with control subjects (Table 1). To our knowledge, this is the first evidence that lower expression of ABCA1 decreases the cholesterol efflux from macrophage, resulting in less HDL formation in T2D patients. In turn, the anti-inflammatory and immunosuppressive functions of ABCA1 are diminished.
Yang et al.,  recently reported that G0S2 expression was significantly reduced in white adipose tissues (WAT) of db/db mice and high-fat-fed wild-type mice. They speculated that decreased insulin sensitivity of adipocytes may lead to the down-regulation of G0S2, but they are not sure whether expression of G0S2 is also decreased in obese humans. Our gene expression profiling provided strong evidence that GOS2 was down-regulated in T2D patients, who also had significant higher BMI compared tohealthy controls.
Type 2 diabetes (T2D) is a complex disease. Although the exact pathogenesis has not been fully understood, family history, physical inactivity, socioeconomic condition, and race/ethnicity contribute to the onset of T2D [13–15]. AAs have the highest incidence and mortality rate of T2D . To better understand the health disparity, we compared the clinical parameters and expression profiles of peripheral WBCs from AA and CAU. As described previously, CAU had a significantly higher level of blood triglycerides (TG) in both the controls and the patients. There was no significant difference in other clinical parameters between two races. However, as shown in Figure 1, AA T2D patients had 2-6 times more altered genes than CAU T2D patients. In AAs, there were 294 down-regulated and 269 up-regulated genes in T2D patients; however, we only saw 134 down-regulated and 47 up-regulated genes in CAUs. The differences were also observed in the top 10 up-regulated and down-regulated genes from the two ethnic groups shown in table 3. None of the up- or down-regulated genes overlapped between AAs and CAUs. In agreement with gene expression differences, pathway enrichment study further suggested that AAs have different pathophysiological mechanisms of T2D from CAU. For example, significant enrichments of top 10 down-regulated pathways in AA patients included "communication between innate and adaptive immune cells" and "primary immunodeficiency signaling", whereas that in CAU patients were "interferon signaling" and "complement system" pathways. The main up-regulated pathways in AA were "inositol metabolism", "Wnt/β-catenin signaling" and "LXR/RXR activation", but these pathways were not significantly altered in CAU. Taken together, differences between AA and CAU participants in expression profiles parallel differences between these groups in pathway enrichment data. Ethnic differences in TG cannot explain ethnic disparities. The possibility that CAU persons have systematically higher levels of TG warrants further study.
In summary, through the analyses of clinical parameters, expression profiles of T2D, and ethnic differences of gene expression and pathway enrichment, we found that the innate immunity was compromised and significantly lower inflammation was associated with T2D. The newly identified ABCA1 and GOS2 in human could be used as genetic markers for diagnosis and pharmaceutical drug design. Our results also indicated that genes and pathways are differentially expressed between AA and CAU T2D patients.
Materials and methods
Subjects and clinical laboratory data
The study was approved by the Institutional Review Board of Tougaloo College. All subjects provided written informed consent for this study. T2D was diagnosed based on American Diabetes Association (ADA) guidelines  and characteristic symptoms of diabetes: a higher BMI and a fasting plasma glucose > 126 mg dl-1 or a 2 h plasma glucose during an oral glucose tolerance test of > 200 mg dl-1. A total of 144 blood samples from healthy controls (n = 60, 32 Caucasians and 28 African Americans), and T2D (n = 84, 40 Caucasians and 44 African Americans) were collected. All subjects were evaluated by age, sex, ethnicity, body mass index (BMI), triacylglycerol (TG), high-density lipoprotein (HDL), low-density lipoprotein (LDL), total cholesterol (TC), and glucose levels.
Total RNA from 8-10 mLs peripheral blood WBCs was obtained using LeukoLock™ Total RNA system (Ambion Inc, Austin, TX) according to the manufacturer's instructions. The quantity and quality of the isolated RNA were determined by Nanodrop spectrophotometry and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).
Gene expression profiling was conducted using Agilent Whole Human Genome1 (4 × 44K) Oligo arrays with ~20,000 genes represented (Agilent Technologies, Palo Alto, CA). Each sample was hybridized against a human universal RNA control (Stratagene, La Jolla, CA). 500 ng of total RNA was amplified and labeled using the Agilent Low RNA Input Fluorescent Linear Amplification Kit, according to manufacturer's protocol. For each two color array, 850 ng of each Cy5- (universal control) and Cy3-labeled (sample) cRNA were mixed and fragmented using the Agilent In Situ Hybridization Kit protocol. Hybridizations were performed for 17 hours in a rotating hybridization oven according to the Agilent 60-mer oligo microarray processing protocol prior to washing and scanning with an Agilent Scanner (G2565AA, Agilent Technologies, Wilmington, DE). Arrays were processed and background corrected with default settings for all parameters with the Agilent Feature Extraction software (v.220.127.116.11).
Microarray data analysis
Microarray data analyses were processed with GeneSpring version 7.0 and 10.0. The sample quality control was based on the Pearson correlation of a sample with other samples in the whole experiment. If the average Pearson correlation with other samples was less than 80%, the sample was excluded for further analysis. If the scanned intensity was less than 5.0 for a probe, it was transformed to 5. A perchip (within) array normalization was performed using 50 percentile values of all the probe values in the array. Per gene (between) array normalization was also applied using the median value of a gene across all samples in the experiment. Probe features were first filtered using flags. A "present" or "absent" flag was defined using the Agilent Feature Extraction 9.5.1 software. Only a probe that had present flags in at least 50% samples of all the arrays was kept for further analysis. Data were subsequently log (base 2) transformed for statistical analysis. To identify differentiated genes between two groups (e g. normal vs. diabetes), an un-paired t-test with cut off p value 0.05 was applied to compare control samples and TNT exposed samples. In addition, 1.4 fold changes were applied to identify more significantly regulated genes.
Gene functional analysis and pathway analysis
Significantly regulated probes were employed for one way hierarchical clustering (only cluster genes) or two-way hierarchical clustering (clustering both genes and samples) using GeneSpring 7.0 and/or 10 (Agilent Technologies, Foster City, CA, USA). A Pearson correlation with average linkage was applied for the clustering. Gene functional categories were classified according to Gene Ontology (GO)  using The Database for Annotation, Visualization and Integrated Discovery (DAVID) [37, 38] and Gofetch  tools. A Gene Ontology functional term enrichment p value less than 0.05 was considered significant. Pathway analysis was performed using the Ingenuity canonical pathways analysis tool. Similar to GO analysis, a pathway with an enrichment p value less than 0.05 was considered to be a significantly regulated pathway (Ingenuity Systems, Inc., Redwood City, CA). Gene networks were constructed based on the Ingenuity knowledge base. A score was assigned to a network according to the fit of the original set of significant genes. This score reflects the negative logarithm of the p value that indicates the likelihood of the focus genes in a network being found together due to random chance .
To validate the microarray data, a real-time PCR assay was performed using Brilliant II SYBR® Green system for five genes. Samples were selected from the original microarray experiments for further RT-PCR testing based on sufficient RNA remaining. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression was used as an internal control and the relative expression of each gene was determined by the ΔΔCt method, comparing expression of test gene to an average GAPDH, and then comparing the T2D group versus healthy control group. The primers used for cDNA amplification are shown in Table 3.
Centers for Disease Control and Prevention: National diabetes fact sheet. 2008, US Department of Health and Human Services, Centers for Disease Control and Prevention: Atlanta, GA
Stumvoll M, Goldstein BJ, van Haeften TW: Type 2 diabetes: principles of pathogenesis and therapy. Lancet. 2005, 365: 1333-1346. 10.1016/S0140-6736(05)61032-X.
Griffith ML, Younk LM, Davis SN: Visceral Adiposity, insulin resistance, and type 2 diabetes. American J Life Med. 2010, 4: 230-243. 10.1177/1559827609360959.
Harris MI: Racial and ethnic differences in health care access and health outcomes for adults with type 2 diabetes. Diabetes Care. 2001, 24: 454-459. 10.2337/diacare.24.3.454.
American Diabetes Association: Diabetes statistics, data from the 2011 national diabetes fact sheet. 2011
Pickup JC, Mattock MB, Chusney GD, Burt D: NIDDM as a disease of the innate immune system: association of acute-phase reactants and interleukin6with metabolic syndrome X. Diabetologia. 1997, 40: 1286-1292. 10.1007/s001250050822.
Spranger J: Inflammatory cytokines and the risk to develop type 2 diabetes: results of the prospective population-based European prospective investigation into cancer and nutrition (EPIC)-potsdam study. Diabetes. 2003, 52: 812-817. 10.2337/diabetes.52.3.812.
Herder C: Inflammation and type 2 diabetes: results from KORA Augsburg. Gesundheitswesen. 2005, 67: S115-S121. 10.1055/s-2005-858252.
Herder C: Elevated levels of the anti-inflammatory interleukin-1 receptor antagonist precede the onset of type 2 diabetes: the Whitehall II study. Diabetes Care. 2009, 32: 421-423. 10.2337/dc08-1161.
Bilan PJ, Samokhvalov V, Koshkina A, Schertzer JD, Samaan MC, Klip A: Direct and macrophage-mediated actions of fatty acids causing insulin resistance in muscle cells. Arch Physiol Biochem. 2009, 115: 176-190. 10.1080/13813450903079314.
Hotamisligil GS: Inflammation and metabolic disorders. Nature. 2006, 444: 860-867. 10.1038/nature05485.
Wellen KE, Hotamisligil GS: Obesity-induced inflammatory changes in adipose tissue. J Clin Invest. 2003, 112: 1785-1788.
Harris KM, Gordon-Larsen P, Chantala K, Udry R: Longitudinal trends in race/ethnic disparities in leading health indicators from adolescence to young adulthood. Arch Pediatr Adolesc Med. 2006, 160: 74-81. 10.1001/archpedi.160.1.74.
Muoio DM, Newgard CB: Mechanisms of disease: molecular and metabolic mechanisms of insulin resistance and beta-cell failure in type 2 diabetes. Nat Rev Mol Cell Biol. 2008, 9 (3): 193-205. 10.1038/nrm2327.
Zhang Q, Wang YF, Huang ES: Changes in racial/ethnic disparities in the prevalence of type 2 diabetes by obesity level among U.S. adults. Ethn Health. 2009, 14 (5): 439-5. 10.1080/13557850802699155.
Wang M, Wang XC, Zhang ZY, Mou B, Hu RM: Impaired mitochondrial oxidative phosphorylation in multiple insulin-sensitive tissues of humans with type 2 diabetes mellitus. J Int Med Res. 2010, 38 (3): 769-81.
Alibegovic AC, Sonne MP, Højbjerre L, Bork-Jensen J, Jacobsen S, Nilsson E, Faerch K, Hiscock N, Mortensen B, Friedrichsen M, Stallknecht B, Dela F, Vaag A: Insulin resistance induced by physical inactivity is associated with multiple transcriptional changes in skeletal muscle in young men. Am J Physiol Endocrinol Metab. 2010, 299 (5): E752-63. 10.1152/ajpendo.00590.2009.
Matsumoto K, Yokoyama S: Gene expression analysis on the liver of cholestyramine-treated type 2 diabetic model mice. Biomed Pharmacother. 2010, 64 (6): 373-8. 10.1016/j.biopha.2010.02.008.
Hayashi Y, Kajimoto K, Iida S, Sato Y, Mizufune S, Kaji N, Kamiya H, Baba Y, Harashima H: DNA microarray analysis of whole blood cells and insulin-sensitive tissues reveals the usefulness of blood RNA profiling as a source of markers for predicting type 2 diabetes. Biol Pharm Bull. 2010, 33 (6): 1033-42. 10.1248/bpb.33.1033.
Wolfs MG, Rensen SS, Bruin-Van Dijk EJ, Verdam FJ, Greve JW, Sanjabi B, Bruinenberg M, Wijmenga C, van Haeften TW, Buurman WA, Franke L, Hofker MH: Co-expressed immune and metabolic genes in visceral and subcutaneous adipose tissue from severely obese individuals are associated with plasma HDL and glucose levels: a microarray study. BMC Med Genomics. 2010, 3: 34-10.1186/1755-8794-3-34.
Permana PA, Nair S, Lee YH: Application of DNA microarray to the study of human adipose tissue/cells. Methods Mol Biol. 2008, 456: 141-54. 10.1007/978-1-59745-245-8_11.
Marselli L, Thorne J, Dahiya S, Sgroi DC, Sharma A, Bonner-Weir S, Marchetti P, Weir GC: Gene expression profiles of Beta-cell enriched tissue obtained by laser capture microdissection from subjects with type 2 diabetes. PLoS One. 2010, 5 (7): e11499-10.1371/journal.pone.0011499.
Elbein SC, Kern PA, Rasouli N, Yao-Borengasser A, Sharma NK, Das SK: Global gene expression profiles of subcutaneous adipose and muscle from glucose-tolerant, insulin-sensitive, and insulin-resistant individuals matched for BMI. Diabetes. 2011, 60 (3): 1019-29. 10.2337/db10-1270.
Grayson BL, Wang L, Aune TM: Peripheral blood gene expression profiles in metabolic syndrome, coronary artery disease and type 2 diabetes 2011. Genes Immun. 2011, 12: 341-351. 10.1038/gene.2011.13.
Yu X, Kong Y, Dore LC, Abdulmalik O, Katein AM, Zhou S, Choi JK, Gell D, Mackay JP, Gow AJ, Weiss MJ: An erythroid chaperone that facilitates folding of alpha-globin subunits for hemoglobin synthesis. J Clin Invest. 2007, 117 (7): 1856-65. 10.1172/JCI31664.
Grünewald E, Kinnell HL, Porteous DJ, Thomson PA: GPR50 interacts with neuronal NOGO-A and affects neurite outgrowth. Mol Cell Neurosci. 2009, 42 (4): 363-37. 10.1016/j.mcn.2009.08.007.
Hynynen R, Suchanek M, Spandl J, Bäck N, Thiele C, Olkkonen VM: OSBP-related protein 2 is a sterol receptor on lipid droplets that regulates the metabolism of neutral lipids. J Lipid Res. 2009, 50 (7): 1305-15. 10.1194/jlr.M800661-JLR200.
Chung S, Tamura K, Furihata M, Uemura M, Daigo Y, Nasu Y, Miki T, Shuin T, Fujioka T, Nakamura Y, Nakagawa H: Overexpression of the potential kinase serine/threonine/tyrosine kinase 1 (STYK 1) in castration-resistant prostate cancer. Cancer Sci. 2009, 100 (11): 2109-14. 10.1111/j.1349-7006.2009.01277.x.
McMahon AR, Merzdorf CS: Expression of the zic1, zic2, zic3, and zic4 genes in early chick embryos. BMC Res Notes. 2010, 16;3: 167-
Bellei E, Rossi E, Lucchi L, Uggeri S, Albertazzi A, Tomasi A, Iannone A: Proteomic analysis of early urinary biomarkers of renal changes in type 2 diabetic patients. Proteomics Clin Appl. 2008, 2 (4): 478-91. 10.1002/prca.200780109.
Kurji KH, Cui JZ, Lin T, Harriman D, Prasad SS, Kojic L, Matsubara JA: Microarray analysis identifies changes in inflammatory gene expression in response amyloid-stimulation of cultured human retinal pigment epithelial cells. Invest Ophthalmol Visual Sci. 2010, 51 (2):
Donath MY, Shoelson SE: Type 2 diabetes as an inflammatory disease. Nat Rev Immunol. 2011, 11 (2): 98-107. 10.1038/nri2925.
Yvan-Charvet L, Wang N, Tall AR: Role of HDL, ABCA1, and ABCG1 transporters in cholesterol efflux and immune responses. Arterioscler Thromb Vasc Biol. 2010, 30: 139-143. 10.1161/ATVBAHA.108.179283.
Yvan-Charvet L, Ranalletta M, Wang N, Han S, Terasaka N, Li R, Welch C, Tall A: Combined deficiency of ABCA1 and ABCG1 promotes foam cell accumulation and accelerates atherosclerosis in mice. J Clin Invest. 2007, 117: 3900-3908.
Yang X, Lu X, Lombès M, Rha GB, Chi YI, Guerin TM, Smart EJ, Liu J: The G(0)/G(1) switch gene 2 regulates adipose lipolysis through association with adipose triglyceride lipase. Cell Metab. 2010, 11: 194-205. 10.1016/j.cmet.2010.02.003.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-9. 10.1038/75556.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: 3-10.1186/gb-2003-4-5-p3.
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-
Pirooznia M, Habib T, Perkins EJ, Deng Y: GOfetcher: a database with complex searching facility for gene ontology. Bioinformatics. 2008, 24: 2561-10.1093/bioinformatics/btn441.
Bredel M, Bredel C, Juric D, Harsh GR, Vogel H: Functional network analysis reveals extended gliomagenesis pathway maps and three novel MYC-interacting genes in human gliomas. Cancer Res. 2005, 65: 8679-10.1158/0008-5472.CAN-05-1204.
This study was supported by NIH/NCMHD/RIMI, P20MD002725. We thank Yan Li from Rush University Medical Center for carefully reading the manuscript.
The authors declare that they have no competing interests.
JM conceived of the study, and contributed to its design and coordination and drafted the manuscript; JA analyzed microarray data and contributed to the writing of the manuscript. XZ made substantial contributions to conception and design, and has been involved in revising it critically for important intellectual content. MB and MS processed the samples, isolated total RNA, and conducted the array hybridizations and part of the data organization. MO provided clinical samples and other laboratory data. JTW participated in clinical data analysis and English proofreading. YD initiated the project, analyzed microarray data and revised the manuscript. All authors approved the final manuscript. XW participated in clinical discussion and helped to improve the manuscript.
About this article
Cite this article
Mao, J., Ai, J., Zhou, X. et al. Transcriptomic profiles of peripheral white blood cells in type II diabetes and racial differences in expression profiles. BMC Genomics 12, S12 (2011). https://doi.org/10.1186/1471-2164-12-S5-S12
- African American
- African American
- Interferon Signaling
- Peripheral White Blood Cell
- Adaptive Immune Cell