Integrative multiomics analysis of human atherosclerosis reveals a serum response factor-driven network associated with intraplaque hemorrhage.

BACKGROUND: While single-omics analyses on human atherosclerotic plaque have been very useful to map stage- or disease-related differences in expression, they only partly capture the array of changes in this tissue and suffer from scale-intrinsic limitations. In order to better identify processes associated with intraplaque hemorrhage and plaque instability, we therefore combined multiple omics into an integrated model. METHODS: In this study, we compared protein and gene makeup of low- versus high-risk atherosclerotic lesion segments from carotid endarterectomy patients, as judged from the absence or presence of intraplaque hemorrhage, respectively. Transcriptomic, proteomic, and peptidomic data of this plaque cohort were aggregated and analyzed by DIABLO, an integrative multivariate classification and feature selection method. RESULTS: We identified a protein-gene associated multiomics model able to segregate stable, nonhemorrhaged from vulnerable, hemorrhaged lesions at high predictive performance (AUC >0.95). The dominant component of this model correlated with αSMA- PDGFRα+ fibroblast-like cell content (p = 2.4E-05) and Arg1+ macrophage content (p = 2.2E-04) and was driven by serum response factor (SRF), possibly in a megakaryoblastic leukemia-1/2 (MKL1/2) dependent manner. Gene set overrepresentation analysis on the selected key features of this model pointed to a clear cardiovascular disease signature, with overrepresentation of extracellular matrix synthesis and organization, focal adhesion, and cholesterol metabolism terms, suggestive of the model's relevance for the plaque vulnerability. Finally, we were able to corroborate the predictive power of the selected features in several independent mRNA and proteomic plaque cohorts. CONCLUSIONS: In conclusion, our integrative omics study has identified an intraplaque hemorrhage-associated cardiovascular signature that provides excellent stratification of low- from high-risk carotid artery plaques in several independent cohorts. Further study revealed suppression of an SRF-regulated disease network, controlling lesion stability, in vulnerable plaque, which can serve as a scaffold for the design of targeted intervention in plaque destabilization.

Hypoxia increases genome-wide bivalent epigenetic marking by specific gain of H3K27me3.

BACKGROUND: Trimethylation at histone H3 lysine 4 (H3K4me3) and lysine 27 (H3K27me3) controls gene activity during development and differentiation. Whether H3K4me3 and H3K27me3 changes dynamically in response to altered microenvironmental conditions, including low-oxygen conditions commonly present in solid tumors, is relatively unknown. Demethylation of H3K4me3 and H3K27me3 is mediated by oxygen and 2-oxoglutarate dioxygenases enzymes, suggesting that oxygen deprivation (hypoxia) may influence histone trimethylation. Using the MCF7 breast epithelial adenocarcinoma cell model, we have determined the relationship between epigenomic and transcriptomic reprogramming as a function of fluctuating oxygen tension. RESULTS: We find that in MCF7, H3K4me3 and H3K27me3 marks rapidly increase at specific locations throughout the genome and are largely reversed upon reoxygenation. Whereas dynamic changes are relatively highest for H3K27me3 marking under hypoxic conditions, H3K4me3 occupation is identified as the defining epigenetic marker of transcriptional control. In agreement with the global increase of H3K27 trimethylation, we provide direct evidence that the histone H3K27me3 demethylase KDM6B/JMJD3 is inactivated by limited oxygen. In situ immunohistochemical analysis confirms a marked rise of histone trimethylation in hypoxic tumor areas. Acquisition of H3K27me3 at H3K4me3-marked loci results in a striking increase in "bivalent" epigenetic marking. Hypoxia-induced bivalency substantially overlaps with embryonal stem cell-associated genic bivalency and is retained at numerous loci upon reoxygenation. Transcriptional activity is selectively and progressively dampened at bivalently marked loci upon repeated exposure to hypoxia, indicating that this subset of genes uniquely maintains the potential for epigenetic regulation by KDM activity. CONCLUSIONS: These data suggest that dynamic regulation of the epigenetic state within the tumor environment may have important consequences for tumor plasticity and biology.

Quantitative analysis of ChIP-seq data uncovers dynamic and sustained H3K4me3 and H3K27me3 modulation in cancer cells under hypoxia.

BACKGROUND: A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understanding the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treatment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets. RESULTS: We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around centromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia-induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells. CONCLUSIONS: In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation.

CyTargetLinker: a cytoscape app to integrate regulatory interactions in network analysis.

INTRODUCTION: The high complexity and dynamic nature of the regulation of gene expression, protein synthesis, and protein activity pose a challenge to fully understand the cellular machinery. By deciphering the role of important players, including transcription factors, microRNAs, or small molecules, a better understanding of key regulatory processes can be obtained. Various databases contain information on the interactions of regulators with their targets for different organisms, data recently being extended with the results of the ENCODE (Encyclopedia of DNA Elements) project. A systems biology approach integrating our understanding on different regulators is essential in interpreting the regulation of molecular biological processes. IMPLEMENTATION: We developed CyTargetLinker (http://projects.bigcat.unimaas.nl/cytargetlinker), a Cytoscape app, for integrating regulatory interactions in network analysis. Recently we released CyTargetLinker as one of the first apps for Cytoscape 3. It provides a user-friendly and flexible interface to extend biological networks with regulatory interactions, such as microRNA-target, transcription factor-target and/or drug-target. Importantly, CyTargetLinker employs identifier mapping to combine various interaction data resources that use different types of identifiers. RESULTS: Three case studies demonstrate the strength and broad applicability of CyTargetLinker, (i) extending a mouse molecular interaction network, containing genes linked to diabetes mellitus, with validated and predicted microRNAs, (ii) enriching a molecular interaction network, containing DNA repair genes, with ENCODE transcription factor and (iii) building a regulatory meta-network in which a biological process is extended with information on transcription factor, microRNA and drug regulation. CONCLUSIONS: CyTargetLinker provides a simple and extensible framework for biologists and bioinformaticians to integrate different regulatory interactions into their network analysis approaches. Visualization options enable biological interpretation of complex regulatory networks in a graphical way. Importantly the incorporation of our tool into the Cytoscape framework allows the application of CyTargetLinker in combination with a wide variety of other apps for state-of-the-art network analysis.

An evaluation of two-channel ChIP-on-chip and DNA methylation microarray normalization strategies.

BACKGROUND: The combination of chromatin immunoprecipitation with two-channel microarray technology enables genome-wide mapping of binding sites of DNA-interacting proteins (ChIP-on-chip) or sites with methylated CpG di-nucleotides (DNA methylation microarray). These powerful tools are the gateway to understanding gene transcription regulation. Since the goals of such studies, the sample preparation procedures, the microarray content and study design are all different from transcriptomics microarrays, the data pre-processing strategies traditionally applied to transcriptomics microarrays may not be appropriate. Particularly, the main challenge of the normalization of "regulation microarrays" is (i) to make the data of individual microarrays quantitatively comparable and (ii) to keep the signals of the enriched probes, representing DNA sequences from the precipitate, as distinguishable as possible from the signals of the un-enriched probes, representing DNA sequences largely absent from the precipitate. RESULTS: We compare several widely used normalization approaches (VSN, LOWESS, quantile, T-quantile, Tukey's biweight scaling, Peng's method) applied to a selection of regulation microarray datasets, ranging from DNA methylation to transcription factor binding and histone modification studies. Through comparison of the data distributions of control probes and gene promoter probes before and after normalization, and assessment of the power to identify known enriched genomic regions after normalization, we demonstrate that there are clear differences in performance between normalization procedures. CONCLUSION: T-quantile normalization applied separately on the channels and Tukey's biweight scaling outperform other methods in terms of the conservation of enriched and un-enriched signal separation, as well as in identification of genomic regions known to be enriched. T-quantile normalization is preferable as it additionally improves comparability between microarrays. In contrast, popular normalization approaches like quantile, LOWESS, Peng's method and VSN normalization alter the data distributions of regulation microarrays to such an extent that using these approaches will impact the reliability of the downstream analysis substantially.

An integrated bioinformatics approach to improve two-color microarray quality-control: impact on biological conclusions.

Omics technology used for large-scale measurements of gene expression is rapidly evolving. This work pointed out the need of an extensive bioinformatics analyses for array quality assessment before and after gene expression clustering and pathway analysis. A study focused on the effect of red wine polyphenols on rat colon mucosa was used to test the impact of quality control and normalisation steps on the biological conclusions. The integration of data visualization, pathway analysis and clustering revealed an artifact problem that was solved with an adapted normalisation. We propose a possible point to point standard analysis procedure, based on a combination of clustering and data visualization for the analysis of microarray data.

Bioinformatics for the NuGO proof of principle study: analysis of gene expression in muscle of ApoE3*Leiden mice on a high-fat diet using PathVisio.

Insulin resistance is a characteristic of type-2 diabetes and its development is associated with an increased fat consumption. Muscle is one of the tissues that becomes insulin resistant after high fat (HF) feeding. The aim of the present study is to identify processes involved in the development of HF-induced insulin resistance in muscle of ApOE3*Leiden mice by using microarrays. These mice are known to become insulin resistant on a HF diet. Differential gene expression was measured in muscle using the Affymetrix mouse plus 2.0 array. To get more insight in the processes, affected pathway analysis was performed with a new tool, PathVisio. PathVisio is a pathway editor customized with plug-ins (1) to visualize microarray data on pathways and (2) to perform statistical analysis to select pathways of interest. The present study demonstrated that with pathway analysis, using PathVisio, a large variety of processes can be investigated. The significantly regulated genes in muscle of ApOE3*Leiden mice after 12 weeks of HF feeding were involved in several biological pathways including fatty acid beta oxidation, fatty acid biosynthesis, insulin signaling, oxidative stress and inflammation.

The public road to high-quality curated biological pathways.

Biological pathways are abstract and functional visual representations of existing biological knowledge. By mapping high-throughput data on these representations, changes and patterns in biological systems on the genetic, metabolic and protein level are instantly assessable. Many public domain repositories exist for storing biological pathways, each applying its own conventions and storage format. A pathway-based content review of these repositories reveals that none of them are comprehensive. To address this issue, we apply a general workflow to create curated biological pathways, in which we combine three content sources: public domain databases, literature and experts. In this workflow all content of a particular biological pathway is manually retrieved from biological pathway databases and literature, after which this content is compared, combined and subsequently curated by experts. From the curated content, new biological pathways can be created for a pathway analysis tool of choice and distributed among its user base. We applied this procedure to construct high-quality curated biological pathways involved in human fatty acid metabolism.

Fasting induces a biphasic adaptive metabolic response in murine small intestine.

BACKGROUND: The gut is a major energy consumer, but a comprehensive overview of the adaptive response to fasting is lacking. Gene-expression profiling, pathway analysis, and immunohistochemistry were therefore carried out on mouse small intestine after 0, 12, 24, and 72 hours of fasting. RESULTS: Intestinal weight declined to 50% of control, but this loss of tissue mass was distributed proportionally among the gut's structural components, so that the microarrays' tissue base remained unaffected. Unsupervised hierarchical clustering of the microarrays revealed that the successive time points separated into distinct branches. Pathway analysis depicted a pronounced, but transient early response that peaked at 12 hours, and a late response that became progressively more pronounced with continued fasting. Early changes in gene expression were compatible with a cellular deficiency in glutamine, and metabolic adaptations directed at glutamine conservation, inhibition of pyruvate oxidation, stimulation of glutamate catabolism via aspartate and phosphoenolpyruvate to lactate, and enhanced fatty-acid oxidation and ketone-body synthesis. In addition, the expression of key genes involved in cell cycling and apoptosis was suppressed. At 24 hours of fasting, many of the early adaptive changes abated. Major changes upon continued fasting implied the production of glucose rather than lactate from carbohydrate backbones, a downregulation of fatty-acid oxidation and a very strong downregulation of the electron-transport chain. Cell cycling and apoptosis remained suppressed. CONCLUSION: The changes in gene expression indicate that the small intestine rapidly looses mass during fasting to generate lactate or glucose and ketone bodies. Meanwhile, intestinal architecture is maintained by downregulation of cell turnover.

Linking microarray reporters with protein functions.

BACKGROUND: The analysis of microarray experiments requires accurate and up-to-date functional annotation of the microarray reporters to optimize the interpretation of the biological processes involved. Pathway visualization tools are used to connect gene expression data with existing biological pathways by using specific database identifiers that link reporters with elements in the pathways. RESULTS: This paper proposes a novel method that aims to improve microarray reporter annotation by BLASTing the original reporter sequences against a species-specific EMBL subset, that was derived from and crosslinked back to the highly curated UniProt database. The resulting alignments were filtered using high quality alignment criteria and further compared with the outcome of a more traditional approach, where reporter sequences were BLASTed against EnsEMBL followed by locating the corresponding protein (UniProt) entry for the high quality hits. Combining the results of both methods resulted in successful annotation of > 58% of all reporter sequences with UniProt IDs on two commercial array platforms, increasing the amount of Incyte reporters that could be coupled to Gene Ontology terms from 32.7% to 58.3% and to a local GenMAPP pathway from 9.6% to 16.7%. For Agilent, 35.3% of the total reporters are now linked towards GO nodes and 7.1% on local pathways. CONCLUSION: Our methods increased the annotation quality of microarray reporter sequences and allowed us to visualize more reporters using pathway visualization tools. Even in cases where the original reporter annotation showed the correct description the new identifiers often allowed improved pathway and Gene Ontology linking. These methods are freely available at http://www.bigcat.unimaas.nl/public/publications/Gaj_Annotation/.

Impact of supervised gene signatures of early hypoxia on patient survival.

BACKGROUND AND PURPOSE: Hypoxia is a common feature of solid tumors associated with therapy resistance, increased malignancy and poor prognosis. Several approaches have been developed with the hope of identifying patients harboring hypoxic tumors including the use of microarray based gene signatures. However, studies to date have largely ignored the strong time dependency of hypoxia-regulated gene expression. We hypothesized that use of time-dependent patterns of gene expression during hypoxia would enable development of superior prognostic expression signatures. MATERIALS AND METHODS: Using published data from the microarray study of Chi et al., we extracted gene signatures correlating with induction during either early or late hypoxic exposure. Gene signatures were derived from in vitro exposed human mammary epithelial cell line (HMEC) under 0% or 2% oxygen. Gene signatures correlating with early and late up-regulation were tested by means of Kaplan-Meier survival, univariate, and multivariate analysis on a patient data set with primary breast cancer treated conventionally (surgery plus on indication radiotherapy and systemic therapy). RESULTS: We found that the two early hypoxia gene signatures extracted from 0% and 2% hypoxia showed significant prognostic power (log-rank test: p=0.004 at 0%, p=0.034 at 2%) in contrast to the late hypoxia signatures. Both early gene signatures were linked to the insulin pathway. From the multivariate Cox-regression analysis, the early hypoxia signature (p=0.254) was found to be the 4th best prognostic factor after lymph node status (p=0.002), tumor size (p=0.016) and Elston grade (p=0.111). On this data set it indeed provided more information than ER status or p53 status. CONCLUSIONS: The hypoxic stress elicits a wide panel of temporal responses corresponding to different biological pathways. Early hypoxia signatures were shown to have a significant prognostic power. These data suggest that gene signatures identified from in vitro experiments could contribute to individualized medicine.

Differential valine metabolism in adipose tissue of low and high fat-oxidizing obese subjects.

Differences in fat metabolism are of importance in relation to energy balance. Low fat-oxidizers (LFO) are thought to be more prone for developing obesity. We studied whether LFO have different fasting adipose tissue (AT) protein profiles than high fat-oxidizers (HFO). Six LFO and six HFO subjects were selected from an obese group (n = 99, body mass index>30 kg/m(2) ) taking part in a multi-center study (Nutrient-Gene interaction in human obesity) based on the postprandial fat oxidation capacity after a high fat load. AT protein profiles were studied by 2-DE. Differential proteins were clustered with MAPPfinder according to their function. Protein profiles of purified blood cells and adipocytes served to confine the comparison to adipocyte-specific proteins in AT profiles of LFO and HFO subjects. LFO had increased mitochondrial ROS scavengers possibly related to long-chain unsaturated fatty acid-induced increases in mitochondrial ROS-production. Carbohydrate oxidation seemed to be reduced since expression of several proteins from the glycolysis pathway was lower in LFO. Up-regulation of the valine catabolism at the level of methylmalonate-semialdehyde dehydrogenase appeared to be (part of) the compensatory mechanism. In conclusion, the fasting AT protein profile of LFO and HFO differ at the level of ROS scavenging, the glycolysis pathway and valine metabolism.

Biologically relevant effects of mRNA amplification on gene expression profiles.

BACKGROUND: Gene expression microarray technology permits the analysis of global gene expression profiles. The amount of sample needed limits the use of small excision biopsies and/or needle biopsies from human or animal tissues. Linear amplification techniques have been developed to increase the amount of sample derived cDNA. These amplified samples can be hybridised on microarrays. However, little information is available whether microarrays based on amplified and unamplified material yield comparable results. In the present study we compared microarray data obtained from amplified mRNA derived from biopsies of rat cardiac left ventricle and non-amplified mRNA derived from the same organ. Biopsies were linearly amplified to acquire enough material for a microarray experiment. Both amplified and unamplified samples were hybridized to the Rat Expression Set 230 Array of Affymetrix. RESULTS: Analysis of the microarray data showed that unamplified material of two different left ventricles had 99.6% identical gene expression. Gene expression patterns of two biopsies obtained from the same parental organ were 96.3% identical. Similarly, gene expression pattern of two biopsies from dissimilar organs were 92.8% identical to each other.Twenty-one percent of reporters called present in parental left ventricular tissue disappeared after amplification in the biopsies. Those reporters were predominantly seen in the low intensity range. Sequence analysis showed that reporters that disappeared after amplification had a GC-content of 53.7+/-4.0%, while reporters called present in biopsy- and whole LV-samples had an average GC content of 47.8+/-5.5% (P <0.001). Those reporters were also predicted to form significantly more (0.76+/-0.07 versus 0.38+/-0.1) and longer (9.4+/-0.3 versus 8.4+/-0.4) hairpins as compared to representative control reporters present before and after amplification. CONCLUSION: This study establishes that the gene expression profile obtained after amplification of mRNA of left ventricular biopsies is representative for the whole left ventricle of the rat heart. However, specific gene transcripts present in parental tissues were undetectable in the minute left ventricular biopsies. Transcripts that were lost due to the amplification process were not randomly distributed, but had higher GC-content and hairpins in the sequence and were mainly found in the lower intensity range which includes many transcription factors from specific signalling pathways.

Up-regulation of CD36/FAT in preadipocytes in familial combined hyperlipidemia.

Familial combined hyperlipidemia (FCHL) shows many features of the metabolic syndrome. The strong genetic component makes it an excellent model to study the genetic background of metabolic syndrome and insulin resistance. Adipose tissue is believed to contribute to, or even underlie, the FCHL phenotype and is an interesting target tissue for gene expression studies. However, interpretation of adipose tissue gene expression experiments is complex since expression differences cannot only arise as a direct consequence of a genetic trait, but may also reflect an adaptation to metabolic influences at the cellular level. In the present study, we measured gene expression levels in cultured primary human preadipocytes from FCHL and control subjects. Since isolated preadipocytes were allowed to replicate for weeks under standardized conditions, the contribution of previous metabolic influences is rather small whereas genetic defects are preserved and expressed in vitro. The main finding was up-regulation of CD36/FAT in FCHL preadipocytes, confirmed in two independent groups of subjects, and a concomitant increase in CD36/FAT-mediated fatty acid uptake. CD36/FAT overexpression has previously been shown to be associated with other insulin-resistant states. The present data suggest that CD36/FAT overexpression in FCHL occurs very early in adipocyte differentiation and may be of genetic origin.

Erythrocyte antioxidant defense response against cigarette smoking in humans--the glutathione S-transferase vulnerability.

Cigarette smoking leads to uptake of a multitude of reactive chemicals including many electrophiles and may also give rise to oxidative stress. Human red blood cells are important targets for electrophilic and oxidant foreign compounds. We investigated the oxidative stress in erythrocytes upon cigarette smoking, and the response of antioxidant defense system against it. With this aim, simultaneous determination of erythrocyte superoxide dismutase (SOD), selenium dependent glutathione peroxidase (Se-GPx), catalase (CAT), glutathione S-transferase (GST) activities and plasma levels of thiobarbituric acid reactive substances (TBARS), and the degree of erythrocyte membrane lipid peroxidation (EMLP) were carried out in blood samples of smokers and their controls. Plasma TBARS levels and EMLP in smokers were significantly higher than the control levels (p < 0.01 and p < 0.005, respectively). SOD activity was diminished in smokers compared to nonsmoker controls (p < 0.005). Erythrocyte Se-GPx activity was also found significantly diminished in smokers (p < 0.005), while plasma Se-GPx activity was not changed. We observed that erythrocyte CAT activity was not different in smokers compared to nonsmoker controls. We found that the erythrocyte GST activity is significantly lower in young adult smokers (3.03 +/- 0.18 U/mg protein; mean +/- SEM; n = 46) than in nonsmoking contemporaries (3.98 +/- 0.26 U/mg protein; mean +/- SEM; n = 41). Together with previously reported data, it can be concluded that the decrease in GST activity leads to extra GST synthesis during erythrocyte proliferation. The same data were also analyzed for the sex differences. The statistically significant differences remained the same between nonsmoker and smoker females. Only EMLP degree and SOD activity were significantly different between nonsmoker and smoker males; however, when compared the parameters between male and female nonsmokers, GST activity was found to be significantly higher in females than that of males.

Gene profiling in atherosclerosis reveals a key role for small inducible cytokines: validation using a novel monocyte chemoattractant protein monoclonal antibody.

BACKGROUND: Pathological aspects of atherosclerosis are well described, but gene profiles during atherosclerotic plaque progression are largely unidentified. METHODS AND RESULTS: Microarray analysis was performed on mRNA of aortic arches of ApoE-/- mice fed normal chow (NC group) or Western-type diet (WD group) for 3, 4.5, and 6 months. Of 10 176 reporters, 387 were differentially (>2x) expressed in at least 1 group compared with a common reference (ApoE-/-, 3- month NC group). The number of differentially expressed genes increased during plaque progression. Time-related expression clustering and functional grouping of differentially expressed genes suggested important functions for genes involved in inflammation (especially the small inducible cytokines monocyte chemoattractant protein [MCP]-1, MCP-5, macrophage inflammatory protein [MIP]-1alpha, MIP-1beta, MIP-2, and fractalkine) and matrix degradation (cathepsin-S, matrix metalloproteinase-2/12). Validation experiments focused on the gene cluster of small inducible cytokines. Real-time polymerase chain reaction revealed a plaque progression-dependent increase in mRNA levels of MCP-1, MCP-5, MIP-1alpha, and MIP-1beta. ELISA for MCP-1 and MCP-5 showed similar results. Immunohistochemistry for MCP-1, MCP-5, and MIP-1alpha located their expression to plaque macrophages. An inhibiting antibody for MCP-1 and MCP-5 (11K2) was designed and administered to ApoE-/- mice for 12 weeks starting at the age of 5 or 17 weeks. 11K2 treatment reduced plaque area and macrophage and CD45+ cell content and increased collagen content, thereby inducing a stable plaque phenotype. CONCLUSIONS: Gene profiling of atherosclerotic plaque progression in ApoE-/- mice revealed upregulation of the gene cluster of small inducible cytokines. Further expression and in vivo validation studies showed that this gene cluster mediates plaque progression and stability.

A systematic review of large scale and heterogeneous gene array data in heart failure.

Microarray analysis has become a widely available tool for the generation of gene expression data on a genomic scale. Since the studies with similar protocols are growing, it has become necessary to systematically revise the large body of literature to decipher the gene expression data. In this review, we analyzed and critically discussed the database presented from 14 published studies that showed the gene expression profile in heart failure (HF) using microarray as a primary tool. After comparing the diverse database from these studies, we explain the protein translational, matri-cellular, immunological and fibrosis-related mechanisms in HF. In addition to previously annotated genes, we analyzed two differentially expressed expressed sequence tags (ESTs) (KIAA0152 and Suppressor of G(Two) allele of the suppressor of kinetochore protein-1, SGT1) in HF and showed how bio-informatic analysis of ESTs can lead to the identification of novel pathways active in HF. We have also discussed the new publicly accessible tools that link the gene expression data to gene ontogeny (GO) and functionality. Finally, we have systematically revised the chromosomal localization of the genes that are specifically up-regulated in HF. We have thus spotted chromosome 1, 2, 11 and 12 as the chromosomal hotspots of HF. This methodical approach will simplify the existing concepts on the evolution and progression of HF and lead us toward the development of newer diagnostic and therapeutic tools. Although modeled to HF, this approach should be of broader scientific interest to elaborate multiple genes and complex pathways.

Thrombospondin-2 is essential for myocardial matrix integrity: increased expression identifies failure-prone cardiac hypertrophy.

Cardiac hypertrophy can lead to heart failure (HF), but it is unpredictable which hypertrophied myocardium will progress to HF. We surmised that apart from hypertrophy-related genes, failure-related genes are expressed before the onset of failure, permitting molecular prediction of HF. Hearts from hypertensive homozygous renin-overexpressing (Ren-2) rats that had progressed to early HF were compared by microarray analysis to Ren-2 rats that had remained compensated. To identify which HF-related genes preceded failure, cardiac biopsy specimens were taken during compensated hypertrophy and we then monitored whether the rat progressed to HF or remained compensated. Among 48 genes overexpressed in failing hearts, we focused on thrombospondin-2 (TSP2). TSP2 was selectively overexpressed only in biopsy specimens from rats that later progressed to HF. Moreover, expression of TSP2 was increased in human hypertrophied hearts with decreased (0.19+/-0.01) versus normal ejection fraction (0.11+/-0.03 [arbitrary units]; P<0.05). Angiotensin II induced fatal cardiac rupture in 70% of TSP2 knockout mice, with cardiac failure in the surviving mice; this was not seen in wild-type mice. In TSP2 knockout mice, angiotensin II increased matrix metalloproteinase (MMP)-2 and MMP-9 activity by 120% and 390% compared with wild-type mice (P<0.05). In conclusion, we identify TSP2 as a crucial regulator of the integrity of the cardiac matrix that is necessary for the myocardium to cope with increased loading and that may function by its regulation of MMP activity. This suggests that expression of TSP2 marks an early-stage molecular program that is activated uniquely in hypertrophied hearts that are prone to fail.

Role of the Rhoa/Rho kinase system in flow-related remodeling of rat mesenteric small arteries in vivo.

In small arteries, a chronic blood flow reduction leads to inward hypotrophic remodeling, while a chronic blood flow elevation induces outward hypertrophic remodeling. The RhoA/Rho kinase system was shown to be modulated by shear stress, and to be involved in other kinds of vascular remodeling. The aim of this study was to investigate the role of RhoA/Rho kinase in flow-related small artery remodeling. Rat mesenteric small arteries were subjected to flow-modifying surgery. After 1, 2, 4, 16, and 32 days, the animals were sacrificed and small arteries were harvested. Messenger RNA was isolated and amplified. Using cDNA microarray analysis, the differential expression of >14,000 genes was analyzed, part of which was confirmed by RT-PCR. In vivo treatment with fasudil (3 mg/kg/day s.c.) was used to test the effect of Rho kinase inhibition. The main findings are that: (1) blood flow alteration modified the expression of approximately 5% of the genes by >2-fold, (2) flow reduction downregulated many RhoA-related cytoskeletal markers of smooth muscle cell phenotype, (3) many RhoA-related genes were rapidly (<1 day) regulated and (4) fasudil treatment potentiated the inward hypotrophic remodeling in response to chronically reduced flow. These results indicate the importance of the RhoA/Rho kinase system in flow-related small artery remodeling.

Effect of vitamin E on glutathione-dependent enzymes.

Reactive oxygen species and various electrophiles are involved in the etiology of diseases varying from cancer to cardiovascular and pulmonary disorders. The human body is protected against damaging effects of these compounds by a wide variety of systems. An important line of defense is formed by antioxidants. Vitamin E (consisting of various forms of tocopherols and tocotrienols) is an important fat-soluble, chain-breaking antioxidant. Besides working as an antioxidant, this compound possesses other functions with possible physiological relevance. The glutathione-dependent enzymes form another line of defense. Two important enzymes in this class are the free radical reductase and glutathione S-transferases (GSTs). The GSTs are a family of phase II detoxification enzymes. They can catalyze glutathione conjugation with various electrophiles. In most cases the electrophiles are detoxified by this conjugation, but in some cases the electrophiles are activated. Antioxidants do not act in isolation but form an intricate network. It is, for instance, known that vitamin E, together with glutathione (GSH) and a membrane-bound heat labile GSH-dependent factor, presumably an enzyme, can prevent damaging effects of reactive oxygen species on polyunsaturated fatty acids in biomembranes (lipid peroxidation). This manuscript reviews the interaction between the two defense systems, vitamin E and glutathione-dependent enzymes. On the simplest level, antioxidants such as vitamin E have protective effects on glutathione-dependent enzymes; however, we will see that reality is somewhat more complicated.