|
|
||||||||
COMPARATIVE AND EVOLUTIONARY PHYSIOLOGY
Department of Biology, Queen's University, Kingston, Ontario, Canada
Submitted 10 March 2004 ; accepted in final form 3 September 2004
| ABSTRACT |
|---|
|
|
|---|
citrate synthase; suborder Scombroidei; nuclear DNA content; posttranscriptional regulation
-coactivator 1
has emerged as a master regulator (32, 54), acting through a network of transcription factors to alter mitochondrial biogenesis in mammalian muscle fibers (13, 20, 33, 53). Although transcriptional regulation of mitochondrial enzymes is clearly important within an individual, it is not yet known whether it contributes to evolutionary variations in mitochondrial abundance that arise between species. In this study, we use scombroid fish (order Perciformes, suborder Scombroidei) to address the origins of variations in mitochondrial content among species that arise over evolutionary time and also among fiber types within individuals. Fish are good models in which to investigate fiber-type differences in mitochondrial properties because, in contrast to mammalian fiber types, their skeletal muscles are divided into distinct red1 (aerobic) and white (glycolytic) muscle fibers with large differences in mitochondrial content (29, 30). In particular, the scombroid fish have many properties that are conducive to the study of evolutionary and developmental variations in mitochondrial abundance. They show marked phenotypic variations in muscle oxidative capacity associated with aerobic activity and thermal physiology (14, 27). The tunas (family Scombridae) have impressive exercise physiology, including high metabolic rates, rapid swimming speeds, long migrations, rapid recovery from burst swimming, endothermy, numerous adaptations for oxygen uptake and delivery, and a unique mode of swimming (9). Also, tunas have a muscle mitochondrial capacity that exceeds that of most fish (39), including the closely related billfishes (family Xiphiidae and Istiophoridae) (6, 14, 27). Billfishes possess a heater organ, which is a modified muscle that warms the brain and eyes to preserve vision in cold water (4). Derived from the superior rectus eye muscle, the heater organ lacks most muscle contractile machinery and contains a very high mitochondrial content (5570% of cell volume) (46, 8). It is not known how this increased mitochondrial content is regulated in the heater tissue, although it is possible that the mechanisms are similar to those found between red and white muscles. This exceptional tissue may be unique to the suborder Scombroidei: an analogous tissue has been found only in the butterfly mackerel (Gasterochisma melampus) of the family Scombridae (7).
In this study, we assess the molecular basis for differences in muscle mitochondrial enzyme content found among fiber types (red muscle, white muscle, cardiac muscle, and heater organ) and also arising over evolutionary time, between clades of scombroid fishes (tunas vs. billfishes). Focusing on citrate synthase (CS) as a marker enzyme for mitochondrial content, we investigate how both genetic (enzyme structure, nuclear content) and regulatory (transcriptional, posttranscriptional) mechanisms influence mitochondrial content in these two contexts. Because CS is thought to be a transcriptionally regulated enzyme in mammalian skeletal muscle (1), we hypothesized that differences among species and tissues would be attributed to differences in CS mRNA levels that are ultimately linked to transcriptional regulators.
| MATERIALS AND METHODS |
|---|
|
|
|---|
RNA isolation and CS cDNA sequence information. We purified total RNA from guanidine thiocyanate extracts using a standard acid phenol protocol (12). Poly(A) RNA was isolated from total RNA using the poly(A) Purist MAG kit (Ambion).
All scombroid CS sequence data were determined from RT-PCR analyses. First-strand cDNA was synthesized from red muscle total RNA and used as the template for subsequent PCR reactions. Degenerate forward (F251) and reverse (R252) primers, based on CS mRNA sequences from human (GenBank accession number NM 004077), mouse (NM 026444), rat (NM 130755), pig (M21197), and zebrafish (BC045362.1), were used to amplify a 590-bp portion of CS cDNA at 50°C using standard PCR conditions (Table 1). Products were sequenced, and this information was used to design specific primers for use in 3' and 5' rapid amplification of cDNA ends (RACE) using the Roche 5'/3' RACE kit. Poly(A) RNA was used as the template for 5' RACE, and gene-specific reverse primers 5'-R1, 5'-R2, 5'-R3, and 5'-R4 were used in cDNA synthesis and subsequent nested amplifications (Table 1). Total RNA was used as the template for 3' RACE and reverse transcribed with an oligo(dT) primer. We performed nested PCR amplifications with this cDNA as the initial template using the specific forward primers 3'-F1, 3'-F2, and 3'-F3, respectively, each with the anchor primer (Table 1). The band from the second PCR for yellowfin and marlin and the third PCR for skipjack were excised and cloned.
|
CS sequence analyses. Scombroid CS sequences were threaded onto the crystal structure determined for pig CS complexed with the reaction products coenzyme A and citrate, listed under the Protein Data Bank (3) identification code 2CTS [PDB] (43). Cn3D 4.1 software was used to visualize these alignments (MMDB, NCBI). Catalytic and ligand binding residues were identified from structural studies (43, 44, 51) and site-directed mutagenesis experiments (16, 28) of pig and chicken CS.
Tissue analyses.
Enzyme activities for CS (EC 4.1.3.7
[EC]
) and cytochrome c-oxidase (COX; EC 1.9.3.1
[EC]
) were measured from whole cell extracts (40) at 25°C, using assays optimized to ensure that substrates were not limiting. A small volume of homogenate (50 µl) was digested in 5 vol of proteinase K digestion buffer (29), which was then assayed for nuclear DNA content by Picogreen (Molecular Probes), with a standard curve constructed from trout genomic DNA. The contribution of mitochondrial DNA to this measure of nuclear DNA content is minimal. Despite its high copy number, mitochondrial DNA represents <2% of total DNA in striated muscles due to its small size (
16,500 bp) relative to the nuclear genome (2, 36).
We measured total RNA levels per gram of tissue directly from the upper (aqueous) phase of a standard RNA extraction protocol using Ribogreen (Molecular Probes). This organic extraction with acid phenol ensures that contaminating DNA levels are low (<2.5% of total nucleic acid content as determined from a Picogreen assay; data not shown). A mock RNA extraction (no tissue added) was performed to obtain an equivalent aqueous phase to use with the standard curve. Both Picogreen and Ribogreen analyses were conducted with a SpectraMAX GEMINI XS spectrofluorometer (Molecular Devices).
Quantitative competitive RT-PCR.
CS mRNA levels were quantified using quantitative competitive (QC) RT-PCR with an RNA competitor template. The competitor template was constructed by using PCR to delete an internal portion of the endogenous template (11). Swordfish CS cDNA was amplified with the forward primer QC-delete (ATGAGGGGGATGAAGGGTCTGGTGTAAGCTTCCAGAGTGTCAGGAGTTGTTGCC), containing a HindIII restriction site, and the reverse primer QC-R1, yielding a product with a 50-bp deletion just after the end of the QC-F1 primer binding site in the CS gene product. The resulting product,
CS, was cloned into pCR 2.1 (Invitrogen), analyzed for orientation by restriction enzyme digestion with HindIII, and sequenced.
The plasmid containing
CS was linearized with BamHI, purified (QIAquick PCR purification kit, Qiagen), and used to generate sense strand
CS RNA using T7 RNA polymerase. The in vitro transcript was DNase treated (1 unit of DNase per µg DNA) and then purified with an RNeasy Mini kit (Qiagen). No contaminating DNA was detectable at the dilution to be used in QC-RT-PCR when tested in negative control PCRs (no RT). The purity and size of the
CS transcript were verified by gel electrophoresis, quantified with Ribogreen (Molecular Probes), diluted to 500 pg/µl (4.05 fmol/µl) in DEPC-treated water, and stored in single-use aliquots at 80°C.
Before the reverse transcription procedure, RNA samples from tissues were quantified by absorbance at 260 nm using a quartz plate treated with 1 M sodium hydroxide to destroy RNases. Aliquots were taken directly from the plate well to ensure accurate quantification. A constant amount of sample RNA (500 ng) and various amounts of competitor RNA transcript were added to each RT tube. Samples were reverse transcribed at 42°C for 1 h. The synthesized cDNA was then amplified by PCR as described previously (31), using strict consensus primers (QC-F1 and QC-R1, Table 1) for all five species at an annealing temperature of 65°C. An aliquot of
-dCTP32 was added to the reaction for the last five cycles. Only reactions with full-length products within 50% of competitor products (once corrected for deoxycytosine content) were used in quantification (47). Controls with no added RT, as well as a negative PCR control, were run in parallel to detect any genomic DNA contamination. None of these negative controls generated PCR products. The amplified products were size fractioned on a 12% polyacrylamide gel, and the signal strength was quantified with a storage phosphor screen (Molecular Dynamics) and analyzed with Imagequant software (Molecular Dynamics). We converted values from nanograms to femtomoles of CS mRNA using the molecular mass of the competitor transcript (8.10 fmol/ng).
Statistical analyses. For all parameters, significant differences (P < 0.05) among tissue types within a species and for a tissue type among species were detected with one-way ANOVAs and identified with the Tukey-Kramer test.
| RESULTS |
|---|
|
|
|---|
6- to 10-fold higher CS activity per gram of tissue than white muscle. The heater organ had the highest activities of all fiber types, with 15- to 30-fold more CS than billfish white muscle.
|
-helix) and are at the periphery of the protein, far from catalytic, ligand-binding, and subunit-interaction sites (Fig. 2). Amino acid hydrophobicity is also conserved in sites with clade-specific differences. Although no crystal structure is available for fish CS, the conservation of structure between mammals (Sus scrofa, Protein Data Bank no. 2CTS) and birds (Gallus gallus, Protein Data Bank no. 3CTS) suggests that the three-dimensional properties of scombroid CS should parallel those of other vertebrates (see Fig. 6).
|
|
|
2-fold higher than mean billfish values in red and white muscle and increased to
3.5-fold higher in tuna cardiac muscle (Fig. 3C). Each tuna species also had a significantly higher CS activity per milligram DNA than each billfish species in all homologous tissues (P < 0.05). Conversely, the large differences in CS activity seen among fiber types are greatly reduced when enzyme levels are expressed relative to nuclear content (Fig. 3B). For example, the
30-fold difference in enzyme activity between the swordfish heater organ and white muscle is reduced to a threefold difference when expressed per milligram DNA (Fig. 1A vs. Fig. 3B). In addition, the 6- to 10-fold difference in CS enzyme activity between the red and cardiac muscle compared with the white muscle is reduced to an approximate threefold difference (Fig. 1A vs. Fig. 3B). These results argue for the large role of constitutive expression in combination with differences in nuclear content in determining fiber-type mitochondrial phenotype, even within the exceptional heater organ. Role of transcriptional and posttranscriptional regulation. To pinpoint the step at which increased enzyme content per milligram DNA is mediated, we investigated the impact of variations in transcriptional control of the CS gene. CS mRNA levels were determined by QC-RT-PCR using primers designed in areas of exact consensus for all five species. CS transcript levels per gram total RNA calculated by QC-RT-PCR (Fig. 4B) were expressed per milligram DNA by correcting for both total DNA (Fig. 3A) and RNA (Fig. 4A) content of the tissue. Within a species, white muscle contained as much or more CS transcript per milligram DNA than the red, cardiac, and heater muscle in 10 of 12 comparisons (Fig. 4C) (significant, P < 0.05). Among clades, mean CS mRNA per milligram DNA values vary <15% between tuna and billfish cardiac muscle, and the billfish white and red muscle are found to have 40% and 3% more CS mRNA than tuna muscles, opposite to CS enzyme trends (Fig. 4D). Each billfish species also contained as much or more CS transcript per milligram DNA than each tuna species in all homologous muscles, with the exception of the swordfish red and cardiac muscles compared with the skipjack (significant, P < 0.05). Therefore, neither intertissue nor interspecies variations in CS enzyme activity (Fig. 3, B and C) are reflected in CS mRNA levels (Fig. 4, C and D), arguing against major differences in transcriptional regulation and suggesting that differences in CS activity found when DNA content is accounted arise from posttranscriptional regulation. Within a species, red and cardiac muscle produce approximately two- to threefold more CS enzyme per femtomole of CS mRNA, and the heater organ produces four- to eightfold more CS enzyme than white muscle (Fig. 5A). Between clades, mean CS enzyme production per femtomole CS mRNA for tunas is about two- to threefold greater than homologous billfish tissues (Fig. 5B). Each tuna species also produces significantly more CS enzyme per femtomole CS mRNA than each billfish species in all homologous muscles, with the exception of the bigeye and skipjack red muscle compared with the swordfish (P < 0.05).
|
|
In this study, we assessed the relative importance of gene dosage and transcriptional and posttranscriptional regulation in determining evolutionary and fiber-type differences in mitochondrial enzyme content in muscles of scombroid fish. Studies on mammalian muscle fiber types suggest that transcriptional regulators play an important role in controlling differences in mitochondrial enzymes between muscles and with physiological challenges (23, 32, 33, 37, 38, 46, 53, 54). However, the genetic bases for differences in mitochondrial content between animals are virtually unknown, although mitochondrial abundance is a heritable trait that can be selected for on population-level time scales (24).
We have compared the mitochondrial capacity of selected tunas and billfishes using CS, a citric acid cycle enzyme located in the mitochondrial matrix, as an index of mitochondrial content. The differences we found in CS activity among fiber types and species are consistent with previous studies (6, 15, 27, 50). Within each clade (i.e., tunas and billfishes), homologous muscles displayed similar CS activities (Fig. 1A). Between clades, tunas had about twofold higher CS activity than billfishes in each homologous muscle (Fig. 1B). Within each species, red muscle and cardiac muscle had similar enzyme activities, both
6- to 10-fold greater than white muscle (Fig. 1A). Billfish heater organ CS activities were about 3-fold greater vs. cardiac and red muscle and 15- to 30-fold higher than activities of white muscle of the same species (Fig. 1A).
CS activity is a good marker for mitochondrial content because it correlates with morphometric measures of mitochondria volume density (42) and milligrams of mitochondrial protein (39). In this study, we also found good correlation between the activities of CS and the inner mitochondrial membrane enzyme COX or complex IV of the electron transport chain (Fig. 1). Although COX is probably a more accurate predictor of maximal mitochondrial respiration, its complex structure and regulation make such functional and regulatory investigations difficult. The genetics and enzymology of CS are well suited to study the evolution of mitochondrial properties. Metazoan CS is a homodimer encoded by a single copy nuclear gene with no known isoforms (17), no allosteric regulators (44), and no known covalent regulators. Thus any difference in CS catalytic activity between tissues also reflects a difference in CS protein content. However, when comparing across species, it is important to assess the potential for genetic differences that could affect the relationship between CS catalytic activity and CS protein content. Changes in enzyme catalytic efficiency could, in principle, occur in CS or other mitochondrial enzymes as an evolutionary mechanism to increase catalytic activity. However, this is the least parsimonious mechanism by which a global increase in mitochondrial capacity could arise because of the complexity of mitochondrial metabolism. For example, many changes in primate COX subunits have structural implications that necessitate coevolution of binding partners (19). We sequenced the full-length coding region of CS from all five species and found no clade-specific mutations that affected functionally relevant sites in the protein. There are no variations in catalytic, substrate binding, or dimerization sites, arguing that the catalytic efficiency of CS is similar between the tunas and billfishes (Figs. 2 and 6). Thus we conclude that CS activity patterns reflect CS content in comparisons between these species as well as among tissues within each species. To assess the regulatory origins of differences in CS content among species and tissues, we considered the importance of CS gene dosage (nuclei per gram tissue), transcriptional regulation (CS mRNA per CS gene), and posttranscriptional regulation (CS activity per CS mRNA).
Importance of muscle DNA content. Myonuclear domain and muscle morphology are important considerations when comparing muscles in terms of gene expression and protein levels. Consider two muscles that differ in myonuclei per gram and therefore CS gene copy number per gram. The CS genes in both tissues could be regulated identically, but a tissue with more nuclei per gram (e.g., red muscle) would be expected to synthesize more CS per gram muscle. Because mitochondrial DNA contributes <2% to total DNA content, total DNA per gram tissue is a valid measure of nuclear DNA per gram tissue (2, 36). Nuclear DNA content should also be a valid measure of the CS gene content of a tissue if all nuclei within an individual are the same. This would not be the case in some tissues, such as polytenic salivary glands, where specific genes are amplified in the genome. This measure should also be valid for our evolutionary comparisons since haploid genome size is comparable among the tunas and billfishes (21, 22). Thus expressing CS per milligram DNA provides a better indication of potential differences in CS gene expression between tissues. Although this may be less important in mammals, where myonuclear domain differs only two- to threefold between skeletal fiber types (49), it might be very important in fish, which display a much greater divergence in metabolic properties (30). The scombroid fish are of particular interest because of the mitochondria-rich heater organ tissue.
In this study, fiber-type differences in CS activity parallel variations in nuclear DNA content. Red muscle and cardiac muscle have many more nuclei and mitochondria than does white muscle (Fig. 3A and Fig. 1A). The heater organ has both the highest CS activity (Fig. 1A) and highest DNA content (Fig. 3A) of all billfish tissues: expressing CS per milligram DNA reduces the 30-fold difference in CS per gram muscle between swordfish heater organ and white muscle to a 3-fold difference (Fig. 1A vs. Fig. 3B). With such a close relationship between nuclear DNA content and CS activity across fiber types, gene regulatory mechanisms (transcriptional, posttranscriptional) need to account for only a fraction of the differences arising between tissues. Instead, our study highlights the important role of constitutive expression in controlling mitochondrial enzyme content across fiber types. Leary et al. (30) found similar evidence for constitutive expression of the mitochondrial enzymes CS and COX in red, cardiac, and white muscles of another teleost fish, Oncorhynchus mykiss (30). Surprisingly, we found that even the mitochondria-rich billfish heater organ is not exceptional with respect to mitochondrial gene expression compared with other billfish muscles. Although our results support the idea that constitutive pathways are primarily responsible for variations in mitochondrial content across scombroid muscles, they do not negate the possibility of major differences in inducible/adaptive gene regulation during physiological challenges, such as exercise and cold acclimation.
Regulation of CS content.
Although gene dosage accounted for most of the differences between fiber types, it did not influence the differences in CS activities between tunas and billfishes (Fig. 3C). Most studies that assess the determinants of mitochondrial gene expression in mammalian striated muscles focus on transcriptional regulators, particularly peroxisome proliferator-activated receptor
-coactivator-1
(33, 37, 38, 46). However, it is important to recognize that many mitochondrial proteins are also regulated by posttranscriptional controls [e.g.,
-subunit of the F1-ATPase complex (26), mitoribosomal protein S12 (35), subunit VIII-L of COX (41), aconitase (18)]. We assessed the relative importance of transcriptional and posttranscriptional regulation in the differences that emerge after DNA content is considered: approximately twofold between clades and threefold among muscles.
After developing a QC-RT-PCR approach that took interspecies differences in CS sequence into consideration, we assessed the relationship between CS transcript levels and CS gene dosage. CS is transcriptionally regulated during the mitochondrial biogenesis that accompanies electric stimulation in mammalian striated muscle (1). Under steady-state conditions, CS mRNA correlates with CS enzyme levels (1). Thus transcriptional regulation is a potential site for evolutionary differences in CS content among species. For example, variations in enzyme levels could arise from mutations in 1) cis-regulatory regions of genes (altering transcription factor binding), 2) cis-regulatory regions of the transcriptional regulators, or 3) transcription factor coding regions, causing structural changes that influence DNA binding or catalytic properties (25, 45). Although transcriptional regulation is increasingly implicated in evolutionary variations, in other nonmitochondrial proteins (25, 45, 52), it is not yet known whether evolutionary differences in mitochondrial protein levels arise through effects on transcriptional regulators. If the differences in CS activity between species and tissues arose by transcriptional regulation, we would also predict a correlation between CS mRNA and CS enzyme levels. However, there was little correlation between CS activity and CS mRNA levels when compared between clades of fish (Fig. 3C vs. Fig. 4D, note differences in scale). Instead, it is the amount of CS enzyme produced per femtomole of CS mRNA, a measure of posttranscriptional regulation, that accounts for the majority of differences in CS enzyme content between tunas and billfishes (Fig. 5B). There is also very little correspondence between CS enzyme and CS mRNA values when compared across fiber types within an individual (Fig. 3B vs. Fig. 4C). Again, it is the amount of CS enzyme produced per femtomole of CS mRNA (Fig. 5A) that accounts for a higher CS enzyme content per milligram DNA in the red, cardiac, and heater muscles (Fig. 3B). These data argue that the differences seen in CS activity among scombroid muscles after DNA content is accounted for likely arise through variations in posttranscriptional processes, rather than transcriptional control.
Posttranscriptional regulation of CS during mitochondrial biogenesis also occurs in the livers of North Sea eelpout (Zoarces viviparous) during cold acclimation (34). The mechanism by which the posttranscriptional regulation of CS occurs in North Sea eelpout liver and in scombroid muscles, both within and among species, is unknown. It is important to note that both studies have measured steady-state levels of CS transcripts and not transcription per se. CS mRNA levels may also be influenced by mRNA stability. However, even if this is a contributing factor, it would not affect the calculations of CS enzyme produced per CS mRNA, and the importance of posttranscriptional controls remains.
As mentioned previously, several mitochondrial enzymes are known to experience characteristic differences in posttranscriptional processing. Furthermore, declines in muscle energetics during aging and disease can also arise through posttranscriptional defects [e.g., protein import (48)]. Disease models, involving loss of function, provide insight into possible sites of posttranscriptional regulation, but they may not be directly applicable to evolutionary or developmental comparisons. Tunas are exceptional athletes among fish; it is likely that they possess unusual regulatory features that enable them to maintain higher CS levels in relation to CS mRNA. The most likely points of differential regulation are enhanced translation initiation, protein stability, or mitochondrial stability.
This is the first work to examine the genetic determinants of evolutionary variations in muscle mitochondrial content. Surprisingly, we found evidence for a substantial role of posttranscriptional differences in determining CS activity patterns among species. Posttranscriptional processes also contribute to fiber-type differences in CS content within scombroid fishes. However, we identify nuclear DNA content as the main determinant of phenotypic differences in mitochondrial content among muscle fiber types. Elucidating the posttranscriptional control mechanisms by which CS is regulated in these two contexts may have broad implications for the regulation of mitochondrial content and will represent an important step in unraveling the possible molecular mechanisms by which physiologically important traits, such as mitochondrial content, may evolve.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 We use the descriptive terms red and white muscle in preference to tetrapod distinctions of slow-oxidative, fast-oxidative/glycolytic, and fast-glycolytic or myosin I, IIa, IIb, or IId/x. The use of the terms red and white avoids implicit links with specific myosin isoform profiles, which in many fish is complicated by early and recent genome duplications. ![]()
| REFERENCES |
|---|
|
|
|---|
) and mitochondrial function by MEF2 and HDAC5. Proc Natl Acad Sci USA 100: 17111716, 2003.
expression in muscle. Proc Natl Acad Sci USA 100: 71117116, 2003.
-F1-ATPase mRNA depends on the regulation of a protein that binds the 3' untranslated region of the mRNA. Mol Cell Biol 17: 52555268, 1997.[Abstract]
drives the formation of slow-twitch muscle fibres. Nature 418: 797801, 2002.[CrossRef][Medline]
This article has been cited by other articles:
![]() |
M. Lucassen, N. Koschnick, L. G. Eckerle, and H.-O. Portner Mitochondrial mechanisms of cold adaptation in cod (Gadus morhua L.) populations from different climatic zones J. Exp. Biol., July 1, 2006; 209(13): 2462 - 2471. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. St-Pierre MASTER OF MITOCHONDRIAL NUMBER J. Exp. Biol., February 1, 2005; 208(3): vii - vii. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |