Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Muscle-specific regulation of right ventricular transcriptional responses to chronic hypoxia-induced hypertrophy by the muscle ring finger-1 (MuRF1) ubiquitin ligase in mice

BMC Medical Genetics volume 19, Article number: 175 (2018) Cite this article

Abstract

Background

We recently identified a role for the muscle-specific ubiquitin ligase MuRF1 in right-sided heart failure secondary to pulmonary hypertension induced by chronic hypoxia (CH). MuRF1−/− mice exposed to CH are resistant to right ventricular (RV) dysfunction whereas MuRF1 Tg + mice exhibit impaired function indicative of heart failure. The present study was undertaken to understand the underlying transcriptional alterations in the RV of MuRF1−/− and MuRF1 Tg + mice.

Methods

Microarray analysis was performed on RNA isolated from the RV of MuRF1−/−, MuRF1 Tg+, and wild-type control mice exposed to CH.

Results

MuRF1−/− RV differentially expressed 590 genes in response to CH. Analysis of the top 66 genes (> 2-fold or < − 2-fold) revealed significant associations with oxidoreductase, transcription regulation, and transmembrane component annotations. The significant genes had promoters enriched for HOXD12, HOXC13, and RREB-1 protein transcription factor binding sites. MuRF1 Tg + RV differentially expressed 150 genes in response to CH. Analysis of the top 45 genes (> 3-fold or < − 3-fold) revealed significant associations with oxidoreductase-metabolic, glycoprotein-transmembrane-integral proteins, and alternative splicing/splice variant annotations. The significant genes were enriched for promoters with ZIC1 protein transcription factor binding sites.

Conclusions

The differentially expressed genes in MuRF1−/− and MuRF1 Tg + RV after CH have common functional annotations related to oxidoreductase (including antioxidant) and transmembrane component functions. Moreover, the functionally-enhanced MuRF1−/− hearts regulate genes related to transcription, homeobox proteins, and kinases/phosphorylation. These studies also reveal potential indirect effects of MuRF1 through regulating Rreb-1, and they reveal mechanisms by which MuRF1 may transcriptionally regulate anti-oxidant systems in the face of right heart failure.

Peer Review reports

Background

Right ventricular heart failure is the most frequent cause of death in patients with pulmonary arterial hypertrophy (PAH) [1, 2]. In general, chronic hypoxia (CH) in the lungs induces an elevation in pulmonary artery pressures, increasing the afterload in the right ventricle. In response, right ventricular (RV) remodeling results, possibly including maladaptive hypertrophy. Recent reviews have outlined the molecular and cellular drivers of right ventricular remodeling in chronic pressure overload, including the induction of cardiomyocyte cell death, increases in collagen synthesis and decreases in collagen degradation, decreases in capillary density and a metabolic shift to glycolysis [3]. Rho-kinase [4], adrenergic signaling [5] and altered angiogenesis have been reported to play roles in the pathogenesis of RV remodeling at the molecular level. Maladaptive RV hypertrophy characteristically exhibits a shift to aerobic glycolysis mediated by FOXO1 [6] and c-Myc [7]. Beyond this, little is currently known about the molecular regulation of these underlying mechanisms that occur during PAH-induced RV heart failure.

The muscle-specific ubiquitin ligase Muscle Ring Finger-1 (MuRF1) regulates both physiological and pathological hypertrophy and metabolism by interacting with specific transcription factors that mediate these processes. MuRF1 interacts with and inhibits serum response factor [8], which is upregulated in pathological pressure overload-induced hypertrophy. As such, MuRF1−/− mice exhibited exaggerated cardiac hypertrophy after transaortic constriction in vivo [8]. MuRF1 interacts with the c-Jun transcription factor in a proteasome-dependent manner to inhibit downstream AP1 activity, implicated in inhibiting IGF-1 signaling [9]. Consistent with this, MuRF1−/− mice voluntarily ran farther and faster without training, improving even more after training compared to wild-type mice [9]. The mechanisms by which cardiomyocyte MuRF1 inhibits the PPARalpha and TRalpha nuclear receptors appear to involve multi-mono-ubiquitination, resulting in changes in metabolism in vivo [9, 10].

Given transcriptional regulation by MuRF1 in cardiomyocytes, we recently hypothesized a role for MuRF1 in the pathophysiology of right ventricle heart failure [11]. We exposed MuRF1−/− and cardiomyocyte-specific MuRF1 transgenic (MuRF1 Tg+) mice to CH, resulting in an elevated pulmonary artery pressure/resistance and measured functional and morphological responses [11]. We found that MuRF1−/− mice exhibited greater RV growth compared to wild-type mice but were resistant to CH-induced changes in ejection fraction. As such, the overall output in MuRF1−/− hearts provided significantly increased perfusion to skeletal muscle compared to sibling wild-type mice [11]. In contrast, MuRF1 Tg + mice exhibited similar patterns in RV weight compared to wild-type mice in response to CH, but there was a reduced ejection fraction, indicative of a maladaptive dilated phenotype (heart failure), consistent with the significant decreases in skeletal muscle perfusion [11]. The present study was undertaken to understand the underlying transcriptional alterations that occur during CH in the right ventricle of MuRF1−/− and MuRF1 Tg + mice using microarray analysis.

Methods

MuRF1−/− and αMHC transgenic MuRF1 mice and the hypoxia-induced pulmonary hypertension model of right sided heart failure

MuRF1−/− and MuRF1 Tg mice and their strain-match wild-type controls (previously described [8, 12]) were shipped to the University of New Mexico. Following a week-long quarantine, mice maintained in standard shoebox cages were moved to either control or normobaric hypoxia chambers with continuous access to food and water ad libitum throughout. The hypoxia chamber was set at 10.0% oxygen (partial pressure of oxygen approximately 65 mmHg in Albuquerque, NM) and monitored both by a digital feedback control system (Biospherix, Parish, NY) as well as by a secondary O2/CO2 monitor (O2Cap, OxiGraf, Inc.; Mountain View, CA), as recently described [11]. Hypoxia exposure lasted 3 weeks with twice-weekly cage changes and standard 12 h light:dark cycle. Procedures were conducted under full isoflurane anesthesia to minimize or eliminate risk of pain and discomfort. Mice were anesthetized with isoflurane and were euthanized by exsanguination, followed by rapid removal of the heart into ice-cold phosphate-buffered saline. The RV free wall was dissected away from the rest of the heart, as previously described [13]. Mice were bred at the University of North Carolina at Chapel Hill and all the hypoxia procedures were conducted at the University of New Mexico with full approval of the UNC-Chapel Hill and University of New Mexico Institutional Animal Care and Use Committees and were carried out in compliance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.

Microarray assay and statistical analysis

Microarrays were performed on RNA isolated from RV tissue from animals challenged with hypoxia. At least three independent samples were assayed from each of the four experimental groups MuRF1−/− (N = 5), its strain matched wild-type MuRF1+/+ (N = 6), MuRF1 Tg + (N = 3), and its strain-matched wild type (N = 3). Gene expression assay was performed using Agilent Whole Mouse 4 × 44 k multiplex format oligo arrays (014868, Agilent Technologies, Santa Clara, CA) following the Agilent 1-color microarray-based gene expression analysis protocol. Data were obtained using Agilent Feature Extraction software (version 12), using the 1-color defaults for all parameters. The Agilent Feature Extraction software performed error remodeling, adjusting for additive and multiplicative noise. To determine differentially expressed genes between groups, an ANOVA with multiple test correction (FDR q-value) was performed using Partek Genomics Suite software (version 6.6). Only one gene for MuRF1−/− and 6 genes for MuRF1 Tg + passed the FDR q-value < 0.05 cutoff; therefore, the set of differentially expressed genes was expanded to include genes from the ANOVA analysis with an unadjusted p-value < 0.01 (see Additional file 1: Table S1 and Additional file 2: Table S2). To minimize inclusion of false positives, genes with fold changes greater than 2 and less than − 2 (for MuRF1−/− vs. MuRF1+/+) or greater than 3 and less than − 3 (for MuRF1 Tg + vs. MuRF1+/+) were selected for downstream bioinformatic analyses. The differential expression of a subset of selected genes was confirmed by RT-PCR. The MAIME-compliant data were submitted to Gene Expression Omnibus and assigned the accession #GSE82345 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=sfavmocqtrupnan&acc=GSE82345, currently embargoed/not publicly available).

Bioinformatics resources and analytic tools used and retrieval of public microarray data

TRANSFAC® FMatch was used to analyze transcription factors (https://portal.biobase-international.com, GeneExplan GmbH). The TRANSFAC MATRIX TABLE, Release 2016.2 was used as the matrix library. To minimize false-positive matches, the matrix file vertebrate_non_redundant_minFP was used for transcription factor binding sites prediction, using cut-off criteria to minimize false-positives, and a p-value threshold set at 0.01. The significant gene sets for the MuRF1−/− (66 unique genes) and MuRF1 Tg + (45 unique genes) were used in the STRING analysis. Duke Gather Literature Net. Genes were analyzed at http://changlab.uth.tmc.edu/gather/ via Literature Net. STRING database. The STRING database (STRING Consortium, Swiss Institute of Bioinformatics, Zurich, Switzerland) of known and predicted protein-protein interactions (https://string-db.org/) includes direct (physical) and indirect (functional) associations aggregated from other (primary) databases [14,15,16]. The significant gene sets for the MuRF1−/− (66 unique genes) and MuRF1 Tg + (45 unique genes) were used in the STRING analysis. The MuRF1−/− significant networks contained 26 nodes (vs. expected 3 nodes), 74 edges, average node degree 5.69, with a clustering coefficient of 0.735. The MuRF1 Tg + significant networks contained 26 nodes, 75 edges (28 expected), average node degree 5.77 and a clustering coefficient of 0.696. The Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 Functional Annotation Bioinformatics database was used to analyze differentially expressed genes (https://david.ncifcrf.gov/) [17, 18]. Analysis of gene lists (gene symbols) were uploaded using the gene accession conversion tool and analyzed using Disease (OMIM_Disease), Functional Categories (COG_Ontology, UP_Keywords, UP_SEQ_FEATURE) Gene Ontology (GOTERM_BP_DIRECT, GOTERM_CC_DIRECT, GOTERM_MF_DIRECT), Pathways (BBID, BIOCARTA, KEGG PATHWAY) and Protein Domains (INTERPRO, PIR_SUPERFAMILY, SMART). Microarray data retrieval and experimental design. The SuHx model of severe angio-obliterative pulmonary artery hypertension and RV failure was used for male Sprague-Dawley rats receiving a single subcutaneous injection of SU5416 (20 mg/kg) prior to 4 weeks of hypoxia [2]. Microarray analysis (Agilent whole rat genome 4 × 44 k microarray slide, Agilent Technologies, Wilmington, DE) was performed on snap-frozen right heart tissues [2]. Study datasets were retrieved from the publicly available NCBI Gene Expression Omnibus (GEO Dataset_GSE42579). Microarray replicates included 12 hypoxia and 8 normoxia from 6 biological replicates per group (hypoxia, normoxia). The GEO2R tool was used to download the ID, adj. p value, p value, t value, B value, log fold change (logFC), gene symbol, gene title, and gene ID of the differentially expressed genes for further analysis and comparison.

RT qPCR analysis of mRNA expression

Total RNA was isolated using RNeasy Mini kits (QIAGEN, Valencia, CA) according to the manufacturer’s protocols. Determination of gene expression was performed using either a two-step RT-qPCR reaction [8, 19] or a one-step RT-qPCR reaction, as previously described [20]. For the two-step reaction, cDNA was made using iScript cDNA Synthesis kit, according to the manufacturer’s protocol (BioRad, Hercules, CA, Cat# 1708890). One microliter of the resulting cDNA product was then amplified in a 20-μl final volume using the TaqMan® Universal PCR Master Mix. Each reaction included 1 μl of mouse specific TaqMan probes for MuRF1 (Mm01185221_m1), Atp2b2 (Mm00437640_m1), Elovl7 (Mm00512434_m1), Gbp3 (Mm00497606_m1), Cxcl9 (Mm00434946), Myl1 (Mm00659043_m1), Casq1 (Mm00486725_m1), GAPDH (glyceraldehyde-3-phosphate dehydrogenase, Probe Mm99999915_g1), or 18S (Hs99999901_s1) in triplicate (ThermoFisher Scientific, Inc., Waltham, MA). The relative expression of mRNA (ΔΔCT algorithm) was determined using GAPDH or 18S as an internal sample loading control.

Statistical analysis

For qRT-PCR analyses, Student’s t-test was performed to determine significant differences between groups. All statistics were performed using GraphPad Prism (La Jolla, CA, Version 7.01) with p values < 0.05 considered statistically significant.

Results

To understand MuRF1’s role in the pathophysiology of CH-induced right heart failure, we performed microarray analysis on right ventricle tissue from MuRF1−/− and MuRF1 Tg + mice (and sibling wild-type controls) after three weeks of CH. Analysis of MuRF1−/− and MuRF1+/+ RV tissue by microarray analysis identified 590 genes to be significantly different from wild-type controls (ANOVA, p < 0.01) (Fig. 1a, Additional file 1: Table S1). We identified the top 35 unique genes increased (> 2-fold) and top 31 unique genes decreased (<− 2-fold) (Fig. 1b) and analyzed them for common transcription factors that may be responsible for the differential regulation of the MuRF1−/− mice using TRANSFAC® FMatch Analysis. We identified three (3) transcription factors that were enriched in significantly altered gene promoters in MuRF1−/− hypoxia-challenged RV compared to wild-type controls (left, Fig. 1c). These transcription factors are listed by the number of the binding sites found in the transcription factor sequences (sites/sequence ratio) reflecting how common these sites were for each transcription factor. The Homeobox D12 (HOXD12), Homeobox C13 (HOXC13), and the Ras Responsive Element Binding Protein 1 (RREB-1) were enriched 10-, 6.5- and 6.5-fold, respectively (Fig. 1c). Of the 66 significant genes used in the FMatch analysis, 15, 8 and 2 of these genes were associated with HOXD12, HOXC13 and RREB-1, respectively (Fig. 1c).

Fig. 1
figure 1

Differential gene expression of MuRF1−/− right ventricle after three weeks CH exposure. a Summary of how differentially expressed genes in hypoxia-challenged MuRF1−/− RV tissue compared to MuRF1+/+ were identified. b Top 35 increased unique genes (> 2-fold) and top 31 decreased unique genes (<− 2-fold) compared to MuRF1+/+ hearts. Asterisk (*) designates same/related sequence for duplicate gene name. c Identification of transcription factors with binding sites of top 35 increased and top 31 decreased genes using TRANSFAC® FMatch. The number of genes in dataset and enrichment (sites/sequence ratio) for each transcription factor is calculated to the right. d STRING analysis of differentially expressed genes for known interactions and e STRING analysis pathway ID statistics. MuRF1−/−: N = 5 biological replicates; MuRF1+/+: N = 6 biological replicates

Full size image

To understand how the differentially expressed genes in the MuRF1−/− right ventricle are known to interact, we performed a STRING analysis to identify relationships that had been reported between these gene products in the literature (Fig. 1d). We identified significantly more interactions than would be expected by chance. Specifically, the network had 26 modes and 74 edges compared to the expected 3 edges due to chance. The functional enrichments identified by STRING Gene Ontology (GO) included the regulation of transcription, gene expression and RNA metabolic processes (Fig. 1e). Overall, HOXD12 had the most interactions (20), while MTF2 and MLL4 had 9 and INTS12 had 6, demonstrating the importance of HOXD12 among these proteins based on our current knowledge of these genes.

To determine the functional significance of the 66 differentially expressed genes from the MuRF1−/− right ventricle in more detail, we performed functional annotations using both STRING and DAVID. STRING analysis of functional categories identified regulation of transcription, gene expression, and RNA metabolism as the significant pathways with low false discovery rates (Fig. 1e). DAVID analysis detailed these functions, identifying five functional clusters across several databases (Fig. 2a). Functional annotation identified 30 genes (Fig. 2a, far left) that exhibited reduced enrichment in five clusters, including Homeobox-related genes (cluster 1), genes involved in oxidoreductase (cluster 2), genes involved in phosphorylation/kinase activity (cluster 3), genes involved in the regulation of transcription (cluster 4), and genes involved in the transmembrane components (cluster 5) (Fig. 2a). We validated the differential expression of several genes in the MuRF1−/− right ventricle after CH, including MuRF1, Atp2b2, Elovl7, Gbp3 and Cxcl9 mRNA by RT-qPCR analysis (Fig. 2b and Additional file 3: Figure S1). Together, these results demonstrate that 3 weeks of CH-induced changes in the MuRF1−/− right ventricle to alter functions associated with homeobox proteins (associated with gene expression), oxidoreductase processes, kinases/phosphorylation, transcriptional regulation and transmembrane components.

Fig. 2
figure 2

Functional annotation of differentially expressed genes in MuRF1−/− right ventricle after three weeks CH exposure. a. Hypercluster analysis of differentially expressed MuRF1−/− genes in the right ventricle after hypoxia challenge using the DAVID v6.8 Functional Annotation Bioinformatics database. Differentially expressed gene symbols are on the left, functional categories along the top (database source in parentheses). The cluster fold enrichment is color-coded (key in lower right). b. Reverse transcriptase quantitative PCR analysis of MuRF1−/− hearts of differentially expressed genes MuRF1, Atp2b2, and Elovl7 mRNA in the right ventricle

Full size image

Analysis of MuRF1 Tg + and MuRF1+/+ RV tissue by microarray analysis identified 150 genes as significantly different from wild-type controls (ANOVA, p < 0.01) (Fig. 3a, Additional file 2: Table S2). We identified the top 27 unique genes increased (> 3-fold) and top 18 unique genes decreased (<− 3-fold) (Fig. 3b) so that we could identify the common transcription factors that may be responsible for the differential regulation of MuRF1Tg + mice using TRANSFAC® FMatch Analysis. Analysis of these 45 genes identified that they were enriched with promoter regions recognized by the ZIC1 (ZInc finger protein of the Cerebellum (Zic) Family Member 1) protein transcription factor in MuRF1Tg + − hypoxia-challenged RV compared to wild-type controls (Fig. 4a). The 2.6-fold enrichment in ZIC1 protein promoter sites (p = 6.76549e-03) was based on the significantly altered genes Jakmip3 (Janus Kinase and Microtubule Interacting Protein 3), Cyp2b13 (Cytochrome P450 Family 2 Subfamily B Member 6), Gucd1 (Guanylyl Cyclase Domain Containing 1), Kansl1l (KAT8 Regulatory NSL Complex Subunit 1 Like) and Rab3gap1 (RAB3 GTPase Activating Protein Catalytic Subunit 1) genes (Fig. 4a).

Fig. 3
figure 3

Differential gene expression of MuRF1 Tg + right ventricle after three weeks CH exposure. a. Summary of how differentially expressed genes in hypoxia-challenged MuRF1Tg + right ventricles were identified. b. Top 27 increased unique genes (> 3-fold) and top 18 decreased unique genes (<− 3-fold) compared to wild-typeMuRF1Tg+ RV tissue. Asterisk (*) designates transcript variant for duplicate gene name. MuRF1 Tg+: N = 3; strain-matched wild type: N = 3

Full size image
Fig. 4
figure 4

Functional annotation of differentially expressed genes in MuRF1 Tg + right ventricles after three weeks CH exposure. a. Identification of ZIC1 protein binding sites based on binding sites in 5 of the differentially expressed genes in MuRF1 Tg + heart (Jakmip3, Cyp2b13, Gucd1, Kansl1l, Rab3gap1). The number of genes in dataset and enrichment (sites/sequence ratio) for each transcription factor is calculated to the right. b. STRING analysis of the top 27 genes increased (> 3-fold) and top 18 genes decreased (<− 3-fold) listed in Fig. 3b. c. Functional annotation of differentially expressed MuRF1 Tg + right ventricle by STRING analysis

Full size image

To identify how the genes that were differentially expressed in the MuRF1 Tg + right ventricle compared to wild-type controls after three weeks of CH, we performed a STRING analysis to identify relationships that had been reported between these gene products in the literature (Fig. 4b). Specifically, the network had 26 modes and 75 edges compared to the expected 28 edges due to chance (PPI enrichment p-value = 3.27e-13). The functional enrichments identified by STRING Gene Ontology (GO) included chromatin modification, histone modification and regulation of gene expression among other related functions (Fig. 4c). Overall, ZIC1 protein had 6 interactions across the STRING network and GL1 and GL13 had 9 and 13 interactions with the network, respectively, demonstrating the importance of ZIC1 within these proteins in the literature.

We analyzed the 45 differentially expressed genes from the MuRF1 Tg + right ventricle after 3 weeks of CH using DAVID analysis to understand the biological relevance of the gene set. We found five clusters, including oxidoreductase-metabolic (cluster 1), glycoprotein-transmembrane-integral proteins (cluster 2), cell membrane/cell junction (cluster 3), signal peptide, glycosylation (cluster 4) and alternative splicing/splice variant (cluster 5) (Fig. 5a). We validated the differential expression of several genes in the MuRF1 Tg + right ventricle after CH, including MuRF1, Casq1 and Myl1 mRNA by RT-qPCR analysis (Fig. 5b). Analysis of the 45 differentially expressed genes from the MuRF1 Tg + right ventricle for their relationships in the literature [21] identified 12 groups, including links with PPAR binding protein (PPARBP), regulators of G-protein signaling (RGS5, RGS8, RGS12, RGS14, RGS20) and v-ros UR2 sarcoma virus oncogene (ROS1) (Additional file 4: Figure S2).

Fig. 5
figure 5

Functional annotation of differentially expressed genes in MuRF1 Tg + hearts compared to wild-type controls. a. Hypercluster analysis of differentially expressed MuRF1−/− genes in the right ventricle after hypoxia challenge using the DAVID v6.8 Functional Annotation Bioinformatics database. Differentially expressed gene symbols are on the left, functional categories along the top (database source in parentheses). The cluster fold enrichment is color-coded (key in lower left). b. Reverse transcriptase quantitative PCR analysis of MuRF1−/− hearts of differentially expressed genes MuRF1, Casq1, and Myl1 mRNA in the right ventricle

Full size image

Comparison of the functional categories identified in the CH-challenged MuRF1−/− right ventricles (vs. wild type) (Fig. 2a) and the MuRF1 Tg + right ventricles (vs. wild type) (Fig. 5a) demonstrates overlapping and unique functional regulation in each model (Table 1). Both the MuRF1−/− and MuRF1 Tg + hearts have differential regulation of oxidoreductase functions and transmembrane component proteins (Table 1). The MuRF1−/− model uniquely regulated genes associated with homeobox proteins, kinase/phosphorylation, and transcriptional regulation (Table 1, left column) while the MuRF1 Tg + model uniquely regulated genes associated with signal peptide/glycosylation and alternative splicing (Table 1, right column).

Table 1 Summary of DAVID functional annotation categories of significant genes. From Figure 3A (MuRF1-/-) and Figure 5A (MuRF1 Tg+) based off differentially expressed genes in Figure 2A and Figure 5A. Matching colors represent related categories found in hearts lacking MuRF1 and hearts with increased cardiomyocyte MuRF1 post-hypoxia challenge compared to wildtype controls
Full size table

To give context to our identification of MuRF1-regulated genes in the right ventricle after 3 weeks of CH, we identified a previous study that performed microarray analysis of the right ventricle after CH [2]. We obtained the publicly available data from NCBI GEO to identify the genes differentially expressed in right heart failure induced by CH (Fig. 6). Of the 25,241 defined genes (Fig. 6a), we identified the top 50 genes increased (> 4.9-fold) and top 50 genes decreased (<− 3.9-fold) compared to controls (Additional file 5: Figure S3) and we analyzed these genes using TRANSFAC® FMatch Analysis to identify transcription factors that bind promoters of the significantly altered hypoxia-challenged RV (Fig. 6b). The most enriched transcription factor protein sites included the CREB binding protein (CPBP, 33.46-fold enriched), Zinc Finger Protein 333 (ZNF333, 20-fold enriched), Sex Determining Region Y (SRY) (18.32-fold enriched), and the Nuclear Factor of Activated T-Cells 2 (NFATc2, 11.66-fold enriched) (Fig. 6b). Analysis of the 100 differentially expressed genes from the CH-challenged right ventricle for their relationships in the literature [21] identified 8 groups, including links with PPAR binding protein (PPARBP), regulators of G-protein signaling (RGS5, RGS8, RGS12, RGS14, RGS18, RGS20) and v-ros UR2 sarcoma virus oncogene (ROS1) (Additional file 6: Figure S4). This analysis provides a context of the transcriptional regulation of the RV remodeling that occurs in response to pulmonary hypertension induced right heart failure, demonstrating the transcriptional control involved (Fig. 6b) and the associated connections with proteins in the literature (Additional file 6: Figure S4). When we compared the TRANSFAC Match Analysis identified in the CH-challenged MuRF1−/− right ventricles (vs. wild type) (Fig. 1c) with the MuRF1 Tg + right ventricles (vs. wild type) (Fig. 4a) and hypoxia-induced right heart failure (Fig. 6b), Rreb-1 was identified as a common transcription factor (Table 2). Finally, considerable overlap was seen in Literature Net analysis of the MuRF1−/− significant genes after CH (vs. wild type) (Additional file 4: Figure S2) and the significant genes found in wild-type rat right hearts after CH (Additional file 6: Figure S4), with almost complete overlap (PPARBP, RGS proteins, and ROS1).

Fig. 6
figure 6

Differentially expressed genes in a rat model of CH-induced pulmonary hypertension and right heart failure. a. Summary of how differentially expressed genes in CH-induced rat heart failure resulted in top 50 increased genes (> 4.9-fold) and top 50 decreased genes (< 3.9-fold) compared to normoxia using data from Drake, et al., 2013 [2]. Genes summarized in Additional file 4: Figure S2. b. Identification of transcription factors with binding sites of top 50 increased and top 50 decreased genes. The number of genes in the dataset and enrichment (sites/sequence ratio) for each transcription factor is calculated to the right

Full size image
Table 2 Summary of transcription factors identified in the MuRF1-/- and MuRF1 Tg+ right ventricles (vs. wildtype) by TRANSFAC® Match Analysis. From Figure 2C (MuRF1-/-) and Figure 4C (MuRF1 Tg+) based off differentially expressed genes in Figure 2B and Figure 4B. Matching colors represent related transcription factors identified in hearts lacking MuRF1 and hearts with increased cardiomyocyte MuRF1 post-hypoxia challenge compared to wildtype controls
Full size table

Discussion

We recently identified that MuRF1−/− mice were resistant to weeks of CH-induced changes in RV ejection fraction and volume (systolic and diastolic) and exhibited significantly increased skeletal muscle perfusion, consistent with better cardiac function compared to wild-type mice [11]. In the present study, we analyzed the MuRF1−/− and MuRF1 Tg + right ventricle transcriptome by microarray analysis to determine the role of cardiomyocyte MuRF1 in response to CH and resulting pulmonary hypertension. The differentially expressed transcripts in the MuRF1−/− hearts (vs. wild type) were enriched by up to 10-fold for HOXD12, HOXC13 and RREB-1 protein promoter regions (Fig. 1c). The role of HOXD12 and HOXC13 in the heart or muscle has not been previously reported to our knowledge; however, their role in transcription and gene expression (as are other homeodomain proteins) is well-annotated in gene ontology (Fig. 1e). Our analysis of a publicly available dataset from a rat model of CH exposure also identified the Rreb-1 (enriched 4.25-fold) as a transcription factor associated with differentially regulated genes, among other transcription factors (e.g., SRY, CRX, MAF, SMAD4) (Fig. 6b). There are several differences in the experimental design of this previous study that used Sprague-Dawley rats (vs. mice in the present study) and pulmonary artery hypertension induced by a single subcutaneous SU5416 (20 mg/kg) injection and 4 weeks of hypoxia (vs. 3 weeks of hypoxia only) [2]. RREB-1 transcription factor binds specifically to the RAS-response elements (RRE) of gene promoters and acts as a coregulatory factor for ninety-one (91) transcription factors in various tissues [22]. In the heart, the RREB-1 protein notably binds with MEF-2, POU3F2, SRF, PPARalpha/RXRalpha, GR and E2F. While the role of RREB-1 has not been investigated directly in the heart, these seven transcription factors have well-defined roles in heart failure. MEF-2 and SRF, in particular, are known to play critical roles in pathological hypertrophy and are part of the “fetal gene” response in pressure overload-induced heart failure [23]. PPARalpha/RXRalpha have critical roles in regulating fatty acid metabolism in cardiomyocytes and in the context of heart failure have been proposed to play roles in the shift in substrate utilization in heart failure [24].

These findings are notable given what is known about cardiac MuRF1 activity according to the literature. In studies of MuRF1−/− mice challenged with pressure overload-induced hypertrophy via transaortic constriction, MuRF1−/− mice experienced enhanced hypertrophy (cardiomyocyte growth) and decreased SRF-regulated gene expression (fetal genes) [8]. These studies found that MuRF1 immunoprecipitated specifically with SRF and directly inhibited cardiomyocyte SRF activity [8]. Similarly, cardiomyocyte MuRF1 has been reported to interact directly with PPARalpha to inhibit its activity, with MuRF1−/− hearts having a five-fold increase in PPARalpha activity [25]. In these studies, increasing MuRF1 expression demonstrated enhancement of PPARalpha nuclear export, dependent upon specific lysines adjacent to nuclear export sequences found to be mono-ubiquitinated [25]. In additional to MuRF1’s regulation of SRF and PPARalpha, cardiac MuRF1 and MuRF2 have been reported to support transcription factor E2F1 in microarray studies, confirmed by chromatin IP studies in the context of developmental cardiac growth [26, 27]. Since pulmonary (RV) hypertension has numerous parallels with left ventricular pressure overload-induced hypertrophy, the present study may suggest that MuRF1’s regulatory effects on SRF may contribute to the phenotype observed. Taken together, the present study suggests that MuRF1’s regulation of SRF, PPARalpha, and E2F1 may be due to MuRF1’s broader regulation of cotranscription factors such as Rreb-1, given Rreb-1’s established connection with SRF, PPARalpha and E2F1 in other tissues [22].

Beyond the evidence linking MuRF1 to the regulation of transcription, the differentially expressed genes in the MuRF1−/− right ventricle after CH were related to oxidoreductase functions. Oxidoreductases are enzymes that catalyze oxidoreduction reactions and have important roles in aerobic and anaerobic metabolism and the utilization of NAD+/NADH. In the heart, catalase is an oxidoreductase that acts on H2O2, breaking it down to harmless hydrogen and water [28, 29]. In the present study, the MuRF1−/− right ventricle had several differentially expressed genes related to the oxidoreductase system (Adh1, Aasdh, Pam, Xdh) (Fig. 2a). Alcohol dehydrogenase 1 (Adh1) has been identified in the heart to have important antioxidant roles in collaboration with catalase [30]. Expression of aminoadipate-semialdehyde dehydrogenase (Aasdh) is highest in the brain, followed by heart and skeletal muscle [31]. In the heart, the membrane bound peptidylglycine alpha-amidating monooxygenase (Pam) has been implicated in the shaping of atrial secretory vesicles related to the secretion of pro-ANP [32]. Cell injury from hyperoxia was associated with increased superoxide radicals, which has been associated with xanthine dehydrogenase (Xdh) activity in heart [33]. MuRF1−/− right ventricles also had differentially expressed genes (Araf, Ip6K3, Nmrk2) associated with kinase activity (Fig. 2a). The role of raf proteins, including A-Raf Proto- Oncogene, Serine/Threonine Kinase (Araf), has been described in cardiac hypertrophy and cardiomyocyte survival, with antiapoptotic activity independent of MKK and ERK activities [34, 35]. Inositol hexakisphosphate kinase 3 (Ip6K3) has been linked to metabolism and lifespan in mice [36], while nicotinamide riboside kinase 2 (Nmrk2) has critical roles in NAD metabolism in mammalian cells [37].

Cardiomyocyte MuRF1 has been indirectly linked to oxidoreductase functions in previous studies. Cardiomyocyte-specific MuRF1 Tg + mice were challenged with an acute cardiac ischemia reperfusion injury both ex vivo and in vitro and demonstrated considerable resistance to injury and functional defects [38]. One mechanism that could explain this protection against acute hypoxia is MuRF1’s inhibition of JNK signaling, mediated by MuRF1’s polyubiquitination of phosphorylated c-Jun (subsequently degraded) [38]. Subsequent studies showed that MuRF1 Tg + heart mitochondria generated significantly less (51–67% less) ROS measured by H2O2 emission at baseline compared to sibling wild-type control mitochondria [39]. How MuRF1 reduces H2O2 emission is not known, but the results from the current study may lend some insights.

We utilized the MuRF1 Tg + mice in present study because of the limited redundancy that is known to occur between MuRF1 and other MuRF family members (MuRF2, MuRF3). This redundancy was initially suggested by overlap of some binding partners (e.g., cardiac troponin I [cTnI]) but not others (e.g., myosin light chain 2) [40]. The partial overlap between cardiac MuRF1, MuRF2 and MuRF3 was also seen in nontargeted metabolomics studies of baseline heart tissues [41]. As seen in previous studies of MuRF1 Tg + mice challenged with pressure overload-induced hypertrophy [12], CH resulted in significant decreases in RV weight and ejection fraction, indicating a maladaptive dilated phenotype (heart failure) [11]. In contrast to the MuRF1−/− mice challenged with CH that resulted in improved skeletal muscle perfusion (vs. wild-type mice), MuRF1 Tg + mice exhibited significant decreases in skeletal muscle perfusion [11]. Identification of transcription factors enriched in the differentially expressed genes in the MuRF1 Tg + right ventricle in the present study identified promoter regions that the ZIC1 protein binds (Fig. 4a). The ZIC1 protein is a transcriptional activator described primarily in the development and maturation of the cerebellum [42]. It has also been reported to have a role in shear flow mechanotransduction in osteocytes [43]. Functionally it retains nuclear GLI1 and GLI3 in the cytoplasm by binding the minimal GLI-consensus sequence 5-TGGGTGGTC-3 [44]. Mutations in Zic gene family members resulted in congenital heart defects in both humans and mice [45]. The ZIC1 protein has also been shown to act as a protein transcription factor upstream Pax3 in the induction of embryonic neural crest [46]. While little is known about the ZIC1 protein in the heart, it has been linked to the homeobox domain protein HOXA1 protein. Specifically, Hoxa1 null embryos revealed their regulation of genes in cardiac NC precursors, including Hnf1b, Foxd3, and Zic1 [47]. Sumoylation regulates the nuclear localization and function of closely related ZIC3 protein family members [48, 49].

Hypoxia-induced RV hypertrophy in cardiomyocyte-specific MuRF1 Tg + mice induced gene expression changes functionally linked to oxidoreductase regulation. Of the five differentially expressed genes in the MuRF1 Tg + post-CH related to oxidoreductase regulation (Fig. 5a, far left column), only one had clear roles in the heart. The Delta 4-Desaturase Sphingolipid 2 (Degs2) is an enzyme involved in de novo ceramide synthesis. One study reported that ethanol exposure upregulated Degs in rats and altered lipids profiles that support beneficial cardioprotective and neuroprotective effects of moderate ethanol consumption [50]. The role of DEGS2 in the heart is unknown. The Cytochrome P450 Family 27 Subfamily B Member 1 (Cyp27b1), also known as 25-Hydroxyvitamin D3 1-alpha-hydroxylase, is responsible for the hydroxylation of calcifediol to calcitriol (the bioactive form of vitamin D). Variants in Diacylglycerol Kinase Kappa (Dgkk) have been described in association with hypospadias [51,52,53,54], but no functions in the heart have been reported. The role of Cyp27b1 in the heart is unknown. Detailed function of the Cytochrome P450 Family 2 Subfamily B Member 13 (Cyp2b13) and Arachidonate 12-Lipoxygenase, 12R Type (Alox12b) were unavailable beyond the functional categorization by DAVID (Fig. 5a).

Hypoxia-induced RV hypertrophy in cardiomyocyte-specific MuRF1 Tg + mice also induced gene expression changes functionally linked to alternative splicing/splice variants. Of the 16 differentially expressed genes in the MuRF1 Tg + post-CH related to alternative splicing (Fig. 5a, far right column), four stood out (Casq1, Myl1, Erbb2ip, Pitx2). Calsequestrin 1 (Casq1) functions as a calcium buffer in the heart with a role in calcium release by interactions with the ryanodine receptor, having important implications for contractility [55,56,57]. Myosin Light Chain 1 (Myl1) is a GATA4 target gene related to IGF-1 signaling in the context of pathological hypertrophy stimuli [58]. The Erbb2 Interacting Protein (Erbb2ip) is present in relatively high levels in the heart and has been found to be down-regulated in heart biopsies from human failing hearts, with a role in attenuated ERK signaling to negatively modulate cardiac hypertrophy experimentally in mice [59]. Mutations in the Paired Like Homeodomain 2 (Pitx2) gene in the heart have been seen in the Tetralogy of Fallot and atrial fibrillation in patients [60,61,62,63,64]. Interestingly, Pitx2 was the only gene in common between the set of 66 differentially expressed genes in the MuRF1−/− hearts and the set of 45 differentially expressed genes in the MuRF1 Tg + hearts. Pitx2 expression was induced in RV in mice lacking MuRF1 and was repressed in response to MuRF1 overexpression, suggesting an exquisite sensitivity of this gene to levels of MuRF1.The Fatty Acid Synthase (FASN aka C93) gene was associated with lipid metabolism reprogramming in arrhythmogenic RV cardiomyopathy [65], while hypoxia has been associated with FASN-dependent free fatty acid production [66,67,68]. Lastly, the RAB3 GTPase Activating Protein Catalytic Subunit 1 (Rab3gap1) has been associated with cardiovascular risk (total cholesterol and HDL) by GWAS association in Framingham Heart Studies [69], with Rab3gap1 loci having associations with an increased risk of SCD [70]. The remaining differentially expressed genes in MuRF1 Tg + hearts after CH exposure associated with alternative splicing (Ptpru, Svopl, Syt16, Ccp110, Jakmip3, Kansl1l) have not been functionally defined in heart and/or cardiac disease to our knowledge.

Conclusions

Analysis of the top differentially expressed transcripts in the MuRF1−/− right ventricle after CH revealed genes involved in the regulation of transcription and oxidoreductase functions. Bioinformatic analysis suggested that the differentially expressed genes in MuRF1−/− right ventricle tissue were linked to HOX (homeobox) transcription factors and Rreb-1. Analysis of the top differentially expressed transcripts in the cardiomyocyte-specific MuRF1 Tg + − right ventricle after CH revealed genes involved in the regulation of oxidoreductase-metabolic, glycoprotein-transmembrane and alternative splicing functions. Together, the MuRF1−/− and MuRF1 Tg + results have common functional annotations related to oxidoreductase (including antioxidant functions) and transmembrane component functions (Table 1). These studies reveal potential indirect effects of MuRF1 on known substrates such as SRF and PPARalpha by regulation of Rreb-1. The results also give insight into potential mechanisms by which MuRF1 regulates antioxidant systems by transcriptionally regulating genes related to oxidoreductase functions.

Abbreviations

ATF-2:

Activating Transcription Factor 2

Atp2b2:

ATPase Plasma Membrane Ca2+ Transporting 2

BEN:

Alpha-helical domain found in diverse proteins including BANP/SMAR1, E5R, and NAC1 protein)

Casq1:

Calsequestrin 1

CDX-2:

Caudal Type Homeobox 2; CPBP, (aka Kruppel Like Factor 6, KLF6)

CRX:

Cone-Rod Homeobox

DR11:

(Aka HLA-DRB1, Major Histocompatibility Complex, Class II, DR Beta 1)

Elovl7:

ELOVL Fatty Acid Elongase 7

FAC1:

Fetal Alz-50 reactive Clone 1 (aka BPTF, Bromodomain PHD Finger Transcription Factor)

GKLF:

Aka KLF4, Kruppel Like Factor 4

GLI:

GLIoma-associated ongogene homolog 1 (Zinc Finger Protein)

HDX:

Highly Divergent Homeobox

HMGIY:

High Mobility Group Protein HMG-I/HMG-Y

HOX:

Homeobox

ING4:

Inhibitor Of Growth Family Member 4

MuRF1 Tg + :

MuRF1 transgenic

MuRF1:

Muscle ring finger-1

MuRF2:

Muscle ring finger-2

Myl1:

Myosin Light Chain 1

NFATc2:

Nuclear Factor of Activated T-cells 2

Rreb-1:

Ras Responsive Element Binding Protein 1

SRT:

Sirtuin 1

Zic1:

Zinc finger protein of the Cerebellum (Zic) Family Member 1

ZNFF333:

Zinc Finger Protein 333

References

  1. Zornoff LA, Skali H, Pfeffer MA, St John Sutton M, Rouleau JL, Lamas GA, Plappert T, Rouleau JR, Moye LA, Lewis SJ, et al. Right ventricular dysfunction and risk of heart failure and mortality after myocardial infarction. J Am Coll Cardiol. 2002;39(9):1450–5.

    Article  Google Scholar 

  2. Drake JI, Gomez-Arroyo J, Dumur CI, Kraskauskas D, Natarajan R, Bogaard HJ, Fawcett P, Voelkel NF. Chronic carvedilol treatment partially reverses the right ventricular failure transcriptional profile in experimental pulmonary hypertension. Physiol Genomics. 2013;45(12):449–61.

    Article  CAS  Google Scholar 

  3. Guihaire J, Noly PE, Schrepfer S, Mercier O. Advancing knowledge of right ventricular pathophysiology in chronic pressure overload: insights from experimental studies. Arch Cardiovasc Dis. 2015;108(10):519–29.

  4. Ikeda S, Satoh K, Kikuchi N, Miyata S, Suzuki K, Omura J, Shimizu T, Kobayashi K, Kobayashi K, Fukumoto Y, et al. Crucial role of rho-kinase in pressure overload-induced right ventricular hypertrophy and dysfunction in mice. Arterioscler Thromb Vasc Biol. 2014;34(6):1260–71.

    Article  CAS  Google Scholar 

  5. Ryan JJ, Archer SL. The right ventricle in pulmonary arterial hypertension: disorders of metabolism, angiogenesis and adrenergic signaling in right ventricular failure. Circ Res. 2014;115(1):176–88.

    Article  CAS  Google Scholar 

  6. Piao L, Sidhu VK, Fang YH, Ryan JJ, Parikh KS, Hong Z, Toth PT, Morrow E, Kutty S, Lopaschuk GD, et al. FOXO1-mediated upregulation of pyruvate dehydrogenase kinase-4 (PDK4) decreases glucose oxidation and impairs right ventricular function in pulmonary hypertension: therapeutic benefits of dichloroacetate. J Mol Med (Berl). 2013;91(3):333–46.

    Article  CAS  Google Scholar 

  7. Piao L, Fang YH, Parikh K, Ryan JJ, Toth PT, Archer SL. Cardiac glutaminolysis: a maladaptive cancer metabolism pathway in the right ventricle in pulmonary hypertension. J Mol Med (Berl). 2013;91(10):1185–97.

    Article  CAS  Google Scholar 

  8. Willis MS, Ike C, Li L, Wang DZ, Glass DJ, Patterson C. Muscle ring finger 1, but not muscle ring finger 2, regulates cardiac hypertrophy in vivo. Circ Res. 2007;100(4):456–9.

    Article  CAS  Google Scholar 

  9. Wadosky KM, Rodriguez JE, Hite RL, Min JN, Walton BL, Willis MS. Muscle RING finger-1 attenuates IGF-I-dependent cardiomyocyte hypertrophy by inhibiting JNK signaling. Am J Physiol Endocrinol Metab. 2014;306(7):E723–39.

    Article  CAS  Google Scholar 

  10. Wadosky KM, Berthiaume JM, Tang W, Zungu M, Portman MA, Gerdes AM, Willis MS. MuRF1 mono-ubiquitinates TRalpha to inhibit T3-induced cardiac hypertrophy in vivo. J Mol Endocrinol. 2016;56(3):273–90.

    Article  CAS  Google Scholar 

  11. Campen MJ, Paffett ML, Colombo ES, Lucas SN, Anderson T, Nysus M, Norenberg JP, Gershman B, Hesterman J, Hoppin J, et al. Muscle RING finger-1 promotes a maladaptive phenotype in chronic hypoxia-induced right ventricular remodeling. PLoS One. 2014;9(5):e97084.

    Article  Google Scholar 

  12. Willis MS, Schisler JC, Li L, Rodriguez JE, Hilliard EG, Charles PC, Patterson C. Cardiac muscle ring finger-1 increases susceptibility to heart failure in vivo. Circ Res. 2009;105(1):80–8.

    Article  CAS  Google Scholar 

  13. Campen MJ, Shimoda LA, O'Donnell CP. Acute and chronic cardiovascular effects of intermittent hypoxia in C57BL/6J mice. J Appl Physiol (1985). 2005;99(5):2028–35.

    Article  CAS  Google Scholar 

  14. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.

    Article  CAS  Google Scholar 

  15. Franceschini A, Lin J, von Mering C, Jensen LJ. SVD-phy: improved prediction of protein functional associations through singular value decomposition of phylogenetic profiles. Bioinformatics. 2016;32(7):1085–7.

    Article  CAS  Google Scholar 

  16. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43(Database issue):D447–52.

    Article  CAS  Google Scholar 

  17. Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  18. Huang d W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.

    Article  Google Scholar 

  19. Willis MS, Rojas M, Li L, Selzman CH, Tang RH, Stansfield WE, Rodriguez JE, Glass DJ, Patterson C. Muscle ring finger 1 mediates cardiac atrophy in vivo. Am J Physiol Heart Circ Physiol. 2009;296(4):H997–H1006.

    Article  CAS  Google Scholar 

  20. Oakley RH, Ren R, Cruz-Topete D, Bird GS, Myers PH, Boyle MC, Schneider MD, Willis MS, Cidlowski JA. Essential role of stress hormone signaling in cardiomyocytes for the prevention of heart disease. Proc Natl Acad Sci U S A. 2013;110(42):17035–40.

    Article  CAS  Google Scholar 

  21. Jenssen TK, Laegreid A, Komorowski J, Hovig E. A literature network of human genes for high-throughput analysis of gene expression. Nat Genet. 2001;28(1):21–8.

  22. Tissue-specific Gene Expression and Regulation (TiGER). http://bioinfo.wilmer.jhu.edu/tiger/db_tf/RREB-1-index.html. Accessed 19 Feb 2018.

  23. Frey N, Olson EN. Cardiac hypertrophy: the good, the bad, and the ugly. Annu Rev Physiol. 2003;65:45–79.

    Article  CAS  Google Scholar 

  24. Neubauer S. The failing heart--an engine out of fuel. N Engl J Med. 2007;356(11):1140–51.

    Article  Google Scholar 

  25. Rodriguez JE, Liao JY, He J, Schisler JC, Newgard CB, Drujan D, Glass DJ, Frederick CB, Yoder BC, Lalush DS, et al. The ubiquitin ligase MuRF1 regulates PPARalpha activity in the heart by enhancing nuclear export via monoubiquitination. Mol Cell Endocrinol. 2015;413:36–48.

    Article  CAS  Google Scholar 

  26. Willis MS, Wadosky KM, Rodriguez JE, Schisler JC, Lockyer P, Hilliard EG, Glass DJ, Patterson C. Muscle ring finger 1 and muscle ring finger 2 are necessary but functionally redundant during developmental cardiac growth and regulate E2F1-mediated gene expression in vivo. Cell Biochem Funct. 2014;32(1):39–50.

    Article  CAS  Google Scholar 

  27. Parry TL, Desai G, Schisler JC, Li L, Quintana MT, Stanley N, Lockyer P, Patterson C, Willis MS. Fenofibrate unexpectedly induces cardiac hypertrophy in mice lacking MuRF1. Cardiovasc Pathol. 2016;25(2):127–40.

    Article  CAS  Google Scholar 

  28. Williams J. The decomposition of hydrogen peroxide by liver catalase. J Gen Physiol. 1928;11(4):309–37.

    Article  CAS  Google Scholar 

  29. Wang X, Tao Y, Huang Y, Zhan K, Xue M, Wang Y, Ruan D, Liang Y, Huang X, Lin J, et al. Catalase ameliorates diabetes-induced cardiac injury through reduced p65/RelA- mediated transcription of BECN1. J Cell Mol Med. 2017;21(12):3420–34.

    Article  CAS  Google Scholar 

  30. Harris C, Wang SW, Lauchu JJ, Hansen JM. Methanol metabolism and embryotoxicity in rat and mouse conceptuses: comparisons of alcohol dehydrogenase (ADH1), formaldehyde dehydrogenase (ADH3), and catalase. Reprod Toxicol. 2003;17(3):349–57.

    Article  CAS  Google Scholar 

  31. Praphanphoj V, Sacksteder KA, Gould SJ, Thomas GH, Geraghty MT. Identification of the alpha-aminoadipic semialdehyde dehydrogenase-phosphopantetheinyl transferase gene, the human ortholog of the yeast LYS5 gene. Mol Genet Metab. 2001;72(4):336–42.

    Article  CAS  Google Scholar 

  32. Labrador V, Brun C, Konig S, Roatti A, Baertschi AJ. Peptidyl-glycine alpha-amidating monooxygenase targeting and shaping of atrial secretory vesicles: inhibition by mutated N-terminal ProANP and PBA. Circ Res. 2004;95(12):e98–109.

    Article  CAS  Google Scholar 

  33. Elsayed NM, Tierney DF. Hyperoxia and xanthine dehydrogenase/oxidase activities in rat lung and heart. Arch Biochem Biophys. 1989;273(2):281–6.

    Article  CAS  Google Scholar 

  34. Guo W, Hao B, Wang Q, Lu Y, Yue J. Requirement of B-Raf, C-Raf, and A-Raf for the growth and survival of mouse embryonic stem cells. Exp Cell Res. 2013;319(18):2801–11.

    Article  CAS  Google Scholar 

  35. Muslin AJ. Role of raf proteins in cardiac hypertrophy and cardiomyocyte survival. Trends Cardiovasc Med. 2005;15(6):225–9.

    Article  CAS  Google Scholar 

  36. Moritoh Y, Oka M, Yasuhara Y, Hozumi H, Iwachidow K, Fuse H, Tozawa R. Inositol Hexakisphosphate kinase 3 regulates metabolism and lifespan in mice. Sci Rep. 2016;6:32072.

    Article  CAS  Google Scholar 

  37. Ratajczak J, Joffraud M, Trammell SA, Ras R, Canela N, Boutant M, Kulkarni SS, Rodrigues M, Redpath P, Migaud ME, et al. NRK1 controls nicotinamide mononucleotide and nicotinamide riboside metabolism in mammalian cells. Nat Commun. 2016;7:13103.

    Article  CAS  Google Scholar 

  38. Li HH, Du J, Fan YN, Zhang ML, Liu DP, Li L, Lockyer P, Kang EY, Patterson C, Willis MS. The ubiquitin ligase MuRF1 protects against cardiac ischemia/reperfusion injury by its proteasome-dependent degradation of phospho-c-Jun. Am J Pathol. 2011;178(3):1043–58.

    Article  CAS  Google Scholar 

  39. Mattox TA, Young ME, Rubel CE, Spaniel C, Rodriguez JE, Grevengoed TJ, Gautel M, Xu Z, Anderson EJ, Willis MS. MuRF1 activity is present in cardiac mitochondria and regulates reactive oxygen species production in vivo. J Bioenerg Biomembr. 2014;46(3):173–87.

    Article  CAS  Google Scholar 

  40. Witt SH, Granzier H, Witt CC, Labeit S. MURF-1 and MURF-2 target a specific subset of myofibrillar proteins redundantly: towards understanding MURF-dependent muscle ubiquitination. J Mol Biol. 2005;350(4):713–22.

    Article  CAS  Google Scholar 

  41. Banerjee R, He J, Spaniel C, Quintana MT, Wang Z, Bain J, Newgard CB, Muehlbauer MJ, Willis MS. Non-targeted metabolomics analysis of cardiac muscle ring Finger-1 (MuRF1), MuRF2, and MuRF3 in vivo reveals novel and redundant metabolic changes. Metabolomics. 2015;11(2):312–22.

    Article  CAS  Google Scholar 

  42. Frank CL, Liu F, Wijayatunge R, Song L, Biegler MT, Yang MG, Vockley CM, Safi A, Gersbach CA, Crawford GE, et al. Regulation of chromatin accessibility and Zic binding at enhancers in the developing cerebellum. Nat Neurosci. 2015;18(5):647–56.

    Article  CAS  Google Scholar 

  43. Kalogeropoulos M, Varanasi SS, Olstad OK, Sanderson P, Gautvik VT, Reppe S, Francis RM, Gautvik KM, Birch MA, Datta HK. Zic1 transcription factor in bone: neural developmental protein regulates mechanotransduction in osteocytes. FASEB J. 2010;24(8):2893–903.

    Article  CAS  Google Scholar 

  44. Mizugishi K, Aruga J, Nakata K, Mikoshiba K. Molecular properties of Zic proteins as transcriptional regulators and their relationship to GLI proteins. J Biol Chem. 2001;276(3):2180–218.

    Article  CAS  Google Scholar 

  45. Srivastava D. Genetic assembly of the heart: implications for congenital heart disease. Annu Rev Physiol. 2001;63:451–69.

    Article  CAS  Google Scholar 

  46. Li B, Kuriyama S, Moreno M, Mayor R. The posteriorizing gene Gbx2 is a direct target of Wnt signalling and the earliest factor in neural crest induction. Dev. 2009;136(19):3267–78.

    Article  CAS  Google Scholar 

  47. Makki N, Capecchi MR. Cardiovascular defects in a mouse model of HOXA1 syndrome. Hum Mol Genet. 2012;21(1):26–31.

    Article  Google Scholar 

  48. Chen L, Ma Y, Qian L, Wang J. Sumoylation regulates nuclear localization and function of zinc finger transcription factor ZIC3. Biochim Biophys Acta. 2013;1833(12):2725–33.

    Article  CAS  Google Scholar 

  49. Bedard JE, Purnell JD, Ware SM. Nuclear import and export signals are essential for proper cellular trafficking and function of ZIC3. Hum Mol Genet. 2007;16(2):187–98.

    Article  CAS  Google Scholar 

  50. Godfrey J, Jeanguenin L, Castro N, Olney JJ, Dudley J, Pipkin J, Walls SM, Wang W, Herr DR, Harris GL, et al. Chronic voluntary ethanol consumption induces favorable ceramide profiles in selectively bred alcohol-preferring (P) rats. PLoS One. 2015;10(9):e0139012.

    Article  Google Scholar 

  51. Xie H, Lin XL, Zhang S, Yu L, Li XX, Huang YC, Lyu YQ, Chen HT, Xu J, Chen F. Asian J Androl. 2018;20(1):85–9. https://doi.org/10.4103/aja.aja_13_17.

  52. Zhuang LK, Wu M, Ye WJ, Liu YD. Single nucleotide polymorphisms of the DGKK gene and hypospadias in Chinese children. Zhonghua Nan Ke Xue. 2014;20(11):991–4.

    CAS  PubMed  Google Scholar 

  53. Ma Q, Tang Y, Lin H, Xu M, Xu G, Fang X, Chen J, Song Z, Li Z, Shi Y, et al. Diacylglycerol kinase kappa (DGKK) variants and hypospadias in Han Chinese: association and meta-analysis. BJU Int. 2015;116(4):634–40.

    Article  CAS  Google Scholar 

  54. van der Zanden LF, van Rooij IA, Feitz WF, Knight J, Donders AR, Renkema KY, Bongers EM, Vermeulen SH, Kiemeney LA, Veltman JA, et al. Common variants in DGKK are strongly associated with risk of hypospadias. Nat Genet. 2011;43(1):48–50.

    Article  CAS  Google Scholar 

  55. Gaburjakova M, Bal NC, Gaburjakova J, Periasamy M. Functional interaction between calsequestrin and ryanodine receptor in the heart. Cell Mol Life Sci. 2013;70(16):2935–45.

    Article  CAS  Google Scholar 

  56. Zarain-Herzberg A, Estrada-Aviles R, Fragoso-Medina J. Regulation of sarco (endo) plasmic reticulum Ca2+-ATPase and calsequestrin gene expression in the heart. Can J Physiol Pharmacol. 2012;90(8):1017–28.

    Article  CAS  Google Scholar 

  57. Beard NA, Wei L, Dulhunty AF. Ca (2+) signaling in striated muscle: the elusive roles of triadin, junctin, and calsequestrin. Eur Biophys J. 2009;39(1):27–36.

    Article  CAS  Google Scholar 

  58. Bisping E, Ikeda S, Sedej M, Wakula P, McMullen JR, Tarnavski O, Sedej S, Izumo S, Pu WT, Pieske B. Transcription factor GATA4 is activated but not required for insulin-like growth factor 1 (IGF1)-induced cardiac hypertrophy. J Biol Chem. 2012;287(13):9827–34.

    Article  CAS  Google Scholar 

  59. Rachmin I, Tshori S, Smith Y, Oppenheim A, Marchetto S, Kay G, Foo RS, Dagan N, Golomb E, Gilon D, et al. Erbin is a negative modulator of cardiac hypertrophy. Proc Natl Acad Sci U S A. 2014;111(16):5902–7.

    Article  CAS  Google Scholar 

  60. Vande Perre P, Zazo Seco C, Patat O, Bouneau L, Vigouroux A, Bourgeois D, El Hout S, Chassaing N, Calvas P. Eur J Med Genet. 2018;61(2):72–8.

  61. Mora C, Serzanti M, Giacomelli A, Beltramone S, Marchina E, Bertini V, Piovani G, Refsgaard L, Olesen MS, Cortellini V, et al. Generation of induced pluripotent stem cells (iPSC) from an atrial fibrillation patient carrying a PITX2 p.M200V mutation. Stem Cell Res. 2017;24:8–11.

    Article  CAS  Google Scholar 

  62. Weng LC, Lunetta KL, Muller-Nurasyid M, Smith AV, Theriault S, Weeke PE, Barnard J, Bis JC, Lyytikainen LP, Kleber ME, et al. Genetic interactions with age, sex, body mass index, and hypertension in relation to atrial fibrillation: the AFGen consortium. Sci Rep. 2017;7(1):11303.

    Article  Google Scholar 

  63. Ocana OH, Coskun H, Minguillon C, Murawala P, Tanaka EM, Galceran J, Munoz-Chapuli R, Nieto MA. A right-handed signalling pathway drives heart looping in vertebrates. Nature. 2017;549(7670):86–90.

    Article  CAS  Google Scholar 

  64. Lee JY, Kim TH, Yang PS, Lim HE, Choi EK, Shim J, Shin E, Uhm JS, Kim JS, Joung B, et al. Korean atrial fibrillation network genome-wide association study for early-onset atrial fibrillation identifies novel susceptibility loci. Eur Heart J. 2017;38(34):2586–94.

    Article  Google Scholar 

  65. Chen L, Yang F, Chen X, Rao M, Zhang NN, Chen K, Deng H, Song JP, Hu SS. Comprehensive myocardial Proteogenomics profiling reveals C/EBPalpha as the key factor in the lipid storage of ARVC. J Proteome Res. 2017;16(8):2863–76.

    Article  CAS  Google Scholar 

  66. Capitanio D, Fania C, Torretta E, Vigano A, Moriggi M, Bravata V, Caretti A, Levett DZH, Grocott MPW, Samaja M, et al. TCA cycle rewiring fosters metabolic adaptation to oxygen restriction in skeletal muscle from rodents and humans. Sci Rep. 2017;7(1):9723.

    Article  Google Scholar 

  67. Sun W, Kato H, Kitajima S, Lee KL, Gradin K, Okamoto T, Poelllinger L. Interaction between von Hippel-Lindau protein and fatty acid synthase modulates hypoxia target gene expression. Sci Rep. 2017;7(1):7190.

    Article  Google Scholar 

  68. Lee HJ, Jung YH, Choi GE, Ko SH, Lee SJ, Lee SH, Han HJ. BNIP3 induction by hypoxia stimulates FASN-dependent free fatty acid production enhancing therapeutic potential of umbilical cord blood-derived human mesenchymal stem cells. Redox Biol. 2017;13:426–43.

    Article  CAS  Google Scholar 

  69. Ma L, Yang J, Runesha HB, Tanaka T, Ferrucci L, Bandinelli S, Da Y. Genome-wide association analysis of total cholesterol and high-density lipoprotein cholesterol levels using the Framingham heart study data. BMC Med Genet. 2010;11:55.

    Article  CAS  Google Scholar 

  70. Huertas-Vazquez A, Nelson CP, Guo X, Reinier K, Uy-Evanado A, Teodorescu C, Ayala J, Jerger K, Chugh H, Wtccc, et al. Novel loci associated with increased risk of sudden cardiac death in the context of coronary artery disease. PLoS One. 2013;8(4):e59905.

Download references

Acknowledgements

We thank Dr. Kevin Gerrish and Liwen Liu of the Molecular Genomics core laboratory (NIEHS) for assistance with microarray analysis.

Funding

This work is supported by the National Institutes of Health (NIEHS to R.H.O. and J.A.C., NIH R01ES014639 to M.J.C. and R01HL104129 to M.S.W.), the Jefferson-Pilot Corporation (Fellowship to M.S.W.), and the Leducq Foundation Transatlantic Networks of Excellence (to M.W.).

Availability of data and materials

The MAIME-compliant datasets analyzed in the current study were submitted to the NCBI Gene Expression Omnibus (GEO) and assigned the accession #GSE82345 (publicly available). A comparison dataset was retrieved from the publicly available NCBI GEO (GEO Dataset_GSE42579). Analysis output including probeset ID, genebank accession #, p-values, FDR q-value, and fold-chage for annotated genes are included in the Additional file 1: Table S1 and Additional file 2: Table S2 for publication with the manuscript.

Author information

Authors and Affiliations

  1. Department of Health and Human Services, Signal Transduction Laboratory, National Institute of Environmental Health Sciences, National Institutes of Health, Research Triangle Park, NC, USA

    Robert H. Oakley & John A. Cidlowski

  2. Department of Pharmaceutical Sciences, College of Pharmacy, University of New Mexico Health Sciences Center, Albuquerque, NM, USA

    Matthew J. Campen & Michael L. Paffett

  3. McAllister Heart Institute, University of North Carolina, Chapel Hill, NC, USA

    Xin Chen, Traci L. Parry, Carolyn Hillhouse & Monte S. Willis

  4. Department of Surgery, University of North Carolina, Chapel Hill, NC, USA

    Zhongjing Wang

  5. Department of Pathology & Laboratory Medicine, Indiana University School of Medicine, 635 Barnhill Drive, Van Nuys MS 5067, Indianapolis, IN, 46202, USA

    Monte S. Willis

  6. Krannert Institute of Cardiology and Division of Cardiology, Department of Internal Medicine, Indiana University School of Medicine, Indianapolis, IN, USA

    Monte S. Willis

Authors
  1. Robert H. Oakley
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Matthew J. Campen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Michael L. Paffett
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Xin Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Zhongjing Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Traci L. Parry
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Carolyn Hillhouse
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. John A. Cidlowski
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Monte S. Willis
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

RHO, MJC, JAC, MSW: Conceived and designed the experiments. RHO, MJC, JAC, and MSW analyzed and interpreted the data. MJC, MLP, CH, TLP, XC, ZW, MSW performed experiments, interpreted results, prepared results for publication, wrote draft methods and interpretations. All authors contributed to the final version of the manuscript and read/approved the final submitted manuscript.

Corresponding author

Correspondence to Monte S. Willis.

Ethics declarations

Ethics approval

Mice were bred at the University of North Carolina at Chapel Hill and all the hypoxia procedures were conducted at the University of New Mexico with full approval of the UNC-Chapel Hill and University of New Mexico Institutional Animal Care and Use Committees and were carried out in compliance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Differentially expressed genes in MuRF1−/− mice challenged with hypoxia (XLSX 59 kb)

Additional file 2:

Table S2. Differentially expressed genes in MuRF1 Tg + mice challenged with hypoxia (XLSX 22 kb)

Additional file 3:

Figure S1. Validation of differential expressed genes in MuRF1−/− hearts. Reverse transcriptase quantitative PCR analysis of MuRF1−/− hearts of differentially expressed genes Gbp3 and Cxcl9 mRNA in the right ventricle. (PDF 26 kb)

Additional file 4:

Figure S2. Literature Net analysis of differentially expressed genes in MuRF1 Tg + right ventricles after chronic hypoxia challenge compared to wildtype controls. Analysis of the top 28 genes increased (> 3 fold) and top 18 genes decreased (<− 3 fold) listed in Fig. 3b using Literature Net on Duke Gather (http://changlab.uth.tmc.edu/gather/gather.py). (PDF 24 kb)

Additional file 5:

Figure S3. Top differentially expressed genes by microarray analysis of right ventricle tissue in rats challenged with pulmonary hypertension-induced right heart failure. A. Top 50 genes increased compared to normoxia controls. B. Top 50 genes decreased compared to normoxia controls. Data downloaded from Gene Expression Omnibus (GEO) GEO2R interface including the 8 normoxia right ventricle replicates (6 biological replicates) and 12 hypoxia right ventricle replicates (6 biological replicates). 43,480 data points were downloaded, with 25.241 named. The top 50 genes here were used in further analysis of related transcirption factors. Data published in Drake, et al., Physiol Genomics 2013 45(12):449–61. (PDF 52 kb)

Additional file 6:

Figure S4. Literature Net analysis of differentially expressed genes in right ventricles after chronic hypoxia challenge compared to normoxic controls using previously published microarray data. Analysis of the top 50 genes increased (> 4.9 fold) and top 50 genes decreased (<− 3.9 fold) listed in Additional file 4: Figure S2 using Literature Net on Duke Gather (http://changlab.uth.tmc.edu/gather/gather.py). Publicly available data obtained from NCBI GEO published in Drake, et al., Physiol Genomics 2013 45(12):449–61, as described in the materials and methods. (PDF 24 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Oakley, R.H., Campen, M.J., Paffett, M.L. et al. Muscle-specific regulation of right ventricular transcriptional responses to chronic hypoxia-induced hypertrophy by the muscle ring finger-1 (MuRF1) ubiquitin ligase in mice. BMC Med Genet 19, 175 (2018). https://doi.org/10.1186/s12881-018-0670-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12881-018-0670-1

Keywords

BMC Medical Genetics

ISSN: 1471-2350

Contact us