|
|
||||||||
National Heart, Lung, and Blood Institute Program in Genomic Applications-HopGene, Philadelphia, Pennsylvania 19104
1 Department of Medicine, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104
2 Department of Medicine, Johns Hopkins Medical Institutions, Baltimore, Maryland 21287
3 Department of Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland 21287
4 Center for Genetic Medicine, Childrens National Medical Center, Washington, District of Columbia 20010
| ABSTRACT |
|---|
|
|
|---|
-skeletal actin and brain natriuretic peptide). Cluster analyses demonstrated that NOS1 deletion caused more pronounced changes in the myocardial transcriptome than did NOS3 deletion, despite similar cardiac phenotypes. These findings suggest that the transcriptional basis for LVH varies depending on the inciting biochemical stimulus. In addition, NOS isoforms appear to play distinct roles in modulating cardiac structure. ventricular remodeling; mice; knockout; gene expression profiling; genomics
| INTRODUCTION |
|---|
|
|
|---|
-skeletal actin, and brain natriuretic peptide (BNP)] is now so widely accepted that it is often included as part of the definition of cardiac hypertrophy. We and others have demonstrated that physiological cardiac nitric oxide (NO) signaling opposes LVH (3, 5, 35). In mice, deletion of nitric oxide synthases (NOS) leads to age-associated LVH as manifested by an increase in both ventricular wall thickness and mass. Notably, deletion of NOS3 (endothelial NOS) causes hypertrophy associated with systemic hypertension, whereas deletion of NOS1 (neuronal NOS) causes LVH with normal systemic blood pressure (3, 11). These findings suggest that NOS1 deficiency induces LVH via load-independent signals and, therefore, may activate a distinct transcriptional pathway.
We hypothesized that deletion of NOS3 or NOS1 would activate different cardiac hypertrophy transcriptional programs. Because a large number of genes have been implicated in LVH pathogenesis, we employed a genomic approach using microarrays and quantitative polymerase chain reaction (qPCR) to analyze multiple genes simultaneously. We tested our hypothesis using three approaches: 1) exploratory analysis of differentially expressed genes in NOS-/- mice, 2) candidate-gene analysis using sets of genes implicated in LVH pathogenesis, and 3) cluster analysis comparing global patterns of gene expression. Our findings demonstrate that hypertrophied NOS3-/- and NOS1-/- hearts have distinct molecular phenotypes, suggesting unique roles of NOS isoforms in the pathogenesis of LVH.
| METHODS |
|---|
|
|
|---|
Microarray hybridization.
Myocardial RNA was isolated from individual NOS3-/-, NOS1-/-, and WT mice (n = 3 in each group; age 20 ± 3 mo) using the Trizol reagent and Qiagen RNeasy columns. Seven of the mice were male, and two were female (one NOS3-/- and one NOS1-/-). Individual cDNAs were prepared from each RNA isolate using reverse transcriptase (GIBCO-BRL SuperScript). Each cDNA was subsequently used as a template to make biotin-labeled cRNA using an in vitro transcription reaction (Enzo), resulting in a single cRNA for each heart. Each cRNA was hybridized with an individual Affymetrix MG-U74A oligonucleotide array, which was subsequently processed and scanned according to the manufacturers instructions. All arrays were hybridized in the same batch to avoid variability in hybridization conditions. Each array quantifies the expression of 12,422 transcripts (including full-length mRNA sequences and expressed sequence tag clusters) derived from build 74 of the UniGene database (http://www.affymetrix.com). Data were saved as raw image files and converted into probe set data (as "* .cel" files) using Microarray Suite (MAS 5.0). All data are available for download at http://www.hopkins-genomics.org.
Analysis of microarray data.
We used robust multi-array average (RMA) to analyze Affymetrix probe set data (13). Software for RMA is available for download (http://www.bioconductor.org) for use in the R-package for statistical computing (http://www.r-project.org). There are four stages to RMA. First, probe set data ("* .cel" files) from all arrays are simultaneously normalized using quantile normalization, which eliminates systematic differences between GeneChips, without significantly altering the relative intensity of probes within a GeneChip. Second, mean optical background level for each array is estimated, and the intensity for each probe is adjusted to remove this. Third, the normalized, background-corrected data is transformed to the log2 scale. Finally, a median-polish procedure is used to combine multiple probes into a single measure of expression for each gene on each array. Three types of analysis were subsequently performed on this normalized expression data as described below.
Differentially expressed genes in NOS-/- hearts compared with WT controls.
To determine regulated genes in NOS3-/- and NOS1-/- myocardium, we applied the following criteria. First, data were filtered to include genes present above background on at least one array according to MAS 5.0. Next, we required at least a 1.8-fold change in expression between WT and NOS-/- in order for a gene to be considered of biological interest. Finally, changes in gene expression were filtered for statistical significance using the following test statistic
![]() |
. Assuming that
follows a t-distribution, this threshold corresponds to a two-sided P < 0.0009 in our data. The identity of each differentially expressed gene was determined using annotation databases (http://www.affymetrix.com/index.affx and http://unchip.org:8080/bio/) or via BLAST searches of the corresponding expressed sequence tag.
Candidate-gene analysis.
We assembled a list of 71 candidate genes previously implicated in the pathogenesis of LVH (see the Supplemental Tables A and B, available as Supplementary Material at the Physiological Genomics web site)1
and identified the corresponding probe sets on the MG-U74A array. In effect, this collection of probe sets defines a hypertrophy-specific subset. Expression levels for each of the candidates were obtained from the normalized expression profiles of NOS3-/-, NOS1-/-, and WT mice (n = 3 for each genotype). For each candidate gene, two-sided unpaired Students t-tests were performed comparing expression levels between NOS-/- mice and their WT controls, and genes that showed a change compared with control at the P < 0.05 level were selected. Expression of these genes in each of the NOS-/- samples was normalized to the mean expression level in WT mice, and the results were displayed in an Eisen plot (7), in which red indicates increased expression and green indicates decreased expression compared with WT.
Cluster analysis.
We compared global patterns of gene expression among the mouse strains using hierarchical clustering. All arrays (NOS3-/-, NOS1-/-, WT) were normalized and converted to measures of absolute gene expression using RMA. Data were filtered to include genes that were present above background on at least one array according to MAS 5.0 and that had a max/min expression ratio across all arrays of at least 1.8. Clusters were constructed using Cluster software (http://rana.lbl.gov) using average linkage clustering and three different measures of similarity: Pearson correlation coefficients, Spearman rank correlation coefficients, and Kendalls tau (7, 8).
Validation.
qPCR was performed to validate changes in seven selected genes using the same myocardial RNA samples analyzed with microarrays. Each RNA sample was treated with DNase I to remove any contaminating genomic DNA and was subsequently used to synthesize cDNA using an oligo(dT)1218 primer and a first-strand synthesis kit (Invitrogen). Primers for qPCR were designed using Primer Express 2.0 software. Each sample was run in duplicate (6 PCR assays for each mouse strain) on a GeneAmp 7900 Sequence Detection System (PE Applied Biosystems) and analyzed using SDS 2.0 software (Applied Biosystems). For each gene of interest, a standard curve was generated using serial dilutions of WT cDNA. The quantity of gene transcript in unknown samples was estimated using this standard curve, with GAPDH as a normalizer. SYBR green reagent (Applied Biosystems) served as a reporter throughout all experiments. Levels of transcript normalized to GAPDH were compared between experimental and control samples using a Wilcoxon rank-sum test, with P < 0.05 as the criterion for statistical significance.
| RESULTS |
|---|
|
|
|---|
Differentially expressed genes in NOS-/- hearts.
Of the 12,422 genes and expressed sequence tags present on the microarrays, 6,409 were detected at levels higher than background in at least one of the strains (NOS3-/-, NOS1-/- or WT). As shown in Fig. 1, 47 genes were differentially expressed in NOS3-/- hearts and 65 were differentially expressed in NOS1-/- hearts compared with WT controls. A complete annotated list of these genes is available in the Supplemental Tables A and B. We categorized each of the differentially expressed genes by known or proposed gene function (Table 1). These genes participate in cellular responses that are universally important for cell proliferation (growth and differentiation, transcription, protein turnover, apoptosis), as well as pathways that may be related specifically to cardiac hypertrophy (structure, metabolism, stress response, redox control). A number of genes related to immunoglobulin synthesis were also differentially expressed, perhaps reflecting different degrees of immune activation in NOS3-/- and NOS1-/- mice compared with WT. Alternatively, they could result from different antigenic exposures in the three mouse strains. Remarkably, the direction of change in expression (increase or decrease) was similar for only 16 genes (Table 2).
|
|
|
We performed a literature search to identify candidate hypertrophy genes, 71 of which were present on the MG-U74A microarray. A complete list of candidates is available in the Supplemental Tables A and B. Of these candidates, 32 (52%) were changed in at least one of the experimental strains compared with WT. These results are displayed in Fig. 2. Unlike classic models of load-induced hypertrophy (24), NOS-/- hearts did not show an isoform switch from
- to ß-MHC or upregulation of
-skeletal actin. (These genes did not meet criteria for differential expression and are not displayed in Fig. 2.) Likewise, there was no upregulation of natriuretic peptides in NOS-/- hearts. In fact,
-skeletal actin and BNP were paradoxically underexpressed in the NOS1-/- strain despite the presence of substantial LVH. We also found a potent downregulation of myocyte-enriched calcineurin interacting protein (MCIP-1), an inhibitor of calcineurin-mediated hypertrophy (25), in the NOS1-/- strain. Transcripts for the Ca2+ handling proteins (L-type Ca2+ channel, SERCA2a, ryanodine receptor) were downregulated to a greater extent in NOS3-/- than in NOS1-/-.
|
|
-skeletal actin in NOS1-/-. Neither qPCR nor microarray detected significant changes in
- or ß-MHC in either strain. In two cases (
-skeletal actin and MCIP-1), the microarrays failed to detect changes that were detected by the more sensitive qPCR technique in the NOS3-/- hearts. There were no disagreements among microarray and qPCR findings in the direction of change in expression compared with control.
|
|
| DISCUSSION |
|---|
|
|
|---|
Some of our findings are counterintuitive, such as the absence of changes in myosin isoforms and paradoxical decreases in BNP and
-skeletal actin. However, they add to a growing literature revealing the diversity of transcriptional changes associated with hypertrophy. Aronow et al. (2) compared expression profiles from mice with hypertrophy induced by four distinct transgenic alterations and found that no single hypertrophy program was responsible. A number of functional and transgenic studies suggest that calcium/calcineurin signaling pathways play a central role in mediating hypertrophy (9, 19). However, pharmacological inhibition of calcineurin as a strategy for preventing LVH has had mixed results (31, 36), indicating that calcineurin signaling may not mediate all forms of hypertrophy. Studies in humans indicate that increased fetal isoform expression, which is regarded as a defining characteristic of hypertrophy in animal models, may not always accompany human heart failure. Rather, fetal isoform expression may remain unchanged, while adult isoform expression decreases (23).
The divergent transcriptional responses in NOS3-/- and NOS1-/- reflect the functional specificity of NOS isoforms in the cardiovascular system (3). NOS1 is localized to the sarcoplasmic reticulum (SR) and modulates the function of the SR calcium release channel (ryanodine receptor). We have shown that deletion of NOS1 impairs basal SR calcium cycling (3). Long-term dysregulation of calcium cycling has been shown to alter cardiac gene transcription and stimulate hypertrophy (9). In contrast, NOS3, which is localized to the sarcolemma, has no effect on basal calcium cycling in the heart, but is a potent cause of systemic hypertension due to effects on the peripheral vasculature (12). These different functional effects, i.e., altered calcium cycling and altered vascular load, may mediate the different downstream changes in cardiac gene transcription that we observed in NOS1-/- and NOS3-/-. The selective downregulation of MCIP-1 in NOS1-/- hearts, which we demonstrated both on microarray and on qPCR (Figs. 2 and 4), offers another link between NOS1 deletion and calcium signaling. MCIP-1 is an endogenous inhibitor of the calcineurin signaling pathway. When transgenically overexpressed, MCIP-1 prevents hypertrophy induced by calcineurin overexpression, ß-adrenergic stimulation, or exercise training (25). Decreased MCIP-1 in NOS1-/- may therefore promote calcineurin signaling and contribute to hypertrophy.
Although we have emphasized the divergent transcriptional changes in NOS-/- mice, some genes showed similar changes in expression in both strains (Table 2). Most notably, a number of HSPs are markedly reduced in both NOS-/- strains. In particular, HSP70 is downregulated, both on microarray and on qPCR (Fig. 4A). HSP70 protects against adverse ventricular remodeling and prevents myocardial infarction in the setting of ischemia/reperfusion injury (4, 20). HSP downregulation in NOS-/- hearts is consistent with previous findings that NO induces HSP70 expression via activation of the transcription factor HSF1 (34). Impaired HSP function in NOS deficiency may contribute to worsening hypertrophy and cardiac dysfunction. Another noteworthy finding is that Bag3, an inhibitor of apoptosis, is underexpressed in both NOS-/- strains (Table 2), consistent with the well-established role of apoptosis in cardiac hypertrophy and remodeling (6).
Among the transgenic causes of hypertrophy explored in model systems, hypertrophy caused by NOS deletion may be especially relevant to human congestive heart failure. Heart failure is characterized by an imbalance between the formation and degradation of reactive oxygen species within the myocardium (14, 18), and the resultant oxidative stress is sufficient to stimulate cardiac hypertrophy (21, 28). Inactivation of NO signaling by reactive oxygen species may be one of the major mechanisms whereby oxidative stress promotes hypertrophy (10). Moreover, treatments that have been shown to attenuate cardiac hypertrophy in heart failure (angiotensin converting enzyme inhibitors, beta-blockers, and hydralazine/nitrates), increase NO bioavailability and promote NO signaling (1, 33). Our findings identify specific transcriptional targets of NO signaling in cardiac hypertrophy, some of which may be targets for pharmacological intervention.
Expression profiling is a powerful technique that offers tremendous promise, but it creates substantial analytic challenges resulting from the analysis of many genes in a small number of samples (22). We addressed these at multiple levels. Of the methods to convert Affymetrix probe set data into measures of gene expression, which include software supplied by the manufacturer, model-based methods (16), and RMA, we chose RMA based on its superiority in the analysis of small data sets (13). For the exploratory analyses, we used strict selection criteria to define differential gene expression to select biologically significant genes and to minimize spurious findings that result from multiple comparisons. Furthermore, we independently validated findings of interest using qPCR. We also used a candidate-gene approach to detect changes in highly relevant genes that could be overlooked in exploratory analyses.
Several limitations warrant mention. Microarray technologies only assess the abundance of transcripts, which may not correlate with the abundance, localization, or function of corresponding proteins. However, many genes (e.g., HSPs) are principally regulated at the level of transcription and therefore may be accurately characterized by microarrays. We chose to use the entire heart to limit the impact of tissue heterogeneity on quantification of transcript abundance; however, this prevented us from discerning chamber-specific alterations in gene transcription. This may be important in light of the tendency of NOS3-/- mice to develop pulmonary hypertension and right ventricular hypertrophy when they are subjected to chronic hypoxemia (29). Although we analyzed 71 candidate hypertrophy genes, our list of candidates was limited to genes present on the microarray, and not all relevant hypertrophy genes could be analyzed. Our data were obtained in murine models of age-associated LVH. For this reason, our findings may not directly parallel the pathogenesis of LVH in humans and may not reflect the early transcriptional changes caused by hypertension or myocardial infarction. However, prior investigations have demonstrated altered NOS signaling in aged hearts, suggesting that age-associated hypertrophy in NOS knockouts is a pathophysiologically relevant phenotype (37).
In summary, our analyses show that deficiency of different NOS isoforms activates divergent transcriptional programs in cardiac hypertrophy. In addition, deficiency of HSPs may play a central role in initiating or potentiating LVH, and downregulation of MCIP-1 may promote LVH in the absence of NOS1. These findings add to a growing literature suggesting that the transcriptional basis of LVH varies depending on the inciting biochemical stimulus. Moreover, we provide new insights into the role of NOS isoforms in maintaining normal cardiac architecture.
| DISCLOSURES |
|---|
|
|
|---|
T. P. Cappola was a recipient of a Pfizer Postdoctoral Fellowship in Cardiovascular Medicine. J. M. Hare is a recipient of a Paul Beeson Physician Faculty Scholars in Aging Research Award.
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: J. M. Hare, Johns Hopkins Hospital, Division of Cardiology, 600 N. Wolfe St., Carnegie 568, Baltimore, MD 21287-6568 (E-mail: jhare{at}mail.jhmi.edu).
10.1152/physiolgenomics.00156.2002.
1 The Supplementary Material for this article (Tables A and B) is available online at http://physiolgenomics.physiology.org/cgi/content/full/00156.2002/DC1. ![]()
| REFERENCES |
|---|
|
|
|---|
in pericardial fluid of patients with heart failure: a potential role for in vivo oxidant stress in ventricular dilatation and progression to heart failure. Circulation 97: 15361539, 1998.This article has been cited by other articles:
![]() |
R. Howden, E. Liu, L. Miller-DeGraff, H. L. Keener, C. Walker, J. A. Clark, P. H. Myers, D. C. Rouse, T. Wiltshire, and S. R. Kleeberger The genetic contribution to heart rate and heart rate variability in quiescent mice Am J Physiol Heart Circ Physiol, July 1, 2008; 295(1): H59 - H68. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. K. Choate, S. M. Murphy, R. Feldman, and C. R. Anderson Sympathetic control of heart rate in nNOS knockout mice Am J Physiol Heart Circ Physiol, January 1, 2008; 294(1): H354 - H361. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Yang and L. A. Barouch Leptin Signaling and Obesity: Cardiovascular Consequences Circ. Res., September 14, 2007; 101(6): 545 - 559. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. P. Saliaris, L. C. Amado, K. M. Minhas, K. H. Schuleri, S. Lehrke, M. St. John, T. Fitton, C. Barreiro, C. Berry, M. Zheng, et al. Chronic allopurinol administration ameliorates maladaptive alterations in Ca2+ cycling proteins and beta-adrenergic hyporesponsiveness in heart failure Am J Physiol Heart Circ Physiol, March 1, 2007; 292(3): H1328 - H1335. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Liang and B. Ventura Physiological genomics in PG and beyond: July to September 2005 Physiol Genomics, October 17, 2005; 23(2): 119 - 124. [Full Text] [PDF] |
||||
![]() |
C. Ojaimi, W. Li, S. Kinugawa, H. Post, A. Csiszar, P. Pacher, G. Kaley, and T. H. Hintze Transcriptional basis for exercise limitation in male eNOS-knockout mice with age: heart failure and the fetal phenotype Am J Physiol Heart Circ Physiol, October 1, 2005; 289(4): H1399 - H1407. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Jacobson, J. W. Barnard, D. N. Grigoryev, S.-F. Ma, R. M. Tuder, and J. G. N. Garcia Simvastatin attenuates vascular leak and inflammation in murine inflammatory lung injury Am J Physiol Lung Cell Mol Physiol, June 1, 2005; 288(6): L1026 - L1032. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Q. Ye, B. A. Simon, J. P. Maloney, A. Zambelli-Weiner, L. Gao, A. Grant, R. B. Easley, B. J. McVerry, R. M. Tuder, T. Standiford, et al. Pre-B-Cell Colony-enhancing Factor as a Potential Novel Biomarker in Acute Lung Injury Am. J. Respir. Crit. Care Med., February 15, 2005; 171(4): 361 - 370. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Khan, K. Lee, K. M. Minhas, D. R. Gonzalez, S. V. Y. Raju, A. D. Tejani, D. Li, D. E. Berkowitz, and J. M. Hare From the Cover: Neuronal nitric oxide synthase negatively regulates xanthine oxidoreductase inhibition of cardiac excitation-contraction coupling PNAS, November 9, 2004; 101(45): 15944 - 15948. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Zheng, Q.-L. Cheng, A.-R. Plati, S. Q. Ye, M. Berho, A. Banerjee, M. Potier, E. A. Jaimes, H. Yu, Y.-F. Guan, et al. The Glomerulosclerosis of Aging in Females: Contribution of the Proinflammatory Mesangial Cell Phenotype to Macrophage Infiltration Am. J. Pathol., November 1, 2004; 165(5): 1789 - 1798. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Jacobson, S. M. Dudek, K. G. Birukov, S. Q. Ye, D. N. Grigoryev, R. E. Girgis, and J. G. N. Garcia Cytoskeletal Activation and Altered Gene Expression in Endothelial Barrier Regulation by Simvastatin Am. J. Respir. Cell Mol. Biol., May 1, 2004; 30(5): 662 - 670. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. B. Pilz and D. E. Casteel Regulation of Gene Expression by Cyclic GMP Circ. Res., November 28, 2003; 93(11): 1034 - 1046. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Hemish, N. Nakaya, V. Mittal, and G. Enikolopov Nitric Oxide Activates Diverse Signaling Pathways to Regulate Gene Expression J. Biol. Chem., October 24, 2003; 278(43): 42321 - 42329. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP |