| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Vascular Biology |
From the Donald W. Reynolds Cardiovascular Clinical Research Center, Division of Cardiovascular Medicine (R. Tabibiazar, R.A.W., J.M.S., E.A.A., P.S.T., T.Q.), Department of Health Research and Policy (B.N., B.E., R. Tibshirani), and the Department of Statistics (B.E., R. Tibshirani), Stanford University School of Medicine, Stanford, Calif; and the Genome Sciences Department (E.M.R.), Lawrence Berkeley National Laboratory, Berkeley, Calif.
Correspondence to Raymond Tabibiazar, Stanford Medical School, Division of Cardiovascular Medicine, 300 Pasteur Drive, Falk CVRC, Stanford, CA 94305 (E-mail rtabibiazar{at}cvmed.stanford.edu) or Thomas Quertermous, Stanford Medical School, Division of Cardiovascular Medicine, 300 Pasteur Drive, Falk CVRC, Stanford, CA 94305 (E-mail tomq1{at}stanford.edu).
| Abstract |
|---|
|
|
|---|
Methods and Results We performed genome-wide transcriptional profiling of aortas from C3H/HeJ and C57Bl/6 mice. Differences in gene expression were identified at baseline as well as during normal aging and longitudinal exposure to high-fat diet. The significance of these genes to the development of atherosclerosis was evaluated by observing their temporal pattern of expression in the well-studied apolipoprotein E model of atherosclerosis.
Conclusion Gene expression differences between the 2 strains suggest that aortas of C57Bl/6 mice have a higher genetic propensity to develop inflammation in response to appropriate atherogenic stimuli. This study expands the repertoire of factors in known disease-related signaling pathways and identifies novel candidate genes for future study.
To gain insights into the molecular pathways that are differentially activated in strains of mice with varied susceptibility to atherosclerosis, we performed comprehensive transcriptional profiling of their vascular wall. Genes identified through these studies expand the repertoire of factors in disease-related signaling pathways and identify novel candidate genes in atherosclerosis.
Key Words: atherosclerosis vascular disease gene expression microarray inflammation oxidative stress
| Introduction |
|---|
|
|
|---|
With the availability of genomic tools that allow genome-wide transcriptional profiling, individual environmental risk factors can be evaluated in a controlled setting to define gene expression differences that are associated with the risk, without a priori knowledge of specific candidate genes. Recent studies have shown that transcriptional profiling can be used to characterize gene expression patterns associated with inherited genetic differences.5 The integration of gene expression data with genetic analyses has been used to identify susceptibility loci for complex diseases, such as diabetes, obesity, and allergic asthma, and to prioritize candidate genes.6,7 Therefore, gene expression analysis of genetically distinct mice with different disease susceptibilities under controlled environmental conditions can potentially identify genes whose expression underlies differences in disease susceptibility.
Prior studies investigating genetic influences on atherosclerosis susceptibility in various mouse strains have provided evidence for involvement of the vascular wall and not peripheral blood cells or systemic cholesterol metabolism in the development of atherosclerosis.8 We thus hypothesized that genes differentially expressed in the vascular wall in genetically resistant C3H versus genetically susceptible C57 mice may play critical roles in the initiation and modulation of atherosclerosis. To identify these genes, we have performed genome-wide transcriptional profiling of mouse aortic tissue from these 2 mouse strains and compared their vascular wall gene expression. Using a comprehensive study design and novel analytical tools, we show that activity of specific biological processes and molecular functions distinguish disease-resistant versus disease-prone mouse strains. Further, the importance of these genes in development of atherosclerosis is validated with studies in the apoE model.
| Methods |
|---|
|
|
|---|
Data Processing
Array image acquisition and feature extraction was performed using the Agilent G2565AA Microarray Scanner and feature extraction software version A.6.1.1. Normalization was performed using a LOWESS algorithm, and Dye-normalized signals were used in calculating log ratios. Features with reference values of <2.5 SDs above background for the negative control features were regarded as missing values. Those features with values in at least two thirds of the experiments and present in at least 1 of the replicates were retained for further analysis. For significance analysis of microarrays (SAM), a K-nearest-neighbor algorithm was applied to impute for missing values.10
Data Analysis
Experimental design and analysis flow chart is depicted in Figure 1. SAM was used to identify genes with statistically different expression between the C3H and C57 mice at baseline.10,1214 For partitioning clustering of the genes with K-means and self-organizing maps (SOMs), we used positive correlation for distance determination and required complete linkage, which uses the greatest distance between genes to ascribe similarity. SOMs and K-means analyses were performed using Expressionist software (GeneData Inc).15 Heatmaps were generated using HeatMap Builder.14,16 For enrichment analysis we used the EASE analysis software, which uses gene ontology (GO) annotation and the Fisher exact test to derive biological themes within particular gene sets.17 For time-course study a new statistical algorithm, the area under curve (AUC) analysis, was devised. For each sequence of 4 triplicate gene expression measurements over time, we first subtracted the measurement at time 0 from all values. We then computed the signed AUC. The area is a natural measure of change over time. These areas were then used to compute an F statistic for comparing C57 and C3H mice across the different diets. A permutation analysis, similar to that used in SAM, was performed to estimate the false discovery rate (Q value or "FDR") for different levels of the F statistic. For ease of presentation, genes which meet our FDR cutoffs will be referred to as "significant" throughout the remainder of this article. All microarray data were submitted to the National Center for Biotechnology Information Gene Expression Omnibus (GEO GSE1560; http://www.ncbi.nlm.nih.gov/geo/).
|
Aortic Lesion Analysis
For select time points within various experimental groups, 5 to 7 female mice were used for histological lesion analysis. Atherosclerosis lesion area was determined as described previously.18
Quantitative Real-Time Reverse TranscriptasePolymerase Chain Reaction
Primers and probes for 10 representative differentially expressed genes were obtained from Applied Biosystems Assays-on-Demand. A Total of 90 reactions were performed from representative RNA samples used for microarray experiments. These included triplicate assay on 3 pools of 5 aortas. cDNA was synthesized and Taqman was performed as described previously.10
| Results |
|---|
|
|
|---|
|
Comparison of C3H and C57 vascular wall gene expression at baseline provided a list of compelling candidate genes, which reflected differences in biological processes such as growth, differentiation, and inflammation, as well as molecular functions, such as cathecholamine synthesis, phosphatase activity, peroxisome function, insulin-like growth factor activity, and antigen presentation (Figure 2). These processes were exemplified by higher expression of genes such as Cdkn1a, Pparbp, protein tyrosine phosphatase-4a2, and Socs5 in C3H mice, compared with genes such as ABCC1, H2-D1, Bat5, IGFBP1, SCD1, and Serpine6b, which demonstrated higher expression in C57 mice (Tables I and II, available online at http://atvb.ahajournals.org). These fundamental baseline gene expression differences may determine disease susceptibility as the mice are exposed to age-related stimuli or dietary challenges.
Age-Related Differences in Gene Expression Patterns Between the Mouse Strains
To further examine the vascular wall gene expression differences between C57 and C3H mice, we performed analyses to identify genes differentially expressed in response to aging (Figure 3). Data were collected at 5 time points over a 40-week period. To identify such genes, we developed the AUC analysis. The AUC analysis relies on a permutation procedure to reduce the number of potential false-positives generated because of multiple testing but still uses the increase in statistical power of time-course experimental design. Comparing C57 versus C3H time-course differences on normal diet with a rigid cutoff (FDR<0.05) did not identify any genes. However, relaxing the AUC stringency (F statistic>10, FDR<0.45) allowed a large number of genes (413) to be included for pathway overrepresentation analysis using GO annotation (Table III, available online at http://atvb.ahajournals.org). Functional annotation and group overrepresentation analysis (Fisher test probability value <0.02) of the resultant differentially expressed genes revealed differences in a number of biological processes, including growth and development as well as a number of molecular functions, such as cell cycle control, regulation of mitosis, and metabolism (Figure 3B). Some of these processes are exemplified by genes with higher expression in C57 mice, such as Aoc1 (prooxidative stress), Bub1 (cell cycle check point), Cyclin B2, and genes with higher expression in C3H, including INHBA and INHBB (Table IV, available online at http://atvb.ahajournals.org).
|
Temporally variable genes identified by AUC analysis were further characterized with K-means clustering to identify dynamic patterns of expression during the aging process (Figure 3C). Clusters 1, 4, and 9 revealed either higher overall expression or temporally increasing levels of expression in C3H mice compared with C57 mice. In contrast, clusters 2, 6, and 14 revealed the opposite pattern. Of the genes which were noted to be differentially expressed in the 2 strains during aging, 51 genes were also differentially expressed at baseline, suggesting that baseline differences of certain genes can further be affected with aging.
Diet-Related Differences in Gene Expression Patterns Between the Mouse Strains
Differential vascular wall response to atherogenic stimuli was determined by comparing temporal gene expression patterns in C57 versus C3H mice on high-fat diet (Figure 4A). Comparing C57 versus C3H time-course differences on high-fat diet with a rigid cutoff (FDR<0.05) identified 35 genes, including Hgfl and Tgfb4, which were downregulated in C57 on high-fat diet. Additional known genes as well as a number of expressed sequence tags were also identified (Table V, available online at http://atvb.ahajournals.org). Using a less stringent AUC cutoff allowed identification of a larger number of genes that could be evaluated with pathway overrepresentation analysis using GO annotation. At this level of stringency (F statistic>10, FDR<0.35), a total of 650 genes with temporally variable expression were identified (Table VI, available online at http://atvb.ahajournals.org). Genes that were also differentially regulated by the aging process (141 of 650 genes) were excluded from further analysis of this group. Thirty-eight of the remaining 509 genes were among those differentially expressed at baseline. Functional annotation and group overrepresentation analysis (Fisher test P<0.02) of these differentially expressed genes revealed differences in biological processes, such as catabolism, oxygen reactive species and superoxide metabolism, and proteo- and peptidolysis as well as molecular functions, such as fatty acid metabolism, oxidoreductase, and methyltransferase activities (Figure 4B). Interestingly, this analysis suggested important differences between the 2 mouse strains with respect to the activity of the peroxisome, microbody, and lysosome. Some of these processes were exemplified by genes with higher expression in C3H mice, such as Ccs, Ephx2, Gpx4, Prdx6 (antioxidants), Sirt3 (transcriptional repressor), PPAR
, and Mcd as well as genes with higher expression in C57 mice, such as Lysyl oxidase and Cdkn1a (Table VII, available online at http://atvb.ahajournals.org). K-means clustering of these genes identified a small number of distinct expression patterns (Figure 4C), with clusters 3 and 9 revealing increased gene expression in C3H mice and clusters 8 and 10 showing the opposite pattern.
|
Evaluation of Strain-Specific Differentially Regulated Genes in the ApoE Model
Using these techniques, we have identified a significant number of genes that are differentially expressed in the atherosclerosis-resistant C3H and susceptible C57 mice, some of which are likely involved in atherogenesis and some of which are likely irrelevant to the process. To further select genes most likely to be involved in atherogenesis, we examined their expression in apoE-deficient mice fed normal or high-fat diet over a period of 40 weeks (Figure 5). We used SOM analysis to visualize the expression profiles of these subsets of genes throughout the development and progression of atherosclerosis in the apoE-deficient mice. The analysis revealed several patterns of gene expression. For example, SOM cluster 8 demonstrated a consistently increasing pattern of expression which correlated with disease progression in the apoE-deficient mice (Figure 5). As evidenced by the pie chart, this cluster is enriched with genes that were identified as more highly expressed in C57 versus C3H mice at baseline (ie, potentially atherogenic). In contrast, clusters 4, 5, and 6 showed decreasing expression with disease progression. The decreased expression of genes in cluster 4 was somewhat attenuated with high-fat challenge of the apoE-deficient mice. This cluster was particularly enriched with genes that had revealed a higher expression in C3H mice (ie, potentially atheroprotective) with atherogenic stimuli and with aging.
|
Given C3H resistance and C57 susceptibility to atherosclerosis, as an initial hypothesis it was postulated that genes with higher expression in C3H mice confer resistance, whereas genes with higher expression in C57 mice may have a proatherogenic role. With this point of reference, gene clusters were further examined (Table X, available online at http://atvb.ahajournals.org). For example, limiting the list of genes in SOM cluster 8 (genes with increased expression with atherosclerosis) to those that also had higher baseline expression in C57 mice yielded an interesting set of genes that may be atherogenic. This group included inflammation related genes, such as H2-D1, Pdgfc, Paf, and Cd47 (Table Xa). Other compelling genes included Agpt2, Mglap, Xdh, Th, and Ctsc. Conversely, limiting the list of genes in clusters 4 and 5 to those with higher expression in C3H mice identified a group of genes with potential athero-protective function. Some of those genes included Ppar
, Pparbp as well as Ptp4a1, and Mcd (Table Xa).
Lesion Analysis in the Genetic Models
To address whether some of the gene expression differences are related to presence of atherosclerotic lesion in C57 mice, we determined the total atherosclerotic burden in the aorta by calculating a percent lesion area in aortas of C57 (n=5) and C3H (n=5) mice. Comparisons were made at time 0 and 40 weeks on normal or high-fat diet. Noncholate-containing high-fat diet was used to prevent caustic effects on the vascular wall. As expected, C57 and C3H mice on either diet did not demonstrate evidence of atherosclerosis throughout the course of the experiment, suggesting that observed gene expression changes cannot be explained by different cellular composition of the vessel wall (Figure XI, available online at http://atvb.ahajournals.org). Although minimal fatty infiltrates were noted on histological evaluation of the aortic root in C57 mice on high-fat diet, there were no obvious changes in inflammatory cell infiltrate.
Quantitative RT-PCR Validation of Expression Differences
To validate the array results with quantitative RT-PCR and assure that our statistical analyses were identifying truly differentially expressed genes, 10 representative genes were assayed by quantitative RT-PCR (Table XII, available online at http://atvb.ahajournals.org). We attempted to recapitulate the essence of the experimental design of the microarray experiment by using several genes from each group of significant genes. As can be seen in Figure XIII (available online at http://atvb.ahajournals.org), there is high degree of correlation between the 2 methodologies (Pearson correlation of 0.86), validating the results of the microarray analyses.
| Discussion |
|---|
|
|
|---|
There have been previous studies using gene expression analysis, including microarray methodology, to identify genes differentially expressed in atherosclerosis.1922 However, because of limitations of sample availability in humans and experimental design in rodent studies, these efforts have tended to identify those genes expressed at high levels at the terminal phase of the disease. In a broader approach to identification of differentially regulated genes, we have used a temporal experimental design and analytical approach that should identify genes that may be only transiently expressed during the disease process. Our AUC algorithm for data addresses the rigorous analysis of temporal microarray data, using an AUC calculation. The area is a natural measure of differential change over time between 2 comparison groups. It is large if there is a consistent rise or fall over time or a sudden change followed by a leveling off. This measure has been coupled with additional steps to help identify significant values and should have broad utility for a variety of applications to longitudinal array data. Additionally, in these studies, we have used statistically rigorous methods to identify terms and pathways that are reflected by the genes identified to be differentially regulated.
By validating the genes differentially regulated between the C3H and C57 mice in the apoE model of atherosclerosis, it has been possible to highlight those genes and pathways most likely to underlie the inherited disease susceptibility differences between these strains. PPAR
and Pparbp, which bind PPAR-
, RXR, and TR-ß-1, were noted to have increased levels in C3H mice and decreased expression during progression of atherosclerosis in apoE mice. The role of peroxisome proliferatoractivated receptors in atherosclerosis, inflammation, and insulin-resistance is becoming increasingly clear.23 Although Ppar
has been shown to play a protective role in atherosclerosis, the role of Pparbp remains poorly understood. Ptp4a1 was noted to have a higher expression in C3H on high-fat diet and lower expression in apoE-deficient mice. Protein tyrosine phosphatases have been implicated in a variety of biological processes including cell-matrix adhesion and cell migration processes. They can also control the formation and maintenance of the actin cytoskeleton and regulate small GTPase molecular switches.24 Nevertheless, their roles have not been clearly defined in atherosclerosis. Mcd (malonyl-coenzyme A [CoA] decarboxylase), which showed a similar pattern of expression to Ptp4a1, catalyzes the conversion of malonyl-CoA to acetyl-CoA in the fatty acid biosynthesis pathway. A recent study has shown that increased levels of Mcd, which decreases circulating free fatty acids and liver triglcyceride content, is associated with improved insulin sensitivity.25 Hmgcs2 (3-hydroxy-3-methylglutaryl [HMG] CoA synthase) was noted to have lower expression in C3H mice as well. This enzyme condenses acetyl-CoA with acetoacetyl-CoA to form HMG-CoA, which is the substrate for HMG-CoA reductase. HMG-CoA reductase inhibitors are widely used to treat hypercholesterolemia and atherosclerosis. The relationship between HMG-CoA reductase inhibition with statin drugs and inflammation is a subject of intense clinical investigation.
Gene expression differences between the 2 strains suggest that C57Bl/6 mice have a higher genetic propensity to develop inflammation in the presence of appropriate atherogenic stimuli. A number of immune modulatory genes had higher expression in C57 versus C3H mice at baseline and were upregulated throughout the progression of atherosclerosis. Cd47 modulates a variety of cell functions, such as adhesion, spreading, and migration. Interaction of Cd47 with thrombospondin-1 has been shown to perpetuate inflammation in rheumatoid synovitis,26 and it may also contribute to platelet adhesion to tumor necrosis factorstimulated vascular endothelial cells.27 Both major histocompatibility complex classes I and II proteins, which function in antigen presenting complexes, are of particular interest given their crucial role in modulation of inflammation. Our data demonstrate the differential expression of a number of these genes in C57 and C3H mice. H2-D1 (Histocompatibility 2-D1 antigen), H2-Eb1, Bat5, and major histocompatibility complex class II antigen A, beta (H2-Ab) demonstrated higher expression levels in the C57 mice. No prior role in atherosclerosis susceptibility has been reported for these genes. However, an intriguing recent study has implicated a potential role of H2-Ab in autoimmune mediated cardiomyopathy.28 Interestingly, hepatic expression of H2-Eb1 and H2Ab1 is upregulated in mice fed high-fat diet, which may suggest a relationship between lipid metabolism and immunity.29 Among the immune related genes with higher expression in C3H, Ifi202b (interferon-activated gene 202b) was noted to have a higher expression in C3H at baseline (data not shown). Ifi202b is a p202 protein whose overexpression is growth inhibitory and which can bind and inhibit the activity of numerous transcription factors including cJun, cFos, NF
B, E2f1, E2f4, Myod, and myogenin.30
Whether these inflammatory pathways or 1 or more of the other fundamental cellular pathways outlined above are primarily responsible for the inherited differences in atherosclerosis susceptibility remains to be elucidated, but this work serves to focus future studies with these mouse models.
| Acknowledgments |
|---|
Received May 10, 2004; accepted November 11, 2004.
| References |
|---|
|
|
|---|
tabibiazar/cross-strain/index.htm.
This article has been cited by other articles:
![]() |
W. S. SHIN and S. G. ROCKSON Animal Models for the Molecular and Mechanistic Study of Lymphatic Biology and Disease Ann. N.Y. Acad. Sci., May 1, 2008; 1131(1): 50 - 74. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Ashley, J. M. Spin, R. Tabibiazar, and T. Quertermous Frontiers in Nephrology: Genomic Approaches to Understanding the Molecular Basis of Atherosclerosis J. Am. Soc. Nephrol., November 1, 2007; 18(11): 2853 - 2862. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Kim, M. Pham, I. Luttrell, D. D. Bannerman, J. Tupper, J. Thaler, T. R. Hawn, E. W. Raines, and M. W. Schwartz Toll-Like Receptor-4 Mediates Vascular Inflammation and Insulin Resistance in Diet-Induced Obesity Circ. Res., June 8, 2007; 100(11): 1589 - 1596. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Romanenko, T. Nakamoto, A. Srivastava, J. E. Melvin, and T. Begenisich Molecular Identification and Physiological Roles of Parotid Acinar Cell Maxi-K Channels J. Biol. Chem., September 22, 2006; 281(38): 27964 - 27972. [Abstract] [Full Text] [PDF] |
||||
![]() |
A.P.J.J. Bijnens, E. Lutgens, T. Ayoubi, J. Kuiper, A.J. Horrevoets, and M.J.A.P. Daemen Genome-Wide Expression Studies of Atherosclerosis: Critical Issues in Methodology, Analysis, Interpretation of Transcriptomics Data Arterioscler. Thromb. Vasc. Biol., June 1, 2006; 26(6): 1226 - 1235. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Tabibiazar, R. A. Wagner, A. Deng, P. S. Tsao, and T. Quertermous Proteomic profiles of serum inflammatory markers accurately predict atherosclerosis in mice Physiol Genomics, April 13, 2006; 25(2): 194 - 202. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Tabibiazar, R. A. Wagner, E. A. Ashley, J. Y. King, R. Ferrara, J. M. Spin, D. A. Sanan, B. Narasimhan, R. Tibshirani, P. S. Tsao, et al. Signature patterns of gene expression in mouse atherosclerosis and their correlation to human coronary disease Physiol Genomics, July 14, 2005; 22(2): 213 - 226. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Lutgens, B. Faber, K. Schapira, C. T.A. Evelo, R. van Haaften, S. Heeneman, K. B.J.M. Cleutjens, A. P. Bijnens, L. Beckers, J. G. Porter, et al. Gene Profiling in Atherosclerosis Reveals a Key Role for Small Inducible Cytokines: Validation Using a Novel Monocyte Chemoattractant Protein Monoclonal Antibody Circulation, June 28, 2005; 111(25): 3443 - 3452. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
ATVB Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2005 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |