<title>Background</title> <p>Peripheral arterial occlusive disease (PAD) is a major manifestation of atherosclerosis and is commonly found in elderly patients. Epidemiological studies have shown that PAD affects 8 to 10 million adults in the United States [##REF##17345122##1##]. Most patients with PAD are asymptomatic. The disease is primarily diagnosed by an ankle brachial index (ABI) &lt; 0.9. The most common symptom of mild-to-moderate PAD is intermittent claudication, which is present in about one third of symptomatic patients [##REF##17345122##1##]. In addition to leg symptoms, patients with PAD are at an increased risk for developing new coronary events and eventually death from cardiovascular disease. Although conventional procedures such as stents, arterectomies, angioplasty, and bypass surgery have been successful in improving clinical symptoms of PAD to a large extent [##REF##17046794##2##], ultimately elimination of the disease may require sophisticated protocols of pharmaceutical interventions, which may depend on better understanding of molecular mechanisms involved in the disease.</p> <p>Previous studies have implicated the involvement of the immune system in atherosclerosis formation and progression. Animal models have been used to test the contributions of components of the immune system [##REF##1411543##3##,##REF##8274468##4##]. Cellular involvement of macrophages was found to be important in the formation and progression of atherosclerosis in animal models [##REF##8274468##4##]. In addition, various immune-related genes have been examined in an atherosclerosis animal model, and genes such as <italic>CXCR6, CXCL10, CXCR3 </italic>and <italic>CXCL16/scavenger receptor </italic>have been shown to be involved in the progression of atherosclerosis in animal models [##REF##16682613##5##, ####REF##12600915##6##, ##REF##17909108##7##, ##REF##16880330##8####16880330##8##]. In humans, many immune cells such as macrophages, lymphocytes, mast cells, and T cells are found in atherosclerosis [##UREF##0##9##]. These findings suggest that the immune system plays important roles in atherogenesis. However, data available to date are primarily derived from studies of atherosclerosis in the coronary or/and the carotid arteries, whereas data derived from clinical samples of PAD appear to be particularly limited.</p> <p>In the past decade, microarray analysis using high-throughput screening technology has emerged as an important tool to study gene expression patterns and to study molecular events in complex diseases [##REF##17145989##10##, ####REF##12912960##11##, ##REF##17634102##12####17634102##12##]. In this study, Affymetrix GeneChips were used to perform gene expression profiling of femoral atherosclerotic lesions to fully characterize the peripheral arterial wall gene expression patterns associated with atherosclerosis. By statistical analysis, hundreds of known and novel genes were identified that differentially express in PAD. Genes involved in immune/inflammatory responses appeared to be significantly enriched in the set of genes up-regulated in different stages of PAD. To further examine the expression patterns of individual genes in the context of particular biological or molecular pathways, gene functional enrichment was performed using Gene Ontology and KEGG database. The results revealed that immune system-related categories and pathways were significantly overrepresented in the progression of the disease, suggesting that up-regulation of immune/inflammatory genes may be critical components of the disease progression expression signature associated with atherosclerosis. These findings may provide new insights and foster a better understanding of the mechanism of PAD.</p>
<title>Methods</title> <title>Tissue Harvest</title> <p>After obtaining informed consent, primary femoral artery specimens containing atherosclerotic lesions were taken from 30 patients undergoing surgical bypass or limb amputation at Shanghai Ninth People's Hospital. The specimens were immediately rinsed once with PBS and cut longitudinally by the surgeon. Three quarters of the samples were stored at once in – 80°C for subsequent total RNA extraction. The remaining samples were embedded in OCT medium and snap frozen for further morphological analysis. Clinical patient parameters were also registered. For controls, five normal femoral arteries were obtained from healthy donors during organ transplantation (male, mean 31.6 years; range 22–45 years). These five samples were without clinical or gross macroscopic signs of atherosclerotic disease. The Local Ethical Committee approved all procedures in this investigation, and proper protocol was followed throughout the entire course of the experiment.</p> <title>Histology</title> <p>For each sample, cryostat sections of 8 um were stained with hematoxylin for 10 min and eosin for 2 min, dehydrated in graded alcohol, and cover-slipped with permanent mounting solution after xylene clearing.</p> <title>RNA Isolation and Quantification</title> <p>Total RNA was isolated from the samples using a Trizol reagent (Invitrogen, Carlsbad, CA) and cleaned up using RNeasy Micro Kit (Qiagen, Valencia, CA) techniques. In brief, for each tissue, at least 100 mg sample was pulverized under liquid nitrogen. After complete disruption of the tissues, the Trizol reagent was added in the amount of 1 ml/100 mg. Total RNA was extracted using the protocol supplied with the Trizol reagent. After isolation, the RNA was cleaned up using the RNeasy Micro Kit. To remove any contaminating genomic DNA, a DNase step was included, following the manufacturer's protocol. The RNA quantity and quality were determined by an Agilent Bioanalyzer 2100 and an Eppendorf Biophotometer. Any RNA samples that showed degradation was excluded from the study.</p> <title>Microarray Experiment</title> <p>One microgram of total RNA was used for generating biotin labeled cRNA. The labeling reaction was performed according to the standard Affymetrix<sup>® </sup>protocol to generate a biotin-labeled cRNA probe. The samples were hybridized to the Affymetrix<sup>® </sup>Human Genome -U133A Genechip, stained, washed and scanned according to the standard Affymetrix<sup>® </sup>protocol. The computer data files to be used in data analysis (*.dat, *.cel, *.chp) were generated with the Affymetrix GeneChip Operating Software (GCOS) Version 1.4 (Affymetrix<sup>®</sup>), using the statistical algorithm provided. All chip samples were scanned using the same instrument and followed the same protocol. Data quality assessment was then performed following the guidance in Affymetrix data analysis fundamentals manual. All quality control results met Affymetrix recommended criteria.</p> <title>Data process and analysis</title> <p>The probe level intensity data were transferred to ArrayAssist<sup>® </sup>Software (StrataGene; La Jolla, CA) for further analysis. For comparison of differential gene expression between different stage groups, the background was removed and data were normalized in accordance to the GC-RMA method [##REF##16108723##43##]. GC-RMA takes into account the GC content of the probe sequences when comparing the expression intensities of the different probesets. Then, the processed gene expression data were transformed into log base 2 and filtered to delete the genes whose detection calls were "absent" in all samples.</p> <p>Microarray data analysis was carried out to identify individual genes that were significantly expressed between classes by the software package SAM (please see Availability &amp; requirements for more information), using Δ = 0.5. Results from the difference analysis were clustered and displayed using the Cluster3.0 and Treeview1.1.0 software (please see Availability &amp; requirements for more information). Each list of differentially expressed genes was analyzed in the context of Gene Ontology (GO) in order to identify groups of genes with similar functions, or processed using MAPPFinder (Gene MicroArray Pathway Profiler; please see Availability &amp; requirements for more information). For each gene ontology term, the probability values were computed based on a hypergeometric distribution test by comparing (a) the number of genes annotated by the gene ontology term in a given list of differentially expressed genes with (b) the expected number of such genes. <italic>Z-score&gt;0 and p</italic>-values &lt; 0.05 were considered significant categories.</p> <p>Similar methods were used to identify curated pathways that were significantly over-represented in the data using KEGG database by using DAVID (please see Availability &amp; requirements for more information). For each pathway, the probability values were computed based on a modified Fisher exact test. EASE <italic>p</italic>-values &lt; 0.05 were considered significant categories. The enriched pathways are not entirely separate from one another. For example, many genes involved in MAPK signaling pathway can also be involved in other pathways, such as NK pathway. The interconnectedness information was manually extracted from the pathway. Because the nature and complexity of these interactions varied from pathway to pathway, a simple line connecting two pathways was used to represent their interaction. The interaction map was generated for the interaction of enriched pathways using CytoScape software.</p> <p>Transcription factor enrichment analysis was also performed. The putative targets of transcription factors from TRANSFAC (v7.4) were discovered by Xie et al [##REF##15735639##44##] and downloaded from the supplementary web site (please see Availability &amp; requirements for more information). All the RefSeq IDs were converted to Entrez Gene ID according to the mapping table downloaded from NCBI web site (please see Availability &amp; requirements for more information). Enrichment of transcription factor targets was performed as described previously [##REF##15920519##45##]. The interaction map was generated for the interaction of enriched transcription factors and their putative target genes using CytoScape software</p> <title>Real-time QPCR Analysis</title> <p>One microgram of total RNA was reverse transcripted using random hexamers and superscript -II reverse transcriptase (Invitrogen, Carlsbad, CA). QPCR was performed by using ABI prism 7900 (ABI, Foster City, CA) and SYBR Green Detection (Toyobo, Japan). Primers were designed by using the Primer Express 2.0 software and verified by using a BLAST search. Sequences of the primers are listed [see Additional file ##SUPPL##11##12##]. The experimental conditions followed the manufacturer's protocol and the data were analyzed with sequence Detection Software 2.0 (ABI, Foster City, CA). Relative expression of mRNA was calculated with the comparative CT method. To standardize the amount of input RNA, the GAPDH gene was included. For each sample, the experiment was performed in triplicate.</p> <title>Western Blotting</title> <p>Proteins were extracted after RNA isolation according to the Introvigen protocol (Invitrogen, Carlsbad, CA) and measured using a Bio-Rad DC protein assay (Bio-Rad, Richmond, CA, USA). Aliquots of protein (100 μg of protein each) were resolved on a 10% SDS-PAGE gel and transferred to a polyvinylidene difluoride membrane (Millipore, Medford, MA, USA). The membrane was incubated with a primary antibody overnight at 4°C and then with a secondary antibody conjugated with alkaline phosphatase (1 h at room temperature), which was detected by a chemiluminescence method. The following polyclonal primary antibodies were used: anti-human TLR7 (1:300, IMGENEX, San Diego, CA), anti-human CTSS (1:400, Abcam Inc), anti-human beta-action (1:10000, Abcam Inc).</p> <title>Statistics</title> <p>The statistical significance of real-time results was examined with the nonparametric Mann-Whitney test, using GraphPad Prism 4. In the experiment, <italic>p </italic>values &lt; 0.05 were considered significantly different between the lesions group and the normal artery group.</p>
<title>Results</title> <title>Patient classification and outcome</title> <p>Histological characterization of 30 collected peripheral artery samples was conducted based on the criteria of the American Heart Association. Of these samples, 15 were classified as grade III (intermediate lesions), one as grade IV and fourteen as grade V (advanced lesions). Among them, 11 intermediate lesions samples (grade III) and 14 advanced lesions samples (grade V) had RNA of sufficient quality and quantity for hybridization. Representative images of the different stages are shown in Figure ##FIG##0##1##. Further details of these 25 samples are listed in Table ##TAB##0##1##. As shown, there was no significant difference between the intermediate lesions and the advanced lesions group except for indications of hypertension. In the intermediate lesions group, 4 patients (36.4%) presented with hypertension, while 9 hypertensive patients (64.3%) were found in the advanced lesions group.</p> <title>Differentially regulated genes in intermediate lesions</title> <p>Identifying differential expression genes was achieved for different stages by using Significance Analysis of Microarrays (SAM) with a false discovery rate (FDR) of 0.5%. Comparative analysis revealed that 366 genes were differentially expressed in intermediate lesions when compared to normal femoral arteries, of which 230 genes were up-regulated and 136 were down-regulated [see Additional file ##SUPPL##0##1##]. The 100 most differentially expressed genes between intermediate lesions and normal femoral arteries are shown in Figure ##FIG##1##2A##. Notably, in the up-regulated genes, up to 85 genes have been reported to be involved in immune response, such as <italic>HLA-DQB1, HLA-DRB1, CCR1, CXCR4, C1QB </italic>and <italic>TLR7 </italic>[see Additional file ##SUPPL##1##2##]. In addition, a large number of genes known to encode proteins crucial for proteolysis (<italic>CTSB, CTSC, CTSD and CTSS</italic>) and cell proliferation (<italic>BTG1, BTG2, CDKN1A, and MCM5</italic>) appeared to be significantly changed. Since BTG1 and BTG2 are known to be involved in anti-proliferation activities, it can be of interest to further investigate their potential roles in PAD in detail. MCM5 is heavily involved in chromosomal stability. Among the down-regulated genes, those involved in calcium signaling (<italic>CAMK2G)</italic>, transport (<italic>SLC22A3, CYP1A1 </italic>and <italic>ATP5H</italic>), metabolism (<italic>GCSH </italic>and <italic>PLA2G4A</italic>), and protein amino acid dephosphorylation (<italic>PTPN20</italic>) were found to be significantly down-regulated. The Gene Ontology functional categories in which intermediate lesions are overrepresented are illustrated in Table ##TAB##1##2## and additional data [see Additional file ##SUPPL##2##3##]. As shown in the table, the most significant biological process categories in the up-regulated genes are immune response, humoral immune response, inflammatory response, and T cell proliferation (Z-score&gt;5). For down-regulated genes, the significant ones mainly represent metabolism and catabolism-related categories.</p> <title>Differentially regulated genes in advanced lesions</title> <p>When advanced lesions were compared to normal femoral arteries, 447 genes were identified, of which 172 genes were up-regulated and 275 were down-regulated [see Additional file ##SUPPL##3##4##]. The list of the 100 most differentially expressed genes is shown in Figure ##FIG##1##2B##. Interestingly, up to 37 genes involved in the immune system response, such as <italic>CCR1, CX3CR1, TLR1 and TLR7</italic>, were found to be up-regulated in advanced lesion [see Additional file ##SUPPL##4##5##], which might suggest that these immune/inflammatory related genes could serve as expression signatures characterizing different stages of PAD. In addition, genes constituting a major portion of the vascular extracellular matrix were significantly up-regulated in advanced lesions, including <italic>COL1A1, COL3A1, COL1A2, COL5A1, COL6A1, COL6A3 and LAMB1</italic>, suggesting that these genes could be involved in the femoral artery occlusion in PAD. GO analysis further confirmed the above findings, by highlighting categories of immune response, humoral immune response, inflammatory response and I-kappaB kinase/NF-kappaB cascades (Z-score&gt;5) (Table ##TAB##2##3## and Additional file ##SUPPL##5##6##). For down-regulated genes, those involved in ion transport (<italic>GRIA2 and SLC22A3</italic>) and protein folding (<italic>DNAJB5</italic>) appeared to be the most significantly down-regulated in advanced lesions. GO analysis showed that the most significant categories for down-regulated genes were response to protein stimulus, RNA metabolism, and protein folding (Table ##TAB##2##3## and Additional file ##SUPPL##5##6##).</p> <p>In parallel, further data analysis revealed that many genes were over-represented in both intermediate and advanced lesions vs. normal controls. Of these genes, 68 were found to be commonly up-regulated and 48 were found to be commonly down-regulated (Figure ##FIG##2##3##). The list of commonly up-regulated genes is available [see Additional file ##SUPPL##6##7##]. Some of these overlapping genes, such as <italic>CTSB, CCR1, ALOX5, and SPP1</italic>, have been previously reported to play important roles in atherogenesis [##REF##16741146##13##, ####REF##16491201##14##, ##REF##16081317##15##, ##REF##15218333##16####15218333##16##]. Accordingly, these commonly regulated genes can therefore be important for the progression of PAD. In contrast, a much larger number of genes appear to be characteristically expressed in either intermediate lesions or advanced lesions, which may therefore serve as stage-specific signatures of PAD.</p> <title>Differential gene expression in disease progression</title> <p>Intermediate lesions and advanced lesions represent different stages in disease progression of PAD. Identification of genes that exhibit characteristic expression patterns in different stages may provide information relevant to the progression of PAD. For this reason, expression profiles of normal arteries, intermediate lesions and advanced lesions were analyzed by the SAM multiclass method. Out of this analysis, 614 genes appeared to be differentially expressed in the progression of PAD with a FDR&lt;0.5% [see Additional file ##SUPPL##7##8##]. Hierarchical clustering analysis suggested that the expression patterns of the genes could be assigned to three major groups (Figure ##FIG##3##4##). The first group represents those genes commonly expressed in both intermediate and advanced lesions (Cluster II). GO terms indicate that these genes are mainly involved in the immune response, inflammatory response, cellular defense and various signaling pathways (Table ##TAB##3##4## and Additional file ##SUPPL##8##9##). These results further support the notion that the immune system may play an important role in the progression of PAD. The second group represents specifically down-regulated genes in advanced lesions (Cluster II). GO terms indicate that these genes are primarily involved in cell cycle, apoptosis, multicellular organism development and protein folding (Table ##TAB##3##4## and Additional file ##SUPPL##8##9##). Genes in the third group are represented by those down-regulated in both intermediate lesions and advanced lesions (Cluster III). GO terms indicate that these genes are mainly involved in neurogenesis, protein modification, RNA splicing, and blood pressure regulation (Table ##TAB##3##4## and see Additional file ##SUPPL##8##9##). In addition, we have performed data analysis restricted to male subjects. Up to 85% genes identified in male subjects are the same as those identified in the total samples (data not shown), which suggests that the potential gender-biases is minimal. Taken together, genes commonly up-regulated in intermediate and advanced stages are typically represented by those involved in immune and inflammatory responses, implicating enhanced immune response activities during the progression of the disease, whereas down-regulated genes in the both disease stages are primarily represented by those involved in various aspects of cell proliferation and differentiation.</p> <title>Validation of gene transcription by real-time PCR</title> <p>Real-time PCR is still the gold standard for quantitative analysis of mRNA. In order to validate the microarray results, RT- PCR was carried out on the same set of samples that were analyzed by the microarray approach. The results were highly correlated with those from the array data. (The correlation coefficient for microarray and RT-PCR was 0.835 ± 0.076). Representative RT-PCR results of 6 genes are shown in Figure ##FIG##4##5##.</p> <title>Transcription factors enrichment analysis</title> <p>Transcription factors appear to play important roles in the development or progression of atherosclerosis [##REF##17084284##17##,##REF##16873729##18##]. To address whether specific transcription factors are involved in the regulation of genes associated with the progression of PAD, we conducted a transcription factor binding site enrichment study by analyzing cross-species conserved binding sites in promoter regions of genes differentially regulated during progression of PAD. Through the Fisher Exact test, binding sites of transcription factor AP-1 and CREB appeared to be significantly enriched (<italic>q-value </italic>&lt; 0.05). AP-1 is a transcription factor known to be involved in various cellular processes. In atherosclerosis, it has been reported in gene regulation of microphages, vascular smooth muscle cells and epithelial cells [##REF##18309110##19##,##REF##18187666##20##]. In disease progression, AP-1 was enriched to regulate expression of 72 genes (Figure ##FIG##5##6A## and Additional file ##SUPPL##9##10##). The enrichment of AP1 binding sites in regulated genes associated with PAD progression may therefore suggest an important role played by this transcription factor in the development of PAD. Through literature mining, indeed, some of the potential targets of AP1 appear to be previously reported as target genes of AP-1 [##REF##15467434##21##, ####REF##9520467##22##, ##REF##12899698##23##, ##REF##18367585##24####18367585##24##]. CREB is a member of the leucine zipper family of DNA binding proteins. This transcription factor binds as a homodimer to the cAMP-responsive element and induces transcription of genes in response to hormonal stimulation [##REF##17993258##25##]. A total of 55 genes were recognized as potential targets of CREB (Figure ##FIG##5##6B## and Additional file ##SUPPL##10##11##). Although further studies are need to elucidate detailed roles played by AP1 and CREB in PAD progression, significantly enriched binding sites and highly correlated with signature genes of PAD progression suggest that these two transcription factors may play critical roles in the development of PAD.</p> <title>Pathways identification by overabundant genes</title> <p>A pathway analysis database, KEGG, was then applied to genes differentially regulated in intermediate and advanced lesions. Several overrepresented pathways were identified, and the enriched pathways appeared not to be independent of one another, many genes involved in one pathway could be also involved in another pathway. This interaction is illustrated in Figure ##FIG##6##7##, and pathway abbreviations can be found in Table ##TAB##4##5##. As demonstrated, many immune-related pathways were significantly over-represented in intermediate and/or advanced lesions including TLR, NK, BCR, FER, APP, CCC and LTEM pathways. These findings, on the one hand, provide evidence supporting previous hypotheses that immune/inflammatory responses play important roles in the development of PAD, and on the other hand, demonstrate that particular components of immune/inflammatory systems can be crucial for the genesis and progression of PAD. For instance, TLR and NK pathways are shown to be particularly overrepresented in both intermediate lesions and advanced lesions, highlighting their functional importance in the disease. The TLR pathway is shown in Figure ##FIG##7##8## with the differentially regulated genes indicated.</p> <title>Protein validation of TLR7 expression</title> <p>Members of the Toll receptor family are key mediators of innate immunity. They respond to various pathogen-associated stimuli and transduce complex signaling responses that are required for inflammation and for the subsequent development of adaptive immunity [##REF##16917510##26##]. In atherosclerosis, TLR-mediated signaling cascades are observed in macrophages, mast cells and endothelial cells [##REF##11861602##27##,##REF##11160251##28##]. Data shown in this setting demonstrate that genes involved in TLR-mediated pathway are significantly up-regulated in intermediate or advanced lesions, including <italic>TLR1, TLR2, TLR7</italic>, and <italic>MyD88</italic>. <italic>TLR1 </italic>and <italic>TLR2 </italic>have been previously reported to be significantly regulated in atherosclerosis and their functional roles have been widely investigated in atherosclerosis [##REF##16211093##29##,##REF##16982924##30##]. However, the expression of <italic>TLR7 </italic>in atherosclerosis has not been reported before. TLR7 mediates innate responses by recognizing oligonucleotide based (RNA-) molecular patterns in endocytic compartments. Our data show that it is significantly up-regulated in both intermediate and advanced lesions. Western-blot analysis was performed to further validate its expression on the protein level (Figure ##FIG##8##9##). The function of TLR7 in atherogenesis is currently under further investigation.</p>
<title>Discussion</title> <p>In the present study, we first examined the gene expression profiles of PAD. Data analysis identified a number of genes that might be significantly correlated with different levels of PAD severity. The list of differentially expressed genes in intermediate and advanced lesions contains many genes which can be important for atherosclerosis. Most of these genes have not been reported to be related to atherosclerosis before. For example, <italic>MAP4K4 </italic>is a member of the serine/threonine protein kinase family. It has been shown to specifically activate MAPK8/JNK and mediate the TNF-alpha signaling pathway [##REF##14966141##31##,##REF##9890973##32##]. In this study, it was significantly and consistently up-regulated in both intermediate and advanced lesions.</p> <p>A large multidisciplinary study is currently underway to comprehensively assess PAD at multiple levels [##REF##17544028##33##], The goal of that study is to investigate 300 symptomatic patients with PAD undergoing medical management with or without vascular intervention by lower extremity angioplasty/stenting or vein graft bypass, and to test the hypothesis that the systemic inflammatory response after vascular intervention influences the local milieu responsible for vascular repair and adaptation [##REF##17544028##33##]. Identification of genes through the work may be significant in the selection of candidate genes that can be investigated through these cases-control genetic epidemiology studies. Our research supports the idea that immune responses play a key role in the development of PAD.</p> <p>In this report, immune related genes were shown to be significantly expressed during the development of PAD. Gene functional analysis further revealed that immune related categories and pathways were significant enriched in the different stages of PAD. In these immune related genes, several genes have been shown to modulate the development of atherosclerosis in mice models. For example, IgG Fc receptors (FcgammaRs) play a role in activating the immune system and in maintaining peripheral tolerance. Previous research suggested that <italic>Fc</italic>γ receptor deficiency protects against atherosclerosis in Apolipoprotein-E knockout mice [##REF##17053192##34##]. The results suggest that broad-range inhibitors of immune and inflammatory responses can be considered as potential targets for the treatment of PAD. However, gene expression patterns of immune related genes can be different in different stages of PAD. For example, in intermediate lesions, MHC class II molecules were significantly up-regulated including <italic>HLA-DMA, HLA-DMB, HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DRB1 and HLA-DRB5</italic>. MHC class II molecules are normally restricted to a subset of antigen presenting dendritic cells, B cells, macrophages, and thymic epithelium cells [##REF##10837054##35##]. These cells can be detected close to CD4+ T cells and present peptides to the T cells. The results suggest that there can be an ongoing immune activation in the intermediate lesions. However, MHC class II molecules were not differentially expressed in advanced lesions, even with a higher false discovery rate, which may suggest that the HLA-mediated immune activation may occur mainly in the progression stages of PAD. In addition, complement molecules were also significantly up-regulated in intermediate lesions, not in advanced lesions. Previous studies have implicated that activation of the complement system is probably associated with the initiation and progression of atherosclerosis [##REF##16516969##36##,##REF##11287977##37##]. Our data thus provide direct evidence from clinical samples demonstrating that complement system mainly play a role in the development stages of PAD. It is therefore conceivable that different and complex immune/inflammatory responses may take place at different stages of PAD.</p> <p>Atherosclerosis is a systemic, multifocal disease leading to various symptoms and clinical events including cardiovascular disease, cerebrovascular disease, and peripheral arterial disease. Our results reveal that many genes identified in the report are also expressed in coronary or carotid atherosclerotic lesions. For example, <italic>C3AR1 </italic>and <italic>C5R1 </italic>are receptors of C3 (C3a) and C5a respectively. A recent study shows that <italic>C3AR1 </italic>and <italic>C5R1 </italic>are expressed in human atherosclerotic coronary plaques [##REF##17234193##38##]. Double immunofluorescence staining has shown that the plaque of cells that express both <italic>C3aR </italic>and <italic>C5aR </italic>are macrophages, T cells, endothelial cells, and sub-endothelial smooth muscle cells. In addition, gene expression changes between atherosclerosis from coronary and carotid artery samples have been measured by microarray technology in recent years. One study using microarray found that 82 genes were differentially expressed in both animal model and human coronary artery atherosclerosis disease [##REF##15870398##39##]. Our data confirmed 29 genes and 18 genes had significantly different expression in intermediate lesions and advanced lesions, respectively. Moreover, these genes had expression trends similar to the ones found in our data, but our data showed higher fold-changes. In these overlapping genes, 14 were reported to be involved in immune response. Another microarray study found that 206 genes were differentially expressed in aortic atherosclerosis samples [##REF##15297278##40##]. Our data confirms 43 genes and 32 genes had significantly different expression in intermediate and advanced lesions (FDR&lt;1%), respectively. Importantly, in these overlapping genes, 15 were reported to be involved in immune response. Taken together, the results suggested that immune response is a common feature in atherosclerosis-related diseases. Our microarray study differs from prior microarray studies in the array type, sample type, sample classification, and analytical techniques. Nevertheless, the high level of overlapping genes suggests that there are similar molecular mechanisms in the development of peripheral arterial disease and other atherosclerosis-related diseases.</p> <p>Several limitations of our approach should be noted. First, hybridization-based microarrays, despite their immense potential, have inherent shortcomings related to deficient standardization of methods employed in normalization, statistical analysis, and so on [##REF##15904488##41##,##REF##17211383##42##]. In this study, we have attempted to limit these shortcomings by selecting subjects who were phenotypically similar to each other except for hypertension. In addition, the initial phases of data analysis, we used different normalization and statistical methods to identify differentially expressed genes. After choosing SAM, we used a rigorous false discovery rate to minimize false positive results. Expression patterns were validated by confirming mRNA expression patterns with conventional molecular techniques. We attempted, based on current literature, to suggest a potential functional role for genes whose expression was markedly altered. Second, atherosclerosis is a slow, progressive disease that may start in childhood; entirely normal arteries can only be obtained from young donors, a factor that can affect gene expression measurements. Although previous research and our data analysis suggest that age had very little effect on genes, further work is needed to identify age-related genes. Third, the relatively small number of patients did not allow us to assess serial changes in the disease development in more detail as would have been possible in animal models [##REF##16516969##36##]. Furthermore, we do not know to what extent the observed changes in gene expression translate into protein synthesis and function, and which genes cause atherosclerosis. Future studies are needed to address these issues.</p>
<title>Conclusion</title> <p>We first examined the gene expression profiles of PAD; the results from this analysis provide an initial step towards a better understanding of molecular mechanisms underlying PAD development. Differences in immune-related responses were observable at the gene expression level. These findings may be significant for understanding the molecular basis of PAD and investigating pharmacological approaches for the prevention and amelioration of atherosclerosis in PAD.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Peripheral arterial disease (PAD), a major manifestation of atherosclerosis, is associated with significant cardiovascular morbidity, limb loss and death. However, mechanisms underlying the genesis and progression of the disease are far from clear. Genome-wide gene expression profiling of clinical samples may represent an effective approach to gain relevant information.</p> <title>Results</title> <p>After histological classification, a total of 30 femoral artery samples, including 11 intermediate lesions, 14 advanced lesions and 5 normal femoral arteries, were profiled using Affymetrix microarray platform. Following real-time RT-PCR validation, different algorithms of gene selection and clustering were applied to identify differentially expressed genes. Under a stringent cutoff, i.e., a false discovery rate (FDR) &lt;0.5%, we found 366 genes were differentially regulated in intermediate lesions and 447 in advanced lesions. Of these, 116 genes were overlapped between intermediate and advanced lesions, including 68 up-regulated genes and 48 down-regulated ones. In these differentially regulated genes, immune/inflammatory genes were significantly up-regulated in different stages of PAD, (85/230 in intermediate lesions, 37/172 in advanced lesions). Through literature mining and pathway analysis using different databases such as Gene Ontology (GO), and the Kyoto Encyclopedia of Gene and Genomics (KEGG), genes involved in immune/inflammatory responses were significantly enriched in up-regulated genes at different stages of PAD(p &lt; 0.05), revealing a significant correlation between immune/inflammatory responses and disease progression. Moreover, immune-related pathways such as Toll-like receptor signaling and natural killer cell mediated cytotoxicity were particularly enriched in intermediate and advanced lesions (P &lt; 0.05), highlighting their pathogenic significance during disease progression.</p> <title>Conclusion</title> <p>Lines of evidence revealed in this study not only support previous hypotheses, primarily based on studies of animal models and other types of arterial disease, that inflammatory responses may influence the development of PAD, but also permit the recognition of a wide spectrum of immune/inflammatory genes that can serve as signatures for disease progression in PAD. Further studies of these signature molecules may eventually allow us to develop more sophisticated protocols for pharmaceutical interventions.</p>
<title>Availability &amp; requirements</title> <p>SAM software package: <ext-link ext-link-type="uri" xlink:href="http://www-stat.stanford.edu/~tibs/SAM/"/></p> <p>Cluster3.0 and Treeview1.1.0 software: <ext-link ext-link-type="uri" xlink:href="http://bonsai.ims.u-tokyo.ac.jp/~mdehoon/software/cluster/software.htm"/></p> <p>MAPPFinder (Gene MicroArray Pathway Profiler): <ext-link ext-link-type="uri" xlink:href="http://www.genmapp.org"/></p> <p>DAVID: <ext-link ext-link-type="uri" xlink:href="http://david.abcc.ncifcrf.gov/"/></p> <p>Xie et al Supplementary Information: <ext-link ext-link-type="uri" xlink:href="http://www.broad.mit.edu/seq/HumanMotifs/"/></p> <p>NCBI web site: <ext-link ext-link-type="ftp" xlink:href="ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/gene2refseq.gz"/></p> <title>Authors' contributions</title> <p>SJF, experiment design and conduction, and manuscript drafting; HGZ, clinical sample processing; JTS, transcription factor enrichment analysis; AA, WPK, KC and LO–M, data analysis and manuscript revising; JQZ, data analysis; YZD, data analysis and literature support; JZ and MEJ, experimental design and manuscript revising; JGJ, experiment design, data analysis and manuscript revising. All the authors read and approved this version of the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported in part by grants from Science and Technology Commission of Shanghai Municipality (project 044319209, 04JC14084 and 06QA14059), Chinese National Key Program for Basic Research (973: 2006CB910405), Chinese National High Tech Program (863: 2007AA02Z335 and 863: 2006AA02Z332), 100-Talent and Knowledge Innovation Programs of Chinese Academy of Science (J. Z), and the Komen Foundation (FAS0703850, WPK and LOM).</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Histological characteristics of various severities of femoral atherosclerotic lesions in PAD patients</bold>. HE stain analysis of histological characteristics of collected femoral non-atherosclerotic arteries and atherosclerotic arteries. 8 um cryostat sections were stained with hematoxylin and eosin, dehydrated in graded alcohol, and cover-slipped with permanent mounting solution after xylene clearing. Three representative samples are listed: normal artery (A, B, C), intermediate lesions (D, E, F), and advanced lesions (G, H, I).</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Heatmap of the 100 most differentially expressed genes in intermediate lesions and advanced lesions, respectively</bold>. SAM analysis reveals genes with differential expression in PAD. This analysis compared plaques from within arteries of either intermediate (n = 11) or advanced lesions (n = 14) to normal control group, respectively. Heatmap representation of the 100 most differentially expressed in intermediate lesions (A) and advanced lesions (B). Samples are displayed in columns and genes in rows. Gene expression is represented as a color, normalized across each row, with brighter red for higher values and brighter green for lower values. Gene symbols are listed to the right. N (Normal control group), Int (intermediate lesions group), Ad(advanced lesions group). The list of differentially expressed genes in intermediate lesions and advanced lesions is provided [see Additional file ##SUPPL##0##1## and Additional file ##SUPPL##3##4##].</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Over-represented genes in both intermediate lesions and advanced lesions</bold>. A. The genes whose expressions were significantly changed in intermediate lesions and advanced lesions, respectively, are shown in a Venn diagram. B, C. Heatmap representation of commonly up-regulated genes (B) and commonly down-regulated genes (C) in overlapping genes, respectively. Samples are displayed in columns and genes in rows. Gene expression is represented as a color, normalized across each row, with brighter red for higher values and brighter green for lower values. Gene symbols are listed to the right. N (Normal control group), Int (intermediate lesions group), Ad (advanced lesions group).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Hierarchical clustering analysis of the differentially expressed genes in disease progression</bold>. The differentially expressed genes analyzed by hierarchical clustering method in disease progression. The genes were classified into three major clusters by visual inspection. Clustering method: Average linking; Similarity measure: Euclidean distance. Samples are displayed in columns and genes in rows. Gene expression is represented as a color, normalized across each row, with brighter red for higher values and brighter green for lower values. N (Normal control group), Int (intermediate lesions group), Ad (advanced lesions group). The list of differentially expressed genes in disease progression is provided [see Additional file ##SUPPL##7##8##].</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p><bold>Real-time PCR and the relative expression level of six genes</bold>. Six genes mRNA in normal femoral artery (N, black round), intermediate lesions (Int, ascending triangle) and advanced lesions (Ad, descending triangle) were determined by real-time PCR and presented as a ratio to GAPDH mRNA. mRNA abundance in intermediate lesions or advanced lesions was differentially expressed (*P &lt; 0.05, and **P &lt; 0.01, respectively) when normal samples were used for comparison.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p><bold>Enriched transcription factors and their putative target genes in disease progression</bold>. Putative targets of transcription factors were curated based on results by Xiaohui Xie. Fisher Exact test showed that two transcription factors (AP-1 and CREB) were significantly enriched in disease progression (q-value &lt;0.05). AP-1 and CREB were enriched to regulate 72 and 55 genes expression, respectively. A, B. The top 20 putative target genes of AP-1 and CREB were listed, respectively. A list of the enriched transcription factor and their putative targets is provided [see Additional file ##SUPPL##9##10## and Additional file ##SUPPL##10##11##].</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p><bold>Interactions of KEGG pathways for differentially expressed genes in different stages of PAD</bold>. Pathways are enriched in intermediate lesions (A) and advanced lesions (B), respectively. Many genes involved in one pathway could also be involved in another pathway. A, B. Networking displayed the interaction of pathways in intermediate lesions and advanced lesions, respectively.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p><bold>Over-representation of Toll-like receptor signaling pathway genes</bold>. Analysis of over-representation of differentially expressed genes in pathway from KEGG. The Toll-like receptor signaling pathway is illustrated with significantly regulated genes highlighted.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p><bold>Significant expression of TLR7 in femoral atherosclerotic lesions</bold>. Western blot analysis of TLR7 in atherosclerotic femoral arteries (As, n = 3) and normal femoral arteries (N, n = 3). The protein level of Cathepsin S, whose expression was previously validated in atherosclerosis, was also examined in femoral atherosclerotic lesions. Beta – actin served as a loading control.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Patient characteristics for the 25 samples in the microarray analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left"><bold>Intermediate lesions group (n = 11)</bold></td><td align="left"><bold>Advanced lesions group (n = 14)</bold></td></tr></thead><tbody><tr><td align="left">Age*</td><td align="left">76.1 ± 4.68</td><td align="left">76.6 ± 3.79</td></tr><tr><td align="left">Male, n (%)</td><td align="left">8(72.7)</td><td align="left">12(85.7)</td></tr><tr><td align="left">Hypertension<sup>†</sup>, n (%)</td><td/><td/></tr><tr><td align="left">Yes</td><td align="left">4(36.4)</td><td align="left">9(64.3)</td></tr><tr><td align="left">no</td><td align="left">7(63.6)</td><td align="left">5(35.7)</td></tr><tr><td align="left">Hypercholesterolemia, n (%)</td><td/><td/></tr><tr><td align="left">Yes</td><td align="left">9(81.8)</td><td align="left">11(78.6)</td></tr><tr><td align="left">no</td><td align="left">2(18.2)</td><td align="left">3(21.4)</td></tr><tr><td align="left">Diabetes, n (%)</td><td/><td/></tr><tr><td align="left">Yes</td><td align="left">7(63.6)</td><td align="left">10(71.4)</td></tr><tr><td align="left">no</td><td align="left">4(36.4)</td><td align="left">4(28.6)</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Biological process categories overrepresented in intermediate lesions</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>GO Name</bold></td><td align="right"><bold>Z-score</bold></td><td align="right"><bold><italic>P</italic>-value</bold></td></tr></thead><tbody><tr><td align="left"><bold>For up-regulated genes </bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">immune response</td><td align="right">14.491</td><td align="right">0.000</td></tr><tr><td align="left">antigen processing, endogenous antigen via MHC class I</td><td align="right">3.821</td><td align="right">0.014</td></tr><tr><td align="left">antigen processing, exogenous antigen via MHC class II</td><td align="right">5.968</td><td align="right">0.001</td></tr><tr><td align="left">cellular defense response</td><td align="right">5.669</td><td align="right">0.000</td></tr><tr><td align="left">cell-mediated immune response</td><td align="right">3.093</td><td align="right">0.038</td></tr><tr><td align="left">T-helper 1 type immune response</td><td align="right">3.214</td><td align="right">0.035</td></tr><tr><td align="left">humoral immune response</td><td align="right">9.451</td><td align="right">0.000</td></tr><tr><td align="left">complement activation</td><td align="right">3.608</td><td align="right">0.016</td></tr><tr><td align="left">inflammatory response</td><td align="right">8.515</td><td align="right">0.000</td></tr><tr><td align="left">neutrophil chemotaxis</td><td align="right">6.629</td><td align="right">0.004</td></tr><tr><td align="left">leukocyte adhesion</td><td align="right">4.767</td><td align="right">0.013</td></tr><tr><td align="left">immune cell activation</td><td align="right">3.268</td><td align="right">0.011</td></tr><tr><td align="left">immune cell migration</td><td align="right">6.25</td><td align="right">0.000</td></tr><tr><td align="left">cell surface receptor linked signal transduction</td><td align="right">2.082</td><td align="right">0.042</td></tr><tr><td align="left">integrin-mediated signaling pathway</td><td align="right">2.633</td><td align="right">0.035</td></tr><tr><td align="left">intracellular signaling cascade</td><td align="right">4.367</td><td align="right">0.000</td></tr><tr><td align="left">MAPKKK cascade</td><td align="right">2.413</td><td align="right">0.037</td></tr><tr><td align="left">myeloid cell differentiation</td><td align="right">2.876</td><td align="right">0.042</td></tr><tr><td align="left">lipid transport</td><td align="right">3.255</td><td align="right">0.011</td></tr><tr><td align="left">Phagocytosis</td><td align="right">4.757</td><td align="right">0.005</td></tr><tr><td align="left">Apoptosis</td><td align="right">2.091</td><td align="right">0.049</td></tr><tr><td align="left">cell proliferation</td><td align="right">4.257</td><td align="right">0.000</td></tr><tr><td align="left">T cell proliferation</td><td align="right">5.489</td><td align="right">0.003</td></tr><tr><td align="left">oxygen transport</td><td align="right">3.255</td><td align="right">0.016</td></tr><tr><td align="left">icosanoid metabolism</td><td align="right">3.524</td><td align="right">0.010</td></tr><tr><td align="left">lipoprotein metabolism</td><td align="right">3.291</td><td align="right">0.017</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><bold>For down-regulated genes</bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">cellular catabolism</td><td align="right">2.708</td><td align="right">0.014</td></tr><tr><td align="left">amine catabolism</td><td align="right">3.599</td><td align="right">0.011</td></tr><tr><td align="left">glycine catabolism</td><td align="right">8.457</td><td align="right">0.003</td></tr><tr><td align="left">glycine metabolism</td><td align="right">5.817</td><td align="right">0.008</td></tr><tr><td align="left">ubiquitin cycle</td><td align="right">2.092</td><td align="right">0.041</td></tr><tr><td align="left">protein folding</td><td align="right">3.088</td><td align="right">0.008</td></tr><tr><td align="left">Endocytosis</td><td align="right">2.521</td><td align="right">0.032</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Biological process categories overrepresented in advanced lesions</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>GO Name</bold></td><td align="right"><bold>Z-score</bold></td><td align="right"><bold><italic>P</italic>-value</bold></td></tr></thead><tbody><tr><td align="left"><bold>For up-regulated genes </bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">immune response</td><td align="right">6.305</td><td align="right">0.000</td></tr><tr><td align="left">cellular defense response</td><td align="right">2.94</td><td align="right">0.026</td></tr><tr><td align="left">humoral immune response</td><td align="right">5.392</td><td align="right">0.000</td></tr><tr><td align="left">inflammatory response</td><td align="right">6.179</td><td align="right">0.000</td></tr><tr><td align="left">immune cell migration</td><td align="right">4.778</td><td align="right">0.014</td></tr><tr><td align="left">intracellular signaling cascade</td><td align="right">4.297</td><td align="right">0.001</td></tr><tr><td align="left">small GTPase mediated signal transduction</td><td align="right">2.372</td><td align="right">0.024</td></tr><tr><td align="left">I-kappaB kinase/NF-kappaB cascade</td><td align="right">6.408</td><td align="right">0.000</td></tr><tr><td align="left">T cell proliferation</td><td align="right">4.186</td><td align="right">0.016</td></tr><tr><td align="left">induction of apoptosis</td><td align="right">3.709</td><td align="right">0.007</td></tr><tr><td align="left">oxygen transport</td><td align="right">4.778</td><td align="right">0.003</td></tr><tr><td align="left">amine transport</td><td align="right">3.051</td><td align="right">0.024</td></tr><tr><td align="left">anion transport</td><td align="right">3.373</td><td align="right">0.006</td></tr><tr><td align="left">phosphate transport</td><td align="right">6.139</td><td align="right">0.000</td></tr><tr><td align="left">organic acid transport</td><td align="right">2.634</td><td align="right">0.039</td></tr><tr><td align="left">nucleoside monophosphate metabolism</td><td align="right">5.853</td><td align="right">0.001</td></tr><tr><td align="left">cell adhesion</td><td align="right">3.958</td><td align="right">0.000</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><bold>For down-regulated genes</bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">response to protein stimulus</td><td align="right">6.79</td><td align="right">0.000</td></tr><tr><td align="left">response to unfolded protein</td><td align="right">6.79</td><td align="right">0.000</td></tr><tr><td align="left">integrin-mediated signaling pathway</td><td align="right">4.36</td><td align="right">0.003</td></tr><tr><td align="left">Apoptosis</td><td align="right">2.838</td><td align="right">0.003</td></tr><tr><td align="left">muscle cell differentiation</td><td align="right">3.771</td><td align="right">0.032</td></tr><tr><td align="left">vasculature development</td><td align="right">3.795</td><td align="right">0.005</td></tr><tr><td align="left">neuron differentiation</td><td align="right">3.993</td><td align="right">0.003</td></tr><tr><td align="left">vitamin transport</td><td align="right">3.771</td><td align="right">0.023</td></tr><tr><td align="left">RNA metabolism</td><td align="right">3.476</td><td align="right">0.000</td></tr><tr><td align="left">RNA splicing</td><td align="right">3.497</td><td align="right">0.005</td></tr><tr><td align="left">glutamine metabolism</td><td align="right">4.901</td><td align="right">0.004</td></tr><tr><td align="left">macromolecule metabolism</td><td align="right">2.214</td><td align="right">0.023</td></tr><tr><td align="left">protein folding</td><td align="right">5.643</td><td align="right">0.000</td></tr><tr><td align="left">cell-matrix adhesion</td><td align="right">2.605</td><td align="right">0.036</td></tr><tr><td align="left">peptide hormone secretion</td><td align="right">5.401</td><td align="right">0.004</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>GO discovered categories for disease progression manner analysis</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>GO Name</bold></td><td align="right"><bold>Z-Score</bold></td><td align="right"><bold><italic>P</italic>-value</bold></td></tr></thead><tbody><tr><td align="left"><bold>For Cluster I</bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">immune response</td><td align="right">8.908</td><td align="right">0.000</td></tr><tr><td align="left">inflammatory response</td><td align="right">7.643</td><td align="right">0.000</td></tr><tr><td align="left">cellular defense response</td><td align="right">5.938</td><td align="right">0.001</td></tr><tr><td align="left">chemotaxis</td><td align="right">7.611</td><td align="right">0.000</td></tr><tr><td align="left">cell surface receptor linked signal transduction</td><td align="right">2.937</td><td align="right">0.009</td></tr><tr><td align="left">cytokine and chemokine mediated signaling pathway</td><td align="right">4.378</td><td align="right">0.004</td></tr><tr><td align="left">JAK-STAT cascade</td><td align="right">3.254</td><td align="right">0.020</td></tr><tr><td align="left">MAPKKK cascade</td><td align="right">3.062</td><td align="right">0.006</td></tr><tr><td align="left">endocytosis</td><td align="right">1.956</td><td align="right">0.049</td></tr><tr><td align="left">angiogenesis</td><td align="right">4.452</td><td align="right">0.000</td></tr><tr><td align="left">cell adhesion</td><td align="right">4.580</td><td align="right">0.000</td></tr><tr><td align="left">ion homeostasis</td><td align="right">3.090</td><td align="right">0.010</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><bold>For Cluster II</bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">apoptosis</td><td align="right">3.871</td><td align="right">0.002</td></tr><tr><td align="left">cell cycle</td><td align="right">4.108</td><td align="right">0.000</td></tr><tr><td align="left">cell growth</td><td align="right">2.285</td><td align="right">0.037</td></tr><tr><td align="left">epidermal cell differentiation</td><td align="right">3.075</td><td align="right">0.044</td></tr><tr><td align="left">protein folding</td><td align="right">3.677</td><td align="right">0.001</td></tr><tr><td align="left">macromolecule metabolic process</td><td align="right">3.022</td><td align="right">0.005</td></tr><tr><td align="left">RNA metabolic process</td><td align="right">5.378</td><td align="right">0.000</td></tr><tr><td align="left">transcription</td><td align="right">4.593</td><td align="right">0.000</td></tr><tr><td align="left">peripheral nervous system development</td><td align="right">2.950</td><td align="right">0.046</td></tr><tr><td align="left">multicellular organismal development</td><td align="right">3.927</td><td align="right">0.000</td></tr><tr><td align="left">vitamin transport</td><td align="right">4.426</td><td align="right">0.008</td></tr><tr><td align="left">response to stimulus</td><td align="right">4.579</td><td align="right">0.000</td></tr><tr><td align="left">circadian rhythm</td><td align="right">5.090</td><td align="right">0.001</td></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left"><bold>For Cluster III</bold></td><td/><td/></tr><tr><td colspan="3"><hr/></td></tr><tr><td align="left">neurogenesis</td><td align="right">2.612</td><td align="right">0.014</td></tr><tr><td align="left">neuron recognition</td><td align="right">3.315</td><td align="right">0.046</td></tr><tr><td align="left">cell motility</td><td align="right">2.093</td><td align="right">0.042</td></tr><tr><td align="left">glycogen metabolic process</td><td align="right">2.999</td><td align="right">0.025</td></tr><tr><td align="left">protein modification process</td><td align="right">3.166</td><td align="right">0.001</td></tr><tr><td align="left">RNA processing</td><td align="right">2.183</td><td align="right">0.042</td></tr><tr><td align="left">RNA splicing</td><td align="right">2.591</td><td align="right">0.021</td></tr><tr><td align="left">phosphate metabolic process</td><td align="right">2.346</td><td align="right">0.024</td></tr><tr><td align="left">blood pressure regulation</td><td align="right">4.028</td><td align="right">0.008</td></tr><tr><td align="left">heart contraction</td><td align="right">2.901</td><td align="right">0.022</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>KEGG biological pathways for differentially expressed genes in different stages</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>KEGG Pathway Name</bold></td><td align="left"><bold>Pathway ID</bold></td><td align="right"><bold>Genes involved </bold></td><td align="right"><bold><italic>P</italic>-Value</bold></td></tr></thead><tbody><tr><td align="left"><bold>For intermediate lesions</bold></td><td/><td/><td/></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left">Type I Diabetes Mellitus(TIDM)</td><td align="left">Hsa04940</td><td align="right">9</td><td align="right">0.002</td></tr><tr><td align="left">Antigen Processing and Presentation(APP)</td><td align="left">Hsa04612</td><td align="right">13</td><td align="right">0.003</td></tr><tr><td align="left">Complement and Coagulation Cascades(CCC)</td><td align="left">Hsa04610</td><td align="right">12</td><td align="right">0.003</td></tr><tr><td align="left">Cell Adhesion Molecules(CAM)</td><td align="left">Hsa04514</td><td align="right">17</td><td align="right">0.005</td></tr><tr><td align="left">Toll-Like Receptor Signaling Pathway(TLR)</td><td align="left">Hsa04620</td><td align="right">13</td><td align="right">0.013</td></tr><tr><td align="left">Natural Killer Cell Mediated Cytotoxicity(NK)</td><td align="left">Hsa04650</td><td align="right">15</td><td align="right">0.029</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><bold>For advanced lesions</bold></td><td/><td/><td/></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left">Focal Adhesion(FA)</td><td align="left">Hsa04510</td><td align="right">25</td><td align="right">0.000</td></tr><tr><td align="left">ECM-Receptor Interaction(ECM)</td><td align="left">Hsa04512</td><td align="right">15</td><td align="right">0.000</td></tr><tr><td align="left">Toll-Like Receptor Signaling Pathway(TLR)</td><td align="left">Hsa04620</td><td align="right">13</td><td align="right">0.002</td></tr><tr><td align="left">Regulation of Actin Cytoskeleton(RAC)</td><td align="left">Hsa04810</td><td align="right">20</td><td align="right">0.005</td></tr><tr><td align="left">Fc Epsilon RI Signaling Pathway(FER)</td><td align="left">Hsa04664</td><td align="right">10</td><td align="right">0.012</td></tr><tr><td align="left">MAPK Signaling Pathway(MAPK)</td><td align="left">Hsa04010</td><td align="right">22</td><td align="right">0.019</td></tr><tr><td align="left">Natural Killer Cell Mediated Cytotoxicity(NK)</td><td align="left">Hsa04650</td><td align="right">13</td><td align="right">0.021</td></tr><tr><td align="left">Long-Term Potentiation(LTP)</td><td align="left">Hsa04720</td><td align="right">8</td><td align="right">0.030</td></tr><tr><td align="left">Cell Communication(CC)</td><td align="left">Hsa01430</td><td align="right">12</td><td align="right">0.032</td></tr><tr><td align="left">B Cell Receptor Signaling Pathway(BCR)</td><td align="left">Hsa04662</td><td align="right">8</td><td align="right">0.042</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional File 1</title><p><bold>Table 1</bold>. 366 differentially expressed genes in intermediate lesions relative to normal femoral arteries. Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Fold Change and q-value(%) which is the lowest FDR are listed in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional File 2</title><p><bold>Table 2</bold>. Immune-related genes in intermediate lesions relative to normal femoral arteries, Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Fold Change and q-value(%) which is the lowest FDR are listed in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional File 3</title><p><bold>Table 3</bold>. Cell component and molecular function categories overrepresented in intermediate lesions, the first half of the table indicates categories highly significant for up-regulated genes; the second half of the table shows categories highly significant for down-regulated genes. The calculated <italic>p</italic>-values and Z-scores for each category are shown.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional File 4</title><p><bold>Table 4</bold>. 447 differentially expressed genes in advanced lesions relative to normal femoral arteries. Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Fold Change and q-value(%) which is the lowest FDR are listed in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional File 5</title><p><bold>Table 5</bold>. Immune-related genes in advanced lesions relative to normal femoral arteries, Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Fold Change and q-value(%) which is the lowest FDR are listed in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S6"><caption><title>Additional File 6</title><p><bold>Table 6</bold>. Cell component and molecular function categories overrepresented in advanced lesions, the first half of the table indicates categories highly significant for up-regulated genes; the second half of the table shows categories highly significant for down-regulated genes. The calculated <italic>p</italic>-values and Z-scores for each category are shown.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S7"><caption><title>Additional File 7</title><p><bold>Table 7</bold>. 68 commonly up-regulated genes in intermediate lesions and advanced lesions, Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Fold Change and q-value(%) which is the lowest FDR are listed in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S8"><caption><title>Additional File 8</title><p><bold>Table 8</bold>. 614 differentially expressed genes in disease progression. Affymetrix Probe Set ID, Cluster ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Contrast which is the standardized mean difference between the gene's expression in that class versus its overall mean expression, and q-value(%)which is the lowest FDR are listed. In the table, contrast1, 2 and 3 represent the standardized mean difference in normal femoral arteries, intermediate lesions, and advanced lesions, respectively.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S9"><caption><title>Additional File 9</title><p><bold>Table 9</bold>. The detail GO overrepresented categories for disease progression in each cluster, the calculated <italic>p</italic>-values and Z-scores for each category are shown in the table.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S10"><caption><title>Additional File 10</title><p><bold>Table 10</bold>. 72 putative AP-1 target genes in disease progression. Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Contrast, and q-value(%) are listed.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S11"><caption><title>Additional File 11</title><p><bold>Table 11</bold>. 55 putative CREB target genes in disease progression. Affymetrix Probe Set ID, Gene Title, Gene Symbol, GO Biological Process, GO Molecular Function, GO Cellular Component, Unigene, Entrez Gene, Ensembl, Chromosome Number, Socre(d), Contrast, and q-value(%) were listed.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S12"><caption><title>Additional File 12</title><p><bold>Table 12</bold>. The primer sequences of selected genes for real-time PCR.</p></caption></supplementary-material>
<table-wrap-foot><p>*as mean ± SD.</p><p>† P &lt; 0.05 (Student t test)</p></table-wrap-foot> <table-wrap-foot><p>GO analysis was applied to differentially expressed genes in intermediate lesions compared with normal controls. Similar significant categories were not included to reduce redundancy. The first half of the table indicates categories highly significant for up-regulated genes; the second half of the table shows categories highly significant for down-regulated genes. Calculated <italic>p</italic>-values and Z-scores for each category are shown. The corresponding cell component categories and molecular function categories are available in additional data [see Additional file ##SUPPL##2##3##]</p></table-wrap-foot> <table-wrap-foot><p>GO analysis was applied to differentially expressed genes in advanced lesions compared with normal controls. Some other selected categories are shown. The first half of the table indicates categories highly significant for up-regulated genes; the second half of the table shows categories highly significant for down-regulated genes. Calculated <italic>p</italic>-values and Z-scores for each category are shown. The corresponding cell component categories and molecular function categories are available in additional data [see Additional file ##SUPPL##5##6##]</p></table-wrap-foot> <table-wrap-foot><p>GO analysis was applied to differentially expressed genes in each cluster group. Some selected categories are shown. The Calculated <italic>p</italic>-values and Z-scores for each category are shown. The detailed GO results are provided in the additional data [see Additional file ##SUPPL##8##9##].</p></table-wrap-foot> <table-wrap-foot><p>Genes identified from SAM under FDR&lt;1% were tested for overrepresentation within the pathways for the KEGG under the modified Fisher Exact test assumption; <italic>p</italic>-Value and genes involved in the pathway are provided for each pathway</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-369-1"/> <graphic xlink:href="1471-2164-9-369-2"/> <graphic xlink:href="1471-2164-9-369-3"/> <graphic xlink:href="1471-2164-9-369-4"/> <graphic xlink:href="1471-2164-9-369-5"/> <graphic xlink:href="1471-2164-9-369-6"/> <graphic xlink:href="1471-2164-9-369-7"/> <graphic xlink:href="1471-2164-9-369-8"/> <graphic xlink:href="1471-2164-9-369-9"/>
<media xlink:href="1471-2164-9-369-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S5.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S6.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S7.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S8.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S9.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S10.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S11.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-369-S12.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Hansson", "Robertson", "S\u00f6derberg-Naucl\u00e9r"], "given-names": ["G\u00f6ranK", "Anna-KarinL", "Cecilia"], "article-title": ["INFLAMMATION AND ATHEROSCLEROSIS"], "source": ["Annual Review of Pathology: Mechanisms of Disease"], "year": ["2006"], "volume": ["1"], "fpage": ["297"], "lpage": ["329"], "pub-id": ["10.1146/annurev.pathol.1.110304.100100"]}]
BMC Genomics. 2008 Aug 1; 9:369