ORIGINAL RESEARCH article

Front. Physiol., 23 August 2016

Sec. Striated Muscle Physiology

Volume 7 - 2016 | https://doi.org/10.3389/fphys.2016.00352

Mitochondrial and Metabolic Gene Expression in the Aged Rat Heart

  • 1. Balke Biodynamics Laboratory, Department of Kinesiology, University of Wisconsin-Madison Madison, WI, USA

  • 2. Departments of Agriculture, Food, and Nutritional Sciences, University of Alberta-Edmonton Edmonton, AB, Canada

Abstract

Aging is associated with a decline in cardiac function. Exercise intervention has been suggested as a way to improve this decrement. Age-related decline in cardiac function is associated with decreases in fatty acid oxidation, mitochondrial function, and AMP-activated protein kinase (AMPK) activity. The molecular mechanisms involved with age-related changes in mitochondrial function and substrate metabolism are poorly understood. We determined gene expression differences in hearts of Young (6 mo), Old (33 mo), and old exercise trained (Old + EXE) (34 mo) FBN rats, using Qiagen PCR arrays for Glucose, Fatty acid, and Mitochondrial metabolism. Old rats demonstrated decreased (p < 0.05) expression for key genes in fatty acid oxidation, mitochondrial function, and AMPK signaling. There were no differences in the expression of genes involved in glucose metabolism with age. These gene expression changes occurred prior to altered protein translation as we found no differences in the protein content of peroxisome proliferator activated receptor gamma, coactivators 1 alpha (PGC-1α), peroxisome proliferator activated receptor alpha (PPARα), and AMPKα2 between young and old hearts. Four months of exercise training did not attenuate the decline in the gene expression in aged hearts. Despite this lack of change in gene expression, exercise-trained rats demonstrated increased exercise capacity compared to their sedentary counterparts. Taken together, our results show that differential expression of genes associated with fatty acid metabolism, AMPK signaling and mitochondrial function decrease in the aging heart which may play a role in age-related declines in fatty acid oxidation, AMPK activity, and mitochondrial function in the heart.

Introduction

Aging is associated with a decline in cardiac function in humans as well as rodents (Lakatta and Sollott, 2002). The decline in cardiac function with age is associated with changes in mitochondrial function and energy metabolism (Abu-Erreish et al., 1977; McMillin et al., 1993; Fannin et al., 1999; Kates et al., 2003; Bhashyam et al., 2007; Chakravarti et al., 2008). Since the myocardium derives nearly all of its energetic needs from the oxidation of pyruvate and fatty acids within the mitochondria, mitochondrial changes have great potential to contribute to cardiac dysfunction with age. Indeed, mitochondrial energetic deficiency with aging has been well-documented (McMillin et al., 1993; Fannin et al., 1999; Phaneuf and Leeuwenburgh, 2002; Chakravarti et al., 2008; Ventura-Clapier et al., 2008). The mechanisms of this mitochondrial dysfunction may include biogenesis that is inadequate to match the increasing demand (Goffart et al., 2004), as well as increased mitochondrial uncoupling and decreased substrate availability (Murray et al., 2004). Several studies have documented age-dependent impairment in the mitochondrial respiratory capacity (Kumaran et al., 2005; Navarro and Boveris, 2007).

Concomitant with mitochondria functional changes there is a change in cardiac substrate utilization during the aging process. At rest, fatty acids are the primary substrate for ATP supply in the myocardium and glucose uptake and oxidation supplies the remainder of the carbon substrates (Wisneski et al., 1985). There is some evidence that this substrate utilization changes with age (Abu-Erreish et al., 1977; McMillin et al., 1993; Sample et al., 2006), with most evidence indicating an age-related reduction in fatty acid oxidation, with the implication of increased reliance on glucose, although this is not universally seen (Sample et al., 2006). Interestingly, it has been determined that glucose utilization itself does not increase (Abu-Erreish et al., 1977; McMillin et al., 1993; Kates et al., 2003).

The molecular mechanisms underlying age-related changes in mitochondrial function or substrate energy metabolism is largely unknown. Decreases in mRNA levels of peroxisome proliferator-activated receptor (PPARα), and some of its downstream targets (i.e., CPT-1, β-hydroxyacyl dehydrogenase) have been observed with age (Iemitsu et al., 2002). In addition, peroxisome proliferator activated receptor γ coactivator-1α (PGC-1α) protein content is reduced in the aging left ventricle (Turdi et al., 2010). PPARα is a transcription factor for genes involved with fatty acid transport and β-oxidation (Huss and Kelly, 2004). PGC-1α, a co-transcription factor is known to stimulate mitochondrial DNA replication and the coding for genes involved with oxidative phosphorylation (Huss and Kelly, 2004). Another mechanism for age-related changes in substrate metabolism may be related to AMP-activated protein kinase (AMPK) activity. Dually activated by AMP and by upstream phosphorylation, this kinase promotes fatty acid oxidation, glucose uptake, and glycogenolysis while it inhibits anabolic processes such as fatty acid synthesis (Munday et al., 1988; Coven et al., 2003; Hawley et al., 2003, 2005; Scott et al., 2004). There is evidence that AMPK activity declines with age (Gonzalez et al., 2004; Turdi et al., 2010). Changes in the expression, protein content, or activity of PPARα, PGC-1α, and AMPK can help explain some of the known changes in substrate metabolism and mitochondrial function with age. However, these processes are the result of complex pathways that require the coordinated expression and function of a large number of genes and proteins, but there are few studies that have examined the effects of aging on the expression of the comprehensive array of genes that are associated with substrate metabolism and mitochondrial function.

Exercise training is known to improve indices of cardiac function in humans (Fortney et al., 1992; Seals et al., 1994) and rodents (Spurgeon et al., 1983; Iemitsu et al., 2002). Specifically, these functional improvements include, increased contractile function (Spurgeon et al., 1983; Fortney et al., 1992; Seals et al., 1994; Iemitsu et al., 2002) and increased maximal oxygen uptake (VO2max) (Ogawa et al., 1992; Seals et al., 1994; Stratton et al., 1994). There is some evidence that exercise can alter the metabolic phenotype of the aging heart. Exercise training in aging Wistar rats increased AMPK activity, PPARα mRNA and protein content, and proteins involved in fatty acid oxidation in the young and aged heart (Iemitsu et al., 2002; Rimbaud et al., 2009; Dobrzyn et al., 2013). Also, lifelong voluntary wheel running in mice increased electron transport chain (ETC) complex IV gene expression (Bronikowski et al., 2003). There is molecular evidence that exercise training can increase the mRNA and protein levels of important proteins involved with fatty acid oxidation and oxidative phosphorylation, these markers suggest that these metabolic pathways may be improved with exercise training (Iemitsu et al., 2002; Bronikowski et al., 2003; Rimbaud et al., 2009; Dobrzyn et al., 2013). However, substrate metabolism and mitochondrial oxidative phosphorylation in the heart are regulated by the transcription, translation, and activity of many genes in order to optimally function.

Therefore, the aims of this study was to: (1) determine the effects of age on the expression of a large number of genes related to the pathways of glucose and fatty acid metabolism, and mitochondrial function; and (2) determine whether exercise training could mitigate age-related changes in the expression of metabolic and mitochondrial genes in the aging rat heart. We hypothesize that expression of genes associated with the pathways of fatty acid metabolism, AMPK signaling, and mitochondrial function will decrease with age and that the addition of exercise training in these aged rats will mitigate this decrease in gene expression.

Materials and methods

Animals

Male Fischer 344 x Brown Norway hybrid rats (FBN), were obtained from the National Institute on Aging colony at Harlan Industries (Indianapolis, IN). The FBN hybrid rat is a long-lived strain with a median life-span of 33 months and a maximum life-span of 40 months. The FBN rat is considered a “healthy aging model” widely used and highly recommended for gerontological research. All rats were confined to standard size rodent cages and housed 2 rats per cage. Rats had access to food and water ad libitum and were acclimated to reverse daylight (12 h dark, 12 h light). Body weights and average food intake were monitored through the course of the study. Rats were randomly assigned to one of three groups: Young (6 month), Old (33–34 month), and Old + Exercise (Old+EXE) (33–34 month). Animal housing and handling was carried out under the guidelines of the University of Wisconsin-Madison Institutional Animal Care and Use Committees and conducted in pathogen-free facilities that are accredited by the American Association of Accreditation of Laboratory Animal Care.

Tissue collection

Old and Old + EXE hearts were removed and flash frozen in liquid nitrogen 72 h from the last maximal exercise test in order to control for transient gene expression changes due to acute exercise (Neufer and Dohm, 1993; Pilegaard et al., 2000). From the Old (n = 9) and Old + EXE (n = 9) rats, 5 hearts were randomly selected for qRT-PCR, Western blot, and citrate synthase experiments, while Young (n = 5) hearts were used for molecular analysis.

Maximal exercise testing

