Pharmacological characterization of ionic currents that regulate the pacemaker rhythm in a weakly electric fish

Dynamics of a neuronal pacemaker in the weakly electric fish Apteronotus

The precise timing of neuronal activity is critical for normal brain function. In weakly electric fish, the medullary pacemaker network (PN) sets the timing for an oscillating electric organ discharge (EOD) used for electric sensing. This network is the most precise biological oscillator known, with sub-microsecond variation in oscillator period. The PN consists of two principle sets of neurons, pacemaker and relay cells, that are connected by gap junctions and normally fire in synchrony, one-to-one with each EOD cycle. However, the degree of gap junctional connectivity between these cells appears insufficient to provide the population averaging required for the observed temporal precision of the EOD. This has led to the hypothesis that individual cells themselves fire with high precision, but little is known about the oscillatory dynamics of these pacemaker cells. As a first step towards testing this hypothesis, we have developed a biophysical model of a pacemaker neuron action potential based on experimental recordings. We validated the model by comparing the changes in oscillatory dynamics produced by different experimental manipulations. Our results suggest that this relatively simple model can capture a large range of channel dynamics exhibited by pacemaker cells, and will thus provide a basis for future work on network synchrony and precision.

Neuronal Dynamics Underlying Communication Signals in a Weakly Electric Fish: Implications for Connectivity in a Pacemaker Network

Neuronal networks can produce stable oscillations and synchrony that are under tight control yet flexible enough to rapidly switch between dynamical states. The pacemaker nucleus in the weakly electric fish comprises a network of electrically coupled neurons that fire synchronously at high frequency. This activity sets the timing for an oscillating electric organ discharge with the lowest cycle-to-cycle variability of all known biological oscillators. Despite this high temporal precision, pacemaker activity is behaviorally modulated on millisecond time-scales for the generation of electrocommunication signals. The network mechanisms that allow for this combination of stability and flexibility are not well understood. In this study, we use an in vitro pacemaker preparation from Apteronotus leptorhynchus to characterize the neural responses elicited by the synaptic inputs underlying electrocommunication. These responses involve a variable increase in firing frequency and a prominent desynchronization of neurons that recovers within 5 oscillation cycles. Using a previously developed computational model of the pacemaker network, we show that the frequency changes and rapid resynchronization observed experimentally are most easily explained when model neurons are interconnected more densely and with higher coupling strengths than suggested by published data. We suggest that the pacemaker network achieves both stability and flexibility by balancing coupling strength with interconnectivity and that variation in these network features may provide a substrate for species-specific evolution of electrocommunication signals.

Dynamic Neuron-Glia Interactions in an Oscillatory Network Controlling Behavioral Plasticity in the Weakly Electric Fish, Apteronotus leptorhynchus

The involvement of glial cells in the regulation of physiological functions is being increasingly recognized, yet their role in plasticity of neural oscillators has remained largely elusive. An excellent model system to address the latter function is the pacemaker nucleus of the weakly electric fish, Apteronotus leptorhynchus. This brainstem oscillator drives the fish's electric organ discharge in a one-to-one fashion, with median frequencies of 880 Hz in males and 740 Hz in females. Morphometric analysis of the pacemaker nucleus has shown that astrocytes outnumber mature neurons seven-fold, and oscillator neurons even 200-fold. A similar dominance of astrocytes occurs among the adult-born cells that differentiate into glia and neurons. The astrocytes form a dense meshwork of cells interconnected by gap junctions. The degree of association of astrocytic fibers with the neural oscillator cells, and the gap-junction coupling between individual astrocytes, exhibit a sexual dimorphism, which parallels the sexual dimorphisms in the output frequency of the pacemaker nucleus, and ultimately in the electric organ discharge of the fish. It is hypothesized that the dynamics in astroglial structure mediate differences in the capacity to buffer potassium, which increases during the generation of action potentials. These differences, in turn, affect the excitability of the neural oscillator cells, and thus the output frequency of the pacemaker nucleus. Comparison of the pacemaker nucleus with other brain oscillators suggests that modulation of the output activity is one of the chief functions of the interaction of glia with the neural oscillator cells.

Genes linked to species diversity in a sexually dimorphic communication signal in electric fish

Sexually dimorphic behaviors are often regulated by androgens and estrogens. Steroid receptors and metabolism are control points for evolutionary changes in sexual dimorphism. Electric communication signals of South American knifefishes are a model for understanding the evolution and physiology of sexually dimorphic behavior. These signals are regulated by gonadal steroids and controlled by a simple neural circuit. Sexual dimorphism of the signals varies across species. We used transcriptomics to examine mechanisms for sex differences in electric organ discharges (EODs) of two closely related species, Apteronotus leptorhynchus and Apteronotus albifrons, with reversed sexual dimorphism in their EODs. The pacemaker nucleus (Pn), which controls EOD frequency (EODf), expressed transcripts for steroid receptors and metabolizing enzymes, including androgen receptors, estrogen receptors, aromatase, and 5α-reductase. The Pn expressed mRNA for ion channels likely to regulate the high-frequency activity of Pn neurons and for neuromodulator and neurotransmitter receptors that may regulate EOD modulations used in aggression and courtship. Expression of several ion channel genes, including those for Kir3.1 inward-rectifying potassium channels and sodium channel β1 subunits, was sex-biased or correlated with EODf in ways consistent with EODf sex differences. Our findings provide a basis for future studies to characterize neurogenomic mechanisms by which sex differences evolve.

The central nervous system transcriptome of the weakly electric brown ghost knifefish (Apteronotus leptorhynchus): de novo assembly, annotation, and proteomics validation

Background

The brown ghost knifefish (Apteronotus leptorhynchus) is a weakly electric teleost fish of particular interest as a versatile model system for a variety of research areas in neuroscience and biology. The comprehensive information available on the neurophysiology and neuroanatomy of this organism has enabled significant advances in such areas as the study of the neural basis of behavior, the development of adult-born neurons in the central nervous system and their involvement in the regeneration of nervous tissue, as well as brain aging and senescence. Despite substantial scientific interest in this species, no genomic resources are currently available.

Results

Here, we report the de novo assembly and annotation of the A. leptorhynchus transcriptome. After evaluating several trimming and transcript reconstruction strategies, de novo assembly using Trinity uncovered 42,459 unique contigs containing at least a partial protein-coding sequence based on alignment to a reference set of known Actinopterygii sequences. As many as 11,847 of these contigs contained full or near-full length protein sequences, providing broad coverage of the proteome. A variety of non-coding RNA sequences were also identified and annotated, including conserved long intergenic non-coding RNA and other long non-coding RNA observed previously to be expressed in adult zebrafish (Danio rerio) brain, as well as a variety of miRNA, snRNA, and snoRNA. Shotgun proteomics confirmed translation of open reading frames from over 2,000 transcripts, including alternative splice variants. Assignment of tandem mass spectra was greatly improved by use of the assembly compared to databases of sequences from closely related organisms. The assembly and raw reads have been deposited at DDBJ/EMBL/GenBank under the accession number GBKR00000000. Tandem mass spectrometry data is available via ProteomeXchange with identifier PXD001285.

Conclusions

Presented here is the first release of an annotated de novo transcriptome assembly from Apteronotus leptorhynchus, providing a broad overview of RNA expressed in central nervous system tissue. The assembly, which includes substantial coverage of a wide variety of both protein coding and non-coding transcripts, will allow the development of better tools to understand the mechanisms underlying unique characteristics of the knifefish model system, such as their tremendous regenerative capacity and negligible brain senescence.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1354-2) contains supplementary material, which is available to authorized users.

Large-scale identification of proteins involved in the development of a sexually dimorphic behavior

Sexually dimorphic behaviors develop under the influence of sex steroids, which induce reversible changes in the underlying neural network of the brain. However, little is known about the proteins that mediate these activational effects of sex steroids. Here, we used a proteomics approach for large-scale identification of proteins involved in the development of a sexually dimorphic behavior, the electric organ discharge of brown ghost knifefish, Apteronotus leptorhynchus. In this weakly electric fish, the discharge frequency is controlled by the medullary pacemaker nucleus and is higher in males than in females. After lowering the discharge frequency by chronic administration of β-estradiol, 2-dimensional difference gel electrophoresis revealed 62 proteins spots in tissue samples from the pacemaker nucleus that exhibited significant changes in abundance of >1.5-fold. The 20 identified protein spots indicated, among others, a potential involvement of astrocytes in the establishment of the behavioral dimorphism. Indeed, immunohistochemical analysis demonstrated higher expression of the astrocytic marker protein GFAP and increased gap-junction coupling between astrocytes in females compared with males. We hypothesize that changes in the size of the glial syncytium, glial coupling, and/or number of glia-specific potassium channels lead to alterations in the firing frequency of the pacemaker nucleus via a mechanism mediating the uptake of extracellular potassium ions from the extracellular space.

Evolution and hormonal regulation of sex differences in the electrocommunication behavior of ghost knifefishes (Apteronotidae)

The ghost knifefishes (family Apteronotidae) are one of the most successful and diverse families of electric fish. Like other weakly electric fish, apteronotids produce electric organ discharges (EODs) that function in electrolocation and communication. This review highlights the diversity in the structure, function and sexual dimorphism of electrocommunication signals within and across apteronotid species. EOD frequency (EODf) and waveform vary as a function of species, sex and/or social rank. Sex differences in EODf are evolutionarily labile; apteronotid species express every pattern of sexual dimorphism in EODf (males>females; males<females; males=females). The direction and magnitude of sex differences in EODf are correlated across species and populations with the responsiveness of EODf to androgens and/or estrogens, which suggests that sex differences evolve through gains and/or losses of hormone sensitivity. During social interactions, apteronotids also modulate their EODs to produce motivational signals known as chirps. Chirp structure differs markedly across species, and many species produce two or more discrete chirp types with potentially different functions. The structure of chirps is sexually dimorphic in all apteronotid species, and chirping is influenced by gonadal steroids and by neuromodulators. Encoding of chirps by the electrosensory system depends on the social context created by the interactions of the EODs of signalers and receivers. Electrosensory systems may thus influence the evolution of signal structure and function, and neuromodulators may coordinately shape the production and reception of electrocommunication signals depending on social context.

The neuroethology of electrocommunication: how signal background influences sensory encoding and behaviour in Apteronotus leptorhynchus

Weakly-electric fish are a well-established model system for neuroethological studies on communication and aggression. Sensory encoding of their electric communication signals, as well as behavioural responses to these signals, have been investigated in great detail under laboratory conditions. In the wave-type brown ghost knifefish, Apteronotus leptorhynchus, transient increases in the frequency of the generated electric field, called chirps, are particularly well-studied, since they can be readily evoked by stimulating a fish with artificial signals mimicking conspecifics. When two fish interact, both their quasi-sinusoidal electric fields (called electric organ discharge, EOD) superimpose, resulting in a beat, an amplitude modulation at the frequency difference between the two EODs. Although chirps themselves are highly stereotyped signals, the shape of the amplitude modulation resulting from a chirp superimposed on a beat background depends on a number of parameters, such as the beat frequency, modulation depth, and beat phase at which the chirp is emitted. Here we review the influence of these beat parameters on chirp encoding in the three primary stages of the electrosensory pathway: electroreceptor afferents, the hindbrain electrosensory lateral line lobe, and midbrain torus semicircularis. We then examine the role of these parameters, which represent specific features of various social contexts, on the behavioural responses of A. leptorhynchus. Some aspects of the behaviour may be explained by the coding properties of early sensory neurons to chirp stimuli. However, the complexity and diversity of behavioural responses to chirps in the context of different background parameters cannot be explained solely on the basis of the sensory responses and thus suggest that critical roles are played by higher processing stages.

Calcium dynamics encode the magnitude of a graded memory underlying sensorimotor adaptation

The role of Ca(2+) in the induction of neural correlates of memory has frequently been described in binary terms despite the fact that many forms of memory are graded in their strength and/or persistence. We find that Ca(2+) dynamics encode the magnitude of sensorimotor adaptation of the electromotor output in a weakly electric fish. The neural correlate of this memory is a synaptically induced Ca(2+)-dependent enhancement of intrinsic excitability of neurons responsible for setting the electromotor output. Changes in Ca(2+) during induction accurately predict the magnitude of this graded memory over a wide range of stimuli. Thus despite operating over a range from seconds to tens of minutes, the encoding of graded memory can be mediated by a relatively simple cellular mechanism.

Phylogenetic comparative analysis of electric communication signals in ghost knifefishes (Gymnotiformes: Apteronotidae)

Electrocommunication signals in electric fish are diverse, easily recorded and have well-characterized neural control. Two signal features, the frequency and waveform of the electric organ discharge (EOD), vary widely across species. Modulations of the EOD (i.e. chirps and gradual frequency rises) also function as active communication signals during social interactions, but they have been studied in relatively few species. We compared the electrocommunication signals of 13 species in the largest gymnotiform family,Apteronotidae. Playback stimuli were used to elicit chirps and rises. We analyzed EOD frequency and waveform and the production and structure of chirps and rises. Species diversity in these signals was characterized with discriminant function analyses, and correlations between signal parameters were tested with phylogenetic comparative methods. Signals varied markedly across species and even between congeners and populations of the same species. Chirps and EODs were particularly evolutionarily labile, whereas rises differed little across species. Although all chirp parameters contributed to species differences in these signals, chirp amplitude modulation, frequency modulation (FM) and duration were particularly diverse. Within this diversity,however, interspecific correlations between chirp parameters suggest that mechanistic trade-offs may shape some aspects of signal evolution. In particular, a consistent trade-off between FM and EOD amplitude during chirps is likely to have influenced the evolution of chirp structure. These patterns suggest that functional or mechanistic linkages between signal parameters(e.g. the inability of electromotor neurons increase their firing rates without a loss of synchrony or amplitude of action potentials) constrain the evolution of signal structure.

Environmental, seasonal, and social modulations of basal activity in a weakly electric fish

The electric organ discharge (EOD) of weakly electric fish encodes information about species, sex, behavioral, and physiological states throughout the lifetime. Its central command is crucial for sensory-motor coordination, and is therefore the target of plastic mechanisms that adapt fish to environmental and social challenges. The EOD waveform of Brachyhypopomus pinnicaudatus is modulated by environmental factors and the neuroendocrine system. In this study we investigate the effects of water temperature and day-night cycle upon EOD rate in this species during the breeding and non-breeding seasons. During the non-breeding season, EOD rate is a linear function of water temperature and exhibits counterclockwise hysteresis. During breeding, a thermal resistance strategy prevents the decrease of EOD rate to cooling. A nocturnal increase of EOD basal rate independent of water temperature and locomotor activity was demonstrated in isolated non-breeding adults and in male-female dyads all year round. An additional increase of nocturnal EOD rate, probably acting as a social courtship signal, was found in breeding dyads. This additional increase of nocturnal EOD rate could not be fully explained by gonadal maturation and was modulated by social stimuli. This study provides novel data on the complex interactions between environment, reproductive cycle, social behavior, and electromotor output in an advantageous model of the vertebrate central nervous system.

Structure and sexual dimorphism of the electrocommunication signals of the weakly electric fish, Adontosternarchus devenanzii

Electrocommunication signals of electric fish vary across species, sexes and individuals. The diversity of these signals and the relative simplicity of the neural circuits controlling them make them a model well-suited for studying the mechanisms, evolution and sexual differentiation of behavior. In most wave-type gymnotiform knifefishes, electric organ discharge (EOD)frequency and EOD modulations known as chirps are sexually dimorphic. In the most speciose gymnotiform family, the Apteronotidae, EOD frequency is higher in males than females in some species, but lower in males than females in others. Sex differences in EOD frequency and chirping, however, have been examined in only three apteronotid species in a single genus, Apteronotus. To understand the diversity of electrocommunication signals, we characterized these behaviors in another genus, Adontosternarchus. Electrocommunication signals of Adontosternarchus devenanzii differed from those of Apteronotus in several ways. Unlike in Apteronotus, EOD frequency was not sexually dimorphic in A. devenanzii. Furthermore,although A. devenanzii chirped in response to playbacks simulating conspecific EODs, the number of chirps did not vary with different stimulus frequencies. A. devenanzii chirps also differed in structure from Apteronotus chirps. Whereas Apteronotus species produce functionally distinct chirp types differing in frequency modulation (FM), A. devenanzii produced only high-frequency chirps that had either single or multiple frequency peaks. Males produced more multi-peaked chirps than females. Thus, the temporal structure of chirps, rather than the amount of FM, delineated chirp types in A. devenanzii. Our results demonstrate that the structure, function and sexual dimorphism of electrocommunication signals are evolutionary labile in apteronotids and may be useful for understanding the diversity of sexually dimorphic behavior.

A "sample-and-hold" pulse-counting integrator as a mechanism for graded memory underlying sensorimotor adaptation

The mechanisms behind the induction of cellular correlates of memory by sensory input and their contribution to meaningful behavioral changes are largely unknown. We previously reported a graded memory in the form of sensorimotor adaptation in the electromotor output of electric fish. Here we show that the mechanism for this adaptation is a synaptically induced long-lasting shift in intrinsic neuronal excitability. This mechanism rapidly integrates hundreds of spikes in a second, or gradually integrates the same number of spikes delivered over tens of minutes. Thus, this mechanism appears immune to frequency-dependent fluctuations in input and operates as a simple pulse counter over a wide range of time scales, enabling it to transduce graded sensory information into a graded memory and a corresponding change in the behavioral output. This adaptation is based on an NMDA receptor-mediated change in intrinsic excitability of the postsynaptic neurons involving the Ca2+-dependent activation of TRP channels.

Regulation and modulation of electric waveforms in gymnotiform electric fish

Weakly electric gymnotiform fish specialize in the regulation and modulation of the action potentials that make up their multi-purpose electric signals. To produce communication signals, gymnotiform fish modulate the waveforms of their electric organ discharges (EODs) over timescales spanning ten orders of magnitude within the animal's life cycle: developmental, reproductive, circadian, and behavioral. Rapid changes lasting milliseconds to seconds are the result of direct neural control of action potential firing in the electric organ. Intermediate-term changes taking minutes to hours result from the action of melanocortin peptides, the pituitary hormones that induce skin darkening and cortisol release in many vertebrates. Long-term changes in the EOD waveform taking days to weeks result from the action of sex steroids on the electrocytes in the electric organ as well as changes in the neural control structures in the brain. These long-term changes in the electric organ seem to be associated with changes in the expression of voltage-gated ion channels in two gene families. Electric organs express multiple voltage-gated sodium channel genes, at least one of which seems to be regulated by androgens. Electric organs also express multiple subunits of the shaker (Kv1) family of voltage-gated potassium channels. Expression of the Kv1 subtype has been found to vary with the duration of the waveform in the electric signal. Our increasing understanding of the mechanisms underlying precise control of electric communication signals may yield significant insights into the diversity of natural mechanisms available for modifying the performance of ion channels in excitable membranes. These mechanisms may lead to better understanding of normal function in a wide range of physiological systems and future application in treatment of disease states involving pathology of excitable membranes.

Distribution and function of potassium channels in the electrosensory lateral line lobe of weakly electric apteronotid fish

Potassium channels are one of the fundamental requirements for the generation of action potentials in the nervous system, and their characteristics shape the output of neurons in response to synaptic input. We review here the distribution and function of a high-threshold potassium channel (Kv3.3) in the electrosensory lateral line lobe of the weakly electric fish Apteronotus leptorhynchus, with particular focus on the pyramidal cells in this brain structure. These cells contain both high-threshold Kv3.3 channels, as well as low-threshold potassium channels of unknown molecular identity. Kv3.3 potassium channels regulate burst discharge in pyramidal cells and enable sustained high frequency firing through their ability to reduce an accumulation of low-threshold potassium current.

Distribution of Kv1-like potassium channels in the electromotor and electrosensory systems of the weakly electric fishApteronotus leptorhynchus

The electromotor and electrosensory systems of the weakly electric fish Apteronotus leptorhynchus are model systems for studying mechanisms of high-frequency motor pattern generation and sensory processing. Voltage-dependent ionic currents, including low-threshold potassium currents, influence excitability of neurons in these circuits and thereby regulate motor output and sensory filtering. Although Kv1-like potassium channels are likely to carry low-threshold potassium currents in electromotor and electrosensory neurons, the distribution of Kv1 alpha subunits in A. leptorhynchus is unknown. In this study, we used immunohistochemistry with six different antibodies raised against specific mammalian Kv1 alpha subunits (Kv1.1-Kv1.6) to characterize the distribution of Kv1-like channels in electromotor and electrosensory structures. Each Kv1 antibody labeled a distinct subset of neurons, fibers, and/or dendrites in electromotor and electrosensory nuclei. Kv1-like immunoreactivity in the electrosensory lateral line lobe (ELL) and pacemaker nucleus are particularly relevant in light of previous studies suggesting that potassium currents carried by Kv1 channels regulate neuronal excitability in these regions. Immunoreactivity of pyramidal cells in the ELL with several Kv1 antibodies is consistent with Kv1 channels carrying low-threshold outward currents that regulate spike waveform in these cells (Fernandez et al., J Neurosci 2005;25:363-371). Similarly, Kv1-like immunoreactivity in the pacemaker nucleus is consistent with a role of Kv1 channels in spontaneous high-frequency firing in pacemaker neurons. Robust Kv1-like immunoreactivity in several other structures, including the dorsal torus semicircularis, tuberous electroreceptors, and the electric organ, indicates that Kv1 channels are broadly expressed and are likely to contribute significantly to generating the electric organ discharge and processing electrosensory inputs.

Pharmacological characterization of ionic currents that regulate high-frequency spontaneous activity of electromotor neurons in the weakly electric fish,Apteronotus leptorhynchus

The neural circuit that controls the electric organ discharge (EOD) of the brown ghost knifefish (Apteronotus leptorhynchus) contains two spontaneous oscillators. Both pacemaker neurons in the medulla and electromotor neurons (EMNs) in the spinal cord fire spontaneously at frequencies of 500-1,000 Hz to control the EOD. These neurons continue to fire in vitro at frequencies that are highly correlated with in vivo EOD frequency. Previous studies used channel blocking drugs to pharmacologically characterize ionic currents that control high-frequency firing in pacemaker neurons. The goal of the present study was to use similar techniques to investigate ionic currents in EMNs, the other type of spontaneously active neuron in the electromotor circuit. As in pacemaker neurons, high-frequency firing of EMNs was regulated primarily by tetrodotoxin-sensitive sodium currents and by potassium currents that were sensitive to 4-aminopyridine and kappaA-conotoxin SIVA, but resistant to tetraethylammonium. EMNs, however, differed from pacemaker neurons in their sensitivity to some channel blocking drugs. Alpha-dendrotoxin, which blocks a subset of Kv1 potassium channels, increased firing rates in EMNs, but not pacemaker neurons; and the sodium channel blocker muO-conotoxin MrVIA, which reduced firing rates of pacemaker neurons, had no effect on EMNs. These results suggest that similar, but not identical, ionic currents regulate high-frequency firing in EMNs and pacemaker neurons. The differences in the ionic currents expressed in pacemaker neurons and EMNs might be related to differences in the morphology, connectivity, or function of these two cell types.

The long-term resetting of a brainstem pacemaker nucleus by synaptic input: a model for sensorimotor adaptation

The cellular mechanisms behind sensorimotor adaptations, such as the adaptation to a sustained change in visual inputs by prism goggles in humans, are not known. Here we present a novel example of long-term sensorimotor adaptation in a well known neuroethological model, the jamming-avoidance response of a weakly electric fish. The adaptation is relatively long lasting, up to 9 hr in vivo, and is likely to be mediated by NMDA receptors. We demonstrate in a brain slice preparation that the pacemaker nucleus is the locus of adaptation and that it responds to long-lasting synaptic stimulation with an increase in the postsynaptic spike frequency persisting for hours after stimulus termination. The mechanism for the neuronal memory behaves as an integrator, and memory duration and strength are quantitatively related to the estimated amount of synaptic stimulation. This finding is contrary to the idea that neurons respond solely to long-lasting synaptic input by turning down their intrinsic excitability. We show that this positive feedback at the cellular level actually contributes to a negative feedback loop at the organismic level if the entire neural circuit and the behavioral link are considered.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than [D[A(BC)]].