SUMMARYIn mammals, the peroxisome proliferator-activated receptor (PPAR) γ coactivator-1 (PGC-1) family members and their binding partners orchestrate remodelling in response to diverse challenges such as diet, temperature and exercise. In this study, we exposed goldfish to three temperatures (4, 20 and 35°C) and to three dietary regimes (food deprivation, low fat and high fat) and examined the changes in mitochondrial enzyme activities and transcript levels for metabolic enzymes and their genetic regulators in red muscle, white muscle, heart and liver. When all tissues and conditions were pooled, there were significant correlations between the mRNA for the PGC-1 coactivators (both α and β) and mitochondrial transcripts (citrate synthase), metabolic gene regulators including PPARα, PPARβ and nuclear respiratory factor-1 (NRF-1). PGC-1β was the better predictor of the NRF-1 axis, whereas PGC-1α was the better predictor of the PPAR axis (PPARα, PPARβ, medium chain acyl CoA dehydrogenase). In contrast to these intertissue/developmental patterns, the response of individual tissues to physiological stressors displayed no correlations between mRNA for PGC-1 family members and either the NRF-1 or PPAR axes. For example, in skeletal muscles, low temperature decreased PGC-1α transcript levels but increased mitochondrial enzyme activities (citrate synthase and cytochrome oxidase) and transcripts for COX IV and NRF-1. These results suggest that in goldfish, as in mammals, there is a regulatory relationship between (i) NRF-1 and mitochondrial gene expression and (ii) PPARs and fatty acid oxidation gene expression. In contrast to mammals, there is a divergence in the roles of the coactivators, with PGC-1α linked to fatty acid oxidation through PPARα, and PGC-1β with a more prominent role in mediating NRF-1-dependent control of mitochondrial gene expression, as well as distinctions between their respective roles in development and physiological responsiveness.

The zebrafish (Danio rerio) is an important organism as a model for understanding vertebrate cardiovascular development. However, little is known about adult ZF cardiac function and how contractile function changes to cope with fluctuations in ambient temperature. The goals of this study were to: 1) determine if high resolution echocardiography (HRE) in the presence of reduced cardiodepressant anesthetics could be used to accurately investigate the structural and functional properties of the ZF heart and 2) if the effect of ambient temperature changes both acutely and chronically could be determined non-invasively using HRE in vivo. Heart rate (HR) appears to be the critical factor in modifying cardiac output (CO) with ambient temperature fluctuation as it increases from 78 ± 5.9 bpm at 18°C to 162 ± 9.7 bpm at 28°C regardless of acclimation state (cold acclimated CA– 18°C; warm acclimated WA– 28°C). Stroke volume (SV) is highest when the ambient temperature matches the acclimation temperature, though this difference did not constitute a significant effect (CA 1.17 ± 0.15 μL at 18°C vs 1.06 ± 0.14 μl at 28°C; WA 1.10 ± 0.13 μL at 18°C vs 1.12 ± 0.12 μl at 28°C). The isovolumetric contraction time (IVCT) was significantly shorter in CA fish at 18°C. The CA group showed improved systolic function at 18°C in comparison to the WA group with significant increases in both ejection fraction and fractional shortening and decreases in IVCT. The decreased early peak (E) velocity and early peak velocity / atrial peak velocity (E/A) ratio in the CA group are likely associated with increased reliance on atrial contraction for ventricular filling.

The teleost-specific whole genome duplication created multiple copies of genes allowing for subfunctionalization of isoforms. In this study, we show that the teleost cardiac Ca2+-binding troponin C (TnC) is the product of two distinct genes: cardiac TnC (cTnC, TnnC1a) and a fish-specific slow skeletal TnC (ssTnC, TnnC1b). The ssTnC gene is novel to teleosts as mammals have a single gene commonly referred as cTnC but which is also expressed in slow skeletal muscle. In teleosts, the data strongly indicate that these are two TnC genes are different paralogs. Because we determined that ssTnC exists across many teleosts but not in basal ray-finned fish (e.g., bichir), we propose that these paralogs are the result of an ancestral tandem gene duplication persisting only in teleosts. Quantification of mRNA levels was used to demonstrate distinct expression localization patterns of the paralogs within the chambers of the heart. In the adult zebrafish acclimated at 28°C, ssTnC mRNA levels are twofold greater than cTnC mRNA levels in the atrium, whereas cTnC mRNA was almost exclusively expressed in the ventricle. Meanwhile, rainbow trout acclimated at 5°C showed cTnC mRNA levels in both chambers significantly greater than ssTnC. Distinct responses to temperature acclimation were also quantified in both adult zebrafish and rainbow trout, with mRNA in both chambers shifting to express higher levels of cTnC in 18°C acclimated zebrafish and 5°C acclimated trout. Possible subfunctionalization of TnC isoforms may provide insight into how teleosts achieve physiological versatility in chamber-specific contractile properties.

