Search

Zebrafish provide a translational model of human cardiac function. Their similar cardiac electrophysiology enables screening of human cardiac repolarization disorders, drug arrhythmogenicity, and novel antiarrhythmic therapeutics. However, while zebrafish cardiac repolarization is driven by delayed rectifier potassium channel current (I Kr ), the relative role of alternate channel transcripts is uncertain. While human ether-a-go-go-related-gene-1a (hERG1a) is the dominant transcript in humans, expression of the functionally distinct alternate transcript, hERG1b, modifies the electrophysiological and pharmacologic I Kr phenotype. Studies of zebrafish I Kr are frequently translated without consideration for the presence and impact of hERG1b in humans. Here, we performed phylogenetic analyses of all available KCNH genes from Actinopterygii (ray-finned fishes). Our findings confirmed zebrafish cardiac zkcnh6a as the paralog of human hERG1a (hKCNH2a), but also revealed evidence of a hERG1b (hKCNH2b)-like N-terminally truncated gene, zkcnh6b, in zebrafish. zkcnh6b is a teleost-specific variant that resulted from the 3R genome duplication. qRT-PCR showed dominant expression of zkcnh6a in zebrafish atrial and ventricular tissue, with low levels of zkcnh6b. Functional evaluation of zkcnh6b in a heterologous system showed no discernable function under the conditions tested, and no influence on zkcnh6a function during the zebrafish ventricular action potential. Our findings provide the first descriptions of the zkcnh6b gene, and show that, unlike in humans, zebrafish cardiac repolarization does not rely upon co-assembly of zERG1a/zERG1b. Given that hERG1b modifies I Kr function and drug binding in humans, our findings highlight the need for consideration when translating hERG variant effects and toxicological screens in zebrafish, which lack a functional hERG1b-equivalent gene., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

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.

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., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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.

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.

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.

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.

In this paper, we take an historical approach to reviewing research into the patterns of metabolic enzymes in muscle in relation to body size, focusing on mitochondrial enzymes. One of the first studies on allometric scaling of muscle enzymes was published in an early issue of this journal (George and Talesara, 1961 Comp. Biochem. Physiol. 3: 267-273). These researchers studied a number of locally available birds and a bat, measuring the activity of the mitochondrial enzyme succinate dehydrogenase in relation to body mass and muscle structure. Though the phenomenon of allometric scaling of metabolism was well recognized even 50 years earlier, this study was one of the first to explore the enzymatic underpinnings of the metabolic patterns in different animals. In this review, we begin by considering the George and Talesara study in the context of this early era in metabolic biochemistry and comparative physiology. We review subsequent studies in the last 50 years that continued the comparative analysis of enzyme patterns in relation to body size in diverse experimental models. This body of work identified a recurrent (though not ubiquitous) reciprocal relationship between oxidative and glycolytic enzymes. In the last 10 years, studies have focused on identifying the molecular mechanisms that determine the muscle metabolic enzyme phenotype., (Copyright 2010 Elsevier Inc. All rights reserved.)

In mammals, the peroxisome proliferator-activated receptor (PPAR) gamma 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 degrees 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 alpha and beta) and mitochondrial transcripts (citrate synthase), metabolic gene regulators including PPARalpha, PPARbeta and nuclear respiratory factor-1 (NRF-1). PGC-1beta was the better predictor of the NRF-1 axis, whereas PGC-1alpha was the better predictor of the PPAR axis (PPARalpha, PPARbeta, 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-1alpha 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-1alpha linked to fatty acid oxidation through PPARalpha, and PGC-1beta 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.