Time course analysis of gene expression identifies multiple genes with differential expression in patients with in-stent restenosis
BMC Medical Genomics volume 4, Article number: 20 (2011)
The vascular disease in-stent restenosis (ISR) is characterized by formation of neointima and adverse inward remodeling of the artery after injury by coronary stent implantation. We hypothesized that the analysis of gene expression in peripheral blood mononuclear cells (PBMCs) would demonstrate differences in transcript expression between individuals who develop ISR and those who do not.
Methods and Results
We determined and investigated PBMC gene expression of 358 patients undergoing an index procedure to treat in de novo coronary artery lesions with bare metallic stents, using a novel time-varying intercept model to optimally assess the time course of gene expression across a time course of blood samples. Validation analyses were conducted in an independent sample of 97 patients with similar time-course blood sampling and gene expression data. We identified 47 probesets with differential expression, of which 36 were validated upon independent replication testing. The genes identified have varied functions, including some related to cellular growth and metabolism, such as the NAB2 and LAMP genes.
In a study of patients undergoing bare metallic stent implantation, we have identified and replicated differential gene expression in peripheral blood mononuclear cells, studied across a time series of blood samples. The genes identified suggest alterations in cellular growth and metabolism pathways, and these results provide the basis for further specific functional hypothesis generation and testing of the mechanisms of ISR.
Cardiovascular disease is the leading cause of death in western countries and a major cause of morbidity and mortality world-wide. Studying coronary atherosclerotic disease (CAD) is challenging for several reasons, since it has substantial environmental and genetic components. Furthermore, despite the nearly universal presence of coronary atherosclerosis, particularly as individuals age, cardiovascular events such as acute coronary syndromes, sudden death and the need for revascularization therapy only occur in some individuals, highlighting the difficulty in precisely defining atherosclerosis phenotypes. In symptomatic patients, revascularization therapy is often required, and percutaneous intervention with balloon angioplasty and stent implantation is a cornerstone of therapy. In-stent restenosis (ISR) is a late complication of stent implantation in which inflammatory and proliferative responses to the vascular injury caused by angioplasty and stenting lead to neointimal hyperplasia within the stent and at its edges over the following weeks and months. Many of the same inflammatory and proliferative processes are activated in the development of atherosclerosis but occur over years or decades. ISR is characterized by proliferative responses to the vascular wound incurred as a result of stent implantation. Therefore, ISR may be viewed as a model phenotype of vascular wound repair for which the mechanisms represent part of the pathologic picture of atherosclerosis, with relatively accelerated wound repair responses operative in the vascular wall and in peripheral blood leukocytes.
In the analysis reported here, we apply a novel method to analyze time-course gene expression data to gene expression profiles of peripheral blood mononuclear cells (PBMCs) of patients enrolled in our study of ISR. The results of the discovery transcriptome analysis of the CardioGene Study were tested for replication in an independent sample of Icelandic patients. We ultimately identified and validated a set of 32 genes that are temporally differentially expressed in the blood of patients who develop ISR after stenting, compared to those who do not develop ISR, highlighting the importance of cellular growth pathways and identifying several biologic candidates for further mechanistic investigation.
Microarray data quality control and filtering
In the CardioGene Study, 312 patients were included (52 with ISR, and 260 who did not develop ISR, as defined by criteria for clinical restenosis) after quality control filtering (Table 1). All 312 patients had baseline gene expression profile data (Figure 1a). Of these, 203 had a blood sample and high quality gene expression profile data at early follow-up (2-4 weeks post-stent), and 166 had high quality gene expression profile data at both early follow-up and late follow-up (6 months post-stent). A total of 681 samples were included in the time course analysis. Box plots showing the distributions of the early and late follow-up times are presented in Figure 1b. From the deCODE sample of patients in Iceland, 97 patients were enrolled and had high quality gene expression profiles at baseline (28 with ISR, 69 without ISR). Of these, 86 patients had a follow-up blood sample and acceptable gene expression data at 6 months post-stent. Thus, a total of 183 samples were included in the replication analysis.
CardioGene discovery analysis of time-course RNA expression in PBMCs
Applying the time varying intercept model, 46 probes, corresponding to 42 distinct genes, were found to be significant at FDR adjusted q-value less than 0.05 out of 12467 probe sets that passed the quality control (Additional File 1). Two patterns of gene expression were observed, with some genes showing divergent gene expression patterns across the time course (Figure 2a) and others showing consistently differential gene expression over the time course (Figure 2b). Descriptions of the genes to which these probe sets map are provided in Additional File 1.
Replication analysis using deCODE time course samples
Among 46 significant probe sets identified in the CardioGene time-course analysis, 36 probe sets had FDR adjusted q-value less than 0.05 in the deCODE replication analysis (Figure 3). Gene expression patterns were largely consistent with the patterns observed in the CardioGene discovery analysis (Figure 2c and 2d). The 36 probe sets mapped to 32 unique genes (Table 2).
Annotation of Genes Identified
Using the DAVID/EASE annotation tool, several categories of genes were identified, falling into multiple gene ontology classifications. No single category achieved a significance score indicating statistically significant over-representation of a specific gene ontology, and several gene ontology classes showed multiple genes (Additional File 2). The categories with the highest representation (i.e. the largest number of genes mapping to these ontologies) included cell growth and/or maintenance, cellular metabolism, catalytic activity, nucleotide and nucleic acid metabolism, transport, protein binding and cellular binding. Genes with prior evidence of expression in the arterial wall after vascular injury included the NAB2 and LAMP2 genes, with data in the prior literature suggestive of functional relevance of these genes, by mediating proliferative responses and recruitment of PBMCs to injured endothelium, respectively[3–13] (Figure 4).
In a sensitivity analysis, we found that the median number of significant probe sets using the same criteria of FDR 0.05 was only 2, the mean was 3.04, and the maximum number was 14 in one iteration, demonstrating a significantly low likelihood that the probesets we identified and replicated were identified due to chance.
Hypothesis-driven analysis of gene expression
We compared our results to published studies examining gene expression in ISR, chiefly in human vascular tissue samples [14, 15]. In one study, peripheral blood total leukocytes were also examined . In the CardioGene baseline blood sample, genes reported by Ashley et al were examined , and none were significantly differentially expressed in the baseline blood sample or in our timecourse results. In another study by Zohlnhofer et al , gene expression data from neointimal vascular tissue and peripheral blood leukocytes were reported. We examined those genes meeting a Bonferroni multiple testing correction from this discovery report and did not find overlap with our gene expression findings.
Using a repeated sampling study design, we investigated PBMC transcript expression in patients undergoing stent implantation, using a novel time-course analysis method. We identified a set of 42 genes with differential temporal expression among patients with and without ISR at one year follow-up in a discovery analysis of the CardioGene Study. Independent replication testing in an Icelandic sample confirmed differential expression of 36 probesets mapped to 32 genes. The gene expression patterns over time may be of interest as well, with consistently expressed genes representing gene expression data that may be able to predict ISR and differentially expressed genes over the time course representing genes with possible direct functional roles in the development of ISR, both of which require further investigation to explore more fully.
Gene expression profiling with DNA microarray technology is a popular tool to monitor the expression level of thousands of genes simultaneously and has been applied in cardiovascular research, to detect patterns of gene expression indicative of underlying disease states[16–21]. Since the data generated represent the temporal abundance of mRNA levels in the sample, measurements of the change of this abundance over the course of disease progression (or any biological process) is therefore both possible and of great scientific interest using this technology. In fact, Yuan and Kendziorski reported that more than one-third of the experiments catalogued in the Gene Expression Omnibus (National Center for Biotechnology Information; http://www.ncbi.nlm.nih.gov/geo/) are from experiments that measure gene expression over time. Early time-course RNA expression studies have focused on identifying clusters of genes with a similar pattern over time[14, 23, 24]. More recently, detecting a differential gene expression pattern over time between several biological groups has become an interesting goal of time-course gene expression data.
To detect differential gene expression pattern over time, we used a time-varying intercept model which can account for differences in sample intervals between patients. The term "differential gene expression pattern over time" can be interpreted in several ways to form suitable questions and, thus, the hypotheses are dependent on the particular experiment under consideration. We considered two related questions where each can roughly correspond to the main effect of the group, and the interaction between group and time points. Consider, for example, a gene with exactly the same expression pattern over time in both groups, but the first group has a higher expression than the second group, consistently over all time points. This is a gene that shows the main effect of the group. On the other hand, consider another gene with similar expression level at two time points, 1 and 2, but this gene's expression increases from time 1 to 2 in one group but decreases in another group. This group-by-time interaction cannot be examined by methods that only test the main effect of the group. The time-varying intercept model we used can detect both the main effect and group-by-time interaction. This method, however, requires a large number of bootstrap resampling to evaluate the significance level of the difference between groups, which can be computationally challenging especially when a large number of genes are tested.
Of the genes we identified, the most extensive prior literature in vascular disease was found for the NAB2 gene, which is also known as EGR1 binding protein 2. Early growth response (EGR) genes, which are transcription factors that are implicated in a wide variety of proliferative and differentiative processes. Nab proteins are necessary for Schwann cells to exit the cell cycle and generate a myelin-specific gene profile and are key regulators of the myelination process of peripheral nerves. NAB2 is expressed in vascular smooth muscle cells in response to injury[3–12, 26] and EGR1 has been identified in a microarray study of in vitro smooth muscle cell proliferation. The LAMP2 gene product protects the lysosomal membrane from proteolytic enzymes within lysosomes and also functions as a receptor for proteins to be imported into lysosomes. Mutations in the LAMP2 gene have been identified in patients with hypertrophic cardiomyopathy, and the gene product mediates adherence of PBMCs to the vascular endothelium. Cellular adhesion to the vascular endothelium has been well-described in animal models and post-mortem human examinations, in atherogenesis and acute vascular injury[30, 31]. In the latter, the extent of leukocyte adhesion is predictive of the degree of subsequent neointimal hyperplasia, which is the key lesion of ISR.
Some genes identified have no apparent role in vascular biology, such as VPS26, VPS41, SRP54, and RAD23B, and comparison to previously published reports in the literature do not show differential gene expression in other studies of restenosis, although these investigations were conducted primarily on vascular tissues or culture vascular smooth muscle cells rather than peripheral blood[17, 18, 27, 32]. The VPS26 and VPS4 genes belongs to a group of vacuolar protein sorting genes and may have a role in lysosome maintenance, and SRP54 is a protein in the signal recognition particle, which directs secretory proteins to membranes as they emerge from the ribosome . Specific vascular or inflammatory cell function has not been described, but derangement of these basic cellular processes may impact vascular and other physiologic functions adversely. RAD23B has a role in DNA (nucleotide excision) repair, and genetic variants in RAD23B have been associated with several solid tumors[35–38]. The association with cancer would suggest a possible link to excess proliferative mechanisms in vascular wound repair, as has been described for many other cell cycle regulatory genes. Genetic variants in other genes we identified are also associated with human diseases. Several genetic variants in ACADM have been associated with medium-chain acyl-CoA dehydrogenase activity, but there is no known vascular implication of this disorder. Variants in PCMT1 and FOLR2 have been associated with neural tube defects[41–43], with no known vascular phenotype in these cases. Overall the findings of this study are hypothesis-generating and can be used to support the rational for investigating the function of specific genes and pathways in adverse vascular remodeling, which is relevant to both ISR and more general CAD phenotypes.
We compared the results of our study to previously published reports of transcriptome analysis in ISR. Our results were negative for replication of these studies which focused primarily on vascular tissue samples, in relatively small sample sets. In one study, peripheral blood total leukocyte gene expression was studied in 10 patients with ISR and atherectomy specimens . While a high degree of correlation between peripheral blood leukocyte and arterial neointima tissue gene expression was identified in a subset of genes, these findings were based upon single measurements in a small sample size and were not replicated in the original report. These prior reports highlight the major difficulty of studying vascular tissues, since access to these tissues in adequately large sample sizes is limited.
Vascular biopsy and atherectomy are performed infrequently as part of routine clinical care and would not support well-powered studies of vascular tissues. Tissue sampling over a time course is not clinically indicated or possible. Additionally, a large degree of intra-individual variability in gene expression was noted in these prior studies of vascular tissues, making replication testing critical, yet this cannot be done without access to additional tissues samples. In our study, we analyzed peripheral blood leukocytes, specifically focusing on the mononuclear fraction which contains primarily B and T lymphocytes and monocytes. Although the analysis of peripheral blood cells would ideally be complemented by similar studies in vascular tissues, studying gene expression profiles in blood leukocytes is biologically relevant due to well-defined interactions with the arterial wall, particularly in the setting of vascular injury and repair as in the setting of ISR. The overlap between vascular and blood gene expression in one prior transcriptome analysis of ISR was supportive of our rationale to study PBMCs. For these reasons, and with the additional prior knowledge that inflammation plays a significant role in the development of ISR, we undertook a study of PBMCs in several hundred patients, with adequately powered replication testing for our top discovery findings. Additionally, we use a time course analysis method that improved our ability to detect gene expression signals between the two comparison groups, overcoming some of the difficulty of substantial variability in single point microarray gene expression data.
To address the possibility of false positives identified with our statistical methods, we conducted replication analyses in the independent sample of deCODE samples and we conducted bootstrap resampling to assess significance of the findings. Through this sensitivity analysis, we demonstrated that the validation of 36 probe sets is not likely to be due to chance. Additional potential limitations of this study of ISR are the use of a clinical restenosis outcome, rather than an angiographic outcome, in which clinically silent ISR may have been missed, and the choice of tissue analyzed, as discussed. The CardioGene and deCODE cohorts differ in the incidence of ISR (16.7% in the CardioGene Study and 28.8% in deCODE) with the patients in the deCODE sample showing overall lower residual percent stenosis in the treated lesion after stent implantation. Also, the proportions with hyperlipidemia and diabetes differ. However, despite the differences in the cohorts, we find replication of a substantial proportion of the discovery findings.
In summary, we have used a method to analyze gene expression in serial blood samples and identified a set of genes that show differential expression in the blood of patients who develop ISR after BMS implantation, compared to those who do not. These gene expression patterns of adverse vascular remodeling suggest possible hypotheses for the mechanisms of injury-induced remodeling observed in both ISR and CAD, since ISR is a niche phenotype occurring in a subset of patients with CAD. Further studies are needed to investigate the functional relevance of these genes and are warranted based upon the findings of this study.
The CardioGene Study was an IRB-approved, prospective cohort study of 358 patients enrolled at the time of bare metal stent (BMS) implantation to treat de novo, previously untreated native coronary artery lesions at William Beaumont Hospital (Royal Oak, Michigan, USA) and the Mayo Clinic (Rochester, Minnesota, USA). Patients were followed for one year to determine ISR outcomes (Figure 1). Enrollment began in February 2002 and was closed in September 2003, prior to the approval and clinical use of drug-eluting stents (DES) in the United States. The protocol was approved by the NHLBI IRB as well as the IRB at each of the clinical enrollment sites. Informed consent was provided by each patient. Standardized case report forms were used to collect baseline clinical data and outcome information in follow-up. Since gene expression discovery studies may suffer from false positives, despite the use of statistical corrections for multiple hypothesis testing, we sought out independent validation of the primary findings of the CardioGene Study in a separate cohort. For independent replication testing, 97 patients undergoing stent implantation with BMS were enrolled in Iceland by collaborators at Landispitali University Hospital and University of Iceland, and RNA samples were provided through a contract with deCODE Genetics, Inc. The study was approved by the National Bioethics Committee and the Data Protection Authority of Iceland, and each patient provided informed consent for participation in the study. Clinical characteristics of both study cohorts are summarized in Table 1.
Consecutive patients presenting to the cardiac catheterization laboratories of the clinical enrollment sites were approached for participation in the study. Follow-up clinical evaluation was performed via patient interview and review of all available medical records at 6 months and 12 months post-stent (Figure 1). ISR was defined as clinical restenosis, which was defined by ischemic symptoms after stent implantation and evidence of flow limitation in the treated vessel by either invasive or non-invasive testing. Follow-up angiography was not specifically performed for the CardioGene Study. Any available angiographic data performed as part of each patient's clinical care was recorded. The Icelandic subset had follow-up angiography of all patients at 6 months. Quantitative coronary angiography (QCA) in the CardioGene Study was performed by a single QCA reader, and by a TIMI core QCA center for the Icelandic subset.
PBMC Microarray Analysis
Whole blood was sampled immediately prior to stenting, 2-4 weeks after stenting and six months post-stent. Blood samples were collected into EDTA-containing tubes at three times from patients enrolled in the CardioGene Study (Figure 1), and two time from patients enrolled in Iceland. Samples from each time point were handled according to the standardized blood handling and RNA isolation protocols, for consistency across time points and clinical enrollment sites. PBMCs were harvested by technicians at each clinical enrollment site following standardized protocols for Ficoll separation of whole blood and osmostic lysis of red blood cells, to avoid technical variation and to minimize globin mRNA contamination. All blood processing was completed within 4 hours of blood draw, at room temperature. Snap-frozen PBMCs samples were shipped to the NIH laboratory, and RNA isolation was performed using the Qiagen RNA isolation protocol, with DNase treatment. RNA labeling and microarray image analysis was according to the Affymetrix U133A protocol. To minimize the impact of batch-to-batch variation, the multiple samples from each patient were processed together through all stages, and the order of RNA labeling and microarray hybridization was randomized. The same methods were used in the CardioGene and deCODE study samples, and microarray analysis was conducted on all samples at the NIH. In both studies, technicians performing the assays were blinded to the clinical status of patients.
Raw image data were normalized using RMA (Robust multichip average) normalization. Several steps of quality control were performed. These included the investigation of the residual plots using affyPLM, gender match using a linear discriminant analysis and shrunken centroid method[45–47]. Probes with values less than 6 in log 2 scale in more than two thirds of samples were excluded, and age and gender were adjusted before the analysis.
In the study design of the time course of blood gene expression profiling, the follow-up times were set at 2 weeks and 6 months for the early and late follow-up. A window was set at each of the time points, with the early follow-up designated as any time 2-4 weeks post-stent and the late follow-up designated as any time between 5-7 months for patients seen at 6 months. If patients did not return for follow-up at the 6 month time point, an attempt was made to have the patient provide the late follow-up blood sample at 12 months post-stent, at the time when final clinical ascertainment for ISR was made. As a result of the use of time windows, we were able to increase ascertainment of follow-up blood samples, but the actual sample collection intervals after stent implantation were not precisely spaced. To analyze such data, we used the time varying intercept model using a B-spline basis of dimension 2 based on only one knot at the median follow-up time (which was 14 days) (Additional File 3). This is because the majority of the early follow-up is centered around 14 days. The model was fit under the null hypothesis using combined samples of ISR and non-ISR patients, and under the alternative hypothesis where separate models were fitted for ISR and non-ISR patients, respectively. The test statistic is the improvement in the fitted model under the alternative hypothesis from that under the null hypothesis, as measured by the relative difference in the residual sum of squares. The bootstrap method was used to estimate the p-value of this test statistic[47, 49]. To obtain accuracy in estimating the p-values, 10 million bootstrap samples were generated. Significant probe sets were selected based on an FDR-adjusted q-value less than 0.05.
The replication analysis was performed in an independent set of samples, using deCODE samples from patients enrolled in Iceland, to validate probe sets identified by the time course analysis of the CardioGene samples. We performed the same analysis using the time varying intercept model. P-values were estimated based on 100,000 bootstrap samples, and FDR adjusted q-value of 0.05 was, again, used to select significant probe sets. Here the adjustment was based on the selected set of genes since it was a replication analysis. As a sensitivity analysis and to evaluate whether any random set of 46 probe sets would result in significant results as we saw in the analysis of deCODE samples, we randomly selected 46 probe sets, and repeated the validation analysis. This procedure was repeated 100 times, and the number of probe sets yielding an FDR value less than 0.05 was observed.
Data analyses were performed based on the R language utilizing libraries freely available via the Bioconductor project.
Box plots were generated using the default R software options. The upper whisker is the data point which is 1.5 times the difference between 75th percentile and median from the upper bound of the box, and lower whisker is 1.5 times the difference between median and 25th percentile from the lower bound of the box. Overall, the whisker represents the data point which is 1.5 times the interquartile range from the box.
Finally, to evaluate genes identified in a prior studies of ISR [23, 24], we analyzed the reported genes in our data, at the baseline timepoint comparing normalized gene expression values between ISR cases and no-ISR controls and also against our timecourse analysis results.
Annotation of probe sets was done using the NetAffx (Affymetrix U133A library set) and the DAVID/EASE software, which provides gene name and functional annotation of probe sets[51, 52]. Analysis of over-representation, using a Fisher's exact test of over-representation of genes from different functional categories based on gene ontology classifications, was performed.
Conflict of Interest Disclosures
The authors declare that they have no competing interests.
bare metal stents
coronary atherosclerotic disease
peripheral blood mononuclear cells
Tanner FC, et al: Expression of cyclin-dependent kinase inhibitors in vascular disease. Circ Res. 1998, 82 (3): 396-403.
Ganesh SK, et al: Rationale and study design of the CardioGene Study: genomics of in-stent restenosis. Pharmacogenomics. 2004, 5 (7): 952-1004. 10.1517/14622422.214.171.1249.
Santiago FS, Atkins DG, Khachigian LM: Vascular smooth muscle cell proliferation and regrowth after mechanical injury in vitro are Egr-1/NGFI-A-dependent. Am J Pathol. 1999, 155 (3): 897-905. 10.1016/S0002-9440(10)65189-9.
Silverman ES, et al: Vascular smooth muscle cells express the transcriptional corepressor NAB2 in response to injury. Am J Pathol. 1999, 155 (4): 1311-7. 10.1016/S0002-9440(10)65233-9.
Goetze S, et al: TNFalpha induces expression of transcription factors c-fos, Egr-1, and Ets-1 in vascular lesions through extracellular signal-regulated kinases 1/2. Atherosclerosis. 2001, 159 (1): 93-101. 10.1016/S0021-9150(01)00497-X.
Hasan RN, Schafer AI: Hemin upregulates Egr-1 expression in vascular smooth muscle cells via reactive oxygen species ERK-1/2-Elk-1 and NF-kappaB. Circ Res. 2008, 102 (1): 42-50. 10.1161/CIRCRESAHA.107.155143.
Morawietz H, et al: Rapid induction and translocation of Egr-1 in response to mechanical strain in vascular smooth muscle cells. Circ Res. 1999, 84 (6): 678-87.
Ling S, et al: Matrix-dependent gene expression of egr-1 and PDGF A regulate angiotensin II-induced proliferation in human vascular smooth muscle cells. Hypertension. 1999, 34 (5): 1141-6.
Sachinidis A, et al: Thromboxane A2 and vascular smooth muscle cell proliferation. Hypertension. 1995, 26 (5): 771-80.
Grote K, et al: Stretch-inducible expression of the angiogenic factor CCN1 in vascular smooth muscle cells is mediated by Egr-1. J Biol Chem. 2004, 279 (53): 55675-81. 10.1074/jbc.M406532200.
McCaffrey TA, et al: High-level expression of Egr-1 and Egr-1-inducible genes in mouse and human atherosclerosis. J Clin Invest. 2000, 105 (5): 653-62. 10.1172/JCI8592.
Santiago FS, et al: New DNA enzyme targeting Egr-1 mRNA inhibits vascular smooth muscle proliferation and regrowth after injury. Nat Med. 1999, 5 (11): 1264-9. 10.1038/15215.
Kannan K, et al: Lysosome-associated membrane proteins h-LAMP1 (CD107a) and h-LAMP2 (CD107b) are activation-dependent cell surface glycoproteins in human peripheral blood mononuclear cells which mediate cell adhesion to vascular endothelium. Cell Immunol. 1996, 171 (1): 10-9. 10.1006/cimm.1996.0167.
Spellman PT, et al: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell. 1998, 9 (12): 3273-97.
Svaren J, et al: NAB2, a corepressor of NGFI-A (Egr-1) and Krox20, is induced by proliferative and differentiative stimuli. Mol Cell Biol. 1996, 16 (7): 3545-53.
Zohlnhofer D, et al: Gene expression profiling of human stent-induced neointima by cDNA array analysis of microscopic specimens retrieved by helix cutter atherectomy: Detection of FK506-binding protein 12 upregulation. Circulation. 2001, 103 (10): 1396-402.
Zohlnhofer D, et al: Transcriptome analysis reveals a role of interferon-gamma in human neointima formation. Mol Cell. 2001, 7 (5): 1059-69. 10.1016/S1097-2765(01)00239-8.
Ashley EA, et al: Network analysis of human in-stent restenosis. Circulation. 2006, 114 (24): 2644-54. 10.1161/CIRCULATIONAHA.106.637025.
Blaschke F, et al: C-reactive protein induces apoptosis in human coronary vascular smooth muscle cells. Circulation. 2004, 110 (5): 579-87. 10.1161/01.CIR.0000136999.77584.A2.
Fleiner M, et al: Arterial neovascularization and inflammation in vulnerable patients: early and late signs of symptomatic atherosclerosis. Circulation. 2004, 110 (18): 2843-50. 10.1161/01.CIR.0000146787.16297.E8.
Heymans S, et al: Increased cardiac expression of tissue inhibitor of metalloproteinase-1 and tissue inhibitor of metalloproteinase-2 is related to cardiac fibrosis and dysfunction in the chronic pressure-overloaded human heart. Circulation. 2005, 112 (8): 1136-44. 10.1161/CIRCULATIONAHA.104.516963.
Yuan M, Kendziorski C: A unified approach for simultaneous gene clustering and differential expression identification. Biometrics. 2006, 62 (4): 1089-98. 10.1111/j.1541-0420.2006.00611.x.
Chu S, et al: The transcriptional program of sporulation in budding yeast. Science. 1998, 282 (5389): 699-705. 10.1126/science.282.5389.699.
DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278 (5338): 680-6. 10.1126/science.278.5338.680.
Le N, et al: Nab proteins are essential for peripheral nervous system myelination. Nat Neurosci. 2005, 8 (7): 932-40.
Lucerna M, et al: NAB2, a corepressor of EGR-1, inhibits vascular endothelial growth factor-mediated gene induction and angiogenic responses of endothelial cells. J Biol Chem. 2003, 278 (13): 11433-40. 10.1074/jbc.M204937200.
Nguyen KT, et al: Molecular responses of vascular smooth muscle cells and phagocytes to curcumin-eluting bioresorbable stent materials. Biomaterials. 2004, 25 (23): 5333-46. 10.1016/j.biomaterials.2003.12.033.
Fukuda M, et al: Cloning of cDNAs encoding human lysosomal membrane glycoproteins, h-lamp-1 and h-lamp-2. Comparison of their deduced amino acid sequences. J Biol Chem. 1988, 263 (35): 18920-8.
Charron P, et al: Danon's disease as a cause of hypertrophic cardiomyopathy: a systematic survey. Heart. 2004, 90 (8): 842-6. 10.1136/hrt.2003.029504.
Libby P: Inflammation in atherosclerosis. Nature. 2002, 420 (6917): 868-74. 10.1038/nature01323.
Ross R: Atherosclerosis--an inflammatory disease. N Engl J Med. 1999, 340 (2): 115-26. 10.1056/NEJM199901143400207.
Nguyen KT, et al: Molecular responses of vascular smooth muscle cells to paclitaxel-eluting bioresorbable stent materials. J Biomed Mater Res A. 2004, 69 (3): 513-24. 10.1002/jbm.a.30020.
Haft CR, et al: Human orthologs of yeast vacuolar protein sorting proteins Vps26, 29, and 35: assembly into multimeric complexes. Mol Biol Cell. 2000, 11 (12): 4105-16.
Lutcke H, et al: The methionine-rich domain of the 54 kDa subunit of signal recognition particle is sufficient for the interaction with signal sequences. Embo J. 1992, 11 (4): 1543-51.
Garcia-Closas M, et al: Genetic variation in the nucleotide excision repair pathway and bladder cancer risk. Cancer Epidemiol Biomarkers Prev. 2006, 15 (3): 536-42. 10.1158/1055-9965.EPI-05-0749.
Mechanic LE, et al: Polymorphisms in nucleotide excision repair genes, smoking and breast cancer in African Americans and whites: a population-based case-control study. Carcinogenesis. 2006, 27 (7): 1377-85. 10.1093/carcin/bgi330.
Huang WY, et al: Nucleotide excision repair gene polymorphisms and risk of advanced colorectal adenoma: XPC polymorphisms modify smoking-related risk. Cancer Epidemiol Biomarkers Prev. 2006, 15 (2): 306-11. 10.1158/1055-9965.EPI-05-0751.
Shen M, et al: Polymorphisms in the DNA nucleotide excision repair genes and lung cancer risk in Xuan Wei, China. Int J Cancer. 2005, 116 (5): 768-73. 10.1002/ijc.21117.
Boehm M, Nabel EG: The cell cycle and cardiovascular diseases. Prog Cell Cycle Res. 2003, 5: 19-30.
Kelly DP, et al: Molecular characterization of inherited medium-chain acyl-CoA dehydrogenase deficiency. Proc Natl Acad Sci USA. 1990, 87 (23): 9236-40. 10.1073/pnas.87.23.9236.
Zhu H, et al: A known functional polymorphism (Ile120Val) of the human PCMT1 gene and risk of spina bifida. Mol Genet Metab. 2006, 87 (1): 66-70. 10.1016/j.ymgme.2005.09.008.
Boyles AL, et al: Neural tube defects and folate pathway genes: family-based association tests of gene-gene and gene-environment interactions. Environ Health Perspect. 2006, 114 (10): 1547-52. 10.1289/ehp.9166.
O'Leary VB, et al: Analysis of the human folate receptor beta gene for an association with neural tube defects. Mol Genet Metab. 2003, 79 (2): 129-33.
Irizarry RA, et al: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-64. 10.1093/biostatistics/4.2.249.
Tibshirani R, et al: Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci USA. 2002, 99 (10): 6567-72. 10.1073/pnas.082099299.
Gentleman RC, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Bolstad BM, et al: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-93. 10.1093/bioinformatics/19.2.185.
Dave SS, et al: Prediction of survival in follicular lymphoma based on molecular features of tumor-infiltrating immune cells. N Engl J Med. 2004, 351 (21): 2159-69. 10.1056/NEJMoa041869.
Huang JZ, W C, Zhou L: Varying-coefficient models and basis function approximations for the analysis of repeated measurements. Biometrika. 2002, 89 (1): 111-128. 10.1093/biomet/89.1.111.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100 (16): 9440-5. 10.1073/pnas.1530509100.
Hosack DA, et al: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Dennis G, et al: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/4/20/prepub
This work was supported by the Division of Intramural Research of the National Heart, Lung and Blood Institute, National Institutes of Health.
Study design (EGN, JJ, KAS, NLG, SKG, WON), CardioGene clinical enrollment and blood and data collection (DRH, KAS, KON, LM, RDS, SKG, WON), deCODE clinical enrollment and data collection (AH, KA, GT, TG), RNA expression data generation (EMB, KON, SKG), Data analysis (JJ, SKG), Manuscript preparation (AH, DRH, EGN, JJ, SKG, WON).
All authors read and approved the final manuscript.
Santhi K Ganesh, Jungnam Joo, Kimberly Skelding contributed equally to this work.
Electronic supplementary material
Additional file 1: Supplementary Table 1. Initial 46 probesets identified in the timecourse analysis of the CardioGene samples. (XLS 38 KB)
Additional file 2: Supplementary Table 2. Gene ontology annotation of the final 32 genes identified and validated. (XLS 30 KB)
Additional file 3: Supporting Materials. Statistical methods and lists of probe IDs discovered and replicated. (DOCX 42 KB)
About this article
Cite this article
Ganesh, S.K., Joo, J., Skelding, K. et al. Time course analysis of gene expression identifies multiple genes with differential expression in patients with in-stent restenosis. BMC Med Genomics 4, 20 (2011). https://doi.org/10.1186/1755-8794-4-20