Zebrafish (Danio rerio) are widely used as vertebrate model in developmental genetics and functional genomics as well as in cardiac structure-function studies. The zebrafish heart has been increasingly used as a model of human cardiac function, in part, due to the similarities in heart rate and action potential duration and morphology with respect to humans. The teleostian zebrafish is in many ways a compelling model of human cardiac function due to the clarity afforded by its ease of genetic manipulation, the wealth of developmental biological information, and inherent suitability to a variety of experimental techniques. However, in addition to the numerous advantages of the zebrafish system are also caveats related to gene duplication (resulting in paralogs not present in human or other mammals) and fundamental differences in how zebrafish hearts function. In this review, we discuss the use of zebrafish as a cardiac function model through the use of techniques such as echocardiography, optical mapping, electrocardiography, molecular investigations of excitation-contraction coupling, and their physiological implications relative to that of the human heart. While some of these techniques (e.g., echocardiography) are particularly challenging in the zebrafish because of diminutive size of the heart (~1.5 mm in diameter) critical information can be derived from these approaches and are discussed in detail in this article.

Zebrafish, as a model for teleost fish, have two paralogous troponin C (TnC) genes that are expressed in the heart differentially in response to temperature acclimation. Upon Ca(2+) binding, TnC changes conformation and exposes a hydrophobic patch that interacts with troponin I and initiates cardiac muscle contraction. Teleost-specific TnC paralogs have not yet been functionally characterized. In this study we have modeled the structures of the paralogs using molecular dynamics simulations at 18°C and 28°C and calculated the different Ca(2+)-binding properties between the teleost cardiac (cTnC or TnC1a) and slow-skeletal (ssTnC or TnC1b) paralogs through potential-of-mean-force calculations. These values are compared with thermodynamic binding properties obtained through isothermal titration calorimetry (ITC). The modeled structures of each of the paralogs are similar at each temperature, with the exception of helix C, which flanks the Ca(2+) binding site; this region is also home to paralog-specific sequence substitutions that we predict have an influence on protein function. The short timescale of the potential-of-mean-force calculation precludes the inclusion of the conformational change on the ΔG of Ca(2+) interaction, whereas the ITC analysis includes the Ca(2+) binding and conformational change of the TnC molecule. ITC analysis has revealed that ssTnC has higher Ca(2+) affinity than cTnC for Ca(2+) overall, whereas each of the paralogs has increased affinity at 28°C compared to 18°C. Microsecond-timescale simulations have calculated that the cTnC paralog transitions from the closed to the open state more readily than the ssTnC paralog, an unfavorable transition that would decrease the ITC-derived Ca(2+) affinity while simultaneously increasing the Ca(2+) sensitivity of the myofilament. We propose that the preferential expression of cTnC at lower temperatures increases myofilament Ca(2+) sensitivity by this mechanism, despite the lower Ca(2+) affinity that we have measured by ITC.

Gene duplication results in extra copies of genes that must coevolve with their interacting partners in multimeric protein complexes. The cardiac troponin (Tn) complex, containing TnC, TnI, and TnT, forms a distinct functional unit critical for the regulation of cardiac muscle contraction. In teleost fish, the function of the Tn complex is modified by the consequences of differential expression of paralogs in response to environmental thermal challenges. In this article, we focus on the interaction between TnI and TnC, coded for by genes that have independent evolutionary origins, but the co-operation of their protein products has necessitated coevolution. In this study, we characterize functional divergence of TnC and TnI paralogs, specifically the interrelated roles of regulatory subfunctionalization and structural subfunctionalization. We determined that differential paralog transcript expression in response to temperature acclimation results in three combinations of TnC and TnI in the zebrafish heart: TnC1a/TnI1.1, TnC1b/TnI1.1, and TnC1a/TnI1.5. Phylogenetic analysis of these highly conserved proteins identified functionally divergent residues in TnI and TnC. The structural and functional effect of these Tn combinations was modeled with molecular dynamics simulation to link divergent sites to changes in interaction strength. Functional divergence in TnI and TnC were not limited to the residues involved with TnC/TnI switch interaction, which emphasizes the complex nature of Tn function. Patterns in domain-specific divergent selection and interaction energies suggest that substitutions in the TnI switch region are crucial to modifying TnI/TnC function to maintain cardiac contraction with temperature changes. This integrative approach introduces Tn as a model of functional divergence that guides the coevolution of interacting proteins.

There is significant interest in the potential utility of small-molecule activator compounds to mitigate cardiac arrhythmia caused by loss of function of hERG1a voltage-gated potassium channels. Zebrafish ( Danio rerio) have been proposed as a cost-effective, high-throughput drug-screening model to identify compounds that cause hERG1a dysfunction. However, there are no reports on the effects of hERG1a activator compounds in zebrafish and consequently on the utility of the model to screen for potential gain-of-function therapeutics. Here, we examined the effects of hERG1a blocker and types 1 and 2 activator compounds on isolated zkcnh6a (zERG3) channels in the Xenopus oocyte expression system as well as action potentials recorded from ex vivo adult zebrafish whole hearts using optical mapping. Our functional data from isolated zkcnh6a channels show that under the conditions tested, these channels are blocked by hERG1a channel blockers (dofetilide and terfenadine), and activated by type 1 (RPR260243) and type 2 (NS1643, PD-118057) hERG1a activators with higher affinity than hKCNH2a channels (except NS1643), with differences accounted for by different biophysical properties in the two channels. In ex vivo zebrafish whole hearts, two of the three hERG1a activators examined caused abbreviation of the action potential duration (APD), whereas hERG1a blockers caused APD prolongation. These data represent, to our knowledge, the first pharmacological characterization of isolated zkcnh6a channels and the first assessment of hERG enhancing therapeutics in zebrafish. Our findings lead us to suggest that the zebrafish ex vivo whole heart model serves as a valuable tool in the screening of hKCNH2a blocker and activator compounds.

No abstract