Open Access

Global gene expression profiling displays a network of dysregulated genes in non-atherosclerotic arterial tissue from patients with type 2 diabetes

  • Vibe Skov1, 2Email author,
  • Steen Knudsen3,
  • Malene Olesen1,
  • Maria L Hansen1, 4 and
  • Lars M Rasmussen1
Cardiovascular Diabetology201211:15

DOI: 10.1186/1475-2840-11-15

Received: 14 November 2011

Accepted: 17 February 2012

Published: 17 February 2012

Abstract

Background

Generalized arterial alterations, such as endothelial dysfunction, medial matrix accumulations, and calcifications are associated with type 2 diabetes (T2D). These changes may render the vessel wall more susceptible to injury; however, the molecular characteristics of such diffuse pre-atherosclerotic changes in diabetes are only superficially known.

Methods

To identify the molecular alterations of the generalized arterial disease in T2D, DNA microarrays were applied to examine gene expression changes in normal-appearing, non-atherosclerotic arterial tissue from 10 diabetic and 11 age-matched non-diabetic men scheduled for a coronary by-pass operation. Gene expression changes were integrated with GO-Elite, GSEA, and Cytoscape to identify significant biological pathways and networks.

Results

Global pathway analysis revealed differential expression of gene-sets representing matrix metabolism, triglyceride synthesis, inflammation, insulin signaling, and apoptosis. The network analysis showed a significant cluster of dysregulated genes coding for both intra- and extra-cellular proteins associated with vascular cell functions together with genes related to insulin signaling and matrix remodeling.

Conclusions

Our results identify pathways and networks involved in the diffuse vasculopathy present in non-atherosclerotic arterial tissue in patients with T2D and confirmed previously observed mRNA-alterations. These abnormalities may play a role for the arterial response to injury and putatively for the accelerated atherogenesis among patients with diabetes.

Keywords

Systems biology Microarray Diabetes mellitus Gene expression Coronary artery disease

Background

Cardiovascular diseases (CVD) in patients with type 2 diabetes are a large and increasing health problem. Increased atherosclerotic lesions are believed to form the basis behind the high frequency of CVD in diabetes; however, epidemiological studies have shown that traditional risk factors, e.g. hypertension and dyslipidemia cannot explain the increased incidence [1]. Unfortunately, the molecular mechanisms leading to atherosclerosis in diabetes are only partially understood [2].

The arterial wall in diabetes harbors not only increased amounts of atherosclerotic plaques, but also diffuse alterations present in non-atherosclerotic parts of the vessel wall. One element of the generalized alterations in the vasculature in diabetes is endothelial dysfunction [3], characterized by increased permeability [4], increased expression of pro-inflammatory molecules [5], and altered vasomotoric responses [6]. Moreover, changes in extracellular matrix components of the tunica media are present in both atherosclerotic and non-atherosclerotic parts of the arterial tree in diabetes. Increased concentrations of collagen type 4 [7], hyaluronic acid [8], osteopontin, osteoprotegerin [9], and metalloproteinases [10] have, for example, been described in conjunction with the presence of high amounts of glucose-derived increased cross-linking of collagens [11]. Decreased amounts of several gene products related to apoptosis have been observed in vascular smooth muscle cells from normal appearing areas of arteries from patients with diabetes [12]. In addition, linear media calcifications occur with increased frequency among patients with glucose intolerance and diabetes and are strong predictors of CVD in these individuals [13, 14]. In accordance, recent studies of non-atherosclerotic arterial alterations in animal models of type 2 diabetes and hyperglycemia demonstrated increased aortic stiffness and upregulation of matrix components [15], increased arterial calcification [16], and accumulation of glycosaminoglycan-rich material [17]. Thus, defects in several molecular pathways seem to be present in the arterial wall in patients with type 2 diabetes. These changes are likely to play important roles in the arterial response to injury and thus in the build-up of atherosclerotic plaques in diabetic patients.

By using DNA microarrays for transcriptional profiling, a large number of genes can be analyzed simultaneously and the resulting data can be integrated with pathways and biological interaction networks to detect coordinated changes in functionally related genes. Recently, we applied gene expression microarrays and identified fibulin-1 (FBLN1) as a molecular marker of cardiovascular disease in non-atherosclerotic tissue from patients with type 2 diabetes [18]. To our knowledge, only one other study has examined global gene expression changes in arterial tissue from diabetic patients; however, the results only reflected different degrees of atherosclerosis and not specific diabetes related alterations [19].

In the present study, we hypothesized that diffuse diabetes-related alterations in non-atherosclerotic arterial tissue are associated with a cluster of dysregulated genes involved in vascular and metabolic cell functions. We applied global pathway and network analysis to unravel the gene expression profile of well-defined, normal-appearing arterial tissue from patients with type 2 diabetes.

Methods

Subjects

Study subjects were recruited from patients undergoing artery by-pass graft surgery at Skejby and Odense University Hospital. Participants included ten men with more than 2 years known duration of type 2 diabetes, defined as diabetes diagnosed after adolescence, and eleven non-diabetic men matched according to age (Table 1). Participants and tissue handling have been described in a study aimed at identifying single genes differentially expressed in type 2 diabetes [18]. Briefly, a small piece of arteria mammaria interna from all patients was dissected into intima-media and surrounding tissue and frozen in liquid nitrogen within a few minutes. The arteries are spare material from the repair vessel from coronary by-pass operations. The investigation was approved by the local Danish ethical committee in Region Mid and South and informed consent was obtained from each subject.
Table 1

Characteristics of diabetic patients and control subjects

 

Diabetic patients

Control subjects

N

10

11

Age (years)

65.6 ± 11.9

65.2 ± 8.5

Systolic blood pressure (mmHg)

143.1 ± 35.6

153.9 ± 16.6

Diastolic blood pressure (mmHg)

76.3 ± 11.6

87.9 ± 8.7

Hemoglobin A1c (fract.)

0.074 ± 0.01

0.060 ± 0.004††

Total cholesterol (mmol/l)

4.1 ± 0.8

4.7 ± 0.5

High Density Lipoprotein (mmol/l)

1.2 ± 0.2

1.4 ± 0.3

Data represent means ± SD; Students T-test for non-paired data used. P < 0.05 and ††P < 0.005 for diabetic patients vs. control subjects

Histology

Tissue for histology was fixed for 24-72 hours in PBS-buffered 4% formaldehyde and subsequently embedded in paraffin. Sections were cut and stained by hematoxylin-eosin, von kossa, and elastin using routine methods [18].

RNA isolation and microarray hybridization

Total RNA was isolated and hybridized to microarrays as previously described [18]. Briefly, total RNA was isolated from intima-media samples by using the TRIzol protocol (Life Technologies, Gaithersburg, MD, USA), and purified RNA was converted to biotin-labeled amplified RNA (aRNA) and hybridized to Affymetrix HG-U133A 2.0 chips (Affymetrix, Santa Clara, CA, USA).

Data processing and statistical analysis

All microarray chips were preprocessed as reported previously [18]. The regularized t-test limma http://www.bioconductor.org was applied to evaluate differences in gene expression between patients with diabetes and control subjects. Genes with an uncorrected P < 0.05 were considered for further analysis.

Global pathway and network analysis

Changes in gene expression at the pathway level between diabetic patients and control subjects were assessed using GO-Elite 1.19 http://www.genmapp.org/go_elite/ and Gene Set Enrichment Analysis (GSEA) 2.0.1 [20]. GO-Elite is a software tool that runs an updated version of MAPPFinder [21] and filters redundant pathways. A total of 113 pathways were analyzed using GO-Elite. A z-score was calculated for each pathway, and 2000 permutations were applied to assess the significance of the z-score http://www.genmapp.org. Applying the MSigDB database http://www.broad.mit.edu/gsea/msigdb in GSEA, we compiled an "artery gene set" consisting of 820 gene sets and pathways. All genes were ranked according to the ttest and an enrichment score (ES) was assigned to each gene set. The statistical significance of the ES was estimated by performing 2000 gene permutations. After gene set enrichment analysis in GSEA, we masked the redundant gene sets by using the leading edge analysis option http://www.broad.mit.edu/gsea.

Gene expression changes between diabetic patients and control subjects were integrated with molecular interaction networks using Cytoscape 2.6.0 [22]. The Human Interactome data set http://www.cytoscape.org/cgi-bin/moin.cgi/Data_Sets, composed of 61,263 protein-protein interactions analogous to 10,000 different proteins, was loaded into the Cytoscape software. To reduce the complexity of the large network, differentially expressed genes from the data set (P < 0.05) were mapped onto the network, and then the filter option in Cytoscape was used to select proteins (nodes) interacting directly with the differentially expressed genes generating a new network. To identify hubs--i.e. interacting nodes with high network degrees assumed to correspond to essential gene products in the cell [23]--with 15 or more interactions, the topology filter option was applied to the new network creating a smaller network. The jActiveModules plugin in Cytoscape was applied to the small network to identify a subnetwork with an overall significant response to the experimental condition. The anneal algorithm and 250.000 iterations were used in the jActiveModules plugin to identify the subnetwork. A z-score was calculated for the subnetwork and the significance of the z-score was assessed by performing 1000 permutations. NetworkAnalyzer is a core plugin in Cytoscape and has been used to calculate topological values for each gene in the subnetwork.

Results

The mean age of the study subjects was 65 years (Table 1). Patients with type 2 diabetes had elevated hemoglobin A1c (P < 0.005) and reduced diastolic blood pressure (P < 0.05). There was no difference between the two groups with respect to systolic blood pressure, HDL, and total cholesterol as previously reported [18].

Histology

Tissue from all subjects appeared normal, i.e. without atherosclerosis, calcified areas, or cellular infiltrations, as revealed by von kossa, hematoxylin-eosin, and elastin staining as previously reported [18]. In a few instances, enough tissue was obtained to make additional histology after dissection to control the separation of intima-media from adventitia. In these cases, dissection was made just outside the lamina elastica externa.

Global gene expression analysis

Of the 22,277 probe sets represented on the array, 477 were downregulated and 169 were upregulated in type 2 diabetic patients versus control subjects as reported previously (P < 0.05) [18]. No probe sets were significantly differentially expressed after controlling for multiple hypothesis testing using the Benjamini-Hochberg method [24].

Analysis of gene sets and biological pathways

To examine the difference in gene expression between the two groups at the pathway level, data were loaded into GSEA 2.0.1 and GO-Elite 1.19. Applying GSEA, caries pulp up (genes upregulated in pulpal tissue from carious teeth) [25], neuroactive ligand receptor interaction, cell adhesion molecules, monocyte dend DN (genes downregulated in monocyte dendritic cells), and cytokine-cytokine receptor interaction were the most significantly upregulated gene sets (FWER < 0.05). The most significantly downregulated gene sets included kidney transplant well UP [26], insulin 2F UP (genes 2 fold upregulated in muscle by insulin) [27], heartfailure atria DN, GH GHRHR KO 24 hrs DN [28], proteasome, and diab neph DN (genes downregulated in the glomeruli of cadaver kidneys from patients with diabetic nephropathy) (FWER < 0.05) (Table 2). Further description of the gene sets can be found in the MsigDB database http://www.broad.mit.edu/gsea/msigdb.
Table 2

Ranking of the top 10 up- and downregulated gene sets using GSEA in intima-media of type 2 diabetic patients vs. control subjects

GENE SET

SIZE

ES

NES

NOM

FDR

FWER

    

p-value

p-value

p-value

Upregulated gene sets

      

Caries pulp UP

202

-0.43

-2.29

< 0.0001

0.0006

0.0005

HSA04080 Neuroactive ligand receptor interaction

227

-0.42

-2.27

< 0.0001

0.0009

0.002

HSA04514 Cell adhesion molecules

114

-0.43

-2.18

< 0.0001

0.002

0.005

Monocyte dend DN

115

-0.43

-2.13

< 0.0001

0.003

0.01

HSA04060 Cytokine-cytokine receptor interaction

225

-0.39

-2.12

< 0.0001

0.003

0.01

Kidney transplant rejection UP

80

-0.45

-2.09

< 0.0001

0.004

0.02

HSA04610 Complement and coagulation cascades

63

-0.45

-2.05

< 0.0001

0.007

0.04

HSA01430 Cell communication

108

-0.41

-2.02

< 0.0001

0.008

0.05

HSA04640 Hematopoietic cell lineage

83

-0.43

-2.01

< 0.0001

0.008

0.06

Osteoclasts sig

37

-0.50

-1.94

< 0.0001

0.02

0.13

Downregulated gene sets

      

Kidney transplant well UP

251

0.56

2.45

< 0.0001

< 0.0001

< 0.0001

Insulin 2F UP

167

0.57

2.41

< 0.0001

< 0.0001

< 0.0001

Heartfailure atria DN

109

0.55

2.20

< 0.0001

0.0002

0.0005

GH GHRHR KO 24 hrs DN

140

0.53

2.18

< 0.0001

0.0001

0.0005

HSA03050 Proteasome

21

0.73

2.16

< 0.0001

0.0002

0.001

Diab neph DN

249

0.49

2.15

< 0.0001

0.0002

0.001

HSA00190 Oxidative phosphorylation

99

0.55

2.15

< 0.0001

0.0001

0.001

Kidney transplant well PBL DN

40

0.61

2.09

< 0.0001

0.0005

0.004

Krebs TCA cycle

29

0.65

2.07

< 0.0001

0.0008

0.008

mRNA processing

41

0.61

2.07

0.0007

0.0008

0.009

Ranking of the gene sets was done using GSEA 2.0.1. ES, enrichment score; NES, normalized enrichment score; NOM, nominal; FDR, false discovery rate; FWER, family wise error rate

The statistical rating of the results in GO-Elite was provided by the z-score, and pathways with a z-score above 2.0 were considered significant [21]. Evaluation of the results revealed significant upregulation of triacylglyceride synthesis, apoptosis, matrix metalloproteinases, and adipogenesis (z > 2.0). Five pathways including cholesterol biosynthesis, TNF-alpha NFkB Netpath 9, Krebs TCA Cycle, non-homologous end joining, and signaling of hepatocyte growth factor receptor were significantly downregulated (z > 2.0) (Table 3).
Table 3

Ranking of up- and downregulated pathways in intima-media of type 2 diabetic patients vs. control subjects using GO-Elite 1.19

MAPP Name

Changed

Measured

On MAPP

Changed

Z

Permute

FWER

 

(n)

(n)

(n)

(%)

score

p-value

p-value

Upregulated pathways

       

Triacylglyceride synthesis

3

22

24

13.64

4.7

0.002

0.23

Apoptosis

4

83

85

4.82

2.5

0.04

0.67

Matrix metalloproteinases

2

30

32

6.67

2.4

0.04

0.67

Adipogenesis

5

129

133

3.88

2.3

0.04

0.67

Downregulated pathways

       

Cholesterol biosynthesis

4

15

15

26.67

4.8

0.004

0.15

TNFalpha NFkB netpath 9

16

187

208

8.56

3.8

0.003

0.15

Krebs TCA cycle

5

32

39

15.63

3.7

0.007

0.18

Non-homologous end joining

2

7

8

28.57

3.6

0.04

0.61

Signaling of hepatocyte growth factor receptor

5

35

36

14.29

3.5

0.02

0.39

Shown are significantly up- and downregulated pathways in intima-media of diabetic patients vs. control subjects. Only non-redundant pathways with a z-score above 2.0 are listed. A foldchange of 1.05 and a p-value < 0.05 were used as filtering criteria. Changed (n): number of changed genes, Measured (n): number of genes on the chip, On MAPP (n): number of genes on the MAPP, Changed (%): percentage of changed genes, FWER: family wise error rate

Computational analysis of the Human Interactome network

To further explore the global gene expression changes associated with arterial disease in type 2 diabetes, data were integrated with the Human interactome network to identify the principal subnetworks i.e. connected regions of the network that show significant changes in expression. 397 differentially expressed genes (P < 0.05) from the microarray data set interacted directly with 2980 proteins (nodes) in the network creating a network composed of 3227 interacting nodes. Filtering of hubs with 15 or more interactions resulted in a network consisting of 1043 nodes. Performing jActiveModules on the 1043 node network unveiled 5 highly overlapping subnetworks, therefore, we decided only to focus on the top network consisting of 74 nodes (z = 9.23) (permute P = 0.01). The statistically significant 74 node subnetwork is depicted in Figure 1. Topological values and number of interactors for all genes are shown in Table 4. The topological coefficient is a relative measure for the extent to which a node shares neighbors with other nodes. The most important cardiometabolic hubs and their interactors are shown in Table 5.
https://static-content.springer.com/image/art%3A10.1186%2F1475-2840-11-15/MediaObjects/12933_2011_Article_452_Fig1_HTML.jpg
Figure 1

Interaction subnetwork of potential arterial disease-associated genes in patients with type 2 diabetes. The subnetwork consists of 74 nodes and 166 edges. Red nodes represent upregulated genes, green nodes represent downregulated genes (P < 0.05), and white nodes represent genes with P > 0.05. The z-score of the network is 9.23 (permute P = 0.01)

Table 4

Topological values and number of interactors for all genes in the network

Gene

TPV

NOI

Gene

TPV

NOI

Gene

TPV

NOI

ACTR3

0.625

2

FURIN

0.52

3

PTK2

0.20

8

ACVR1

0.7

2

FYN

0.27

6

PTMA

0

1

ARF1

0.5

2

GGA2

0

1

PTPN12

0.51

3

BCR

0.26

6

GRB10

0.39

4

PTPRF

0.36

3

BIRC4

0.33

3

IGF2

0.5

2

RABAC1

0.25

4

C1orf103

0

1

IGFBP3

0.33

3

RAP1A

0.33

3

CANX

0.27

5

INSR

0.17

10

RASA1

0.22

5

CASP10

0.33

3

IRS2

0.16

8

RHOA

0.33

3

CASP4

0

1

ITGAV

0.39

3

RIF1

0.5

2

CASP7

0.13

8

JAK2

0.24

8

RNF11

0

1

CD82

0.71

2

KRAS

0.59

2

ROCK1

0.5

2

CDC42

0.16

8

LRP1

0.32

5

RPS6KB1

0.57

2

DAB2

0

1

LYN

0.31

5

RRAS2

0.3

4

DCN

0.5

2

MAP3K7IP2

0

1

SMAD4

0.23

5

DVL2

0

1

MAX

0.25

4

SMURF2

0.24

5

EIF1B

0.31

4

MBIP

0

1

STK11

0.63

2

ELN

0.33

3

MMP14

0.41

3

STK24

0.5

2

EPOR

0.42

3

MMP2

0.32

4

SYNCRIP

0.5

2

F2

0.41

4

NCK1

0.39

3

TBK1

0.5

2

F8

0.31

5

NRAS

0.25

5

UBTF

0

1

FAM107A

0

1

OAT

0.63

2

UNC119

0.33

4

FAS

0.25

4

PARD3

0.24

5

VWF

0.57

3

FBLN1

0.5

2

PLEKHF2

0.25

4

YWHAQ

0.55

2

FBLN2

0

1

PPP1CB

0.5

2

YWHAZ

0.36

3

FRS2

0.39

3

PRKAR1A

0

1

   

TPV: topological value; NOI: number of interactors. The TPV is a number between 0-1. Nodes that have one neighbor are assigned a topological value of 0

Table 5

Important cardiometabolic hubs and their interactors in non-atherosclerotic diabetic arterial tissue

Hub

Interactor

Hub

Interactor

Hub

Interactor

Hub

Interactor

INSR

ARF1

CDC42

ACTR3

IRS2

EPOR

SMAD4

FBLN1

 

FRS2

 

C1ORF103

 

INSR

 

MAX

 

GRB10

 

CASP7

 

KRAS

 

NRAS

 

IGF2

 

OAT

 

NRAS

 

PARD3

 

IRS2

 

PARD3

 

PTPRF

 

RRAS2

 

JAK2

 

RIF1

 

UBTF

MMP2

DCN

 

PTK2

 

RPS6KB1

 

YWHAQ

 

IGFBP3

 

PTPN12

 

UNC119

 

YWHAZ

 

ITGAV

 

RASA1

     

MMP14

 

SYNCRIP

      

Discussion

In the present study, we used high-density oligonucleotide arrays to elucidate global gene expression patterns in non-atherosclerotic, non-calcified, normal-appearing arterial tissue from patients with type 2 diabetes. Previously, we used these data in a study aimed at identifying new molecular markers of arterial disease in type 2 diabetes. Examinations were done at the single gene level with microarrays, and we were able to validate several genes with q-RT-PCR including FBLN1 as significantly upregulated in the vascular wall both at the mRNA and protein level [18]. In the present bioinformatic study, we use our mRNA expression data to explore biological pathways and networks that are dysregulated in the arterial wall in diabetes. Although our data indicate that differences in gene expression between arterial tissue from patients with diabetes and non-diabetic subjects are modest at the single gene level [18], we demonstrate that clusters of highly interconnected genes are significantly dysregulated at the transcriptional level in arterial tissue in diabetes. Providing credibility to our results, a number of gene products of these pathways has previously been found dysregulated in the vessel wall in diabetic vasculopathy, as discussed below.

The results from GSEA demonstrate significant upregulation of gene sets and pathways related to cytokine-growth factor- and hormonal actions in arterial tissue from diabetic individuals. This include the gene sets caries pulp up, cell adhesion molecules, monocyte dend DN, and cytokine-cytokine receptor interaction with caries pulp up as the most significantly upregulated gene set. Interestingly, there is an increasing evidence supporting an association between periodontitis and diabetic complications [29, 30] suggesting that genes related to periodontitis could play a role for the development of diabetic vasculopathy. Because the tissue in our study is without atherosclerosis and cellular infiltration, it seems that the vascular smooth muscle cells themselves may express an inflammatory phenotype. In agreement, previous observations showed that vascular smooth muscle cells produce a range of peptide factors that have been suggested to play a role in diabetic vasculopathy [9]. In addition, this finding is in line with previous observations that vascular cells in diabetes may display increased inflammatory capacity, at least partly related to effects of toll-like receptors [31, 32].

Over the past decades there has been much debate regarding the relative importance of hyperinsulinism and insulin resistance in vascular cells from diabetic patients [3336]. Interestingly, the insulin 2F-pathway, representing genes upregulated 2-fold by insulin stimulation of skeletal muscle in healthy subjects [27], was significantly downregulated in our study, supporting the idea that at least some effects of insulin in vascular tissue are influenced by insulin resistance and not by hyperinsulinism. This is in agreement with several human and experimental observations [3739].

Surprisingly, a triglyceride pathway was the most significantly upregulated pathway using GO-Elite. Several organs, such as muscle and liver, are known to accumulate triglycerides in type 2 diabetes [40]; however, no information concerning arterial tissue is currently available. Our data are compatible with the hypothesis that accumulation of triglycerides could also take place in the arterial wall.

Since the apoptotic process has been found to be regulated in vascular smooth muscle cells from arterial tissue from individuals with type 2 diabetes [12], it is of interest that our data demonstrate an increased expression of the apoptotic pathway using GO-Elite. The apoptotic pathway may be regulated in arterial tissue in diabetes; however, further studies are warranted to assess the precise implication of apoptosis in the pathogenesis of diabetic arteriopathy. Of interest, also a pathway defined as genes downregulated in glomeruli from patients with diabetic nephropathy was downregulated in the arterial wall in diabetes [41]. This observation may be an indication of a common set of dysregulated genes present in both micro- and macro-angiopathy in diabetes.

Although many detailed unique processes are suggested to be implicated in the development of arterial disease in diabetes, an understanding of the disease mechanisms as an integrated whole may be warranted. Network analysis provides a broad insight into biology in the context of known functional interrelationships among proteins. Our network analysis demonstrates a statistically significant cluster of genes that are dysregulated in the arterial wall in diabetes.

Looking for important hubs in the network, it is interesting to note that the insulin receptor (INSR) appears, although it is not itself regulated. It seems that many genes, interacting with the receptor, are indeed regulated in the arterial wall in diabetes, which is in accordance with the idea that dysfunctional effects of insulin may play a role in diabetic arterial disease (35-37). Another relevant hub is SMAD4, which seems to have an important place in the network, interacting both with intracellularly expressed genes like RRAS2 and PARD3, as well as an extracellular matrix gene, fibulin-1 (FBLN1). SMAD4 is an important intracellular signaling molecule in the TGF-beta system which is in agreement with other observations showing that this system is involved in the development of matrix accumulations in diabetic complications [42, 43]. A recent transcriptome analysis of human diabetic kidney disease pointed towards CDC42 signaling as an important dysregulated pathway [44], which is in accordance with the presence of CDC42 as an important hub in the network we present here. MMP2 is an important hub in the network with connections to other relevant genes including MMP14, DCN and ITGAV. This part of the network is in agreement with other previously reported data showing that MMP2 is dysregulated in the arterial wall in diabetes [45] and seems to indicate that matrix remodeling may be an important feature of the non-atherosclerotic arterial disease seen in diabetes.

Conclusions

We used transcriptional profiling on well-defined non-atherosclerotic arterial samples from diabetic individuals. Using pathway and network analysis, our data display a statistically significant cluster of dysregulated genes in the arteries of diabetic patients, which is in accordance with the presence of a diffuse diabetic macroangiopathy, similar to the diabetic microangiopathy. Our approach has not previously been used, but point towards dysregulated pathways related to matrix metabolism, triglyceride synthesis, inflammation, as well as insulin signaling and apoptosis. Dysregulated gene interactions and pathways in the cells of the arterial wall in diabetes may play important roles in the arterial response to injury and atherosclerosis.

Data deposition footnote

Data are available from GEO (http://www.ncbi.nlm.nih.gov/geo, Accession No. GSE13760).

Abbreviations

GSEA: 

Gene Set Enrichment Analysis

ES: 

Enrichment score

VWF: 

von Willebrand factor

LEP: 

Leptin

MMP2: 

Matrix metallopeptidase 2

FBLN1: 

Fibulin1

FBLN2: 

Fibulin2

EGR3: 

Early growth factor 3

HBA2: 

Hemoglobin alpha 2

FAS: 

TNF receptor superfamily, member 6

TNFAIP6: 

Tumor necrosis factor, alpha-induced protein 6

ELN: 

Elastin

Declarations

Acknowledgements

The technical assistance of Lone Sørensen and Anne-Marie Jakobsen is greatly appreciated. This study was supported by the NOVO Nordisk Foundation, the Danish Medical Research Council, the Danish Diabetes Association, the Augustinus Foundation, Carl and Ellen Hertz's Foundation, Engineer K.A. Rohde and wives Foundation, the A.P. Møller Foundation, Th. Maigaards Eftf. Mrs. Lily Benthine Lunds Fond of the 1 June 1978, and the Bagger-Sorensen Foundation.

Authors’ Affiliations

(1)
Department of Clinical Biochemistry and Pharmacology, Odense University Hospital
(2)
Department of Clinical Genetics, Odense University Hospital
(3)
Medical Prognosis Institute A/S
(4)
Department of Cardiothoracic and Vascular Surgery, Odense University Hospital

References

  1. Kannel WB, McGee DL: Diabetes and glucose tolerance as risk factors for cardiovascular disease: the Framingham study. Diabetes Care. 1979, 2: 120-126. 10.2337/diacare.2.2.120.View ArticlePubMed
  2. Goldberg IJ: Why does diabetes increase atherosclerosis? I don't know!. J Clin Invest. 2004, 114: 613-615.PubMed CentralView ArticlePubMed
  3. De Vriese AS, Verbeuren TJ, van d V, Lameire NH, Vanhoutte PM: Endothelial dysfunction in diabetes. Br J Pharmacol. 2000, 130: 963-974. 10.1038/sj.bjp.0703393.PubMed CentralView ArticlePubMed
  4. Feldt-Rasmussen B: Microalbuminuria, endothelial dysfunction and cardiovascular risk. Diabetes Metab. 2000, 26 (Suppl 4): 64-66.PubMed
  5. Ribau JC, Hadcock SJ, Teoh K, DeReske M, Richardson M: Endothelial adhesion molecule expression is enhanced in the aorta and internal mammary artery of diabetic patients. J Surg Res. 1999, 85: 225-233. 10.1006/jsre.1999.5682.View ArticlePubMed
  6. Rask-Madsen C, Ihlemann N, Krarup T, Christiansen E, Kober L, Nervil KC, Torp-Pedersen C: Insulin therapy improves insulin-stimulated endothelial function in patients with type 2 diabetes and ischemic heart disease. Diabetes. 2001, 50: 2611-2618. 10.2337/diabetes.50.11.2611.View ArticlePubMed
  7. Rasmussen LM, Ledet T: Aortic collagen alterations in human diabetes mellitus. Changes in basement membrane collagen content and in the susceptibility of total collagen to cyanogen bromide solubilisation. Diabetologia. 1993, 36: 445-453. 10.1007/BF00402282.View ArticlePubMed
  8. Heickendorff L, Ledet T, Rasmussen LM: Glycosaminoglycans in the human aorta in diabetes mellitus: a study of tunica media from areas with and without atherosclerotic plaque. Diabetologia. 1994, 37: 286-292. 10.1007/BF00398056.View ArticlePubMed
  9. Olesen P, Ledet T, Rasmussen LM: Arterial osteoprotegerin: increased amounts in diabetes and modifiable synthesis from vascular smooth muscle cells by insulin and TNF-alpha. Diabetologia. 2005, 48: 561-568. 10.1007/s00125-004-1652-8.View ArticlePubMed
  10. Chung AW, Booth AD, Rose C, Thompson CR, Levin A, van BC: Increased matrix metalloproteinase 2 activity in the human internal mammary artery is associated with ageing, hypertension, diabetes and kidney dysfunction. J Vasc Res. 2008, 45: 357-362. 10.1159/000119755.View ArticlePubMed
  11. Sims TJ, Rasmussen LM, Oxlund H, Bailey AJ: The role of glycation cross-links in diabetic vascular stiffening. Diabetologia. 1996, 39: 946-951. 10.1007/BF00403914.View ArticlePubMed
  12. Ruiz E, Gordillo-Moscoso A, Padilla E, Redondo S, Rodriguez E, Reguillo F, Briones AM, van BC, Okon E, Tejerina T: Human vascular smooth muscle cells from diabetic patients are resistant to induced apoptosis due to high Bcl-2 expression. Diabetes. 2006, 55: 1243-1251. 10.2337/db05-0949.View ArticlePubMed
  13. Neubauer B: A quantitative study of peripheral arterial calcification and glucose tolerance in elderly diabetics and non-diabetics. Diabetologia. 1971, 7: 409-413. 10.1007/BF01212055.View ArticlePubMed
  14. Lehto S, Niskanen L, Suhonen M, Ronnemaa T, Laakso M: Medial artery calcification. A neglected harbinger of cardiovascular complications in non-insulin-dependent diabetes mellitus. Arterioscler Thromb Vasc Biol. 1996, 16: 978-983. 10.1161/01.ATV.16.8.978.View ArticlePubMed
  15. Sista AK, O'Connell MK, Hinohara T, Oommen SS, Fenster BE, Glassford AJ, Schwartz EA, Taylor CA, Reaven GM, Tsao PS: Increased aortic stiffness in the insulin-resistant Zucker fa/fa rat. Am J Physiol Heart Circ Physiol. 2005, 289: H845-H851. 10.1152/ajpheart.00134.2005.View ArticlePubMed
  16. Heinonen SE, Leppanen P, Kholova I, Lumivuori H, Hakkinen SK, Bosch F, Laakso M, Yla-Herttuala S: Increased atherosclerotic lesion calcification in a novel mouse model combining insulin resistance, hyperglycemia, and hypercholesterolemia. Circ Res. 2007, 101: 1058-1067. 10.1161/CIRCRESAHA.107.154401.View ArticlePubMed
  17. Renard CB, Kramer F, Johansson F, Lamharzi N, Tannock LR, von Herrath MG, Chait A, Bornfeldt KE: Diabetes and diabetes-associated lipid abnormalities have distinct effects on initiation and progression of atherosclerotic lesions. J Clin Invest. 2004, 114: 659-668.PubMed CentralView ArticlePubMed
  18. Cangemi C, Skov V, Poulsen MK, Funder J, Twal WO, Gall MA, Hjortdal V, Jespersen ML, Kruse TA, Aagard J, Parving HH, Knudsen S, Hoilund-Carlsen PF, Rossing P, Henriksen JE, Argraves WS, Rasmussen LM: Fibulin-1 is a Marker for Arterial Extracellular Matrix Alterations in Type 2 Diabetes. Clin Chem. 2011, 57: 1556-1565. 10.1373/clinchem.2011.162966.View ArticlePubMed
  19. Fassot C, Briet M, Rostagno P, Barbry P, Perret C, Laude D, Boutouyrie P, Bozec E, Bruneval P, Latremouille C, Laurent S: Accelerated arterial stiffening and gene expression profile of the aorta in patients with coronary artery disease. J Hypertens. 2008, 26: 747-757. 10.1097/HJH.0b013e3282f4b3d0.View ArticlePubMed
  20. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.PubMed CentralView ArticlePubMed
  21. Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR: MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol. 2003, 4: R7-10.1186/gb-2003-4-1-r7.PubMed CentralView ArticlePubMed
  22. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMed CentralView ArticlePubMed
  23. Barabasi AL, Oltvai ZN: Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004, 5: 101-113. 10.1038/nrg1272.View ArticlePubMed
  24. Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995, 57: 289-300.
  25. McLachlan JL, Smith AJ, Bujalska IJ, Cooper PR: Gene expression profiling of pulpal tissue reveals the molecular complexity of dental caries. Biochim Biophys Acta. 2005, 1741: 271-281.View ArticlePubMed
  26. Flechner SM, Kurian SM, Head SR, Sharp SM, Whisenant TC, Zhang J, Chismar JD, Horvath S, Mondala T, Gilmartin T, Cook DJ, Kay SA, Walker JR, Salomon DR: Kidney transplant rejection and tissue injury by gene profiling of biopsies and peripheral blood lymphocytes. Am J Transplant. 2004, 4: 1475-1489. 10.1111/j.1600-6143.2004.00526.x.PubMed CentralView ArticlePubMed
  27. Rome S, Clement K, Rabasa-Lhoret R, Loizon E, Poitou C, Barsh GS, Riou JP, Laville M, Vidal H: Microarray profiling of human skeletal muscle reveals that insulin regulates approximately 800 genes during a hyperinsulinemic clamp. J Biol Chem. 2003, 278: 18063-18068. 10.1074/jbc.M300293200.View ArticlePubMed
  28. Govoni KE, Lee SK, Chadwick RB, Yu H, Kasukawa Y, Baylink DJ, Mohan S: Whole genome microarray analysis of growth hormone-induced gene expression in bone: T-box3, a novel transcription factor, regulates osteoblast proliferation. Am J Physiol Endocrinol Metab. 2006, 291: E128-E136. 10.1152/ajpendo.00592.2005.PubMed CentralView ArticlePubMed
  29. Lazenby MG, Crook MA: The innate immune system and diabetes mellitus: the relevance of periodontitis? A hypothesis. Clin Sci (Lond). 2010, 119: 423-429. 10.1042/CS20100098.View Article
  30. Campus G, Salem A, Uzzau S, Baldoni E, Tonolo G: Diabetes and periodontal disease: a case-control study. J Periodontol. 2005, 76: 418-425. 10.1902/jop.2005.76.3.418.View ArticlePubMed
  31. Li J, Jin C, Cleveland JC, Ao L, Xu D, Fullerton DA, Meng X: Enhanced inflammatory responses to toll-like receptor 2/4 stimulation in type 1 diabetic coronary artery endothelial cells: the effect of insulin. Cardiovasc Diabetol. 2010, 9: 90-10.1186/1475-2840-9-90.PubMed CentralView ArticlePubMed
  32. Erridge C, Samani NJ: Saturated fatty acids do not directly stimulate Toll-like receptor signaling. Arterioscler Thromb Vasc Biol. 2009, 29: 1944-1949. 10.1161/ATVBAHA.109.194050.View ArticlePubMed
  33. Erikstrup C, Pedersen LM, Heickendorff L, Ledet T, Rasmussen LM: Production of hyaluronan and chondroitin sulphate proteoglycans from human arterial smooth muscle--the effect of glucose, insulin, IGF-I or growth hormone. Eur J Endocrinol. 2001, 145: 193-198. 10.1530/eje.0.1450193.View ArticlePubMed
  34. Stout RW: Development of vascular lesions in insulin-treated animals fed a normal diet. Br Med J. 1970, 3: 685-687. 10.1136/bmj.3.5724.685.PubMed CentralView ArticlePubMed
  35. Stout RW: Insulin and atherosclerosis. Lancet. 1969, 2: 327-328.View ArticlePubMed
  36. Bansilal S, Farkouh ME, Fuster V: Role of insulin resistance and hyperglycemia in the development of atherosclerosis. Am J Cardiol. 2007, 99: 6B-14B.View ArticlePubMed
  37. Rask-Madsen C, Ihlemann N, Krarup T, Christiansen E, Kober L, Nervil KC, Torp-Pedersen C: Insulin therapy improves insulin-stimulated endothelial function in patients with type 2 diabetes and ischemic heart disease. Diabetes. 2001, 50: 2611-2618. 10.2337/diabetes.50.11.2611.View ArticlePubMed
  38. Jiang ZY, Lin YW, Clemont A, Feener EP, Hein KD, Igarashi M, Yamauchi T, White MF, King GL: Characterization of selective resistance to insulin signaling in the vasculature of obese Zucker (fa/fa) rats. J Clin Invest. 1999, 104: 447-457. 10.1172/JCI5971.PubMed CentralView ArticlePubMed
  39. Rask-Madsen C, Li Q, Freund B, Feather D, Abramov R, Wu IH, Chen K, Yamamoto-Hiraoka J, Goldenbogen J, Sotiropoulos KB, Clermont A, Geraldes P, Dall'Osso C, Wagers AJ, Huang PL, Rekhter M, Scalia R, Kahn CR, King GL: Loss of insulin signaling in vascular endothelial cells accelerates atherosclerosis in apolipoprotein E null mice. Cell Metab. 2010, 11: 379-389. 10.1016/j.cmet.2010.03.013.PubMed CentralView ArticlePubMed
  40. Reusch JE, Draznin BB: Atherosclerosis in diabetes and insulin resistance. Diabetes Obes Metab. 2007, 9: 455-463. 10.1111/j.1463-1326.2006.00620.x.View ArticlePubMed
  41. Baelde HJ, Eikmans M, Doran PP, Lappin DW, de HE, Bruijn JA: Gene expression profiling in glomeruli from human kidneys with diabetic nephropathy. Am J Kidney Dis. 2004, 43: 636-650. 10.1053/j.ajkd.2003.12.028.View ArticlePubMed
  42. Wogensen L, Nielsen CB, Hjorth P, Rasmussen LM, Nielsen AH, Gross K, Sarvetnick N, Ledet T: Under control of the Ren-1c promoter, locally produced transforming growth factor-beta1 induces accumulation of glomerular extracellular matrix in transgenic mice. Diabetes. 1999, 48: 182-192. 10.2337/diabetes.48.1.182.View ArticlePubMed
  43. Ziyadeh FN: Mediators of diabetic renal disease: the case for tgf-Beta as the major mediator. J Am Soc Nephrol. 2004, 15 (Suppl 1): S55-S57.View ArticlePubMed
  44. Woroniecka KI, Park AS, Mohtat D, Thomas DB, Pullman JM, Susztak K: Transcriptome analysis of human diabetic kidney disease. Diabetes. 2011, 60: 2354-2369. 10.2337/db10-1181.PubMed CentralView ArticlePubMed
  45. Chung AW, Yang HH, Sigrist MK, Brin G, Chum E, Gourlay WA, Levin A: Matrix metalloproteinase-2 and -9 exacerbate arterial stiffening and angiogenesis in diabetes and chronic kidney disease. Cardiovasc Res. 2009, 84: 494-504. 10.1093/cvr/cvp242.View ArticlePubMed

Copyright

© Skov et al; licensee BioMed Central Ltd. 2012

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.