Skip to content

Advertisement

You're viewing the new version of our site. Please leave us feedback.

Learn more

BMC Genomics

Open Access

Whole genome expression analysis within the angiotensin II-apolipoprotein E deficient mouse model of abdominal aortic aneurysm

  • Catherine Rush1,
  • Moses Nyara1,
  • Joseph V Moxon1,
  • Alexandra Trollope1,
  • Bradford Cullen1 and
  • Jonathan Golledge1Email author
BMC Genomics200910:298

https://doi.org/10.1186/1471-2164-10-298

Received: 04 September 2008

Accepted: 06 July 2009

Published: 06 July 2009

Abstract

Background

An animal model commonly used to investigate pathways and potential therapeutic interventions relevant to abdominal aortic aneurysm (AAA) involves subcutaneous infusion of angiotensin II within the apolipoprotein E deficient mouse. The aim of this study was to investigate genes differentially expressed in aneurysms forming within this mouse model in order to assess the relevance of this model to human AAA.

Results

Using microarrays we identified genes relevant to aneurysm formation within apolipoprotein E deficient mice. Firstly we investigated genes differentially expressed in the aneurysm prone segment of the suprarenal aorta in these mice. Secondly we investigated genes that were differentially expressed in the aortas of mice developing aneurysms relative to those that did not develop aneurysms in response to angiotensin II infusion. Our findings suggest that a host of inflammation and extracellular matrix remodelling pathways are upregulated within the aorta in mice developing aneurysms. Kyoto Encyclopedia of Genes and Genome categories enriched in the aortas of mice with aneurysms included cytokine-cytokine receptor interaction, leukocyte transendothelial migration, natural killer cell mediated cytotoxicity and hematopoietic cell lineage. Genes associated with extracellular matrix remodelling, such as a range of matrix metalloproteinases were also differentially expressed in relation to aneurysm formation.

Conclusion

This study is the first report describing whole genome expression arrays in the apolipoprotein E deficient mice in relation to aneurysm formation. The findings suggest that the pathways believed to be critical in human AAA are also relevant to aneurysm formation in this mouse model. The findings therefore support the value of this model to investigate interventions and mechanisms of human AAA.

Background

Abdominal aortic aneurysm (AAA) affects approximately 5% of men and 1% of women aged >60 years [1]. The principal concern with the condition is aortic rupture which is frequently fatal. Currently the only treatment option for AAA is surgery and approximately 25,000 aortic repairs are performed annually in the USA [2]. AAA accounts for approximately 15,000 deaths annually in the USA despite the increasing numbers of elective aortic repairs [2, 3]. There is increasing interest in using animal models of AAA to investigate mechanisms important in aneurysm development and progression in order to develop new non-interventional treatments and better ways of monitoring disease progression [4]. Currently rodent models of AAA are most commonly employed, particularly through infusion of angiotensin II in hyperlipidaemic mice, such as apolipoprotein E deficient (ApoE-/-) animals [57]. This mouse model has a number of similarities to human AAA such as the preponderance of aneurysm formation in males and the focal nature of aortic dilatation [5, 6]. The model also displays a number of differences from human AAA, for example aneurysms commonly involve the suprarenal aorta, dissection is an important histological finding and the infrarenal aorta is rarely if ever affected [7]. In contrast the infrarenal aorta is the most common site of AAA in humans. These and other disparities between this mouse model and human AAA stimulated us to investigate the gene expression profile of aneurysms in the ApoE-/- mouse model. The aim of the current study was to identify genes and pathways associated with aneurysm formation in the angiotensin II infused ApoE-/- mouse model. We used three approaches. Firstly we aimed to identify genes which might underlie the protection of the infrarenal aorta from aneurysm formation in ApoE-/- mice. Secondly we aimed to identify genes which were differentially expressed within the aortas of mice developing aneurysms in comparison to those that did not. Finally we selectively examined whether the differences in gene expression associated with aneurysm formation translated into similar changes in protein expression. The design of these studies is illustrated in Table 1. Our findings suggest the importance of chemokines, pro-inflammatory cytokines, leukocyte transendothelial migration mechanisms, a number of different signalling pathways (such as the Janus kinase pathway) and proteolytic mechanisms in aneurysm formation in these mice supporting the relevance of this model to human AAA [810].
Table 1

Summary of the included studies.

Study

Comparison

Design

Method

Mice

Genes identified

1

Paired supra and infrarenal aortas

Pooled samples

Codelink mouse whole genome bioarray (36227 genes)

10

Differentially expressed between AAA prone supra and infrarenal aorta

2

Whole aortas of mice exposed to angiotensin II that did and did not develop aneurysms and saline infused controls

Individual mice

Illumina mouse sentrix 6 microarray (46,628 genes).

18

Differentially expressed between the 3 groups revealing aneurysm-associated and "protective" genes

3

Suprarenal aortas of mice exposed to angiotensin II that did and did not develop aneurysms

Individual mice

ELISA and IHC

18*

Validate findings from the microarray in study 2 for Tnfrsf11b, Tgfb1, and inflammatory cell markers.

* 10 further mice died and were excluded from this study.

Results

Comparison of gene expression in aneurysm prone suprarenal compared to aneurysm resistant infrarenal aorta

A total of 26,522 transcripts (73% of the reference list) were expressed above background and compared between segments. A total of 304 transcripts were differentially expressed (1.5 fold, p < 0.05) between supra and infrarenal aortic segments (53 upregulated within the suprarenal aorta and 251 upregulated within the infrarenal aorta). A full list of transcripts is given in Additional File 1 and examples of potentially relevant genes shown in Table 2. A number of these genes may be relevant to the predilection of the suprarenal aorta to aneurysm formation based on knowledge of their functional importance or their previous identification in human AAA association studies [9]. Ucn3 for example, which was downregulated in the suprarenal aorta, has been demonstrated to antagonize angiotensin induced production of oxygen derived free radicals in vitro and play a vascular protective role in animal models [11, 12]. Pik3ap1, Dapk1 and Cd59 have previously been identified to be upregulated within human AAA biopsies [9]. Since vascular smooth muscle apoptosis and inflammation have been demonstrated to be important mechanisms in human AAA, upregulation of Dapk1 and Cd59 may be relevant to the predilection of the suprarenal aorta to aneurysm formation in ApoE-/- mice [8]. Perhaps more important to the predilection of the suprarenal aorta to aneurysm formation are the downregulation of a number of genes coding for structural extracellular matrix proteins including Lama3, Myo18b and the cytoskeletal protein Tuba8 (Table 2). Lama3 has been identified to be downregulated within human AAA biopsies [9].
Table 2

Examples of genes differentially expressed between paired supra and infrarenal aortas.

Symbol

Gene ID

Gene Name

Mean expression ratio*

P value

Adprhl1

234072

ADP-ribosylhydrolase like 1

3.11

0.015

Lgtn

16865

Ligatin

2.91

0.041

Fabp3

14077

Fatty acid binding protein 3

1.86

0.043

Pik3ap1†

83490

phosphoinositide-3-kinase adaptor protein 1

1.85

0.030

Dapk1†

69635

Death associated protein kinase 1 (Dapk1)

1.73

0.028

Tuba8

53857

tubulin, alpha 8

1.68

0.045

Cd59a†

12509

CD59a antigen

1.68

0.007

Cxcl13

55985

chemokine (C-X-C motif) ligand 13

1.56

0.021

Tgfbr3†

21814

Transforming growth factor, beta receptor III

1.53

0.030

Myo18b

74376

myosin XVIIIb

0.59

0.006

Ucn3

83428

urocortin 3

0.57

0.022

Srd5a2

94224

steroid 5 alpha-reductase 2

0.56

0.029

Dsc3

13507

desmocollin 3

0.52

0.037

Lama3†

16774

laminin, alpha 3

0.52

0.025

Trat1†

77647

T cell receptor associated transmembrane adaptor 1

0.51

0.008

Shbg

20415

sex hormone binding globulin

0.42

0.034

Akr1b7

11997

aldo-keto reductase family 1, member B7

0.24

< 0.001

Shown are a subset of the 304 gene transcripts that were differentially expressed in suprarenal compared to infrarenal aortic segments (1.5-fold, uncorrected p value < 0.05 calculated with paired t-test; for a full list see Additional File 1). The genes were selected from the full list based on the potential relevance to AAA. Genes highlighted with † have previously been identified as associated with human AAA in a whole genome expression study [9]. *Ratio of relative gene expression in supra compared to infrarenal aortic segments from 10 paired samples (5 pairs of arrays due to pooling of 2 aortic segments for each analysis).

Genes expressed in aortas from mice with aneurysms compared to those without aneurysms and controls

Maximum aortic diameters measured in the aortic arch, thoracic, suprarenal and infrarenal aortas in saline controls and mice exposed to angiotensin II that did and did not develop aneurysms are shown in Table 3 for animals used in the microarray experiment (study 2). Aortic diameters were similar between non-aneurysmal angiotensin II perfused mice and saline controls. We used Volcano Plots to identify differentially expressed genes between groups using our pre-defined criteria, i.e. 2-fold differential expression and uncorrected p value < 0.05 by unpaired t-test. The numbers of transcripts differentially expressed in comparison between the 3 groups of mice are shown in Additional File 2. 531 transcripts were differentially expressed (two fold change, p < 0.05) between aortas of mice exposed to angiotensin II which developed aneurysms and those that did not develop aneurysms (aneurysm resistant) (446 up and 85 downregulated, a full list is given in Additional File 3). 1196 transcripts were differentially expressed (two fold change, p < 0.05) between aortas of mice exposed to angiotensin II which developed aneurysms and saline controls (749 up and 447 downregulated, Additional File 4). 654 transcripts were differentially expressed (two fold change, p < 0.05) between aortas of mice exposed to angiotensin II which did not develop aneurysms and saline controls (128 up and 526 downregulated, Additional File 5). To further analyse these results, we concentrated on the transcripts differentially expressed between the aortas of mice that had been exposed to angiotensin II and had or had not developed AAAs (n = 531 listed in Additional File 3). We considered that genes identified by this analysis would be most likely to represent those relevant to aneurysm formation or resistance rather than simply reflecting the effects of angiotensin II. We also wanted to incorporate some of the findings relevant to the saline control group as they enabled us to have a measure of the baseline level of expression of different genes. This group allowed us to better postulate genes upregulated in the AAA resistant aortas, and thus of potential protective value against AAA. We carried out hierachical clustering analysis incorporating expression levels from all three groups of mice for the 531 transcripts that were differentially expressed in angiotensin II treated mice with and without aneurysms. This analysis is illustrated in terms of a heatmap, with a gene tree showing clustering of individual transcripts, in Figure 1. Each line in the heatmap represents a transcript with mean upregulated (red) or downregulated (green) in expression within the three groups. No change in expression is indicated in black. Hierachical clustering analysis indicated four main patterns of transcript expression (Figure 1, patterns 1–4). Transcripts grouped in pattern 1 (104 transcripts) showed decreased expression in aneurysm resistant aortas relative to other groups and included genes such as Ccl5, and Mmp9 which have been noted to be upregulated in human AAA biopsies (Additional File 6) [9]. Pattern 2 transcripts (325 transcripts) showed increased expression in aneurysms and included many genes associated with inflammation including matrix metalloproteinase encoding genes (e.g. Mmp2), genes encoding markers of inflammatory cells (Cd14, Cd68) and many cytokine and chemokine genes (e.g. Il1b, Il6, Ccl4, Ccl7, Ccl8, Ccl19, MMP2) (Table 4 shows a selected list and Additional File 7 a full list). Many of these genes, particularly those related to leukocyte chemotaxis, have previously been identified as upregulated within human AAA biopsies (highlighted in Table 4) [9, 10]. Pattern 3 transcripts (17 transcripts) were downregulated in the aortas of mice resistant to aneurysm formation. The difference between pattern 1 and pattern 3 transcripts is that in pattern 1, genes are most highly expressed in the AAA group, whereas pattern 3 genes are, in general, more highly expressed in the saline group. This can be clearly seen on the heatmap (Figure 1) where pattern 1 genes in the AAA group are mainly coloured red compared to red/black in saline group, whereas pattern 3 genes are red/black in the AAA group and red in the saline group. Although the clustering algorithim has separated these groups of genes, we believe that both pattern 1 and pattern 3 genes can be considered "pathogenic". Additional File 8 shows transcripts (out of those in pattern 3) that were associated with known genes and were significantly downregulated (2-fold, p < 0.05) with respect to both saline controls and mice that develop aneurysms. Two of these "pathogenic" pattern 3 genes, namely Apoc1 (Apolipoprotein C-I) and Stxbp2 (syntaxin binding protein 2), have been previously demonstrated to be upregulated within biopsies of human AAA [9]. Pattern 4 transcripts (85 transcripts) were upregulated in mice resistant to aneurysm formation (examples shown in Table 5 and a full list in Additional File 9) and we considered that upregulation of these may be "protective" for aneurysm formation. It is interesting to note that a number of the genes identified in pattern 4 have been previously shown to be downregulated within human AAA biopsies supporting their potential protective role (see Table 5 for examples) [9].
Table 3

Maximum aortic diameters of mice exposed to saline or angiotensin II that did or did not develop macroscopic aneurysms used in the microarray experiment (study 2).

Group

Diameter (mm)

 

Arch

Thoracic

Suprarenal

Infrarenal

Average Maximum#

Angiotensin II and aortic aneurysm (n = 5)

1.61 ± 0.18

1.11 ± 0.15

2.06 ± 0.48

0.72 ± 0.18

1.38 ± 0.17

Angiotensin II and no aortic aneurysm (n = 7)

1.34 ± 0.10

0.9 ± 0.10

1.13 ± 0.19

0.65 ± 0.09

1.0 ± 0.07

Saline controls (n = 6)

1.35 ± 0.11

1.16 ± 0.11

0.97 ± 0.08

0.67 ± 0.06

1.0 ± 0.06

Data represent mean values ± standard deviation; #calculated by averaging maximum diameters for all 4 aorta regions in each mouse.

Table 4

Examples of genes upregulated in the aortas of mice with aneurysms.

Gene symbol

Gene ID

Description; synonym

Fold increase cf no AAA

P value

Fold increase cf saline

P value

Chemokines/receptors

      

Cxcl10

15945

chemokine (C-X-C motif) ligand 10; IP-10

2.766

0.0284

2.487

0.0374

Cxcl12

20315

chemokine (C-X-C motif) ligand 12; SDF-1

3.022

0.0121

3.03

0.00358

Cxcl14

57266

chemokine (C-X-C motif) ligand 14; MIP-2g

3.812

0.0092

2.801

0.0373

Ccl2

20296

chemokine (C-C motif) ligand 2; MCP-1

2.263

0.021

3.325

0.00493

Ccl4

20303

chemokine (C-C motif) ligand 4; MIP-1B

4.232

0.0154

7.864

0.00302

Ccl7

20306

chemokine (C-C motif) ligand 7; MCP-3

3.183

0.00247

5.009

0.000346

Ccl8

20307

chemokine (C-C motif) ligand 8; MCP-2

3.12

0.00645

3.911

0.000657

Ccl19

24047

chemokine (C-C motif) ligand 19; ELC

2.119

0.00758

2.629

0.00228

Ccl21a

65956

chemokine (C-C motif) ligand 21; SLC

3.577

0.000615

2.66

0.0052

Cx3cr1

13051

Cx3cr1

2.602

0.00131

3.99

0.000356

Ccr5

12774

CD195

2.224

0.00187

3.358

0.000184

MMPs

      

Mmp2

17390

matrix metalloproteinase 2

2.513

0.000969

3.123

0.000609

Mmp12

17381

matrix metalloproteinase 12

2.172

0.0232

7.993

0.00031

Mmp13

17386

matrix metalloproteinase 13

6.428

0.02

13.62

0.00657

Mmp14

17387

matrix metalloproteinase 14

2.345

0.00432

3.99

0.000764

Cell lineage

      

Cd14

12475

CD14 antigen

2.068

0.00045

2.549

4.29E-05

Cd68

12514

CD68 antigen

2.907

0.00865

5.776

0.00162

Csf3r

12986

colony stimulating factor 3 receptor (granulocyte); CD114.

2.627

9.99E-05

2.53

0.00056

FceR1g

14127

Fc receptor, IgE, high affinity I, gamma polypeptide; CD23

2.835

0.0275

3.691

0.0125

Anpep

16790

alanyl (membrane) aminopeptidase; CD13

2.216

0.000594

3.312

5.46E-05

Cytokines and receptors

      

Il6

16193

interleukin 6

2.61

0.0411

3.92

0.0136

Il1b

16176

interleukin 1 beta

3.966

0.0158

5.89

0.00501

Antigen processing and presentation

      

Cd74

16149

Ia-associated invariant chain

2.376

0.00037

2.488

7.03E-05

H2-Ab1

14961

histocompatibility 2, class II antigen A, beta 1

2.799

0.000852

3.108

5.10E-06

Shown are genes which were significantly up-regulated (2-fold, uncorrected p value < 0.05 calculated with unpaired t-test) within the aortas of mice with AAAs by comparison to both mice that received angiotensin II and did not develop AAAs and also saline control mice. Genes were selected from the full list of differentially expressed genes (shown in Additional Files 3 and 4) based on their upregulation in both comparisons, likely functional importance and previous association with human AAA. Transcripts associated with no recognised gene are not included. Genes highlighted with † have previously been identified as associated with human AAA in a whole genome expression study or other human AAA studies [9, 10]. Genes highlighted in bold were selected for validation by real time PCR (see Table 6). Cf: compared to; ns: not significant

Table 5

Examples of genes upregulated in the aortas of mice exposed to angiotensin II which did not develop aneurysms. These genes are potentially protective against AAA.

Gene symbol

Gene ID

Description

Fold increase cf AAA

P value

Fold increase cf saline

P value

Cald1

109624

caldesmon 1

2.053388

0.000625

2.379

5.64E-06

Acta2

11475

actin, alpha 2, smooth muscle, aorta

2.178649

0.0284

2.536

1.84E-06

Sost

74499

sclerostin

2.659574

0.0118

2.167

9.01E-06

Kcnmb1

16533

potassium large conductance calcium-activated channel, subfamily M, beta member 1.

2.004008

0.0129

2.262

7.45E-07

Dstn

56431

destrin

2.169197

0.00498

2.242

6.70E-06

Timp4

110595

tissue inhibitor of metalloproteinase 4

2

1.99E-05

ns

ns

Pdlim3

53318

PDZ and LIM domain 3

2.12766

0.00386

ns

ns

Rgs17

56533

regulator of G-protein signaling 17

2.439024

0.00106

ns

ns

Hspa1l

15482

heat shock protein 1-like

2.020202

0.00844

ns

ns

Gabra3

14396

gamma-aminobutyric acid (GABA-A) receptor, subunit alpha 3

2.028398

0.00482

ns

ns

Btc

12223

betacellulin, epidermal growth factor family member

2.314815

0.0412

ns

ns

Fbxo30

71865

F-box protein 30

2.325581

0.000911

ns

ns

Hspa1a

193740

heat shock protein 1A

2.262443

0.00944

ns

ns

Mtap1b

17755

microtubule-associated protein 1B

2.252252

0.0363

ns

ns

Xirp1

22437

cardiomyopathy associated 1

2.155172

0.0165

ns

ns

Klk10

69540

Kallikrein 10

2.136752

7.42E-05

ns

ns

Slc22a1

20517

solute carrier family 22 (organic cation transporter), member 1

2.873563

0.000794

ns

ns

Shown are examples of genes which were significantly up-regulated (2-fold, uncorrected p value < 0.05 calculated with unpaired t-test) within the aortas of mice that received angiotensin II and did not develop AAAs by comparison to mice which developed AAAs. Also shown is comparison with findings from the aortas of saline control mice. Genes were selected from the full list of differentially expressed genes (shown in Additional File 9) based on their likely functional importance and previous association with human AAA. Genes highlighted with † have previously been identified as down-regulated within human AAA in a whole genome expression study [9]. Cf: compared to; ns: not significant.

Figure 1

Heatmap of 531 genes differentially expressed in the aortas of angiotensin II treated mice that developed aneurysms (n = 5) versus those that did not develop aneurysms (n = 7). A complete list of the genes used for this analysis is provided in Additional File 3. Expression in the saline controls (n = 6) for each of these genes is also shown. Each line in the heatmap represents a gene with mean upregulated (red) or downregulated (green) expression for the three groups. Black colour denotes no change in expression. Different patterns of gene expression are indicated as Patterns 1–4 and associated gene lists are provided in Additional Files 6, 7, 8 and 9 and Table 4. A Volcano plot (two fold change, p < 0.05) comparing the aneurysm group with aneurysm resistant group was used to select genes for this analysis.

Based on previous studies demonstrating the importance of inflammation in AAA we selected a number of inflammation associated markers to validate which were differentially expressed between aortic aneurysms and control aortas [9]. Initially we confirmed the up-regulation of Il6, Ccl4, Ccl8 and MMP2 within aortic aneurysms using real time PCR (Table 6). These genes were highly upregulated in aortic aneurysms by comparison to aortas from both control groups in microarray studies and real time PCR (see Tables 4 and 6). We also confirmed the differential expression of the Tnfrsf11b (osteoprotegerin) and Tgfb1 (TGF-b1) genes across the three groups translated to differences in protein expression using samples from study 3. We measured the concentration of OPG and TGFb-1 in the suprarenal segments of angiotensin II treated mice that developed AAA (mean maximum diameter 2.67 ± 0.24, n = 9) and those that did not (mean maximum diameter 1.26 ± 0.25, n = 9, p < 0.01) and found that both proteins were expressed at higher concentrations in aortas of mice with aneurysms, p < 0.05 (Figure 2A). Using immunohistochemistry we demonstrated that OPG was found within the medial layers and remodeled adventitia of aortas with aneurysm, whereas little was detected in aneurysm-resistant aortas or saline controls (Figure 2B).
Table 6

Validation of genes up-regulated in AAA by real time PCR.

Gene symbol

Gene ID

Description; synonym

Relative expression in AAA

Relative expression in no AAA

Relative expression in saline control

P value AAA v no AAA

P value AAA v saline control

Ccl4

20303

chemokine (C-C motif) ligand 4; MIP-1B

52.71 (19.80–133.40)

5.55 (4.87–9.91)

1.40 (0.84–2.49)

0.009

0.004

Ccl8

20307

chemokine (C-C motif) ligand 8; MCP-2

8.84 (5.94–12.25)

1.66 (1.27–2.73)

0.84 (0.55–1.18)

0.009

0.004

Mmp2

17390

matrix metalloproteinase 2

4.58 (3.80–11.00)

2.30 (1.50–3.17)

0.97 (0.81–1.06)

0.004

0.004

Il6

16193

interleukin 6

70.87 (14.40–184.59)

9.49 (2.54–16.83)

0.46 (0.17–1.26)

0.030

0.004

Shown are median (inter-quartile range) of relative gene expression (compared to house-keeping gene) in whole aortas used in the micro-array studies removed from mice exposed to angiotensin II that developed AAAs (n = 5) or no AAA (n = 7) or saline controls (n = 6).

Figure 2

Verification of microarray data for two genes in which aortic expression was related to aneurysm formation. A) ELISA for aortic osteoprotegerin and transforming growth factor beta-1 (right panel) confirmed microarray expression data (left panel). Gene expression data are presented as boxplots of normalised intensity for Tnfsf11b (OPG) and Tgfb1 (TGFb-1) for each group (aneurysmal mice n = 5, non-aneurysmal mice n = 7, saline controls n = 6). ELISA data from suprarenal AAAs (n = 9) and suprarenal aortas without AAA (n = 9) of ApoE-/- mice infused with angiotensin II. B) Immunohistochemistry confirms OPG localisation to aortic media (inset) and thrombus in suprarenal aortas with aneurysm. Aneurysm resistant and saline control aortas showed minimal OPG staining. Scale bars represent 500 μm, except inset where bar represents 50 μm.

Investigation of known functional KEGG pathways related to the differentially expressed genes

KEGG pathway analysis of transcripts upregulated in the AAA group compared to the no AAA group using Webgestalt, revealed that of the 302 known genes in the input list, 116 of these could be assigned to one or more KEGG pathways (Additional File 10). This pathway analysis revealed a significant enrichment for genes involved in leukocyte recruitment, transendothelial migration and B lymphocyte receptor signalling and identified a number of cell lineage specific genes associated with aneurysm (Table 7). Chemokine genes upregulated in aneurysms included Ccl21 and Ccl19 (which attract T cells and dendritic cells), Ccl7 (attracts monocytes), Cxcl9 (attracts T cells), Cxcl12 (attracts lymphocytes) and Cxcl5 (attracts neutrophils), suggesting that all these inflammatory cell types are recruited to mouse aortic aneurysms. Genes that can be considered as markers for the following cell types were upregulated in aneurysms (Additional Files 5 and 6): neutrophils (Csf3r); T lymphocytes (Cd3d) and B lymphocytes (Blnk, Cd72, various immunoglobulin genes); and monocyte/macrophages (Cd14, Cd68). We used IHC to confirm the recruitment of neutrophils, T and B lymphocytes, dendritic cells and monocyte/macrophages to aneurysmal but not aneurysm-resistant aortas (Figure 3).
Table 7

Top 10 KEGG pathways enriched in mouse aortic aneurysms.

Pathway (KEGG linked)

Total genes in pathway

Number genes upregulated in pathway in AAA v no AAA

Representative Genes (Entrez Gene ID)

P value*

Cytokine-cytokine receptor interaction

249

28

Cxcl5 (20311), Cxcl9 (17329), Cxcl10 (15945), Cxcl12 (20315), Cxcl13 (55985), Cxcl16 (66102), Cxcl14 (57266), Xcl1(16963), Cx3cr1 (13051), Cxcr4 (12767), Ccl19(24047), Ccl21a (20298), Ccl21b (18829), Ccl2 (20296), Ccl12 (20293), Ccl4 (20303), Ccl7 (20306), Ccl5 (20304), Ccl8 (20307), Ccr5 (12774), Il6 (16193), Csf3r (12986), Csf2rb2 (12984), Il7r (16197), Il2rg (16186), Il10ra (16154), Tnfrsf13b (57916), Il1b (16176)

8.70e-26

Leukocyte transendothelial migration

115

15

Cldn11 (18417), Mmp2 (17390), Mmp9 (17395), Ptpn11, Cybb (13058), Ncf2 (17970), Ncf4 (17972), Thy1 (21838), Cxcl12 (20315), Cxcr4 (12767), Rac2 (19354), Vav1 (22324), Ptk2b(19229), Cxcl13 (55985), Cxcl14 (57266)

4.20e-15

Focal adhesion

194

11

Col3a1 (12825), Col1a1 (12842), Comp (12845), Reln (19699), Thbs1 (21825), Thbs2 (21826), Igf1 (16000), Igf1r (16000), Vav1 (22324), Rac2 (19354), Itga11(319480), Itga6 (16403)

1.22e-7

B cell receptor signaling pathway

69

10

Cd72(12517), Cd79b (15985), Fcgr2b (14130), Inpp5d(16331), Btk (12229), Blnk (17060), Vav1(22324), Rac2(19354), Rasgrp3(240168), Ifitm1(68713)

3.75e-11

Natural killer cell mediated cytotoxicity

129

10

Hcst (23900), Klrd1 (16643), Fcgr3(14131), Fcer1g (14127), Tyrobp (22177, Ptk2b (19229), Vav1 (22324), Lcp2(16822), Cd48 (12506), Rac2 (19354)

8.70e-9

Hematopoietic cell lineage

86

10

H2-Eb1(14969), Il7r (16197), Cd3d (12500), Fcgr1(14129), Anpep (16790), Il6 (16193), Il1b(16176), Csf3r (12986), Itga6(16403), Cd14(12475)

3.24e-10

Jak-STAT signaling pathway

155

8

Il6(16193), Csfrb2(12984), Csf3r(12986), Il10ra(16154), Il2rg(16186), Il7r(16197), Jak3 (16453), Socs3(12702)

1.18e-6

ECM-receptor interaction

81

8

Reln(19699), Itga11(319480), Itga6(16403), Col3a1(12825), Col1a1(12842), Thbs1(21825), Thbs2(21826), Sdc3(20970)

1.43e-7

Toll-like receptor signaling pathway

101

7

Cd14(12475), Il6 (16193), Il1b(16176), Ccl5(20304), Ccl4(20303), Cxcl10(15945), Cxcl9(17329),

3.52e-6

Cell adhesion molecules (CAMs)

158

7

H2-Ab1(14961), H2-Eb1(14969), H2-DMb1(14999), Selpl (20345), Sdc3(20970), Itga6(16403), Cldn11(18417)

4.42e-5

Pathways were identified using WebGestalt by comparison of differentially expressed genes in the aortas of mice that did and did not develop aneurysms following angiotensin II infusion. Transcripts with >2 fold increase in expression in AAA compared to no AAA and p < 0.05 were submitted for pathway analysis. Pathways are linked to the KEGG website http://www.genome.jp and genes are linked to gene information. P value indicates the uncorrected significance of enrichment calculated from the hypergeometric test.

Figure 3

Immunohistochemistry demonstrating inflammatory cells in the suprarenal aorta of mice with aneurysms demonstrating functional confirmation of microarray expression data. Boxplots show normalised expression of genes that are markers of various haematopoietic cell lineages; Csf3r (neutrophils), Cd3d (T lymphocytes), Blnk (B lymphocytes), Itgax (dendritic cells) and Cd68 (monocyte/macrophages) from aneurysmal (n = 5), non-aneurysmal (n = 7) and saline control (n = 6) mice. Shown are sections stained for (A) neutrophils, (B) T lymphocytes (C) B lymphocytes (D) dendritic cells (E), macrophages in mice with aneurysms (i) and those without aortic dilatation (ii). Scale bars represent 100 μm.

Discussion

Mouse models of human diseases are potentially important tools with which to investigate mechanisms involved in the pathology and identify new treatments. Currently a number of animal models of human AAA are available [4]. Many of these models require significant interventions such as exposing the aorta and subsequently infusing or painting elastolytic solution on a segment of the artery to induce weakening and inflammation. Currently the angiotensin II infusion model is most commonly used, possibly due to the relatively ease with which aneurysms are produced and the lack of requirement for what might be considered artificial manipulation of the aorta. To our knowledge this is the first study to examine whole genome expression in relation to aneurysm formation within this model. Our findings demonstrate that aneurysms forming in ApoE-/- mice have marked influx by a range of inflammatory cells and upregulation of cytokines which have all been previously demonstrated in human AAA [8, 9, 1317]. These findings confirm those from more selective previous studies within this model [57, 1822]. Importantly they suggest the value of this model for investigating the inflammatory and cytokine aspects of human AAA. We also demonstrated the upregulation of a range of chemokines, cytokines and proteolytic enzymes, such as Ccl4, Ccl8, Il6 and Mmp2 which have been previously implicated in human AAA [9, 10, 23, 24]. The upregulation of the latter genes in AAAs was validated using real time PCR. We identified genes up or downregulated in the aortas of mice resistant to aneurysm formation suggesting potentially protective and pathological roles of these genes respectively in AAA. Sclerostin expression is decreased within human AAA biopsies compared to controls, however the significance of this finding is unclear [9]. We found increased expression of the sclerostin gene (Sost) in the aortas of mice resistant to aneurysm formation compared to both mice with aneurysms and saline controls, suggesting that sclerostin may play a role in inhibiting aortic dilatation. The ability of sclerostin to antagonise transforming growth factor beta, which has been linked to aortic aneurysm development, could be of significance in an aortic protection role for this gene [25]. Acta2 (vascular smooth muscle actin) expression was increased in aortas resistant to aneurysm formation compared to other groups. Acta2 encodes the single most abundant protein in vascular smooth muscle cells, and mutations in the human gene are associated with ascending thoracic aortic aneurysms and dissections [26]. Increased expression of Acta2 and other genes encoding smooth muscle cell proteins (Cald1, Dstn), in aortas protected from aneurysm, highlight the importance of vascular smooth muscle cells in maintaining vascular integrity. We observed downregulation of the genes encoding kininogen (kng1), Apolipoprotein CI (Apoc1) and the neutrophil-associated leucine-rich alpha-2-glycoprotein 1 (Lrg1) amongst others, in aortas of mice resistant to aneurysm suggesting that they may be pathological genes. The role of these in aneurysm formation warrants further investigation.

We also investigated genes that might underlie the predilection for suprarenal aneurysm formation within this model, highlighting downregulation of a number of genes which may explain the susceptibility of the suprarenal aorta to aneurysm formation. Urocortin has been shown in vitro and within animal models to inhibit the effects of the renin-angiotensin system, thus downregulation of Ucn3 (urocortin 3) within the suprarenal aorta of ApoE-/- mice may be relevant to its predilection for aneurysm formation following angiotensin II infusion [11, 12]. We also demonstrated downregulation of Lama3 (laminin alpha 3) within the suprarenal aorta. This extracellular matrix protein has been demonstrated to be present in reduced concentrations within some human AAAs and therefore may also be relevant to the preponderance of the suprarenal aorta to aneurysm formation [27, 28]. In vitro laminin also plays important roles in modulating inflammation and MMP production suggesting its relevance to aneurysm formation [29, 30]. Further studies will be required to investigate these candidate genes for example employing deficient mice.

A small number of previous studies have used expression arrays to examine gene transcripts associated with AAA in rodent models or patient samples [9, 28, 3134] (Table 8). There are a number of issues to consider when interpreting the results of the present study and those previously conducted, such as sample sizes, controls used, model employed (in the case of animal studies) and adjustment for multiple testing. Only one of the previous studies was a whole genome expression assessment, similar to the current study, with other studies examining much smaller numbers of transcripts. The sample sizes employed in all studies means that the power to detect differentially expressed genes is limited; therefore all current studies would have been expected to generate a large number of false negatives. We have therefore concentrated on identifying genes which were identified in the current study and also in these previous investigations and highlighted these throughout our tables (a summary of our findings in comparison with previous studies is provided in table 8). A number of chemokines were identified in all studies underlining the findings of the current investigation. Other common mechanisms identified to be associated with AAA included pro-inflammatory cytokines and matrix proteolysis. The design of the current experiment has some advantages over previous studies in the use of controls employed. We utilised two control sample types within study 2, i.e. aortas from mice receiving saline infusion but importantly also from mice that received angiotensin II and did not develop aneurysms. This has enabled us to try and isolate effects associated with AAA rather than just those due to angiotensin II. Findings in previous rodent studies may have related in part simply due to the effects of elastase [32, 33]. Controls for human expression array studies are problematic for a number of reasons, including the large amount of patient to patient variations, and the effects of factors that are impossible to control for in small samples such as differences in medication taken and co-morbidities [9, 28, 34].
Table 8

Comparison of results from this study and those from other expression arrays for AAA.

Ref.

Transcripts examined

Species

Number of samples/controls

Tissue of interest

Control tissue

Similar genes upregulated

9

18057 oligonucleotides microarray

Human

6/7

AAA

Normal#

Chemokines (e.g. CCL2, CCL4, CCL8, CCR5, CXCL5, CXCR4, TNFRSF13B)

Pro-inflammatory cytokines (e.g. IL1B, IL2RG)

Matrix metalloproteinase (e.g. MMP9)

Cell lineage markers (e.g. CD53, CD68, CD72)

32‡

1181 cDNA clones

Mice

3/3**

Elastase induced AAA

Inactivated elastase infused aortas

Pro-inflammatory cytokines (e.g. IL1B, IL6, INF, IGF1)

Chemokines (e.g. CCL7)

JAK-STAT signalling (e.g. SOC3)

31‡

375 cDNA clones

Human

1/1*

AAA body

AAA neck

Chemokines (e.g. CXCR2)

28‡

265 cDNA clones

Human

7/5/5†

AAA

AOD and normal†

Matrix metalloproteinase (e.g. MMP9, MMP12)

Pro-inflammatory cytokines (e.g. IL1B, IL6)

Chemokines (e.g. CCL4, CCR5)

Fibrinolysis (e.g. uPA)

33‡

8799 cDNA clones microarray

Rat

1/1‡‡

Elastase induced AAA

Saline infused aortas

Oxidative stress (Heme oxygenase, lipoxygenase)

34‡

1176 cDNA clones

Human

4/4

AAA

normal††

Cathepsins (e.g. cathepsin H)

Matrix metalloproteinase (e.g. MMP9)

Chemokines (e.g. CXCR4, CCL5)

*Results provided limited; †5 controls from aortic occlusive disease (AOD) from patients undergoing bypass and 5 controls from organ donors which were not age-matched; ‡These studies did not adjust for multiple testing; **pooling of samples was performed with 3 pools of interest compared with 3 controls for 4 different time periods after elastase infusion; ††Controls were from organ donor who were not matched; ‡‡Appears to have been comparison of one pool of samples from elastase induced aaas and one pool from controls examined on two microarrays each; #controls from post-mortem samples. Ref. = reference.

This study is the first to carry out a whole genome expression analysis within the angiotensin II- Apo E-/- model. We used an explorative approach using relatively small numbers of arrays (28 in total were included in the array studies). We were not adequately powered to adjust for multiple testing given the number of transcripts examined. Unlike some other microarray studies we have used real time PCR, IHC and ELISA to confirm some of our findings rather than relying on RNA assessment alone [9]. Numerous studies have now demonstrated the validity of microarray platforms therefore we felt it would be useful to demonstrate the functional relevance of altered RNA expression [9, 35]. Our findings suggest the value of this mouse model as one to investigate the role of inflammation, cell recruitment and proteolysis in AAA.

Conclusion

In conclusion our study supports the value of the angiotensin II infused ApoE-/- mouse model for investigating mechanisms and interventions relevant to human AAA.

Methods

Study design

Infusion of angiotensin II induces aortic dilatation particularly affecting the suprarenal aorta in ApoE-/- mice [5, 6, 1821]. Based on studies carried out in our and other laboratories the response to angiotensin II is variable, with some mice developing large aneurysms but other animals appearing resistant to aneurysm formation with aortic diameters similar to that of saline controls (Additional file 12). The infrarenal aorta is protected from aneurysm formation in these mice. To assess the likely signalling pathways relevant to aneurysm development and progression within this mouse model we carried out three studies (Table 1): 1) We compared RNA expression within segments of supra and infrarenal aortas from 13 week old male ApoE-/- mice unexposed to angiotensin II (n = 10). 2) We compared RNA expression from whole aortas of 17 week old male ApoE-/- mice which had been exposed to angiotensin II (1.44 μg/kg/min) for 4 weeks where there was clear evidence of aortic aneurysm formation (n = 5) with that of mice failing to develop aneurysms (n = 7) and those exposed to saline infusion (n = 6). 3) We selected 2 genes and 5 cellular pathways identified as upregulated within the aortas of mice with aneurysms for validation using protein assessments including ELISAs and immunohistochemistry (IHC). The selection of genes was based on their association with pathways we had identified from pathway analysis of the data and recognised ability to be able to assess them using ELISAs or IHC. For study 3, twenty-eight additional 17 week old male ApoE-/- mice were infused with angiotensin II for 4 weeks. Ten mice in study 3 died prematurely due to aortic rupture and were excluded from further analysis in order to avoid post-mortem or post rupture effects confounding assessments, leaving 18 animals in this group.

Mice

This investigation conformed to the Guide for the Care and Use of Laboratory Animals published by the US National Institutes of Health (NIH Publication No. 85-23, revised 1996). Ethical approval was obtained from the local institutional committee prior to commencement of the study. Male ApoE-/- mice were obtained from Animal Resources Centre, Canning Vale, Western Australia aged 6–8 weeks. Mice were fed standard chow, and at 13 weeks either their aortas were harvested (study 1, n = 10) or angiotensin II (n = 12 study 2; n = 28 study 3) or saline (n = 6 study 2) infusion commenced. Aortic harvesting and osmotic minipump placement (Model 2004, ALZET, Durect Corporation, Cupertino, California, USA) was carried out under ketamine (150 mg/kg i.p.) and xylazine (10 mg/kg i.p.) anaesthesia as previously described [22]. The 28 day infusion period and dose of angiotensin II were selected as these have been routinely used in other studies within this model [5, 6, 1821]. In studies 1 and 2, mice were euthanized by carbon dioxide asphyxiation, the aortas perfused with RNAlater® (Qiagen, Doncaster, Victoria, Australia), and harvested from arch to iliac bifurcation and stored at -80°C for later analysis. For study 3 mice were euthanized, aortas were perfused with PBS and the suprarenal aortic segments divided into 2 and stored separately at -20°C for later IHC and cytokine assessment.

Aorta Morphometry

Aortas were placed on a black background and digitally photographed (Coolpix 4500, Nikon). Maximum diameters of the aortic arch, thoracic, suprarenal and infrarenal aorta were determined from the images using computer-aided analysis (Scion Image, Scion Corporation). Preliminary studies (n = 27) established that these measurements could be repeated with good intra-observer reproducibility (Coefficient of repeatability 0.98, 95% confidence intervals 0.975–0.982, and coefficient of variation 4%).

Microarrays

Total RNA was extracted from homogenised aortas using the TRI reagent (Sigma) and RNeasy Mini kits (Qiagen) and analysed on the HP 2100 Bioanalyzer (Agilent Technologies) for integrity using Eukaryotic total RNA nano chips (Agilent) (Additional File 11). All samples had RNA integrity scores between 6.9 and 8.1. For study 1, RNA extracted from the pooled suprarenal aortas of two mice was compared with that obtained from the infrarenal aortic segments of the same mice. A total of 10 mice (i.e. 5 pairs of pooled segments) were included in study 1. For study 2, RNA extracted from the whole aortas of mice with aneurysms (n = 5, Additional File 12A) was compared with that extracted from the whole of aortas of saline controls (n = 6, Additional File 12C) and mice exposed to angiotensin II which did not develop aneurysms (n = 7, Additional File 12B). No sample pooling was carried out in study 2. RNA hybridization was performed, and gene expression profiles determined, using CodeLink Mouse Whole Genome Bioarray chips (GE healthcare, Amersham, Bioscience) (study 1, n = 10) and Illumina Mouse Sentrix 6 version 1.1 Beadchips (study 2, n = 18). For study 1 total RNA from infrarenal and suprarenal aortas was supplied to GenUS BioSystems Inc. (Northbrook, Illinois USA) who conducted cDNA synthesis, hybridization and scanning using standard protocols. In brief, 10 μg of total RNA was used for each 'sample versus sample' comparison and first and second strand cDNA were prepared using the Codelink™ iExpress Expression assay Reagent kit (GE healthcare Amersham). Biotinylated cRNA target was prepared from the cDNA template by linear amplification using biotin-dNTPs and verified on a Bioanalyzer 2100. The cRNA was fragmented to uniform size and once again verified on the Bioanalyzer 2100. Each fragmented, biotin-labeled cRNA was added to a Codelink Mouse Whole Genome array which contained 34,957 probes targeting unique mouse transcripts. The bioarrays were washed, exposed to Cy5 streptavidin and scanned using GenePix 4000B laser scanner. Scanned image files were examined using CodeLink image and data analysis software (GE healthcare, Amersham, Bioscience). The image information was converted into spot intensity values using CodeLink™ Expression Analysis software (GE healthcare). The generated values were exported to GeneSpring GX 7.3.1 software (Silicon Genetic, USA) for further analysis. Study 2 labelling, hybridisation and scanning was done at the SRC Microarray Facility, University of Queensland, Brisbane, Australia according to the manufacturer's instructions. For study 2 total RNA was amplified in a single-round of in vitro transcription amplification that allowed incorporation of biotin-labeled nucleotides using the Illumina TotalPrep RNA amplification kit (Ambion, Inc., Austin, TX). mRNA samples were assessed for integrity and purity prior to hybridization using the Bioanalyser 2100 with mRNA Nano chips (Agilent Technologies). cRNA of each sample was hybridized to an Illumina Mouse WG-6 V1.1 BeadChip followed by washing, blocking, and streptavidin-Cy3 staining steps, and scanning with a high-resolution Illumina BeadArray reader scanner. The data extraction was performed by using Illumina Bead Studio V2.3.41 software with the output being raw, non-normalized bead summary values. The BeadStudio matrix contains the summarised expression values (Avg_Signal), standard error of the bead replicates (BEADSTDEV), number of beads used (Avg_NBEADS) and a detection score, which estimates the confidence limit of detection of a gene. The performance of the built-in controls that accompany each Illumina beadchip experiment was assessed as part of the Bead Studio V2 experiment performance report and was found to be satisfactory. Controls included Housekeeping controls for intactness of the biological specimens, negative controls to establish gene expression detection limits and hybridization controls including low and high stringency controls and biotin signal generation controls. The raw data matrix extracted from Beadstudio was uploaded into GeneSpring GX 7.3.1 (Silicon Genetics, Redwood City, CA) software for downstream analysis. Details of GeneSpring analysis including data transformations, per chip and per gene normalisations are described below in Analysis and Design.

Real time PCR

Using RNA obtained for the mice employed in the micro-array investigations in study 2 we validated findings for 4 genes (Ccl4, Ccl8, Il6 and Mmp2) using real time PCR. RNA samples from mice exposed to angiotensin II that developed aneurysms (n = 5), mice exposed to angiotensin II that did not develop aneurysms (n = 7) and saline control mice (n = 6) were included in the analysis. The QuantiTect SYBR Green one-step RT-PCR Kit (Qiagen) and Quantitect Primer Assays (Qiagen) (Il6, QT00098875; Ccl4, QT00154616; Ccl8, QT00128548; Mmp2, QT00116116) were used according to the manufacturer's instructions with 25 ng of total RNA as template. Primers for mouse Gapdh were used to amplify the housekeeping gene (Qiagen, QT01658692). Standard curves for each gene including the housekeeping gene were constructed using duplicate sets of five (5) 10-fold serial dilutions of equal volumes of the pooled saline control RNAs. Negative and "minus-RT" controls were also included. All reactions were independently repeated in duplicate to ensure the reproducibility of the results. Cycling parameters were as follows: 50°C 30 min for RT; 95°C 15 min; 40 cycles of 94°C 15 sec, 55°C 30 sec, 72°C 30 sec. Data were viewed and analysed using the Rotor-Gene's real-time analysis software (Rotor-Gene 6000; Corbett Life Science, Sydney). The relative expression of the gene of interest in each sample was calculated by the Rotor-Gene software using the concentration-Ct-standard curve method and normalised using the average expression of Gapdh for each sample. SSPS statistical software was used to calculate median and interquartile ranges for expression of each gene of interest in each group.

IHC

Serial cryostat sections 7 μm thick were cut from suprarenal aortas with and without aneurysms prior to staining for inflammatory cells. Serial frozen sections were air-dried, fixed in acetone for 10 min at -20°C, air dried and rehydrated with PBS before being incubated in 3% H2O2/0.1% sodium azide/PBS to block endogenous peroxidase. For macrophage detection, sections were blocked in 2% normal goat serum in PBS followed by staining using pan-macrophage antibody (clone MOMA-2, Abcam AB33451), and goat anti-rat HRP (Millipore AP136P). Rat IgG (Sigma I4131) was used as isotype control. Detection of osteoprotegerin (OPG) was carried out using biotinylated goat anti-mouse OPG (R&D BAF459) diluted to 1 μg/ml in blocking buffer followed by Vectastain Elite ABC-HRP (Vector Laboratories). Biotinylated goat-IgG was used as isotype control (Vector Laboratories VEBI1001). T cells (CD3e FITC clone 145-2C11, BD Biosciences 553061), B cells (CD45R/B220 biotin clone RA3-6B2, BD Biosciences 550385), dendritic cells (CD11c biotin clone HL3, BD Biosciences 553800) and neutrophils (Ly6G/Ly6C(Gr-1) FITC clone RB6-8C5, BD Biosciences 553127) were demonstrated using either anti-FITC HRP (Invitrogen), biotinyl-tyramide (Perkin Elmer), SA-HRP (Perkin Elmer) for FITC-conjugated antibodies or SA-HRP, biotinyl-tyramide, SA-HRP for biotinylated antibodies. Appropriate isotype controls were added to other sections (all BD Biosciences). Slides were incubated in the peroxidase substrate 3, 3'-diamminobenzidine (ImmPACT DAB, Vector), counterstained in Mayer's Haematoxylin, dehydrated, cleared in xylene and mounted in Depex mounting medium. Sections were photographed using a Nikon Eclipse 50i microscope, Digital Sight camera and NIS-elements software.

Cytokine studies

Protein was extracted from individual frozen suprarenal aortic segments by homogenising in buffer (10 mM cacodylic acid, 60 mM L-arginine, 0.25% triton x-100 in PBS, pH 7.2) and centrifuging at 18,000 × g at 4°C for 20 min. Supernatant protein was quantified by the Bradford technique (Protein Assay, Bio-Rad, Hercules, California, USA). Concentrations of OPG and transforming growth factor beta-1 (TGFb-1) were assessed using commercial ELISAs (Quantikine, R&D Systems, MOP00 for OPG, MB100B for TGFb-1) and expressed as pg/mg of protein. Initial tests were performed to determine an appropriate amount of extracted aortic protein to load in each well to suit the ELISA standard curves (10 μg for OPG ELISA, 20 μg for TGFb-1 ELISA). Serum and extracted aortic proteins from OPG deficient mice were confirmed as having zero readings using the OPG ELISA (data not shown). Batch analysis of the samples was carried out to facilitate comparison between groups. We have previously reported excellent reproducibility of similar assays [36].

Analysis and design

We designed the current study as an exploratory analysis expecting a minimal level of false discovery inevitable in analysis of large microarray gene probes. Sample size calculations were conducted using SPCalc [37, 38]. For study 1 (which included paired samples) we aimed to detect a 1.5-fold difference at a power of 0.8 (alpha uncorrected 0.05), while for study 2 (which included unpaired samples) we aimed to detect a 2-fold difference. Based on the number of gene probes to be examined (26,522) we wished to control the false positives to a mean of 2 [17, 18]. We estimated the total number of arrays needed for study 1 and 2 were 10 and 18 respectively. Microarray data were imported into Genespring GX 7.3.1 (Agilent) for analysis. The data from the individual bioarrays have been deposited in NCBIs Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO accession number GSE7006 (GEO Accession viewer) for study 1 and GSE12591 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=rduxbsiqmqswefi&acc=GSE12591 for study 2. For study 1 which used paired samples, normalisations were as follows: each chip was normalised to the median intensity of the array; each gene was normalised relative to the suprarenal sample of each pair. Thus the expression of each SRA sample was normalised to 1. Volcano plots were used to identify differentially expressed genes with >1.5 fold difference in expression and p < 0.05. Statistical analysis was carried out using a paired t-test (unequal variances) comparing the mean infrarenal/suprarenal ratio to a baseline value of 1 for each sample pair. In study 2 we followed the standard normalization procedures recommended for the GeneSpring GX 7.3.1 software for one-colour array data. In brief, data transformation was corrected for a low signal, with values recorded at <0.01 increased to the minimum (0.01). Default settings were used for experiment normalisation which included chip normalisation to 50th percentile and gene normalisation to the median. The cross gene error model was used to further filter low quality data. We sought to identify genes with a 2-fold differential expression within the aortas of mice with and without aneurysms based on an uncorrected p value of < 0.05 included in our sample size estimate. We initially compared findings from mice that received angiotensin II and developed AAAs (n = 5) with both other groups of mice (n = 7 and 6). Genes showing a greater than 2 fold difference in expression between groups (unpaired t-test (unequal variance) p < 0.05), were considered to be differentially expressed. We reasoned that the most significant findings would relate to genes differentially expressed between angiotensin II perfused mice, that did (n = 5) and did not (n = 7) develop aneurysms. Volcano plots were used to identify differentially expressed genes between the aortas of these two groups of mice. This group of 531 genes was further investigated by hierarchical clustering in GeneSpring, in which a tree of transcripts or genes is built by successively finding the two most similar gene expression patterns from the full transcript set. Individual samples were then arranged in a condition tree according to overall similarity. We also included the saline group in this clustering analysis in order to show how these 531 genes are expressed in control aortas. Genes with similar expression patterns were grouped as hierarchical clusters with distances between samples computed using Pearson correlations for similarity measures and average linkage as the clustering algorithm. Hierachical clustering revealed several major nodes in the gene tree structure which indicated at least 4 different patterns of gene expression. Lists of differentially expressed genes were examined for biologically relevant associations using Gene Ontologies and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analysis in the web-based software Webgestalt http://bioinfo.vanderbilt.edu/webgestalt. Quantitative real time PCR outcomes in study 2 and concentrations of cytokines measured in study 3 were compared between groups using Mann Whitney U test.

Declarations

Acknowledgements

Grants from the National Institute of Health, USA (RO1 HL080010) and NHMRC (project grant 540403) supported this work. JG is supported by a Practitioner Fellowships from the NHMRC, Australia (431503). These funding sources took no part in the design of this study, data collection or analysis, and the decision to submit the work for publication. Study 2 was supported by the Australian Research Council's, Special Research Centre for Functional and Applied Genomics (Institute for Molecular Bioscience) Microarray Facility.

Authors’ Affiliations

(1)
The Vascular Biology Unit, School of Medicine, James Cook University

References

  1. Fleming C, Whitlock EP, Beil TL, Lederle FA: Screening for abdominal aortic aneurysm: a best-evidence systematic review for the U.S. Preventive Services Task Force. Ann Intern Med. 2005, 142: 203-211.View ArticlePubMedGoogle Scholar
  2. Timaran CH, Veith FJ, Rosero EB, Modrall JG, Arko FR, Clagett GP, Valentine RJ: Endovascular aortic aneurysm repair in patients with the highest risk and in-hospital mortality in the United States. Arch Surg. 2007, 142: 520-524. 10.1001/archsurg.142.6.520.View ArticlePubMedGoogle Scholar
  3. Miniño A, Heron M, Smith B, Kochanek K: Deaths: Final data for 2004. National vital statistics reports. Hyattsville, MD. [http://www.cdc.gov/nchs/products/pubs/pubd/hestats/finaldeaths04/finaldeaths04.htm]
  4. Daugherty A, Cassis LA: Mouse models of abdominal aortic aneurysms. Arterioscler Thromb Vasc Biol. 2004, 24: 429-434. 10.1161/01.ATV.0000118013.72016.ea.View ArticlePubMedGoogle Scholar
  5. Manning MW, Cassis LA, Huang J, Szilvassy SJ, Daugherty A: Abdominal aortic aneurysms: fresh insights from a novel animal model of the disease. Vasc Med. 2002, 7: 45-54. 10.1191/1358863x02vm413ra.View ArticlePubMedGoogle Scholar
  6. Daugherty A, Manning MW, Cassis LA: Angiotensin II promotes atherosclerotic lesions and aneurysms in apolipoprotein E-deficient mice. J Clin Invest. 2000, 105: 1605-1612. 10.1172/JCI7818.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Saraff K, Babamusta F, Cassis LA, Daugherty A: Aortic dissection precedes formation of aneurysms and atherosclerosis in angiotensin II-infused, apolipoprotein E-deficient mice. Arterioscler Thromb Vasc Biol. 2003, 23: 1621-1626. 10.1161/01.ATV.0000085631.76095.64.View ArticlePubMedGoogle Scholar
  8. Schönbeck U, Sukhova GK, Gerdes N, Libby P: T(H)2 predominant immune responses prevail in human abdominal aortic aneurysm. Am J Pathol. 2002, 161: 499-506.PubMed CentralView ArticlePubMedGoogle Scholar
  9. Lenk GM, Tromp G, Weinsheimer S, Gatalica Z, Berguer R, Kuivaniemi H: Whole genome expression profiling reveals a significant role for immune function in human abdominal aortic aneurysms. BMC Genomics. 2007, 8: 237-10.1186/1471-2164-8-237.PubMed CentralView ArticlePubMedGoogle Scholar
  10. Knox JB, Sukhova GK, Whittemore AD, Libby P: Evidence for altered balance between matrix metalloproteinases and their inhibitors in human aortic diseases. Circulation. 1997, 95: 205-212.View ArticlePubMedGoogle Scholar
  11. Honjo T, Inoue N, Shiraki R, Kobayashi S, Otsui K, Takahashi M, Hirata K, Kawashima S, Yokozaki H, Yokoyama M: Endothelial urocortin has potent antioxidative properties and is upregulated by inflammatory cytokines and pitavastatin. J Vasc Res. 2006, 43: 131-138. 10.1159/000090132.View ArticlePubMedGoogle Scholar
  12. Rademaker MT, Charles CJ, Espiner EA, Frampton CM, Lainchbury JG, Richards AM: Endogenous urocortins reduce vascular tone and renin-aldosterone/endothelin activity in experimental heart failure. Eur Heart J. 2005, 26: 2046-2054. 10.1093/eurheartj/ehi227.View ArticlePubMedGoogle Scholar
  13. Bobryshev YV, Lord RS: Vascular-associated lymphoid tissue (VALT) involvement in aortic aneurysm. Atherosclerosis. 2001, 154: 15-21. 10.1016/S0021-9150(00)00441-X.View ArticlePubMedGoogle Scholar
  14. Galle C, Schandené L, Stordeur P, Peignois Y, Ferreira J, Wautrecht JC, Dereume JP, Goldman M: Predominance of type 1 CD4+ T cells in human abdominal aortic aneurysm. Clin Exp Immunol. 2005, 142: 519-527.PubMed CentralPubMedGoogle Scholar
  15. Forester ND, Cruickshank SM, Scott DJ, Carding SR: Functional characterization of T cells in abdominal aortic aneurysms. Immunology. 2005, 115: 262-270. 10.1111/j.1365-2567.2005.02157.x.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Folkesson M, Kazi M, Zhu C, Silveira A, Hemdahl AL, Hamsten A, Hedin U, Swedenborg J, Eriksson P: Presence of NGAL/MMP-9 complexes in human abdominal aortic aneurysms. Thromb Haemost. 2007, 98: 427-433.PubMedGoogle Scholar
  17. Middleton RK, Lloyd GM, Bown MJ, Cooper NJ, London NJ, Sayers RD: The pro-inflammatory and chemotactic cytokine microenvironment of the abdominal aortic aneurysm wall: a protein array study. J Vasc Surg. 2007, 45: 574-580. 10.1016/j.jvs.2006.11.020.View ArticlePubMedGoogle Scholar
  18. Manning MW, Cassis LA, Daugherty A: Differential effects of doxycycline, a broad-spectrum matrix metalloproteinase inhibitor, on angiotensin II-induced atherosclerosis and abdominal aortic aneurysms. Arterioscler Thromb Vasc Biol. 2003, 23: 483-488. 10.1161/01.ATV.0000058404.92759.32.View ArticlePubMedGoogle Scholar
  19. Wang YX, Martin-McNulty B, da Cunha V, Vincelette J, Lu X, Feng Q, Halks-Miller M, Mahmoudi M, Schroeder M, Subramanyam B, Tseng JL, Deng GD, Schirm S, Johns A, Kauser K, Dole WP, Light DR: Fasudil, a Rho-kinase inhibitor, attenuates angiotensin II-induced abdominal aortic aneurysm in apolipoprotein E-deficient mice by inhibiting apoptosis and proteolysis. Circulation. 2005, 111: 2219-2226. 10.1161/01.CIR.0000163544.17221.BE.View ArticlePubMedGoogle Scholar
  20. King VL, Trivedi D, Gitlin JM, Loftin CD: Selective cyclooxygenase-2 inhibition with celecoxib decreases angiotensin II-induced abdominal aortic aneurysm formation in mice. Arterioscler Thromb Vasc Biol. 2006, 26: 1137-1143. 10.1161/01.ATV.0000216119.79008.ac.View ArticlePubMedGoogle Scholar
  21. Martin-McNulty B, Tham DM, da Cunha V, Ho JJ, Wilson DW, Rutledge JC, Deng GG, Vergona R, Sullivan ME, Wang YX: 17 Beta-estradiol attenuates development of angiotensin II-induced aortic abdominal aneurysm in apolipoprotein E-deficient mice. Arterioscler Thromb Vasc Biol. 2003, 23: 1627-1632. 10.1161/01.ATV.0000085842.20866.6A.View ArticlePubMedGoogle Scholar
  22. Moran CS, Cullen B, Campbell JH, Golledge J: Interaction between angiotensin II, osteoprotegerin and peroxisome proliferator-activated receptorã in abdominal aortic aneurysm. J Vasc Res. 2009, 46: 209-217. 10.1159/000163019.View ArticlePubMedGoogle Scholar
  23. McMillan WD, Patterson BK, Keen RR, Pearce WH: In situ localization and quantification of seventy-two-kilodalton type IV collagenase in aneurysmal, occlusive, and normal aorta. J Vasc Surg. 1995, 22: 295-305. 10.1016/S0741-5214(95)70144-3.View ArticlePubMedGoogle Scholar
  24. McMillan WD, Patterson BK, Keen RR, Shively VP, Cipollone M, Pearce WH: In situ localization and quantification of mRNA for 92-kD type IV collagenase and its inhibitor in aneurysmal, occlusive, and normal aorta. terioscler Thromb Vasc Biol. 1995, 15: Ar1139-1144.View ArticleGoogle Scholar
  25. Habashi JP, Judge DP, Holm TM, Cohn RD, Loeys BL, Cooper TK, Myers L, Klein EC, Liu G, Calvi C, Podowski M, Neptune ER, Halushka MK, Bedja D, Gabrielson K, Rifkin DB, Carta L, Ramirez F, Huso DL, Dietz HC: Losartan, an AT1 antagonist, prevents aortic aneurysm in a mouse model of Marfan syndrome. Science. 2006, 312: 117-121. 10.1126/science.1124287.PubMed CentralView ArticlePubMedGoogle Scholar
  26. Guo DC, Pannu H, Tran-Fadulu V, Papke CL, Yu RK, Avidan N, Bourgeois S, Estrera AL, Safi HJ, Sparks E, Amor D, Ades L, McConnell V, Willoughby CE, Abuelo D, Willing M, Lewis RA, Kim DH, Scherer S, Tung PP, Ahn C, Buja LM, Raman CS, Shete SS, Milewicz DM: Mutations in smooth muscle alpha-actin (ACTA2) lead to thoracic aortic aneurysms and dissections. Nat Genet. 2007, 3912: 1488-1493. 10.1038/ng.2007.6.View ArticleGoogle Scholar
  27. Cotrufo M, De Santo L, Della Corte A, Di Meglio F, Guerra G, Quarto C, Vitale S, Castaldo C, Montagnani S: Basal lamina structural alterations in human asymmetric aneurismatic aorta. Eur J Histochem. 2005, 49: 363-370.View ArticlePubMedGoogle Scholar
  28. Armstrong PJ, Johanning JM, Calton WC, Delatore JR, Franklin DP, Han DC, Carey DJ, Elmore JR: Differential gene expression in human abdominal aorta: aneurysmal versus occlusive disease. J Vasc Surg. 2002, 35: 346-355. 10.1067/mva.2002.121071.View ArticlePubMedGoogle Scholar
  29. Adair-Kirk TL, Atkinson JJ, Kelley DG, Arch RH, Miner JH, Senior RM: A chemotactic peptide from laminin alpha 5 functions as a regulator of inflammatory immune responses via TNF alpha-mediated signaling. J Immunol. 2005, 174: 1621-1629.View ArticlePubMedGoogle Scholar
  30. Adair-Kirk TL, Atkinson JJ, Broekelmann TJ, Doi M, Tryggvason K, Miner JH, Mecham RP, Senior RM: A site on laminin alpha 5, AQARSAASKVKVSMKF, induces inflammatory cell production of matrix metalloproteinase-9 and chemotaxis. J Immunol. 2003, 171: 398-406.View ArticlePubMedGoogle Scholar
  31. Yamagishi M, Higashikata T, Ishibashi-Ueda H, Sasaki H, Ogino H, Iihara K, Miyamoto S, Nagaya N, Tomoike H, Sakamoto A: Sustained upregulation of inflammatory chemokine and its receptor in aneurysmal and occlusive atherosclerotic disease: results form tissue analysis with cDNA macroarray and real-time reverse transcriptional polymerase chain reaction methods. Circ J. 2005, 69: 1490-1495. 10.1253/circj.69.1490.View ArticlePubMedGoogle Scholar
  32. Van Vickle-Chavez SJ, Tung WS, Absi TS, Ennis TL, Mao D, Cobb JP, Thompson RW: Temporal changes in mouse aortic wall gene expression during the development of elastase-induced abdominal aortic aneurysms. J Vasc Surg. 2006, 43: 1010-1020. 10.1016/j.jvs.2006.01.004.View ArticlePubMedGoogle Scholar
  33. Yajima N, Masuda M, Miyazaki M, Nakajima N, Chien S, Shyy JY: Oxidative stress is involved in the development of experimental abdominal aortic aneurysm: a study of the transcription profile with complementary DNA microarray. J Vasc Surg. 2002, 36: 379-85. 10.1067/mva.2002.124366.View ArticlePubMedGoogle Scholar
  34. Tung WS, Lee JK, Thompson RW: Simultaneous analysis of 1176 gene products in normal human aorta and abdominal aortic aneurysms using a membrane-based complementary DNA expression array. J Vasc Surg. 2001, 34: 143-50. 10.1067/mva.2001.113310.View ArticlePubMedGoogle Scholar
  35. Barnes M, Freudenberg J, Thompson S, Aronow B, Pavlidis P: Experimental comparison and cross-validation of the Affymetrix and Illumina gene expression analysis platforms. Nucleic Acids Res. 2005, 33: 5914-5923. 10.1093/nar/gki890.PubMed CentralView ArticlePubMedGoogle Scholar
  36. Clancy P, Oliver L, Jayalath R, Buttner P, Golledge J: Assessment of a serum assay for quantification of abdominal aortic calcification. Arterioscler Thromb Vasc Biol. 2006, 26: 2574-2576. 10.1161/01.ATV.0000242799.81434.7d.View ArticlePubMedGoogle Scholar
  37. Qiu W, Lee ML: SPCalc: A web-based calculator for sample size and power calculations in micro-array studies. Bioinformation. 2006, 1: 251-2.PubMed CentralView ArticlePubMedGoogle Scholar
  38. Lee ML, Whitmore GA: Power and sample size for DNA microarray studies. Stat Med. 2002, 21: 3543-70. 10.1002/sim.1335.View ArticlePubMedGoogle Scholar

Copyright

© Rush et al; licensee BioMed Central Ltd. 2009

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.

Advertisement