Sensing and responding to energetic stress: The role of the AMPK-PGC1α-NRF1 axis in control of mitochondrial biogenesis in fish.

Remodeling the muscle metabolic machinery in mammals in response to energetic challenges depends on the energy sensor AMP-activated protein kinase (AMPK) and its ability to phosphorylate PPAR γ coactivator 1 α (PGC1α), which in turn coactivates metabolic genes through direct and indirect association with DNA-binding proteins such as the nuclear respiratory factor 1 (NRF1) (Wu et al., 1999). The integrity of this axis in fish is uncertain because PGC1α i) lacks the critical Thr177 targeted by AMPK and ii) has mutations that may preclude binding NRF1. In this study we found no evidence that AMPK regulates mitochondrial gene expression through PGC1α in zebrafish and goldfish. AICAR treatment of zebrafish blastula cells increased phosphorylation of AMPK and led to changes in transcript levels of the AMPK targets mTOR and hexokinase 2. However, we saw no increases in mRNA levels for genes associated with mitochondrial biogenesis, including PGC1α, NRF1, and COX7C, a cytochrome c oxidase subunit. Further, AMPK phosphorylated mammalian peptides of PGC1α but not the corresponding region of zebrafish or goldfish in vitro. In vivo cold acclimation of goldfish caused an increase in mitochondrial enzymes, AMP and ADP levels, however AMPK phosphorylation decreased. In fish, the NRF1-PGC1α axis may be disrupted due to insertions in fish PGC1α orthologs within the region that serves as NRF1 binding domain in mammals. Immunocopurification showed that recombinant NRF1 protein binds mammalian but not fish PGC1α. Collectively, our studies suggest that fish have a disruption in the AMPK-PGC1α-NRF1 pathway due to structural differences between fish and mammalian PGC1α.

Evolution of the oxygen sensitivity of cytochrome c oxidase subunit 4.

Vertebrates possess two paralogs of cytochrome c oxidase (COX) subunit 4: a ubiquitous COX4-1 and a hypoxia-linked COX4-2. Mammalian COX4-2 is thought to have a role in relation to fine-tuning metabolism in low oxygen levels, conferred through both structural differences in the subunit protein structure and regulatory differences in the gene. We sought to elucidate the pervasiveness of this feature across vertebrates. The ratio of COX4-2/4-1 mRNA is generally low in mammals, but this ratio was higher in fish and reptiles, particularly turtles. The COX4-2 gene appeared unresponsive to low oxygen in nonmammalian models (zebrafish, goldfish, tilapia, anoles, and turtles) and fish cell lines. Reporter genes constructed from the amphibian and reptile homologues of the mammalian oxygen-responsive elements and hypoxia-responsive elements did not respond to low oxygen. Unlike the rodent ortholog, the promoter of goldfish COX4-2 did not respond to hypoxia or anoxia. The protein sequences of the COX4-2 peptide showed that the disulfide bridge seen in human and rodent orthologs would be precluded in other mammalian lineages and lower vertebrates, all of which lack the requisite pair of cysteines. The coordinating ligands of the ATP-binding site are largely conserved across mammals and reptiles, but in Xenopus and fish, sequence variations may disrupt the ability of the protein to bind ATP at this site. Collectively, these results suggest that many of the genetic and structural features of COX4-2 that impart responsiveness and benefits in hypoxia may be restricted to the Euarchontoglires lineage that includes primates, lagomorphs, and rodents.

Origins of interspecies variation in mammalian muscle metabolic enzymes.

Do the transcriptional mechanisms that control an individual's mitochondrial content, PGC1α (peroxisome proliferator-activated receptor γ coactivator-1α) and NRF1 (nuclear respiratory factor-1), also cause differences between species? We explored the determinants of cytochrome c oxidase (COX) activities in muscles from 12 rodents differing 1,000-fold in mass. Hindlimb muscles differed in scaling patterns from isometric (soleus, gastrocnemius) to allometric (tibialis anterior, scaling coefficient = -0.16). Consideration of myonuclear domain reduced the differences within species, but interspecies differences remained. For tibialis anterior, there was no significant scaling relationship in mRNA/g for COX4-1, PGC1α, or NRF1, yet COX4-1 mRNA/g was a good predictor of COX activity (r(2) = 0.55), PGC1α and NRF1 mRNA correlated with each other (r(2) = 0.42), and both could predict COX4-1 mRNA (r(2) = 0.48 and 0.52) and COX activity (r(2) = 0.55 and 0.49). This paradox was resolved by multivariate analysis, which explained 90% of interspecies variation, about equally partitioned between mass effects and PGC1α (or NRF1) mRNA levels, independent of mass. To explore the determinants of PGC1α mRNA, we analyzed 52 mammalian PGC1α proximal promoters and found no size dependence in regulatory element distribution. Likewise, the activity of PGC1α promoter reporter genes from 30 mammals showed no significant relationship with body mass. Collectively, these studies suggest that not all muscles scale equivalently, but for those that show allometric scaling, transcriptional regulation of the master regulators, PGC1α and NRF1, does not account for scaling patterns, though it does contribute to interspecies differences in COX activities independent of mass.

Evolution of the nuclear-encoded cytochrome oxidase subunits in vertebrates.

Vertebrate mitochondrial cytochrome c oxidase (COX) possesses 10 nuclear-encoded subunits. Six subunits have paralogs in mammals, but the origins and distribution of isoforms among vertebrates have not been analyzed. We used Bayesian phylogenetic analysis to interpret the origins of each subunit, inferring the roles of gene and genome duplications. The paralogous ancestries of five genes were identical throughout the major vertebrate taxa: no paralogs of COX6c and COX7c, two paralogs of COX4 and COX6a, and three paralogs of COX7a. Two genes had an extra copy in teleosts (COX5a, COX5b), and three genes had additional copies in mammals (COX6b, COX7b, COX8). Focusing on early vertebrates, we examined structural divergence and explored transcriptional profiles across zebrafish tissues. Quantitative transcript profiles revealed dramatic differences in transcript abundance for different subunits. COX7b and COX4 transcripts were typically present at very low levels, whereas COX5a and COX8 were in vast excess in all tissues. For genes with paralogs, two general patterns emerged. For COX5a and COX8, there was ubiquitous expression of one paralog, with the other paralog in lower abundance in all tissues. COX4 and COX6a shared a distinct expression pattern, with one paralog dominant in brain and gills and the other in muscles. The isoform profiles in combination with phylogenetic analyses show that vertebrate COX isoform patterns are consistent with the hypothesis that early whole genome duplications in basal vertebrates governed the isoform repertoire in modern fish and tetrapods, though more recent lineage-specific gene/genome duplications also play a role in select subunits.