| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Clinical Research |
From the Departments of Medicine I (A.S.B., E.A., L.Z., P.K., S.K., M.H., H.K., G.S., M.N.) and Cardiac Surgery (E.K.), University Hospital Grosshadern, and the Department of Medical Informatics, Biometrics, and Epidemiology (S.M., M.D.), Ludwig-Maximilians-University, Munich, Germany; Aventis Pharma Deutschland GmbH (M.G., K.S., M.B.), Frankfurt am Main, Germany; and the Institute of Physiology (M.B.), Christian- Albrechts-University, Kiel, Germany.
Correspondence to Dr med Andreas S. Barth, Department of Medicine I, University Hospital Grosshadern, Marchioninistrasse 15, 81377 Munich, Germany. E-mail andreas.barth{at}med.uni-muenchen.de
| Abstract |
|---|
|
|
|---|
Key Words: atrial fibrillation gene expression atrial myocardium functional genomics
| Introduction |
|---|
|
|
|---|
Moreover, atrial fibrillation has been shown to be associated with an atrial-to-ventricular switch by expressing ventricular myosin isoforms, which is regarded as part of a dedifferentiation process of fibrillating atrial myocardium.5 To explore the hypothesis that pmAF leads to upregulation of ventricular transcripts other than those encoding for contractile proteins, we used a genome-wide approach to compare atrial mRNA expression of patients with pmAF to atrial mRNA expression of patients with sinus rhythm (SR) as well as to ventricular (LV) gene expression. Our findings indicate that pmAF leads to repression of atrial transcriptional activity and induces a ventricular-like genomic signature in atrial myocardium.
| Materials and Methods |
|---|
|
|
|---|
mRNA Preparation and Hybridization to Affymetrix U133A+B Oligonucleotide Microarrays
Total RNA was purified from homogenized deep-frozen tissue samples following the TRIZOL standard protocol as described by the manufacturer (GibcoBRL) and quantified by photometry. Then, 150 ng total RNA was used for quality control using the RNA 6000 Nano LabChip kit and Agilent Bioanalyzer 2100 (Agilent Technologies).
Synthesis of double-stranded cDNA and of biotinylated cRNA probes as well as subsequent steps of hybridization, washing, and staining followed instructions provided by Affymetrix.
Data Extraction and Statistical Analysis
One set of U133A+B chips (together representing 44 928 probe sets) was used for each human heart sample. The original data files have been deposited in Gene Expression Omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo) and are accessible through GEO Series accession number GSE2240. All U133 GeneChips passed quality control to eliminate scans with abnormal characteristics, ie, abnormal low or high dynamic range, high perfect match saturation, high pixel noise, grid misalignment problems, and low mean signal-to-noise ratio. Low level and statistical analysis was done using R 1.8.1. Microarray data of all atrial and ventricular samples was normalized simultaneously using robust multiarray average (RMA)6 implemented in Bioconductors "affy"-package. Probe sets with "absent" calls in more than 50% of tissue samples in either group (SR, pmAF, or LV) were filtered out and omitted from further analysis. To determine differentially expressed genes between the groups, Welch t test was applied. Considering multiple hypothesis testing, differences in expression were regarded as statistically significant if a false discovery rate (FDR) of q <0.03 was achieved.7 Duplicate probe sets within the gene lists were eliminated by means of the "MatchMiner" software tool.8 Hierarchical clustering was performed by using an average linkage algorithm with Euclidean distance measure comprised in the Genesis software.9 Functional classification was based on the "cellular component," "biological process," and "molecular function" of GO using "FatiGO."10 After accounting for uncharacterized transcripts, up to 70% of all known genes were annotated across all GO levels. As many genes are not annotated at a particular level of ontology, three adjacent ontology levels are presented to increase the amount of information.
Validation by Semiquantitative Real-Time Reverse Transcription PCR
Expression patterns of 13 genes, selected with respect to the central role in the proposed working model of AF (Figure 3) were validated by real-time SYBR-Green PCR. A detailed list of the genes examined by RT-PCR including the design of the primers and a description of the protocols used for RT-PCR are given in the expanded Materials and Methods section in the online data supplement.
| Results |
|---|
|
|
|---|
and NFATc3, and concomitant upregulation of NFATc3-inhibitory transcripts, GSK3B and DSCR1. Moreover, deregulation of other calcium-dependent pathways, including CaMK- and MAPK-pathways, was observed, as 7 of 9 transcripts encoding for genes of MAPK pathways with altered expression in pmAF were downregulated (online Table 3). In addition, downregulation of several additional transcription factors like NF-
B, CREB1, STAT3, and STAT5 emphasized the general repression of transcriptional processes in pmAF.
Transcriptional Basis of Distinct Remodeling Pathways in Fibrillating Atrium
Based on GO, distinct functional classes could be identified for differentially expressed genes (online Tables 4 and 5), supporting the notion of several distinct areas of remodeling: consistent with prominent fibrosis and structural remodeling in pmAF, we noted changes in the expression of genes related to extracellular matrix composition and turnover with upregulation of collagen transcripts and extracellular matrix-cleaving proteases, like MMP9 and ADAM19. Downregulation of the corresponding tissue inhibitors of matrix metalloproteinases (eg, TIMP3, showing the most significant change of all genes deregulated in pmAF) also adds to adverse matrix remodeling. In addition, transforming growth factor (TGF)-ß1, a particularly important profibrotic cytokine, was upregulated. Downregulation of the TGF-ß1inhibitory transcript decorin (DCN) and upregulation of the insulin-like growth factor 2 receptor (IGF2R), which controls the extracellular bioavailability of TGF-ß is expected to augment TGF-ß1 levels. Deregulation of TGF-ß signaling is also exemplified by downregulation of signaling molecules downstream of TGF-ß1, including TGF-ß receptor 2, fibrillin 1, and SMAD4. Besides TGF-ß, we found evidence for transcriptional activation of additional profibrotic growth factors (upregulation of platelet-derived growth factor, PDGF, and concomitant downregulation of the PDGF receptor; downregulation of the fibroblast growth factor receptors 1 and 2, FGFR1+2) and neurohumoral signaling pathways (downregulation of endothelin receptors A and B, EDNRA and EDNRB; online Tables 3 and 5, "cellular component" and "molecular function," level 5).
There is also evidence for contractile remodeling in pmAF, as troponin I1 (TNNI1) was downregulated. Moreover, deregulation of several genes linked to a cardiomyopathic phenotype was observed (downregulation of ErbB2, ErbB3, and ErbB4; downregulation of lamin A/C (LMNA) and several genes involved in the formation of the dystrophin glycoprotein complex, including dystrophin, dystroglycan 1,
-sarcoglycan, and ß2-syntrophin; online Table 3).
Regarding electrophysiological remodeling, 5 transcripts encoding for subunits of L-type calcium channels were found to exhibit downregulation (online Table 3). The repolarization process was also disturbed, reflected in deregulation of several potassium channel transcripts (KCNIP2, KCNJ3, KCNJ5, KCNJ2, and KCNK3). In addition to changes of sarcolemmal ion channels, sarcoplasmic Ca2+ handling is expected to be altered by upregulation of RYR3 (a dominant-negative splice variant of the ryanodine receptor) and downregulation of ATP2A2 (sarcoplasmic Ca2+-ATPase). Changes in expression levels of KCNK3 and 12 additional transcripts were validated using real-time RT-PCR assays, confirming expression changes for 10 of 13 selected transcripts (Table 1).
|
Energy processes emerged as further area of pronounced changes, as 5 of 12 GO classes for biological processes found to be deregulated in pmAF serve metabolic functions (online Table 4, "biological process," level 5). Except for fatty acid and glycoprotein metabolism, there was a trend toward upregulation for all metabolic classes, as identified by GO (Figure 1 and online Table 4, "biological process," level 5). Upregulation of genes involved in the formation of the mitochondrial inner membrane equally demonstrates adaptation to an increased metabolic demand in pmAF (online Table 4, "cellular component," level 4).
|
Permanent AF Leads to Expression of a Ventricular-Like Signature in Atrium
Unsupervised clustering of genes deregulated in pmAF for all 30 atrial (pmAF+SR) and 5 LV samples revealed that atrial pmAF samples clustered with ventricular samples, rather than atrial SR samples (Figure 1 in the online data supplement). This is unexpected, as we identified 4 times more transcripts to be differentially expressed between atrium and ventricle than between pmAF and SR using identical technology and statistical methods: 3300 and 2974 transcripts showed higher expression values in atria and ventricles, respectively.11 Comparison of atrial- and ventricular-predominant transcripts with the genes deregulated in pmAF revealed that the ventricular-type transcription pattern of fibrillating atrial myocardium is mainly the result of downregulation of atrial-specific transcripts and only to a lesser degree of upregulation of ventricular-specific transcripts (Table 2).
|
Expression of a ventricular signature is also reflected in GO classes, as functional classes previously shown to be characteristic of atrial myocardium, such as those involved in signal transduction and cell communication,11 were downregulated in pmAF, whereas functional classes overrepresented in ventricular myocardium (eg, metabolic processes11) were found to be upregulated in pmAF (Figure 1 and online Table 4, "biological process").
| Discussion |
|---|
|
|
|---|
Transcriptional Basis of Distinct but Interconnected Remodeling Pathways in Fibrillating Atrium
Remodeling in AF has been recognized to consist of 3 major processes, affecting structural, contractile, and electrophysiological properties.1 Analysis based on GO identified major changes in several genes connecting these seemingly distinct areas of remodeling. For instance, downregulation of TIMP3 represents a possible link between extracellular matrix remodeling and depressed contractile function, as TIMP3-deficient mice have recently been shown to develop a cardiomyopathic phenotype with chamber dilatation and contractile dysfunction.12 Likewise, conditional mutations of the ErbB2 receptor in myocardium lead to dilated cardiomyopathy, and prolonged exposure of cultured human and rat cardiomyocytes with anti-ErbB2 antibody causes loss of contractile function.13 Thus, downregulation of the ErbB2 receptor might not only provide an explanation for extracellular matrix remodeling with dilatation in pmAF, but may also be linked to the contractile dysfunction that persists for some time after cardioversion of AF. Similarly, downregulation of several genes involved in the formation of the dystrophin glycoprotein complex, which is essential in providing the mechanical link from the intracellular cytoskeleton to the extracellular matrix, would also promote a cardiomyopathic phenotype.14
Other mechanisms, like endothelin activation, initiated to compensate for the contractile insufficiency in pmAF, may also be maladaptive in the long run by promoting cardiac fibrosis through activation of TGF-ß. Alterations of TGF-ß signaling have been connected to a number of cardiovascular and fibrotic diseases.15 Loss-of-function mutations in TGF-ß receptor 2 and fibrillin 1 have recently been linked to Marfan syndrome, an extracellular matrix disorder associated with atrial dilatation and atrial fibrillation.16 This suggests that the observed downregulation of both TGF-ß receptor 2 (TGFBR2) and fibrillin 1 (FBN1) in pmAF might provide the link to extracellular matrix remodeling. Interestingly, one of our AF patients had been diagnosed with Marfan syndrome based on the clinical phenotype. Figure 2 demonstrates individual measurements for FBN1 and TGFBR2 in all tissue specimens, illustrating that downregulation of these genes was not restricted to the patient with Marfan syndrome but common to all AF patients. Considering further genes for which mutations are associated with hereditary atrial fibrillation, it is important to note that Sébillon et al17 identified a lamin A/C missense mutation with early-onset AF as presenting symptom, preceding the development of dilated cardiomyopathy. Based on the observation that the same genes that have been shown to lead to hereditary forms of dilated atrial chambers and contractile dysfunction are also downregulated in our patients with acquired, valvular AF (eg, LMNA, TGFBR2, FBN1, TIMP3, and members of the dystrophin complex), it is tempting to speculate that these genes also contribute significantly to atrial dilatation and contractile dysfunction in acquired AF. Although transcriptional changes observed in the present study are the consequence rather than the cause of pmAF, it appears that some of the changes that occur initially to provide a substrate for pmAF are similar to alterations in the course of AF supporting the concept that AF promotes its own substrate ("AF begets AF").18
|
In addition to the well-established remodeling processes affecting contractile, electrophysiological, and structural properties, we identified prominent deregulation of genes belonging to metabolic classes, suggesting changes in atrial bioenergetics with an increased energy demand in pmAF. In line with this finding, cellular hypertrophy is among the most consistent histological findings in fibrillating atrial tissue, suggested to reflect an adaptation to altered environmental conditions.19 Induction of AF has been shown to increase atrial oxygen consumption and coronary flow almost 3-fold,20 possibly the result of high wall stress and the dramatic increase in atrial depolarization frequency with enhanced energy demand for ion transport processes including calcium-handling to maintain physiological cytosolic ion levels. In a goat model of AF, phosphocreatine decreased by 60% in atrial myocytes after 1 week of AF, indicating enhanced demand for high-energy phosphates.21 As phosphocreatine levels increased thereafter and returned to baseline levels after 16 weeks of AF,21 compensatory mechanisms are assumed to be activated to restore the metabolic balance between energy demand and supply in pmAF. Our finding of prominent upregulation of metabolic processes as identified by GO strongly supports this concept (Figure 1). The well-established atrial-to-ventricular switch of myosin isoforms might also be related to the increase in metabolic demand in pmAF: contraction has been shown to be 5-fold more economical in ventricular than atrial myocardium, mainly because of the greater relative contribution of the slow ventricular ß-myosin heavy chain (MHC) isoform compared with the fast atrial
-MHC transcript in ventricular muscle.22 Thus, increased expression of the ventricular-predominant slow ß-MHC isoform in fibrillating atrial myocardium may represent an attempt to improve economy of contraction.
Generalized Transcriptional Downregulation in Permanent AF, Likely Mediated by Altered Calcium Homeostasis
In pmAF, we observed prominent downregulation of transcriptional processes with more than twice as many genes being downregulated than upregulated. The predominant pattern of transcriptional repression is reflected in downregulation of functional classes linked to transcriptional processes as identified by GO. In addition, our data are in agreement with the first microarray study in an animal model of AF5 and a recent study in human myocardium,23 also suggesting predominant downregulation of transcription in AF.
Even if transcript levels of several myocyte-specific markers were unchanged between SR and pmAF samples (cardiac myosin heavy polypeptide 7, cardiac troponin T2, and calsequestrin 2; data not shown), suggesting that the cardiomyocyte compartment has not been diminished in pmAF, a reduction in cardiomyocyte number might still have occurred, as myolysis and cellular hypertrophy are well-established features of pmAF.24 Therefore, because of simultaneous reduction of myocyte number and increase in the expression of cardiac-specific genes per single cell no change in overall cardiac-specific gene expression might be apparent in pmAF samples. As laser microdissection has not been performed, we cannot determine transcription on a myocyte level. However, the degree of AF-induced myolysis is expected to be modest24 and cannot fully account for the pronounced reduction of transcriptional processes.
As to the mechanisms regulating gene transcription in myocardial cells, calcium homeostasis, contractile function, and stretch have been shown to play key roles.25,26 Because all these factors are altered after prolonged AF,2729 they must be considered potential initiators of altered gene transcription. Specifically, calcium is a potent activator of cardiac gene transcription. Although Ca2+ overload has clearly been established as an early event after initiation of AF,1 there is no direct experimental evidence to support the notion of a chronic calcium overload state in persistent and pmAF. In fact, after an initial increase, steady-state calcium levels begin to decrease as early as 30 minutes after initiation of AF and normalize within the first 48 hours.27 Moreover, Ausma and coworkers28 demonstrated decreased steady-state calcium levels after 16 weeks of AF in goats. Likewise, isolated canine atrial myocytes displayed decreased resting intracellular calcium levels and blunted calcium transients after a 4-week period of rapid atrial pacing.29 The idea that high intracellular calcium levels may no longer be present after prolonged AF is strongly supported by the experimental and clinical observation that calcium-channel blockers like verapamil are useful in preventing calcium-overload and electrical remodeling only when given very early or during paroxysms of AF but fail to have beneficial effects in persistent and permanent AF.30 Given our observation of transcriptional repression of several calcium-dependent pathways including calcineurin-NFATc3, MAPK, and CaMK in pmAF and the central, well-established role of calcium as a regulator of gene expression, we propose that the decrease in transcription during prolonged episodes of AF results from alterations in calcium homeostasis with alterations in the spatiotemporal properties of the calcium transients (Figure 3).
|
Permanent AF Leads to a Ventricular- and Fetal-Like Genomic Signature in Atrial Myocardium
Atrial and ventricular myocardium exhibit fundamental differences on a morphological, ultrastructural, and functional basis, because ventricular myocardium, above all, functions as contractile pump, thereby contrasting to the prominent neuroendocrine and regulatory functions of atrial myocardium.11,31
In AF, an atrial-to-ventricular switch has been shown for some individual genes in atrial samples of hypertrophied32 and cardiomyopathic hearts,33 eg, contractile protein isoforms like myosin heavy chains. It has been hypothesized that the atrial-to-ventricular switch of MHCs observed in these pathophysiological conditions is the result of the rise in intra-atrial pressure, similar to hemodynamic changes observed during postnatal development.32 Our approach of microarray analysis with simultaneous examination of thousands of genes extends these observations in single genes to a general pattern of regulation toward a ventricular gene program in AF. Central in the changes are genes related to energy metabolism, also reflected in ultrastructural changes occurring during long-term AF that resemble those seen in fetal cardiomyocytes and in hibernating myocardium. This led Ausma et al19 to put forward the concept of "cellular adaptation through dedifferentiation" in AF, where a fetal phenotype is thought to facilitate cardiomyocyte survival during stress conditions including stretch and ischemia.34 The coordinated transcriptional downregulation of enzymes controlling fatty acid oxidation and concomitant upregulation of enzymes involved in glucose utilization reported in our study is consistent with reinduction of a fetal energy metabolic program35 (Figure 1 and online Table 4). Equally, ultrastructural characteristics in AF like the perinuclear accumulation of glycogen are seen in fetal cardiomyocytes and are suggestive of a metabolic shift,36 as the glycogen storage is thought to act as supply for high-energydemanding processes and to enhance mechanical cell stability.37 Thus, the observed expression of a ventricular and fetal gene program with upregulation of metabolic genes and downregulation of atrial-specific genes appears to be an integral part of this adaptive response to pmAF. Given the high metabolic demand of fibrillating atrial myocardium, one might also speculate that the extensive suppression of transcriptional processes and hence protein synthesis is an attempt to save energy by restricting synthesis to a smaller number of proteins, possibly similar to chronic hibernation.38
Limitation of the Present Study
Analysis of human myocardial tissue is complicated by biological variability, which is expected to be greater than in well-controlled animal models of AF. Different etiologies and duration of AF, differences in age, gender, and medications, as well as individual course of the disease of patients from which tissue specimens were obtained contribute to variability in gene expression when examining human tissue samples. In the present study, however, we were able to detect robust functional patterns extending beyond this biological variability. The validity of our data-set is corroborated by expression patterns of a wide variety of genes known to be regulated in pmAF (including NPPB,39 several transcripts of L-type calcium currents,3 KCNJ2,3 KCNJ3,3 KCNJ5,3 KCNIP2,40 ATP2A2,3 MMP9,2 vWF,41 EDNRA,42 EDNRB,42 and mortalin43). Still, many transcripts reported to be regulated in pmAF were not found to be differentially expressed between SR and pmAF in the present study (eg, different myosin isoforms5). In addition to the biological heterogeneity inherent to any study of human tissue specimens as outlined here, this is partly because of the incomplete annotational status of the human genome at the time the Affymetrix U133A+B arrays were designed. Even if we seek to present a comprehensive portrait of the transcriptional changes in pmAF, this picture can by no means be complete. Moreover, it has to be kept in mind that mRNA levels do not always directly translate into changes in protein or functional levels, as additional regulatory mechanisms on a posttranslational level exist.
AF is usually initiated and maintained in the left atrium, because the cycle length is shorter in left than in right atrium during AF. Accordingly, a recent microarray study of a pig model of AF suggested that the genomic changes are more pronounced in left than right atrium with considerable overlap in the genomic response.5 As only right atrial tissue was available for the present study, future studies will be needed to account for regional heterogeneity in AF at the atrial level.
| Acknowledgments |
|---|
| Footnotes |
|---|
This manuscript contains part of Elisabeth Arnoldis doctoral thesis, in preparation at the Medical Faculty of the Ludwig-Maximilians-University, Munich, Germany.
Original received November 18, 2004; revision received March 24, 2005; accepted March 28, 2005.
| References |
|---|
|
|
|---|
Related Article:
This article has been cited by other articles:
![]() |
B. Burstein and S. Nattel Atrial Fibrosis: Mechanisms and Clinical Relevance in Atrial Fibrillation J. Am. Coll. Cardiol., February 26, 2008; 51(8): 802 - 809. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mayr, S. Yusuf, G. Weir, Y.-L. Chung, U. Mayr, X. Yin, C. Ladroue, B. Madhu, N. Roberts, A. De Souza, et al. Combined metabolomic and proteomic analysis of human atrial fibrillation. J. Am. Coll. Cardiol., February 5, 2008; 51(5): 585 - 594. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Adam, G. Frost, F. Custodis, M. A. Sussman, H.-J. Schafers, M. Bohm, and U. Laufs Role of Rac1 GTPase Activation in Atrial Fibrillation J. Am. Coll. Cardiol., July 24, 2007; 50(4): 359 - 367. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Cardin, E. Libby, P. Pelletier, S. Le Bouter, A. Shiroshita-Takeshita, N. Le Meur, J. Leger, S. Demolombe, A. Ponton, L. Glass, et al. Contrasting Gene Expression Profiles in Two Canine Models of Atrial Fibrillation Circ. Res., February 16, 2007; 100(3): 425 - 433. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. S. Barth, R. Kuner, A. Buness, M. Ruschhaupt, S. Merk, L. Zwermann, S. Kaab, E. Kreuzer, G. Steinbeck, U. Mansmann, et al. Identification of a Common Gene Expression Signature in Dilated Cardiomyopathy Across Independent Microarray Studies J. Am. Coll. Cardiol., October 17, 2006; 48(8): 1610 - 1617. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. El-Armouche, P. Boknik, T. Eschenhagen, L. Carrier, M. Knaut, U. Ravens, and D. Dobrev Molecular Determinants of Altered Ca2+ Handling in Human Chronic Atrial Fibrillation Circulation, August 15, 2006; 114(7): 670 - 680. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. V. Podgoreanu and J. P. Mathew Prophylaxis Against Postoperative Atrial Fibrillation: Current Progress and Future Directions JAMA, December 28, 2005; 294(24): 3140 - 3142. [Full Text] [PDF] |
||||
![]() |
M. V. Podgoreanu and D. A. Schwinn New Paradigms in Cardiovascular Medicine: Emerging Technologies and Practices: Perioperative Genomics J. Am. Coll. Cardiol., December 6, 2005; 46(11): 1965 - 1977. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. F. Tomaselli Ventricularization of Atrial Gene Expression in the Fibrillating Heart? Circ. Res., May 13, 2005; 96(9): 923 - 924. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||