Assessment of peak exercise capacity occurred each month beginning at 30 months until 34 months of age. The exercise test started with a treadmill speed at 5 m/min and progressively increased by 3 m/min at each exercise stage. Animals ran at each new treadmill speed for 3 min to assure steady-state values. To encourage the rats to run, the treadmill was equipped with an electric shock grid at the rear of the treadmill. The shock grid was set to deliver a 0.2 mA current, which gives an uncomfortable shock but does not physically harm or injure the rat. The test was terminated when animals were no longer able to maintain position on the treadmill, and the highest speed was recorded as peak exercise capacity. Blood lactate levels were measured during the graded exercise tests. Blood lactate levels were used to quantify relative exercise intensity. At the end of each 3 min interval, animals were briefly removed from the treadmill, immobilized with the tail extended, the lateral tail vein pierced, and a drop (~25 μL) of blood was collected on a lactate strip inserted into a lactate meter (Lactate Plus Meter, Nova Biomedical). Blood lactate values were recorded at each exercise intensity level during the exercise test.

Training protocol

Trained animals were exercised on a motor driven treadmill during their 12 h dark cycle. Before training began, animals were treadmill acclimatized for 2 weeks at 5 m/min for 5 min during the first week and at 5 m/min for 10 min the second week. After acclimatization, maximal exercise tests were performed and blood lactate measurements were taken at each successive speed on the treadmill during the test. The endurance trained rats began exercise training at 30 months of age with a training program that following 2 weeks of treadmill acclimation consisted of 30 min/day, 5 days per week at a speed that corresponded to each animal's lactate threshold. Thus, the training speed was adjusted each month based on the results of the maximal exercise tests (above). All exercise training sessions included a 3 min warm-up period at 5 m/min. No negative stimuli (electric shock) were used during the daily exercise training of the animals to minimize stress involved in exercise for these aging animals.

RNA isolation and cDNA preparation

Rats were terminated using isoflurane/pneumothorax euthanasia and hearts were removed and flash frozen in liquid nitrogen. Total RNA from each rat was isolated from the left ventricular free wall using RNeasy® Microarray Tissue Mini Kit (Qiagen) according to the manufacturer instruction. Immediately after elution RNA concentration and purity was measured spectrophotometrically (Beckman-Coulter). For each sample 700 ng of total RNA was reversed transcribed using the RT2 First Strand Kit (Qiagen). The reaction was performed at 42°C for 15 min followed by a termination step at 95°C for 5 min. cDNA was stored at −20°C until qRT-PCR.

qRT-PCR data analysis

A total of 251 genes of interest across three PCR arrays from Qiagen were assayed for this study (Glucose Metabolism RT2 Profiler PCR Array, Fatty Acid Metabolism RT2 Profiler PCR Array and RT2 Custom PCR Array for mitochondrial function). A list of all of the genes contained in each array can be found in Tables 13. PCR was performed on Step-One Plus PCR system (Applied Biosystems), according to the manufacturer's instructions. For data analysis, the ΔΔCt method was used with the aid of a Microsoft Excel spreadsheet containing algorithms and a student's t-test was used to analyze differences in fold-changes in gene expression provided by the manufacturer (Hassmann-Poznańska et al., 2013; Wu et al., 2013; Okada et al., 2014). The ΔΔCt method for calculating differences in gene expression are as follows: (ΔCt = CtGOI − CtAVG HKG, where GOI is the gene of interest and HKG is the housekeeping gene selected. ΔΔCt = (ΔCt (Experimental group) −ΔCt (Control)). The housekeeping gene selected for data analysis was Ribosomal protein, large P1 (Rplp1). Fold-changes were then calculated and expressed as log-normalized ratios of values from Old/Young, Old + EXE/Young and Old + EXE/Old heart tissues.

Table 1

SymbolDescription
AclyATP citrate lyase
Aco1Aconitase 1, soluble
Aco2Aconitase 2, mitochondrial
AglAmylo-1,6-glucosidase, 4-alpha-glucanotransferase
AldoaAldolase A, fructose-bisphosphate
AldobAldolase B, fructose-bisphosphate
AldocAldolase C, fructose-bisphosphate
Bpgm2,3-bisphosphoglycerate mutase
CsCitrate synthase
DlatDihydrolipoamide S-acetyltransferase
DldDihydrolipoamide dehydrogenase
DlstDihydrolipoamide S-succinyltransferase (E2 component of 2-oxo-glutarate complex)
Eno1Enolase 1, (alpha)
Eno2Enolase 2, gamma, neuronal
Eno3Enolase 3, beta, muscle
Fbp1Fructose-1,6-bisphosphatase 1
Fbp2Fructose-1,6-bisphosphatase 2
FhFumarate hydratase 1
G6pcGlucose-6-phosphatase, catalytic subunit
G6pc3Glucose 6 phosphatase, catalytic, 3
G6pdGlucose-6-phosphate dehydrogenase
GalmGalactose mutarotase (aldose 1-epimerase)
GapdhGlyceraldehyde-3-phosphate dehydrogenase
GapdhsGlyceraldehyde-3-phosphate dehydrogenase, spermatogenic
GckGlucokinase
GpiGlucose phosphate isomerase
Gsk3aGlycogen synthase kinase 3 alpha
Gsk3bGlycogen synthase kinase 3 beta
Gys1Glycogen synthase 1, muscle
Gys2Glycogen synthase 2
H6pdHexose-6-phosphate dehydrogenase (glucose 1-dehydrogenase)
Hk2Hexokinase 2
Hk3Hexokinase 3 (white cell)
Idh1Isocitrate dehydrogenase 1 (NADP+), soluble
Idh2Isocitrate dehydrogenase 2 (NADP+), mitochondrial
Idh3aIsocitrate dehydrogenase 3 (NAD+) alpha
Idh3bIsocitrate dehydrogenase 3 (NAD+) beta
Idh3gIsocitrate dehydrogenase 3 (NAD), gamma
Mdh1Malate dehydrogenase 1, NAD (soluble)
Mdh1bMalate dehydrogenase 1B, NAD (soluble)
Mdh2Malate dehydrogenase 2, NAD (mitochondrial)
OgdhlOxoglutarate dehydrogenase-like
PcPyruvate carboxylase
Pck1Phosphoenolpyruvate carboxykinase 1 (soluble)
Pck2Phosphoenolpyruvate carboxykinase 2 (mitochondrial)
Pdha2Pyruvate dehydrogenase (lipoamide) alpha 2
PdhbPyruvate dehydrogenase (lipoamide) beta
PdhxPyruvate dehydrogenase complex, component X
Pdk1Pyruvate dehydrogenase kinase, isozyme 1
Pdk2Pyruvate dehydrogenase kinase, isozyme 2
Pdk3Pyruvate dehydrogenase kinase, isozyme 3
Pdk4Pyruvate dehydrogenase kinase, isozyme 4
Pdp2Pyruvate dehyrogenase phosphatase catalytic subunit 2
PdprPyruvate dehydrogenase phosphatase regulatory subunit
PfklPhosphofructokinase, liver
Pgam2Phosphoglycerate mutase 2 (muscle)
Pgk1Phosphoglycerate kinase 1
Pgk2Phosphoglycerate kinase 2
Pgls6-phosphogluconolactonase
Pgm1Phosphoglucomutase 1
Pgm2Phosphoglucomutase 2
Pgm3Phosphoglucomutase 3
Phka1Phosphorylase kinase, alpha 1
PhkbPhosphorylase kinase, beta
Phkg1Phosphorylase kinase, gamma 1
Phkg2Phosphorylase kinase, gamma 2 (testis)
PklrPyruvate kinase, liver, and RBC
Prps1Phosphoribosyl pyrophosphate synthetase 1
Prps1l1Phosphoribosyl pyrophosphate synthetase 1-like 1
PyglPhosphorylase, glycogen, liver
PygmPhosphorylase, glycogen, muscle
RbksRibokinase
RpiaRibose 5-phosphate isomerase A
SdhaSuccinate dehydrogenase complex, subunit A, flavoprotein (Fp)
SdhbSuccinate dehydrogenase complex, subunit B, iron sulfur (Ip)
SdhcSuccinate dehydrogenase complex, subunit C, integral membrane protein
SdhdSuccinate dehydrogenase complex, subunit D, integral membrane protein
Sucla2Succinate-CoA ligase, ADP-forming, beta subunit
Suclg1Succinate-CoA ligase, alpha subunit
Suclg2Succinate-CoA ligase, GDP-forming, beta subunit
Taldo1Transaldolase 1
TktTransketolase
Tpi1Triosephosphate isomerase 1
Ugp2UDP-glucose pyrophosphorylase 2

Gene list Glucose metabolism array.

Table 2

SymbolDescription
Acaa1aAcetyl-Coenzyme A acyltransferase 1A
Acaa2Acetyl-Coenzyme A acyltransferase 2
Acad10Acyl-Coenzyme A dehydrogenase family, member 10
Acad11Acyl-Coenzyme A dehydrogenase family, member 11
Acad9Acyl-Coenzyme A dehydrogenase family, member 9
AcadlAcyl-Coenzyme A dehydrogenase, long-chain
AcadmAcyl-Coenzyme A dehydrogenase, C-4 to C-12 straight chain
AcadsAcyl-Coenzyme A dehydrogenase, C-2 to C-3 short chain
AcadsbAcyl-Coenzyme A dehydrogenase, short/branched chain
AcadvlAcyl-Coenzyme A dehydrogenase, very long chain
Acat1Acetyl-coenzyme A acetyltransferase 1
Acat2Acetyl-Coenzyme A acetyltransferase 3
Acot12Acyl-CoA thioesterase 12
Acot2Acyl-CoA thioesterase 2
Acot3Acyl-CoA thioesterase 3
Acot7Acyl-CoA thioesterase 7
Acot8Acyl-CoA thioesterase 8
Acot9Acyl-CoA thioesterase 9
Acox1Acyl-Coenzyme A oxidase 1, palmitoyl
Acox2Acyl-Coenzyme A oxidase 2, branched chain
Acox3Acyl-Coenzyme A oxidase 3, pristanoyl
Acsbg1Acyl-CoA synthetase bubblegum family member 1
Acsbg2Acyl-CoA synthetase bubblegum family member 2
Acsl1Acyl-CoA synthetase long-chain family member 1
Acsl3Acyl-CoA synthetase long-chain family member 3
Acsl4Acyl-CoA synthetase long-chain family member 4
Acsl5Acyl-CoA synthetase long-chain family member 5
Acsl6Acyl-CoA synthetase long-chain family member 6
Acsm2aAcyl-CoA synthetase medium-chain family member 2
Acsm3Acyl-CoA synthetase medium-chain family member 3
Acsm4Acyl-CoA synthetase medium-chain family member 4
Acsm5Acyl-CoA synthetase medium-chain family member 5
Aldh2Aldehyde dehydrogenase 2 family (mitochondrial)
Bdh13-hydroxybutyrate dehydrogenase, type 1
Bdh23-hydroxybutyrate dehydrogenase, type 2
Cpt1aCarnitine palmitoyltransferase 1a, liver
Cpt1bCarnitine palmitoyltransferase 1b, muscle
Cpt1cCarnitine palmitoyltransferase 1c
Cpt2Carnitine palmitoyltransferase 2
CratCarnitine acetyltransferase
CrotCarnitine O-octanoyltransferase
Decr12,4-dienoyl CoA reductase 1, mitochondrial
Decr22,4-dienoyl CoA reductase 2, peroxisomal
Echs1Enoyl Coenzyme A hydratase, short chain, 1, mitochondrial
Eci2Enoyl-Coenzyme A delta isomerase 2
EhhadhEnoyl-Coenzyme A, hydratase/3-hydroxyacyl Coenzyme A dehydrogenase
Fabp1Fatty acid binding protein 1, liver
Fabp2Fatty acid binding protein 2, intestinal
Fabp3Fatty acid binding protein 3, muscle, and heart
Fabp4Fatty acid binding protein 4, adipocyte
Fabp5Fatty acid binding protein 5, epidermal
Fabp6Fatty acid binding protein 6, ileal
Fabp7Fatty acid binding protein 7, brain
GcdhGlutaryl-Coenzyme A dehydrogenase
GkGlycerol kinase
Gk2Glycerol kinase 2
Gpd1Glycerol-3-phosphate dehydrogenase 1 (soluble)
Gpd2Glycerol-3-phosphate dehydrogenase 2, mitochondrial
HadhaHydroxyacyl-Coenzyme A dehydrogenasea
Hmgcl3-hydroxymethyl-3-methylglutaryl-Coenzyme A lyase
Hmgcs13-hydroxy-3-methylglutaryl-Coenzyme A synthase 1 (soluble)
Hmgcs23-hydroxy-3-methylglutaryl-Coenzyme A synthase 2 (mitochondrial)
LipeLipase, hormone sensitive
LplLipoprotein lipase
MceeMethylmalonyl CoA epimerase
MutMethylmalonyl-Coenzyme A mutase
Oxct2a3-oxoacid CoA transferase 2A
PecrPeroxisomal trans-2-enoyl-CoA reductase
Ppa1Pyrophosphatase (inorganic) 1
Prkaa1Protein kinase, AMP-activated, alpha 1 catalytic subunit
Prkaa2Protein kinase, AMP-activated, alpha 2 catalytic subunit
Prkab1Protein kinase, AMP-activated, beta 1 non-catalytic subunit
Prkab2Protein kinase, AMP-activated, beta 2 non-catalytic subunit
PrkacaProtein kinase, cAMP-dependent, catalytic, alpha
PrkacbProtein kinase, cAMP dependent, catalytic, beta
Prkag1Protein kinase, AMP-activated, gamma 1 non-catalytic subunit
Prkag2Protein kinase, AMP-activated, gamma 2 non-catalytic subunit
Prkag3Protein kinase, AMP-activated, gamma 3 non-catalytic subunit
Slc27a1Solute carrier family 27 (fatty acid transporter), member 1
Slc27a2Solute carrier family 27 (fatty acid transporter), member 2
Slc27a3Solute carrier family 27 (fatty acid transporter), member 3
Slc27a4Solute carrier family 27 (fatty acid transporter), member 4
Slc27a5Solute carrier family 27 (fatty acid transporter), member 5
Slc27a6Solute carrier family 27 (fatty acid transporter), member 6

Gene list for Fatty Acid metabolism array.

a

Hydroxyacyl-Coenzyme A dehydrogenase 3-ketoacyl-Coenzyme A thiolase/enoyl-Coenzyme A hydratase (trifunctional protein), alpha subunit.

Table 3

SymbolDescription
Ppargc1Peroxisome proliferator-activated receptor gamma, coactivator 1 alpha
Nrf1Nuclear respiratory factor 1
TfamTranscription factor A, mitochondrial
Sir2Sirtuin (Silent mating type information regulation 2 homolog) 2
EsrraEstrogen related receptor, alpha
Prkaa2Protein kinase, AMP-activated, alpha 2 catalytic subunit
CREB1cAMP response element binding protein 1
CaMK4Calcium/calmodulin-dependent protein kinase IV
SOD2Superoxide dismutase 2, mitochondrial
UCP2Uncoupling protein 2 (mitochondrial proton carrier)
UCP3Uncoupling protein 3 (mitochondrial proton carrier)
BaxBcl2–associated X protein
Bcl2B-cell CLL/lymphoma 2
Casp3Caspase 3
GabpaGA binding protein transcription factor (alpha subunit) (i.e., NRF-2a)
UCP1Uncoupling protein 1 (mitochondrial proton carrier)
PparaPeroxisome proliferator activated receptor alpha
Slc2a4Solute carrier family 2 (facilitated glucose transporter) member 4
CD36CD 36 molecule (thrombospondin receptor)
AcacbAcetyl-Coenzyme A carboxylase beta
MlycdMalonyl-CoA decarboxylase
Ndufa1NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 1
Ndufa10NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 10
Ndufa11NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 11
Ndufa2NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 2
Ndufa5NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 5
Ndufa6NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 6
Ndufa7NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 7
Ndufa8NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 8
Ndufa9NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 9
Ndufab1NADH dehydrogenase (ubiquinone) 1 alpha/beta subcomplex, 1
Ndufb2NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 2
Ndufb3NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 3
Ndufb5NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 5
Ndufb6NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 6
Ndufb7NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 7
Ndufb8NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 8
Ndufb9NADH dehydrogenase (ubiquinone) 1 beta subcomplex, 9
Ndufc2NADH dehydrogenase (ubiquinone) 1 subcomplex unknown, 2
Ndufs1NADH dehyrogenase (ubiquinone) Fe-S protein 1
Ndufs2NADH dehyrogenase (ubiquinone) Fe-S protein 2
Ndufs3NADH dehyrogenase (ubiquinone) Fe-S protein 3
Ndufs4NADH dehyrogenase (ubiquinone) Fe-S protein 4
Ndufs6NADH dehyrogenase (ubiquinone) Fe-S protein 6
Ndufs7NADH dehyrogenase (ubiquinone) Fe-S protein 7
Ndufs8NADH dehyrogenase (ubiquinone) Fe-S protein 8
Ndufv1NADH dehydrogenase (ubiquinone) flavoprotein 1
Ndufv2NADH dehyrogenase (ubiquinone) flavoprotein 2
Bcs1lBCS1-like yeast
Cyc1Cytochrome c-1
UqcrbUbiquinol-cytochrome c reductase binding protein
Uqcrc1Ubiquinol-cytochrome c reductase core protein 1
Uqcrc2Ubiquinol-cytochrome c reductase core protein 2
Uqcrfs1Ubiquinol-cytochrome c reductase, Rieske iron-sulfur polypeptide 1
UqcrhUbiquinol-cytochrome c reductase hinge protein
UqcrqUbiquinol-cytochrome c reductase, complex III subunit VII
Cox15COX 15 homolog, cytochrome c oxidase assembly protein (yeast)
Cox17COX 17 cytochrome c oxidase assembly homolog (S. cerevisiae)
Cox4i1Cytochrome c oxidase subunit IV isoform 1
Cox4i2Cytochrome c oxidase subunit IV isoform 2
Cox5aCytochrome c oxidase subunit Va
Cox5bCytochrome c oxidase subunit Vb
Cox6a1Cytochrome c oxidase subunit VIa polypeptide 1
Cox6a2Cytochrome c oxidase subunit VIa polypeptide 2
Cox6cCytochrome c oxidase subunit VIc
Cox7a2Cytochrome c oxidase subunit VIIa polypeptide 2
Cox7a2lCytochrome c oxidase subunit VIIa polypeptide 2 like
Cox7bCytochrome c oxidase subunit VIIb polypeptide
Cox8aCytochrome c oxidase subunit VIIIa
Cox8cCytochrome c oxidase subunit VIIIc
mfn1Mitofusin 1
mfn2Mitofusin 2
fis1Fission 1 (mitochondrial outer membrane) homolog (S. cerevisiae)
lonp1Lon protease
Aifm2Apoptosis-inducing factor, mitochondrial-associated 2
Bcl2l1Bcl2-like 1
ClpbClpB caseinolytic peptidase B homolog (E. coli) (i.e. HSP 78)
pnpt1Polyribonucleotide nucleotidyltransferase 1
Me1Malic enzyme 1, NADP(+)-dependent, cytosolic
Foxo3Forkhead box O3
Camkk2Calcium/calmodulin-dependent protein kinase kinase 2,beta
Stk11Serine/threonine kinase 11 (i.e., LKB1)
Ppargc1bPeroxisome proliferator-activated receptor gamma, coactivator 1 beta
Tp53Tumor protein p53

Gene list for the Custom array.

Western blotting

Whole cell lysate from isolated left ventricle was prepared using CelLytic™ MT Cell Lysis Reagent and 1:100 dilution of Protease Inhibitor Cocktail (Sigma, St. Louis, MO). The protein concentration was determined using Bio-Rad Protein Assay (BioRad, Hercules, CA). Thirty micrograms of whole cellular protein per lane was separated by SDS-PAGE with a 4–12% Bis-Tris Criterion™ XT gel (XT MOPS running buffer) and blotted onto a nitrocellulose membrane. The membrane was incubated with blocking buffer (5% non-fat dry milk/TBS/0.1% Tween 20) at room temperature for 1 h. The membrane was then probed with the primary antibodies diluted in blocking buffer overnight at 4°C. Subsequently, membranes were incubated with horseradish peroxidase-conjugated anti-rabbit antibody diluted 1:2000 in blocking buffer. Blots were developed with Clarity™ Western ECL substrate (BioRad, Hercules, CA) and imaged with GE ImageQuant LAS-4000 (GE Healthcare Bio-Sciences, Pittsburgh, PA). The image of the blots was uploaded and densitometry analysis was done with Image Studio Lite (LI-COR Biosciences, Lincoln, NE). Protein content was measured from the densitometry units from PGC-1α, AMPKα2,PPARα, and normalized to vinculin.

Antibodies

Rabbit anti-PGC-1α, AMPKα2, and Vinculin antibodies were acquired from Cell Signaling Technology, Inc., Beverly, MA. Rabbit anti-PPARα antibody was purchased from Santa Cruz Biotechnology, Inc., Santa Cruz, CA.

Citrate synthase activity

Maximal citrate synthase (CS) activity was determined in left ventricular homogenates using (Citrate Synthase Assay Kit, Sigma-Aldrich, St Louis, MO, USA) according to the manufacturer's protocol with absorbance kinetically measured at 412 nm at baseline and after addition of oxaloacetate (Sigma-Aldrich). CS activity was normalized to protein content with tissue protein determined using the Bio-Rad protein assay.

Statistical analysis

A two-tailed Student's t-test was used for analyzing differences in gene expression and exercise performance. A one-way ANOVA was utilized to analyze differences in protein content and citrate synthase activity. Post-hoc analysis was performed when the one-way ANOVA was significant using Tukey's LSD to assess between group differences. Significance was determined at p < 0.05.

Results

Effects of exercise training on exercise performance

The effect of 4 months of exercise training on exercise performance in 34 month old trained (Old + EXE) FBN rats compared to sedentary (Old SED) rats was significant (Figure 1). We measured exercise capacity in two different ways; (1) treadmill speed at lactate threshold (LT) and (2) peak exercise capacity (the final treadmill speed achieved during the maximal exercise test). Exercise training significantly improved lactate threshold (p < 0.01) in the Old + EXE rats compared to Old SED rats (Figure 1A). Peak exercise capacity was also greater (p = 0.01) in the Old + EXE rats compared to Old SED rats (Figure 1B).

Figure 1

Gene expression changes with age

To determine the effects of age on the expression of genes involved in cardiac glucose and fatty acid metabolism as well as mitochondrial function we performed qRT-PCR using three different PCR arrays. Table 4 shows all genes in the energetic pathways for which expression was significantly different between Old and Young rats. Between the three arrays there were a total of 44 genes that were differentially expressed (p < 0.05) in Old compared to Young. Of the 44 genes differentially expressed, 42 of these genes had decreased expression in Old compared to Young. Notably, the expression profile of three important pathways of energy production were altered with age; fatty acid oxidation (FAO), mitochondrial biogenesis, and AMPK signaling. There were no significant changes in glucose metabolism gene expression with age. Figure 2 shows changes in expression with age in key genes involved in the pathways of FAO, mitochondrial biogenesis, and AMPK signaling.

Table 4

Gene nameGene symbolFold Regulationp-value
FA METABOLISM
Acyl-CoA thioesterase 12Acot122.880.026
Acyl-CoA thioesterase 7Acot7–1.860.017
Acyl-CoA thioesterase 8Acot8–1.440.02
Acyl-Coenzyme A oxidase 1, palmitoylAcox1–2.10.032
Acyl-CoA synthetase bubblegum family member 1Acsbg12.090.034
Acyl-CoA synthetase long-chain family member 6Acsl6–1.980.029
3-hydroxybutyrate dehydrogenase, type 1Bdh1–2.550.036
Carnitine palmitoyltransferase 1cCpt1c–1.490.027
Carnitine palmitoyltransferase 2Cpt2–2.020.001
Carnitine acetyltransferaseCrat–2.690.011
2,4-dienoyl CoA reductase 1, mitochondrialDecr1–2.050.03
Enoyl Coenzyme A hydratase, short chain, 1, mitochondrialEchs1–2.390.004
Enoyl-Coenzyme A delta isomerase 2Eci2–1.790.049
Enoyl-Coenzyme A, hydratase/3-hydroxyacyl Coenzyme A dehydrogenaseEhhadh–3.290.034
Fatty acid binding protein 3, muscle and heartFabp3–2.840.017
Glutaryl-Coenzyme A dehydrogenaseGcdh–1.970.013
Glycerol-3-phosphate dehydrogenase 2, mitochondrialGpd2–2.640.028
Hydroxyacyl-Coenzyme A dehydrogenase/3-ketoacyl-Coenzyme A thiolase/enoyl-Coenzyme A hydratase (trifunctional protein), alpha subunitHadha–2.230.035
3-hydroxymethyl-3-methylglutaryl-Coenzyme A lyaseHmgcl–1.510.03
Lipase, hormone sensitiveLipe–2.140.03
Lipoprotein lipaseLpl–2.310.025
Peroxisomal trans-2-enoyl-CoA reductasePecr–1.670.013
Pyrophosphatase (inorganic) 1Ppa1–1.780.03
Protein kinase, AMP-activated, alpha 2 catalytic subunitPrkaa2–2.090.034
Protein kinase, AMP-activated, beta 1 non-catalytic subunitPrkab1–3.960.041
Protein kinase, AMP-activated, beta 2 non-catalytic subunitPrkab2–1.70.041
Protein kinase, cAMP-dependent, catalytic, alphaPrkaca–2.030.005
Protein kinase, cAMP dependent, catalytic, betaPrkacb–1.650.031
Protein kinase, AMP-activated, gamma 1 non-catalytic subunitPrkag1–1.770.002
Solute carrier family 27 (fatty acid transporter), member 1Slc27a1–2.020.03
MITOCHONDRIAL FUNCTION
Peroxisome proliferator—activated receptor gamma, coactivator 1 alphaPpargc1–2.140.03
NADH dehydrogenase (ubiquinone) Fe-S protein 7Ndufs7–2.240.015
Ubiquinol - cytochrome c reductase, complex III subunit VIIUqcrq–1.480.05
COX17 cytochrome c oxidase assembly homolog (S. cerevisiae)Cox17–1.680.011
Cytochrome c oxidase subunit IV isoform 1Cox4i1–1.720.02
Cytochrome c oxidase subunit VIIa polypeptide 2 likeCox7a2l–1.350.049
Mitofusin 1mfn1–2.410.01
Mitofusin 2mfn2–3.210.028
Clpb caseinolytic peptidase B homolog (E. coli) (i.e. HSP78)clpb–2.30.039
Calcium/calmodulin-dependent protein kinase kinase 2, betaCamkk2–1.770.017
Serine/threonine kinase 11 (i.e. LKB1)Stk11–1.860.029
Peroxisome proliferator - activated receptor gamma, coactivator 1 betaPpargc1b–2.690.012
Tumor protein p53Tp53–2.190.025

Age-related gene expression changes in the left ventricle.

Figure 2

Effects of exercise training on gene expression in aged hearts

Differences in gene expression were observed between Old and Old + EXE rat hearts (Table 5). 70 genes associated with glucose metabolism, FA metabolism and mitochondrial function were altered with exercise training in cardiac tissue of old rats compared to that of old sedentary rats. Of these 70 genes, only three were upregulated (Prkaa1, Slc27a1, Slc27a4), while 67 genes were downregulated with exercise training and involved genes from the Kreb's cycle, glucose transport (Glut4), fatty acid oxidation, and the mitochondrial electron transport chain.

Table 5

Gene nameGene symbolFold Regulationp-value
GLUCOSE METABOLISM
Amylo-1,6-glucosidase, 4-alpha-glucanotransferaseAgl–6.150.024
Aldolase A, fructose-bisphosphateAldoa–5.280.036
Enolase 1, (alpha)Eno1–4.10.042
Enolase 2, gamma, neuronalEno2–3.270.019
Glucose phosphate isomeraseGpi–5.590.005
Glycogen synthase kinase 3 alphaGsk3a–3.170.046
Isocitrate dehydrogenase 3 (NAD+) betaIdh3b–7.70.04
Isocitrate dehydrogenase 3 (NAD), gammaIdh3g–5.410.027
Malate dehydrogenase 2, NAD (mitochondrial)Mdh2–6.310.009
Pyruvate dehydrogenase kinase, isozyme 1Pdk1–9.290.004
Pyruvate dehydrogenase kinase, isozyme 2Pdk2–6.480.006
Phosphoglucomutase 3Pgm3–1.70.037
Phosphorylase kinase, gamma 1Phkg1–2.780.011
Ribose 5-phosphate isomerase ARpia–2.380.018
Succinate dehydrogenase complex, subunit A, flavoprotein (Fp)Sdha–3.930.038
Succinate-CoA ligase, ADP-forming, beta subunitSucla2–5.640.021
Solute carrier family 2 (facilitated glucose transporter), member 4Slc2a4–10.450.025
FA METABOLISM
Acetyl-Coenzyme A acyltransferase 1AAcaa1a–9.790.001
Acetyl-Coenzyme A acyltransferase 2Acaa2–11.70.029
Acyl-Coenzyme A dehydrogenase family, member 10Acad10–7.020.002
Acyl-Coenzyme A dehydrogenase family, member 11Acad11–14.380.002
Acyl-Coenzyme A dehydrogenase family, member 9Acad9–8.810.008
Acyl-Coenzyme A dehydrogenase, long-chainAcadl–21.180.002
Acyl-Coenzyme A dehydrogenase, C-4 to C-12 straight chainAcadm–15.170.009
Acyl-Coenzyme A dehydrogenase, C-2 to C-3 short chainAcads–17.090.001
Acyl-Coenzyme A dehydrogenase, short/branched chainAcadsb–7.750.007
Acyl-Coenzyme A dehydrogenase, very long chainAcadvl–9.490.002
Acetyl-coenzyme A acetyltransferase 1Acat1–14.790.004
Acetyl-Coenzyme A acetyltransferase 3Acat2–6.250.003
Acyl-CoA thioesterase 2Acot2–18.92< 0.001
Acyl-CoA thioesterase 3Acot3–3.520.024
Acyl-CoA thioesterase 7Acot7–4.430.027
Acyl-CoA thioesterase 8Acot8–2.680.002
Acyl-CoA thioesterase 9Acot9–5.410.001
Acyl-Coenzyme A oxidase 1, palmitoylAcox1–5.160.007
Acyl-Coenzyme A oxidase 3, pristanoylAcox3–4.470.015
Acyl-CoA synthetase long-chain family member 1Acsl1–18.06< 0.001
Acyl-CoA synthetase long-chain family member 3Acsl3–3.60.008
Acyl-CoA synthetase long-chain family member 4Acsl4–2.460.011
Acyl-CoA synthetase long-chain family member 5Acsl5–2.230.001
Acyl-CoA synthetase long-chain family member 6Acsl6–5.540.01
3-hydroxybutyrate dehydrogenase, type 1Bdh1–10.480.04
3-hydroxybutyrate dehydrogenase, type 2Bdh2–5.030.008
Carnitine palmitoyltransferase 1a, liverCpt1a–3.040.033
Carnitine palmitoyltransferase 1b, muscleCpt1b–7.930.003
Carnitine palmitoyltransferase 2Cpt2–7.020.004
Carnitine acetyltransferaseCrat–4.270.018
Fatty acid binding protein 3, muscle and heartFabp3–7.310.034
Protein kinase, AMP-activated, alpha 1 catalytic subunitPrkaa11.830.042
Protein kinase, cAMP-dependent, catalytic, alphaPrkaca–2.890.041
Protein kinase, cAMP dependent, catalytic, betaPrkacb–1.430.047
Solute carrier family 27 (fatty acid transporter), member 1Slc27a11.730.033
Solute carrier family 27 (fatty acid transporter), member 4Slc27a43.150.002
CD36 molecule (thrombospondin receptor)CD36–2.970.001
Peroxisome proliferator activated receptor alphaPpara–9.280.029
MITOCHONDRIAL FUNCTION
Sirtuin (silent mating type information regulation 2 homolog) 1 (S. cerevisiae)Sirt1–3.270.036
NADH dehyrogenase (ubiquinone) 1 alpha subcomplex, 1Ndufa1–2.240.04
NADH dehyrogenase (ubiquinone) 1 alpha subcomplex, 6 (B14)Ndufa6–4.70.041
NADH dehyrogenase (ubiquinone) 1 alpha subcomplex, 7Ndufa7–5.670.033
NADH dehyrogenase (ubiquinone) 1 alpha subcomplex, 8Ndufa8–10.410.005
NADH dehyrogenase (ubiquinone) 1 alpha subcomplex, 9Ndufa9–5.710.018
NADH dehyrogenase (ubiquinone) 1, alpha/beta subcomplex, 1Ndufab1–5.670.001
NADH dehyrogenase (ubiquinone) 1 beta subcomplex, 2Ndufb2–4.490.006
NADH dehyrogenase (ubiquinone) 1 beta subcomplex, 5Ndufb5–4.720.032
NADH dehyrogenase (ubiquinone) 1 beta subcomplex, 8Ndufb8–4.920.01
NADH dehyrogenase (ubiquinone) Fe-S protein 6Ndufs6–2.940.047
Cytochrome c-1Cyc1–5.470.033
Ubiquinol-cytochrome c reductase core protein 2Uqcrc2–5.380.016
Ubiquinol-cytochrome c reductase, Rieske iron-sulfur polypeptide 1Uqcrfs1–5.110.045

Exercise-related changes in gene expression in Old FBN rats.

Western blot analysis

In order to determine if the gene expression changes with age and exercise training were associated with altered protein content, we selected one protein from fatty acid oxidation (PPARα), AMPK signaling (AMPKα2), and mitochondrial biogenesis/function (PGC-1α) categories shown in Figure 2. Each of these genes had decreased mRNA expression and decreased protein content between Young and Old, although not significantly for PPARα (Figure 3). PGC-1α and AMPKα2 protein content was decreased in the Old + EXE compared to Young and not different than Old, but PPARα protein content was greater in the Old + EXE rat hearts compared to Old rat hearts.

Figure 3

Citrate synthase activity

To determine mitochondrial function and volume we assayed citrate synthase activity in left ventricular homogenates. We found that Old rat hearts had increased citrate synthase activity compared to Young hearts, and the Old + EXE hearts demonstrated no differences in citrate synthase activity compared to Young or Old hearts (Figure 4).

Figure 4

Discussion

To our knowledge this is the first study to determine the effects of age on the expression of a comprehensive group of genes related to cardiac substrate metabolism and mitochondrial function. The other primary aim of this study was to determine whether exercise training in aged rats could alter the age-related gene expression phenotype. We hypothesized that genes associated with fatty acid oxidation, AMPK signaling, and mitochondrial biogenesis/function would be decreased with age and that exercise training would mitigate these changes in aged rat hearts. We found that aging results in the decreased expression of many genes involved with energy metabolism and mitochondrial function but found that exercise training did not improve the downregulation of these genes. In fact, exercise training in aged rats resulted in the downregulation of 67 genes associated with energy metabolism and mitochondrial function compared to aged sedentary rat hearts. Genes associated with glucose metabolism were unaffected by age. The declines in metabolic gene expression were profound, as we observed decreased expression in 42 genes involved with fatty acid metabolism and mitochondrial function in aged hearts compared to young hearts. In aged sedentary hearts altered protein content of PGC-1α and AMPKα2 corresponded to declines in gene expression compared to young hearts. We observed that cardiac muscle citrate synthase activity; a common measure of mitochondrial volume in skeletal muscle (Larsen et al., 2012) was increased in the aged heart compared to the young heart, while exercise training in aged hearts showed no differences between young and old sedentary rat hearts. Overall, these results suggest a significant change in the expression of genes associated with cardiac metabolic pathways with age that is not improved with exercise training.

We determined that a large number of genes involved with fatty acid metabolism and mitochondrial energy metabolism/biogenesis were downregulated with age, consistent with previous reports suggesting age-related reductions in fatty acid oxidation (Abu-Erreish et al., 1977; McMillin et al., 1993; Kates et al., 2003), AMPK activity (Gonzalez et al., 2004; Turdi et al., 2010; Zhao et al., 2014), and mitochondrial function (Fannin et al., 1999; Wanagat et al., 2002; Kumaran et al., 2005; Bhashyam et al., 2007; Preston et al., 2008; Jian et al., 2011). Specifically, we found genes associated with FAO (i.e., CPT-2, HADHA), AMPK signaling (i.e., AMPKα2, CaMKK2, LKB1), and mitochondrial biogenesis (PGC-1α, PGC-1β) with age-related decrements in expression in the myocardium. Cardiac PGC-1α expression has previously been shown to decline with age (Preston et al., 2008; Turdi et al., 2010), although others have observed no age-related change in cardiac PGC-1α gene expression (LeMoine et al., 2006). Our results demonstrated that the genes associated with fatty acid oxidation that were downregulated with age were primarily PPARα regulated transcripts. As in previous studies, we found that expression of FAO genes downstream of PPARα were decreased (Iemitsu et al., 2002; LeMoine et al., 2006). However, our finding of a decline in PPARα regulated genes associated with FAO in aged hearts occurred despite no change in expression of the PPARα gene itself between young and old hearts. Although, this decline in genes associated with FAO is likely explained by reduced PPARα protein content in old hearts compared to young. This finding suggests that in the aged heart, changes in PPARα is likely to occur downstream of gene expression. Lastly, we demonstrated that the gene expression of AMPKα2 decreased with age, which is the major catalytic subunit in the heart and responsible for the phosphorylation of downstream proteins that confer AMPK's effect on energy metabolism in the heart (Dolinsky and Dyck, 2006). Also, we found decreased gene expression in the aging heart of two protein kinases (CaMKK2 & LKB1) that have been found to phosphorylate AMPK. These findings suggest that declining transcripts of AMPK and two upstream activators (CaMKK2 & LKB1), may play a role in the reduced AMPK activity that has been previously observed in aged hearts (Gonzalez et al., 2004; Turdi et al., 2010).

There is little information with regard to gene expression and protein content changes in exercise-trained aged rat hearts with regard to fatty acid oxidation, AMPK signaling, or mitochondrial function, but functional studies led us to hypothesize that exercise training would improve the age-related declines in metabolic gene expression. Surprisingly, we found that exercise training did not attenuate the age-related downregulation in the expression of genes involving fatty acid oxidation, AMPK signaling, and mitochondrial function. Specifically, in old exercise-trained hearts compared to old hearts we found a significant decline in the expression of the PPARα gene. PPARα regulated genes that were downregulated by exercise were Acyl CoA dehydrogenases (Acad), CD36, CPT1b, and CPT2. There were no differences in AMPK signaling genes (Figure 2B) between Old + EXE and Old rats; however, there was an upregulation of AMPKα1 in the exercise trained rats. AMPKα1 is ubiquitously expressed in cells and has lower levels of expression in the myocardium compared to AMPKα2 (Dolinsky and Dyck, 2006). We also found that compared to old rat hearts, exercise trained rat hearts demonstrated further downregulation of many genes involved with glucose transport (Glut4), Kreb's cycle and mitochondrial function (complex I and III in the electron transport chain). This downregulation of genes we observed in exercise trained aged hearts will require further work in order to better understand how exercise training induced downregulation of these genes affects cardiac function.

In order to determine if changes in gene expression with age were associated with altered protein content we measured the protein content of PGC-1α, PPARα, and AMPKα2. We found that the protein levels of PGC-1α and AMPKα2 were decreased in Old hearts compared to Young hearts and PPARα protein content trended toward a significant decline which was similar to the mRNA expression for these genes. In a previous report, mitochondrial oxygen consumption and expression of genes associated with mitogenesis and mitochondrial energy metabolism were both decreased, but mitochondrial number was increased and the mitochondrial marker, citrate synthase was not different in the senescent hearts compared to young (Preston et al., 2008). We observed a similar trend, of an increase in citrate synthase activity in Old hearts compared to Young hearts despite a decrease in gene expression of a large number of genes associated with mitochondrial biogenesis.

Western blot analysis determined that in the old-exercise trained rats PGC-1α and AMPKα2 protein content decreased compared to young hearts, similar to our findings of decreased gene expression of these two genes. However, we found that PGC-1α mRNA expression was reduced in Old + EXE hearts compared to Old hearts, yet the protein content of PGC-1α was similar. Also, we found that the PPARα mRNA expression was significantly reduced in Old + EXE hearts compared to Young and Old hearts, yet PPARα protein content was significantly greater in Old + EXE hearts compared to Old hearts and similar compared to Young hearts. These results are somewhat conflicted to a report demonstrating that exercise training attenuated the age-related decrement in PPARα mRNA expression and protein content in aged rat myocardium (Iemitsu et al., 2002). This study indicated that decreased mRNA expression did result in similar declines in protein content. However, functional changes do not always coincide with coordinate changes in gene expression or protein content (Burelle et al., 2004; Zhu et al., 2013; Lai et al., 2014). Changes in gene expression due to exercise training in aged hearts may or may not coincide with changes in that gene's protein product. Specifically, one study demonstrated that exercise training elicited adaptations that lead to an increase in exercise capacity along with augmented glucose and fatty acid oxidation in the myocardium despite no changes in the protein content of several proteins involved with glucose and fatty acid metabolism (Burelle et al., 2004). Our results demonstrated despite many genes involved with FAO in the exercise trained aged heart was reduced to aged sedentary hearts, PPARα protein content was greater in exercise trained aged hearts compared to aged sedentary hearts. Future studies looking at the protein content of major proteins involved with β-oxidation and fatty acid transport would enable us to see if exercise training in aged hearts tends to uncouple gene expression from protein content.

Citrate synthase is a common marker of mitochondrial volume used to indicate endurance exercise adaptations in skeletal muscle (Larsen et al., 2012). We examined citrate synthase activity in cardiac muscle and found, surprisingly, that aging resulted in an increase in citrate synthase activity and exercise training diminished this effect. Citrate synthase activity was not different in old-exercise trained hearts compared to young or old hearts. This is consistent with previous studies demonstrating that exercise training in rats did not elicit increases in cardiac citrate synthase activity (Oscai et al., 1971; Murakami et al., 1995; Zonderland et al., 1999; Siu et al., 2003; Rimbaud et al., 2009). We found that exercise training improved functional exercise capacity (Figure 4) despite either no changes in gene expression compared to Old hearts or in some cases, a further reduction in the expression of genes associated with energy metabolism and mitochondrial function in the heart. These data suggest that exercise training may impact myocardial energy metabolism and mitochondrial function downstream of gene expression. Also, exercise is known to induce adaptations in skeletal muscle (Hall et al., 1994; Bengtsson et al., 2001; Betik et al., 2008; Kang et al., 2013), which may have been responsible for the increased exercise capacity in our old exercise-trained rats.

Limitations

One limitation to this study is that we did not determine whether exercise training in young rats leads to a similar downregulation of the expression of these cardiac genes that we found in the exercise trained aged hearts. Cardiac gene expression changes due to exercise training in young rats have been well-studied. These studies showed that mitochondrial or metabolic gene expression in the young rat heart to either increase (Hall et al., 1994; Rimbaud et al., 2009; Dobrzyn et al., 2013; Wadley et al., 2016) or not change (Murakami et al., 1995; Iemitsu et al., 2003; Alessio et al., 2014) with exercise training compared to young sedentary rats. Specifically, young hearts respond to exercise training by increasing the expression of genes associated with glucose transport (Hall et al., 1994; Rimbaud et al., 2009), fatty acid oxidation (Rimbaud et al., 2009; Dobrzyn et al., 2013), and mitochondrial biogenesis (i.e., PGC-1α and Cox4il; Rinaldi et al., 2013; Wadley et al., 2016). The results of our study (further decreases in gene expression with exercise training compared to sedentary aging) compared to these previous studies suggest that gene expression changes due to exercise training may be different in the hearts of aged exercise-trained rats compared to young hearts. Future works looking at post-translational modifications and protein activity in genes associated with fatty acid oxidation and mitochondrial function may elucidate molecular mechanisms involved in potential differential exercise training responses between young and old rat hearts.

Another limitation to this study is that our primary endpoint measure was the expression of genes associated with metabolic signaling pathways, substrate energy metabolism and mitochondrial function. Along with our data, previous reports have indicated that changes in tissue mRNA levels are not necessarily predictive of changes in protein levels (Hall et al., 1994; Zhu et al., 2013; Lai et al., 2014). Interestingly, other studies have reported that roughly two-thirds of the variance in protein content may be explained by mRNA concentration (Lu et al., 2007; Shankavaram et al., 2007). However, due to an increasing appreciation of miRNA interactions, post-translational modifications effecting protein activity, and the unpredictable protein response to stimuli (i.e., exercise training), mRNA is not always a good predictor of protein abundance (Burelle et al., 2004; Lu et al., 2007; Shankavaram et al., 2007). Future studies looking at expression changes in cardiac metabolic and mitochondrial genes, may need to include miRNA expression changes and post-translational protein modifications which will provide a more thorough understanding of metabolic and mitochondrial function in the aged and exercise trained rat heart.

Conclusion

This study was the first to utilize a comprehensive approach in the study of age and exercise effects in aged hearts on substrate metabolism and mitochondrial function by an integration of gene expression, protein content, and protein activity. We found that cardiac aging results in the downregulation of a large number of genes associated with fatty acid metabolism, AMPK signaling and mitochondrial function, and exercise in aged-rats did not attenuate these changes and resulted in a further downregulation of expression in many metabolic and mitochondrial genes compared to aged sedentary rat hearts. We also found that gene expression changes may or may not coincide with protein expression. Taken together, our results demonstrated extensive age-related molecular changes in fatty acid metabolism, AMPK signaling and mitochondrial function. These molecular changes are far-reaching and cannot be described by changes in a single gene or a given gene's protein product.

Funding

This study was supported by National Institutes of Health grant AG030423 (JA).

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Statements

Author contributions

GB, GD, and JA designed the study. GB, SM, JS, and GD performed the experiments. GB, GD, JS, and SM analyzed the data. GB and GD drafted the manuscript and guarantors of the paper.

Acknowledgments

We would like to thank Morgan Wochinski and Mac Weninger for assisting with the daily exercise regimen of the animals.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    Abu-ErreishG. M.NeelyJ. R.WhitmerJ. T.WhitmanV.SanadiD. R. (1977). Fatty acid oxidation by isolated perfused working hearts of aged rats. Am. J. Physiol.232, E258E262.

  • 2

    AlessioH. M.AnsinelliH.ThreadgillC.HagermanA. E. (2014). Comparison of gene and protein expressions in rats residing in standard cages with those having access to an exercise wheel. Biomed Res. Int.2014:950516. 10.1155/2014/950516

  • 3

    BengtssonJ.GustafssonT.WidegrenU.JanssonE.SundbergC. J. (2001). Mitochondrial transcription factor, A., and respiratory complex IV increase in response to exercise training in humans. Pflügers Arch.443, 6166. 10.1007/s004240100628

  • 4

    BetikA. C.BakerD. J.KrauseD. J.McConkeyM. J.HeppleR. T. (2008). Exercise training in late middle-aged male Fischer 344 X brown norway F1-hybrid rats improves skeletal muscle aerobic function. Exp. Physiol.93, 863871. 10.1113/expphysiol.2008.042069

  • 5

    BhashyamS.ParikhP.BolukogluH.ShannonA. H.PorterJ. H.ShenY. T.et al. (2007). Aging is associated with myocardial insulin resistance and mitochondrial dysfunction. Am. J. Physiol. Heart Circ. Physiol.293, H3063H3071. 10.1152/ajpheart.00163.2007

  • 6

    BronikowskiA. M.CarterP. A.MorganT. J.GarlandT.Jr.UngN.PughT. D.et al. (2003). Lifelong voluntary exercise in the mouse prevents age-related alterations in gene expression in the heart. Physiol. Genomics12, 129138. 10.1152/physiolgenomics.00082.2002

  • 7

    BurelleY.WamboltR. B.GristM.ParsonsH. L.ChowJ. C.AntlerC.et al. (2004). Regular exercise is associated with a protective metabolic phenotype in the rat heart. Am. J. Physiol. Heart Circ. Physiol.287, H1055H1063. 10.1152/ajpheart.00925.2003

  • 8

    ChakravartiB.OsegueraM.DalalN.FathyP.MallikB.RavalA.et al. (2008). Proteomic profiling of aging in the mouse heart: altered expression of mitochondrial proteins. Arch. Biochem. Biophys.474, 2231. 10.1016/j.abb.2008.02.001

  • 9

    CovenD. L.HuX.CongL.BergeronR.ShulmanG. I.HardieD. G.et al. (2003). Physiological role of AMP-Activated Protein Kinase in the heart: graded activation during exercise. Am. J. Physiol. Endocrinol. Metab.285, E629E636. 10.1152/ajpendo.00171.2003

  • 10

    DobrzynP.PyrkowskaA.DudaM. K.BednarskiT.MaczewskiM.LangfortJ.et al. (2013). Expression of lipogenic genes is upregulated in the heart with exercise training-induced but not pressure overload-induced left ventricular hypertrophy. Am. J. Physiol. Endocrinol. Metab.304, E1348E1358. 10.1152/ajpendo.00603.2012

  • 11

    DolinskyV. W.DyckJ. R. (2006). Role of AMP-Activated Protein Kinase in healthy and diseased hearts. Am. J. Physiol. Heart Circ. Physiol.291, H2557H2569. 10.1152/ajpheart.00329.2006

  • 12

    FanninS. W.LesnefskyE. J.SlabeT. J.HassanM. O.HoppelC. L. (1999). Aging selectively decreases oxidative capacity in rat heart interfibrillar mitochondria. Arch. Biochem. Biophys.372, 399407. 10.1006/abbi.1999.1508

  • 13

    FortneyS.TankersleyC.LightfootJ. T.DrinkwaterD.ClulowJ.GerstenblithG.et al. (1992). Cardiovascular responses to lower body negative pressure in trained and untrained older men. J. Appl. Physiol. (Bethesda, Md. : 1985) 73, 26932700.

  • 14

    GoffartS.von Kleist-RetzowJ. C.WiesnerR. J. (2004). Regulation of mitochondrial proliferation in the heart: power-plant failure contributes to cardiac failure in hypertrophy. Cardiovasc. Res.64, 198207. 10.1016/j.cardiores.2004.06.030

  • 15

    GonzalezA. A.KumarR.MulliganJ. D.DavisA. J.SaupeK. W. (2004). Effects of aging on cardiac and skeletal muscle AMPK activity: basal activity, allosteric activation, and response to in vivo hypoxemia in mice. Am. J. Physiol. Regul. Integr. Comp. Physiol.287, R1270R1275. 10.1152/ajpregu.00409.2004

  • 16

    HallJ. L.MazzeoR. S.PodolinD. A.CarteeG. D.StanleyW. C. (1994). Exercise training does not compensate for age-related decrease in myocardial GLUT-4 content. J. Appl. Physiol. (Bethesda, Md. : 1985) 76, 328332.

  • 17

    Hassmann-PoznańskaE.TarantaA.BialukI.PoznańskaM.ZajączkiewiczH.WinnickaM. M. (2013). Analysis of gene expression profiles in tympanic membrane following perforation using PCR array in rats–preliminary investigation. Int. J. Pediatr. Otorhinolaryngol.77, 17531759. 10.1016/j.ijporl.2013.08.009

  • 18

    HawleyS. A.BoudeauJ.ReidJ. L.MustardK. J.UddL.MäkeläT. P.et al. (2003). Complexes between the LKB1 tumor suppressor, STRAD alpha/beta and MO25 alpha/beta are upstream kinases in the AMP-activated protein kinase cascade. J. Biol.2:28. 10.1186/1475-4924-2-28

  • 19

    HawleyS. A.PanD. A.MustardK. J.RossL.BainJ.EdelmanA. M.et al. (2005). Calmodulin-dependent protein kinase kinase-beta is an alternative upstream kinase for AMP-activated protein kinase. Cell Metab.2, 919. 10.1016/j.cmet.2005.05.009

  • 20

    HussJ. M.KellyD. P. (2004). Nuclear receptor signaling and cardiac energetics. Circ. Res.95, 568578. 10.1161/01.RES.0000141774.29937.e3

  • 21

    IemitsuM.MiyauchiT.MaedaS.SakaiS.FujiiN.MiyazakiH.et al. (2003). Cardiac hypertrophy by hypertension and exercise training exhibits different gene expression of enzymes in energy metabolism. Hypertens. Res.26, 829837. 10.1291/hypres.26.829

  • 22

    IemitsuM.MiyauchiT.MaedaS.TanabeT.TakanashiM.Irukayama-TomobeY.et al. (2002). Aging-induced decrease in the PPAR-alpha level in hearts is improved by exercise training. Am. J. Physiol. Heart Circ. Physiol.283, H1750H1760. 10.1152/ajpheart.01051.2001

  • 23

    JianB.YangS.ChenD.ZouL.ChathamJ. C.ChaudryI.et al. (2011). Aging influences cardiac mitochondrial gene expression and cardiovascular function following hemorrhage injury. Mol. Med. (Cambridge, Mass.) 17, 542549. 10.2119/molmed.2010.00195

  • 24

    KangC.ChungE.DiffeeG.JiL. L. (2013). Exercise training attenuates aging-associated mitochondrial dysfunction in rat skeletal muscle: role of PGC-1α. Exp. Gerontol.48, 13431350. 10.1016/j.exger.2013.08.004

  • 25

    KatesA. M.HerreroP.DenceC.SotoP.SrinivasanM.DelanoD. G.et al. (2003). Impact of aging on substrate metabolism by the human heart. J. Am. Coll. Cardiol.41, 293299. 10.1016/S0735-1097(02)02714-6

  • 26

    KumaranS.SubathraM.BaluM.PanneerselvamC. (2005). Supplementation of L-Carnitine improves mitochondrial enzymes in heart and skeletal muscle of aged rats. Exp. Aging Res.31, 5567. 10.1080/03610730590882846

  • 27

    LaiL.LeoneT. C.KellerM. P.MartinO. J.BromanA. T.NigroJ.et al. (2014). Energy metabolic reprogramming in the hypertrophied and early stage failing heart: a multisystems approach. Circul. Heart Fail.7, 10221031. 10.1161/CIRCHEARTFAILURE.114.001469

  • 28

    LakattaE. G.SollottS. J. (2002). Perspectives on mammalian cardiovascular aging: humans to molecules. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol.132, 699721. 10.1016/S1095-6433(02)00124-1

  • 29

    LarsenS.NielsenJ.HansenC. N.NielsenL. B.WibrandF.StrideN.et al. (2012). Biomarkers of mitochondrial content in skeletal muscle of healthy young human subjects. J. Physiol.590(Pt 14), 33493360. 10.1113/jphysiol.2012.230185

  • 30

    LeMoineC. M.McClellandG. B.LyonsC. N.Mathieu-CostelloO.MoyesC. D. (2006). Control of mitochondrial gene expression in the aging rat myocardium. Biochem. Cell Biol.84, 191198. 10.1139/o05-169

  • 31

    LuP.VogelC.WangR.YaoX.MarcotteE. M. (2007). Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat. Biotechnol.25, 117124. 10.1038/nbt1270

  • 32

    McMillinJ. B.TaffetG. E.TaegtmeyerH.HudsonE. K.TateC. A. (1993). Mitochondrial metabolism and substrate competition in the aging Fischer rat heart. Cardiovasc. Res. 27, 22222228.

  • 33

    MundayM. R.CampbellD. G.CarlingD.HardieD. G. (1988). Identification by amino acid sequencing of three major regulatory phosphorylation sites on rat acetyl-CoA carboxylase. Eur. J. Biochem.175, 331338. 10.1111/j.1432-1033.1988.tb14201.x

  • 34

    MurakamiT.ShimomuraY.FujitsukaN.SugiyamaS. (1995). Differential adaptation to endurance training between heart and gastrocnemius muscle mitochondria in rats. Biochem. Mol. Biol. Int.36, 285290.

  • 35

    MurrayA. J.AndersonR. E.WatsonG. C.RaddaG. K.ClarkeK. (2004). Uncoupling proteins in human heart. Lancet364, 17861788. 10.1016/S0140-6736(04)17402-3

  • 36

    NavarroA.BoverisA. (2007). The mitochondrial energy transduction system and the aging process. Am. J. Physiol. Cell Physiol.292, C670C686. 10.1152/ajpcell.00213.2006

  • 37

    NeuferP. D.DohmG. L. (1993). Exercise induces a transient increase in transcription of the GLUT-4 gene in skeletal muscle. Am. J. Physiol.265(6 Pt 1), C1597C1603.

  • 38

    OgawaT.SpinaR. J.MartinW. H.IIIKohrtW. M.SchechtmanK. B.HolloszyJ. O.et al. (1992). Effects of aging, sex, and physical training on cardiovascular responses to exercise. Circulation86, 494503. 10.1161/01.CIR.86.2.494

  • 39

    OkadaM.IshkitievN.YaegakiK.ImaiT.TanakaT.FukudaM.et al. (2014). Hydrogen sulfide increases hepatic differentiation of human tooth-pulp stem cells compared with human bone-marrow stem cells. Int. Endod. J.47, 11421150. 10.1111/iej.12262

  • 40

    OscaiL. B.MoléP. A.HolloszyJ. O. (1971). Effects of exercise on cardiac weight and mitochondria in male and female rats. Am. J. Physiol.220, 19441948. 10.1249/00005768-197100310-00057

  • 41

    PhaneufS.LeeuwenburghC. (2002). Cytochrome c release from mitochondria in the aging heart: a possible mechanism for apoptosis with age. Am. J. Physiol. Regul. Integr. Comp. Physiol.282, R423R430. 10.1152/ajpregu.00296.2001

  • 42

    PilegaardH.OrdwayG. A.SaltinB.NeuferP. D. (2000). Transcriptional regulation of gene expression in human skeletal muscle during recovery from exercise. Am. J. Physiol. Endocrinol. Metab.279, E806E814.

  • 43

    PrestonC. C.OberlinA. S.HolmuhamedovE. L.GuptaA.SagarS.SyedR. H.et al. (2008). Aging-induced alterations in gene transcripts and functional activity of mitochondrial oxidative phosphorylation complexes in the heart. Mech. Ageing Dev.129, 304312. 10.1016/j.mad.2008.02.010

  • 44

    RimbaudS.SanchezH.GarnierA.FortinD.BigardX.VekslerV.et al. (2009). Stimulus specific changes of energy metabolism in hypertrophied heart. J. Mol. Cell. Cardiol.46, 952959. 10.1016/j.yjmcc.2009.01.013

  • 45

    RinaldiB.DonniacuoM.SodanoL.GrittiG.SignorielloS.ParrettaE.et al. (2013). Effects of sildenafil on the gastrocnemius and cardiac muscles of rats in a model of prolonged moderate exercise training. PLoS ONE8:e69954. 10.1371/journal.pone.0069954

  • 46

    SampleJ.ClelandJ. G.SeymourA. M. (2006). Metabolic remodeling in the aging heart. J. Mol. Cell. Cardiol.40, 5663. 10.1016/j.yjmcc.2005.09.018

  • 47

    ScottJ. W.HawleyS. A.GreenK. A.AnisM.StewartG.ScullionG. A.et al. (2004). CBS domains form energy-sensing modules whose binding of adenosine ligands is disrupted by disease mutations. J. Clin. Invest.113, 274284. 10.1172/JCI19874

  • 48

    SealsD. R.HagbergJ. M.SpinaR. J.RogersM. A.SchechtmanK. B.EhsaniA. A. (1994). Enhanced left ventricular performance in endurance trained older men. Circulation89, 198205. 10.1161/01.CIR.89.1.198

  • 49

    ShankavaramU. T.ReinholdW. C.NishizukaS.MajorS.MoritaD.CharyK. K.et al. (2007). Transcript and protein expression profiles of the NCI-60 cancer cell panel: an integromic microarray study. Mol. Cancer Ther.6, 820832. 10.1158/1535-7163.MCT-06-0650

  • 50

    SiuP. M.DonleyD. A.BrynerR. W.AlwayS. E. (2003). Citrate synthase expression and enzyme activity after endurance training in cardiac and skeletal muscles. J. Appl. Physiol. (Bethesda, Md. : 1985) 94, 555560. 10.1152/japplphysiol.00821.2002

  • 51

    SpurgeonH. A.SteinbachM. F.LakattaE. G. (1983). Chronic exercise prevents characteristic age-related changes in rat cardiac contraction. Am. J. Physiol.244, H513H518.

  • 52

    StrattonJ. R.LevyW. C.CerqueiraM. D.SchwartzR. S.AbrassI. B. (1994). Cardiovascular responses to exercise. effects of aging and exercise training in healthy men. Circulation89, 16481655. 10.1161/01.CIR.89.4.1648

  • 53

    TurdiS.FanX.LiJ.ZhaoJ.HuffA. F.DuM.et al. (2010). AMP-activated protein kinase deficiency exacerbates aging-induced myocardial contractile dysfunction. Aging Cell9, 592606. 10.1111/j.1474-9726.2010.00586.x

  • 54

    Ventura-ClapierR.GarnierA.VekslerV. (2008). Transcriptional control of mitochondrial biogenesis: the central role of PGC-1alpha. Cardiovasc. Res.79, 208217. 10.1093/cvr/cvn098

  • 55

    WadleyG. D.LakerR. C.McConellG. K.WlodekM. E. (2016). Endurance training in early life results in long-term programming of heart mass in rats. Physiol. Rep.4:e12720. 10.14814/phy2.12720

  • 56

    WanagatJ.WolffM. R.AikenJ. M. (2002). Age-associated changes in function, structure and mitochondrial genetic and enzymatic abnormalities in the Fischer 344 X Brown Norway F(1) hybrid rat heart. J. Mol. Cell. Cardiol.34, 1728. 10.1006/jmcc.2001.1483

  • 57

    WisneskiJ. A.GertzE. W.NeeseR. A.GruenkeL. D.CraigJ. C. (1985). Dual carbon-labeled isotope experiments using D-[6-14C] glucose and L-[1,2,3-13C3] lactate: a new approach for investigating human myocardial metabolism during ischemia. J. Am. Coll. Cardiol.5, 11381146. 10.1016/S0735-1097(85)80016-4

  • 58

    WuD.LeeY. C.LiuH. C.YuanR. Y.ChiouH. Y.HungC. H.et al. (2013). Identification of TLR downstream pathways in stroke patients. Clin. Biochem.46, 10581064. 10.1016/j.clinbiochem.2013.05.059

  • 59

    ZhaoL.ZouX.FengZ.LuoC.LiuJ.LiH.et al. (2014). Evidence for association of mitochondrial metabolism alteration with lipid accumulation in aging rats. Exp. Gerontol.56, 312. 10.1016/j.exger.2014.02.001

  • 60

    ZhuY.SotoJ.AndersonB.RiehleC.ZhangY. C.WendeA. R.et al. (2013). Regulation of fatty acid metabolism by mTOR in adult murine hearts occurs independently of changes in PGC-1α. Am. J. Physiol. Heart Circul. Physiol.305, H41H51. 10.1152/ajpheart.00877.2012

  • 61

    ZonderlandM. L.BärP. R.ReijneveldJ. C.SpruijtB. M.KeizerH. A.GlatzJ. F. (1999). Different metabolic adaptation of heart and skeletal muscles to moderate-intensity treadmill training in the rat. Eur. J. Appl. Physiol. Occup. Physiol.79, 391396. 10.1007/s004210050527

Summary

Keywords

mitochondrial biogenesis, gene expression, exercise training, fatty acid metabolism, glucose metabolism

Citation

Barton GP, Sepe JJ, McKiernan SH, Aiken JM and Diffee GM (2016) Mitochondrial and Metabolic Gene Expression in the Aged Rat Heart. Front. Physiol. 7:352. doi: 10.3389/fphys.2016.00352

Received

16 March 2016

Accepted

02 August 2016

Published

23 August 2016

Volume

7 - 2016

Edited by

Leonardo Roever, Federal University of Uberlândia, Brazil

Reviewed by

Jop Van Berlo, University of Minnesota, USA; Min Zhang, King's College London, UK

Updates

Copyright

*Correspondence: Gary M. Diffee

This article was submitted to Striated Muscle Physiology, a section of the journal Frontiers in Physiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics