The FGF14(F145S) mutation disrupts the interaction of FGF14 with voltage-gated Na+ channels and impairs neuronal excitability

The physiology of the axon initial segment

The action potential generally begins in the axon initial segment (AIS), a principle confirmed by 60 years of research; however, the most recent advances have shown that a very rich biology underlies this simple observation. The AIS has a remarkably complex molecular composition, with a wide variety of ion channels and attendant mechanisms for channel localization, and may feature membrane domains each with distinct roles in excitation. Its function may be regulated in the short term through the action of neurotransmitters, in the long term through activity- and Ca(2+)-dependent processes. Thus, the AIS is not merely the beginning of the axon, but rather a key site in the control of neuronal excitability.

Bioluminescence methodology for the detection of protein-protein interactions within the voltage-gated sodium channel macromolecular complex

Protein-protein interactions are critical molecular determinants of ion channel function and emerging targets for pharmacological interventions. Yet, current methodologies for the rapid detection of ion channel macromolecular complexes are still lacking. In this study we have adapted a split-luciferase complementation assay (LCA) for detecting the assembly of the voltage-gated Na+ (Nav) channel C-tail and the intracellular fibroblast growth factor 14 (FGF14), a functionally relevant component of the Nav channelosome that controls gating and targeting of Nav channels through direct interaction with the channel C-tail. In the LCA, two complementary N-terminus and C-terminus fragments of the firefly luciferase were fused, respectively, to a chimera of the CD4 transmembrane segment and the C-tail of Nav1.6 channel (CD4-Nav1.6-NLuc) or FGF14 (CLuc-FGF14). Co-expression of CLuc-FGF14 and CD4-Nav1.6-NLuc in live cells led to a robust assembly of the FGF14:Nav1.6 C-tail complex, which was attenuated by introducing single-point mutations at the predicted FGF14:Nav channel interface. To evaluate the dynamic regulation of the FGF14:Nav1.6 C-tail complex by signaling pathways, we investigated the effect of kinase inhibitors on the complex formation. Through a platform of counter screenings, we show that the p38/MAPK inhibitor, PD169316, and the IκB kinase inhibitor, BAY 11-7082, reduce the FGF14:Nav1.6 C-tail complementation, highlighting a potential role of the p38MAPK and the IκB/NFκB pathways in controlling neuronal excitability through protein-protein interactions. We envision the methodology presented here as a new valuable tool to allow functional evaluations of protein-channel complexes toward probe development and drug discovery targeting ion channels implicated in human disorders.

Fibroblast growth factor homologous factor 1 interacts with NEMO to regulate NF-κB signaling in neurons

Neuronal survival and plasticity critically depend on constitutive activity of the transcription factor nuclear factor-κB (NF-κB). We here describe a role for a small intracellular fibroblast growth factor homologue, the fibroblast growth factor homologous factor 1 (FHF1/FGF12) in the regulation of NF-κB activity in mature neurons. FHF's have previously been described to control neuronal excitability, and mutations in FHF isoforms give rise to a form of progressive spinocerebellar ataxia. Using a protein-array approach, we identified FHF1b as a novel interactor of the canonical NF–κB modulator IKKγ/NEMO. Co-immunoprecipitation, pull-down and GAL4-reporter experiments, as well as proximity ligation assays confirmed the interaction of FHF1 and NEMO, and demonstrated that a major site of interaction occurred within the axon initial segment. Fhf1 gene silencing strongly activated neuronal NF-κB activity and increased neurite lengths, branching patterns and spine counts in mature cortical neurons. The effects of FHF1 on neuronal NF-κB activity and morphology required the presence of NEMO. Our results imply that FHF1 negatively regulates the constitutive NF-κB activity in neurons.

From cradle to grave: the multiple roles of fibroblast growth factors in neural development

The generation of a functional nervous system involves a multitude of steps that are controlled by just a few families of extracellular signaling molecules. Among these, the fibroblast growth factor (FGF) family is particularly prominent for the remarkable diversity of its functions. FGFs are best known for their roles in the early steps of patterning of the neural primordium and proliferation of neural progenitors. However, other equally important functions have emerged more recently, including in the later steps of neuronal migration, axon navigation, and synaptogenesis. We review here these diverse functions and discuss the mechanisms that account for this unusual range of activities. FGFs are essential components of most protocols devised to generate therapeutically important neuronal populations in vitro or to stimulate neuronal repair in vivo. How FGFs promote the development of the nervous system and maintain its integrity will thus remain an important focus of research in the future.

Dysregulated phosphorylation of Ca(2+) /calmodulin-dependent protein kinase II-α in the hippocampus of subjects with mild cognitive impairment and Alzheimer's disease

Alzheimer's disease (AD) is a progressive, neurodegenerative disorder and the most prevalent senile dementia. The early symptom of memory dysfunction involves synaptic loss, thought to be mediated by soluble amyloid-beta (Aβ) oligomers. These aggregate species target excitatory synapses and their levels correlate with disease severity. Studies in cell culture and rodents have shown that oligomers increase intracellular calcium (Ca(2+)), impairing synaptic plasticity. Yet, the molecular mechanism mediating Aβ oligomers' toxicity in the aged brain remains unclear. Here, we apply quantitative immunofluorescence in human brain tissue from clinically diagnosed mild cognitive impaired (MCI) and AD patients to investigate the distribution of phosphorylated (active) Ca(2+) /calmodulin-dependent protein kinase-α (p(Thr286)CaMKII), a critical enzyme for activity-dependent synaptic remodeling associated with cognitive function. We show that p(Thr286)CaMKII immunoreactivity is redistributed from dendritic arborizations to neural perikarya of both MCI and AD hippocampi. This finding correlates with cognitive assessment scores, suggesting that it may be a molecular read-out of the functional deficits in early AD. Treatment with oligomeric Aβ replicated the observed phenotype in mice and resulted in a loss of p(Thr286)CaMKII from synaptic spines of primary hippocampal neurons. Both outcomes were prevented by inhibiting the phosphatase calcineurin (CaN). Collectively, our results support a model in which the synaptotoxicity of Aβ oligomers in human brain involves the CaN-dependent subcellular redistribution of p(Thr286)CaMKII. Therapies designed to normalize the homeostatic imbalance of neuronal phosphatases and downstream dephosphorylation of synaptic p(Thr286)CaMKII should be considered to prevent and treat early AD.

Fibroblast growth factor homologous factor 13 regulates Na+ channels and conduction velocity in murine hearts

RATIONALE

Fibroblast growth factor homologous factors (FHFs), a subfamily of fibroblast growth factors (FGFs) that are incapable of functioning as growth factors, are intracellular modulators of Na(+) channels and have been linked to neurodegenerative diseases. Although certain FHFs have been found in embryonic heart, they have not been reported in adult heart, and they have not been shown to regulate endogenous cardiac Na(+) channels or to participate in cardiac pathophysiology.

OBJECTIVE

We tested whether FHFs regulate Na(+) channels in murine heart.

METHODS AND RESULTS

We demonstrated that isoforms of FGF13 are the predominant FHFs in adult mouse ventricular myocytes. FGF13 binds directly to, and colocalizes with, the Na(V)1.5 Na(+) channel in the sarcolemma of adult mouse ventricular myocytes. Knockdown of FGF13 in adult mouse ventricular myocytes revealed a loss of function of Na(V)1.5-reduced Na(+) current density, decreased Na(+) channel availability, and slowed Na(V)1.5-reduced Na(+) current recovery from inactivation. Cell surface biotinylation experiments showed ≈45% reduction in Na(V)1.5 protein at the sarcolemma after FGF13 knockdown, whereas no changes in whole-cell Na(V)1.5 protein or in mRNA level were observed. Optical imaging in neonatal rat ventricular myocyte monolayers demonstrated slowed conduction velocity and a reduced maximum capture rate after FGF13 knockdown.

CONCLUSION

These findings show that FHFs are potent regulators of Na(+) channels in adult ventricular myocytes and suggest that loss-of-function mutations in FHFs may underlie a similar set of cardiac arrhythmias and cardiomyopathies that result from Na(V)1.5 loss-of-function mutations.

Identification of novel interaction sites that determine specificity between fibroblast growth factor homologous factors and voltage-gated sodium channels

Fibroblast growth factor homologous factors (FHFs, FGF11-14) bind to the C termini (CTs) of specific voltage-gated sodium channels (VGSC) and thereby regulate their function. The effect of an individual FHF on a specific VGSC varies greatly depending upon the individual FHF isoform. How individual FHFs impart distinctive effects on specific VGSCs is not known and the specificity of these pairwise interactions is not understood. Using several biochemical approaches combined with functional analysis, we mapped the interaction site for FGF12B on the Na(V)1.5 C terminus and discovered previously unknown determinants necessary for FGF12 interaction. Also, we demonstrated that FGF12B binds to some, but not all Na(V)1 CTs, suggesting specificity of interaction. Exploiting a human single nucleotide polymorphism in the core domain of FGF12 (P149Q), we identified a surface proline that contributes a part of this pairwise specificity. This proline is conserved among all FHFs, and mutation of the homologous residue in FGF13 also leads to loss of interaction with a specific VGSC CT (Na(V)1.1) and loss of modulation of the resultant Na(+) channel function. We hypothesized that some of the specificity mediated by this proline may result from differences in the affinity of the binding partners. Consistent with this hypothesis, surface plasmon resonance data showed that the P149Q mutation decreased the binding affinity between FHFs and VGSC CTs. Moreover, immunocytochemistry revealed that the mutation prevented proper subcellular targeting of FGF12 to the axon initial segment in neurons. Together, these results give new insights into details of the interactions between FHFs and Na(V)1.x CTs, and the consequent regulation of Na(+) channels.

Autosomal dominant cerebellar ataxias

Cerebellar ataxias with autosomal dominant transmission (ADCA) are far rarer than sporadic cases of cerebellar ataxia. The identification of genes involved in dominant forms has confirmed the genetic heterogeneity of these conditions and of the underlying mechanisms and pathways. To date, at least 28 genetic loci and, among them, 20 genes have been identified. In many instances, the phenotype is not restricted to cerebellar dysfunction but includes more complex multisystemic neurological deficits. Seven ADCA (SCA1, 2, 3, 6, 7, 17, and dentatorubro-pallido-luysian atrophy) are caused by repeat expansions in the corresponding proteins; phenotype-genotype correlations have shown that repeat size influences the progression of the disease, its severity and clinical differences among patients, including the phenomenon of anticipation between generations. All other ADCA are caused either by non-coding repeat expansions, conventional mutations or large rearrangements in genes with different functions. This review will focus on the genetic features of ADCA and on the clinical differences among the different forms.

Morphine modulation of thrombospondin levels in astrocytes and its implications for neurite outgrowth and synapse formation

Opioid receptor signaling via EGF receptor (EGFR) transactivation and ERK/MAPK phosphorylation initiates diverse cellular responses that are cell type-dependent. In astrocytes, multiple μ opioid receptor-mediated mechanisms of ERK activation exist that are temporally distinctive and feature different outcomes. Upon discovering that chronic opiate treatment of rats down-regulates thrombospondin 1 (TSP1) expression in the nucleus accumbens and cortex, we investigated the mechanism of action of this modulation in astrocytes. TSP1 is synthesized in astrocytes and is released into the extracellular matrix where it is known to play a role in synapse formation and neurite outgrowth. Acute morphine (hours) reduced TSP1 levels in astrocytes. Chronic (days) opioids repressed TSP1 gene expression and reduced its protein levels by μ opioid receptor and ERK-dependent mechanisms in astrocytes. Morphine also depleted TSP1 levels stimulated by TGFβ1 and abolished ERK activation induced by this factor. Chronic morphine treatment of astrocyte-neuron co-cultures reduced neurite outgrowth and synapse formation. Therefore, inhibitory actions of morphine were detected after both acute and chronic treatments. An acute mechanism of morphine signaling to ERK that entails depletion of TSP1 levels was suggested by inhibition of morphine activation of ERK by a function-blocking TSP1 antibody. This raises the novel possibility that acute morphine uses TSP1 as a source of EGF-like ligands to activate EGFR. Chronic morphine inhibition of TSP1 is reminiscent of the negative effect of μ opioids on EGFR-induced astrocyte proliferation via a phospho-ERK feedback inhibition mechanism. Both of these variations of classical EGFR transactivation may enable opiates to diminish neurite outgrowth and synapse formation.

Channelrhodopsin-2 localised to the axon initial segment

The light-gated cation channel Channelrhodopsin-2 (ChR2) is a powerful and versatile tool for controlling neuronal activity. Currently available versions of ChR2 either distribute uniformly throughout the plasma membrane or are localised specifically to somatodendritic or synaptic domains. Localising ChR2 instead to the axon initial segment (AIS) could prove an extremely useful addition to the optogenetic repertoire, targeting the channel directly to the site of action potential initiation, and limiting depolarisation and associated calcium entry elsewhere in the neuron. Here, we describe a ChR2 construct that we localised specifically to the AIS by adding the ankyrinG-binding loop of voltage-gated sodium channels (Na_vII-III) to its intracellular terminus. Expression of ChR2-YFP-Na_vII-III did not significantly affect the passive or active electrical properties of cultured rat hippocampal neurons. However, the tiny ChR2 currents and small membrane depolarisations resulting from AIS targeting meant that optogenetic control of action potential firing with ChR2-YFP-Na_vII-III was unsuccessful in baseline conditions. We did succeed in stimulating action potentials with light in some ChR2-YFP-Na_vII-III-expressing neurons, but only when blocking KCNQ voltage-gated potassium channels. We discuss possible alternative approaches to obtaining precise control of neuronal spiking with AIS-targeted optogenetic constructs and propose potential uses for our ChR2-YFP-Na_vII-III probe where subthreshold modulation of action potential initiation is desirable.

Building and maintaining the axon initial segment

The axon initial segment is a unique neuronal subregion involved in the initiation of action potentials and in the control of axonal identity. Recent work has helped our understanding of how this specialised structure develops, not least in identifying possible mechanisms leading to the localisation of the AIS's master organiser protein, ankyrin-G. The most exciting current work, however, focuses on later aspects of AIS function and plasticity. Recent studies have shown that the AIS is subdivided into distinct structural and functional domains, have demonstrated how the AIS acts as a cytoplasmic barrier for axonal transport, and have discovered that the AIS can be surprisingly plastic in its responses to alterations in neuronal activity.

Genome-wide association identifies candidate genes that influence the human electroencephalogram

Complex psychiatric disorders are resistant to whole-genome analysis due to genetic and etiological heterogeneity. Variation in resting electroencephalogram (EEG) is associated with common, complex psychiatric diseases including alcoholism, schizophrenia, and anxiety disorders, although not diagnostic for any of them. EEG traits for an individual are stable, variable between individuals, and moderately to highly heritable. Such intermediate phenotypes appear to be closer to underlying molecular processes than are clinical symptoms, and represent an alternative approach for the identification of genetic variation that underlies complex psychiatric disorders. We performed a whole-genome association study on alpha (alpha), beta (beta), and theta (theta) EEG power in a Native American cohort of 322 individuals to take advantage of the genetic and environmental homogeneity of this population isolate. We identified three genes (SGIP1, ST6GALNAC3, and UGDH) with nominal association to variability of theta or alpha power. SGIP1 was estimated to account for 8.8% of variance in power, and this association was replicated in US Caucasians, where it accounted for 3.5% of the variance. Bayesian analysis of prior probability of association based upon earlier linkage to chromosome 1 and enrichment for vesicle-related transport proteins indicates that the association of SGIP1 with theta power is genuine. We also found association of SGIP1 with alcoholism, an effect that may be mediated via the same brain mechanisms accessed by theta EEG, and which also provides validation of the use of EEG as an endophenotype for alcoholism.

Molecular pathology of the fibroblast growth factor family

The human fibroblast growth factor (FGF) family contains 22 proteins that regulate a plethora of physiological processes in both developing and adult organism. The mutations in the FGF genes were not known to play role in human disease until the year 2000, when mutations in FGF23 were found to cause hypophosphatemic rickets. Nine years later, seven FGFs have been associated with human disorders. These include FGF3 in Michel aplasia; FGF8 in cleft lip/palate and in hypogonadotropic hypogonadism; FGF9 in carcinoma; FGF10 in the lacrimal/salivary glands aplasia, and lacrimo-auriculo-dento-digital syndrome; FGF14 in spinocerebellar ataxia; FGF20 in Parkinson disease; and FGF23 in tumoral calcinosis and hypophosphatemic rickets. The heterogeneity in the functional consequences of FGF mutations, the modes of inheritance, pattern of involved tissues/organs, and effects in different developmental stages provide fascinating insights into the physiology of the FGF signaling system. We review the current knowledge about the molecular pathology of the FGF family.

Physiologic alterations in ataxia: channeling changes into novel therapies

The ataxias constitute a heterogeneous group of diseases in which cerebellar dysfunction typically underlies the major neurologic manifestations. It is increasingly clear that ataxia can result directly from mutations in ion channels or from perturbations in ion channel physiology in the absence of a primary channel defect. Neuronal dysfunction stemming from perturbed channel activity likely explains some motor deficits in episodic and degenerative ataxias. Understanding these pathophysiologic changes may reveal novel therapeutic targets for symptomatic treatment of ataxia.

SCA27 caused by a chromosome translocation: further delineation of the phenotype

We report of a spinocerebellar ataxia (SCA)27 in a daughter and her mother whose karyotype is 46, XX t(5;13)(q31.2;q33.1). The translocation breakpoint is identical in both patients, disrupting the gene-encoding fibroblast growth factor 14 isoform b (FGF14-1b). Clinically, both show signs of SCA, although the daughter is the most affected with early onset cerebellar ataxia, microcephaly, and severe mental retardation. FGF14-1b is the predominant isoform in brain, where it interacts with the voltage gated Na channel. Fgf14(-/-) mice develop ataxia and paroxysmal dyskinesia and have cognitive deficits. One missense and one non-sense mutation in FGF14 have previously been linked to SCA27. Truncation of one allele in our patients suggests that haploinsuffiency of FGF14 can cause SCA27.

FGF14 N-terminal splice variants differentially modulate Nav1.2 and Nav1.6-encoded sodium channels

The Intracellular Fibroblast Growth Factor (iFGF) subfamily includes four members (FGFs 11-14) of the structurally related FGF superfamily. Previous studies showed that the iFGFs interact directly with the pore-forming (alpha) subunits of voltage-gated sodium (Nav) channels and regulate the functional properties of sodium channel currents. Sequence heterogeneity among the iFGFs is thought to confer specificity to this regulation. Here, we demonstrate that the two N-terminal alternatively spliced FGF14 variants, FGF14-1a and FGF14-1b, differentially regulate currents produced by Nav1.2 and Nav1.6 channels. FGF14-1b, but not FGF14-1a, attenuates both Nav1.2 and Nav1.6 current densities. In contrast, co-expression of an FGF14 mutant, lacking the N-terminus, increased Nav1.6 current densities. In neurons, both FGF14-1a and FGF14-1b localized at the axonal initial segment, and deletion of the N-terminus abolished this localization. Thus, the FGF14 N-terminus is required for targeting and functional regulation of Nav channels, suggesting an important function for FGF14 alternative splicing in regulating neuronal excitability.

Crystal structure of a fibroblast growth factor homologous factor (FHF) defines a conserved surface on FHFs for binding and modulation of voltage-gated sodium channels

Voltage-gated sodium channels (Nav) produce sodium currents that underlie the initiation and propagation of action potentials in nerve and muscle cells. Fibroblast growth factor homologous factors (FHFs) bind to the intracellular C-terminal region of the Nav alpha subunit to modulate fast inactivation of the channel. In this study we solved the crystal structure of a 149-residue-long fragment of human FHF2A which unveils the structural features of the homology core domain of all 10 human FHF isoforms. Through analysis of crystal packing contacts and site-directed mutagenesis experiments we identified a conserved surface on the FHF core domain that mediates channel binding in vitro and in vivo. Mutations at this channel binding surface impaired the ability of FHFs to co-localize with Navs at the axon initial segment of hippocampal neurons. The mutations also disabled FHF modulation of voltage-dependent fast inactivation of sodium channels in neuronal cells. Based on our data, we propose that FHFs constitute auxiliary subunits for Navs.

The functional organization and assembly of the axon initial segment

Action potential initiation, modulation, and duration in neurons depend on a variety of Na+ and K+ channels that are highly enriched at the axon initial segment (AIS). The AIS also has high densities of cell adhesion molecules (CAMs), modulatory proteins, and a unique extracellular matrix (ECM). In contrast to other functional domains of axons (e.g. the nodes of Ranvier and axon terminals) whose development depends on the interactions with different cells (e.g. myelinating glia and postsynaptic cells), the recruitment and retention of AIS proteins is intrinsically specified through axonal cytoskeletal and scaffolding proteins. We speculate that the AIS has previously unappreciated forms of plasticity that influence neuronal excitability, and that AIS plasticity is regulated by the developmental or activity-dependent modulation of scaffolding protein levels rather than directly altering ion channel expression.

The BTB/kelch protein, KRIP6, modulates the interaction of PICK1 with GluR6 kainate receptors

Neuronal proteins of the BTB/kelch and PDZ domain families interact with different regions of the cytoplasmic C-terminal domain of the GluR6 kainate receptor subunit. The BTB/kelch protein KRIP6 binds within a 58 amino acid segment of GluR6 proximal to the plasma membrane. In contrast, PDZ domain proteins, such as PICK1 and PSD95, interact with the last 4 residues of the GluR6 C-terminus. KRIP6 reduces peak currents mediated by recombinant GluR6 receptors and by native kainate receptors in neurons, whereas PICK1 stabilizes kainate receptors at synapses. Thus, protein-protein interactions at the C-terminal domain of GluR6 are important for regulating kainate receptor physiology. Here, we show by co-clustering and co-immunoprecipitation that KRIP6 interacts with PICK1 in heterologous cells. In addition, we demonstrate a novel modulation of GluR6 receptors by PICK1 resulting in increased peak current and relative desensitization of GluR6-mediated currents, phenotypes opposite to those produced by KRIP6. Importantly, these effects cancel out when KRIP6 and PICK1 are co-expressed together with GluR6. KRIP6 and PICK1 strongly co-cluster and co-immunoprecipitate regardless of the presence of GluR6. Immunofluorescence analysis reveals that GluR6 can either join the KRIP6-PICK1 clusters or remain separate; however, co-expression of KRIP6 reduces the fraction of PICK1 that co-immunoprecipitates with GluR6. Taken together, these results indicate that, in addition to a previously demonstrated direct interaction with the GluR6 C-terminal domain, KRIP6 regulates kainate receptors by inhibiting PICK1 modulation via competition or a mutual blocking effect.

Voltage-gated Na+ channels: potential for beta subunits as therapeutic targets

BACKGROUND

Voltage gated Na(+) channels (VGSCs) contain a pore-forming alpha subunit and one or more beta subunits. VGSCs are involved in a wide variety of pathophysiologies, including epilepsy, cardiac arrhythmia, multiple sclerosis, periodic paralysis, migraine, neuropathic and inflammatory pain, Huntington's disease and cancer. Increasing evidence implicates the beta subunits as key players in these disorders.

OBJECTIVE

To review the recent literature describing the multifunctional roles of VGSC beta subunits in the context of their role(s) in disease.

METHODS

An extensive review of the literature on beta subunits.

RESULTS/CONCLUSION

beta subunits are multifunctional. As components of VGSC complexes, beta subunits mediate signaling processes regulating electrical excitability, adhesion, migration, pathfinding and transcription. beta subunits may prove useful in disease diagnosis and therapy.

FGF14 regulates the intrinsic excitability of cerebellar Purkinje neurons

A missense mutation in the fibroblast growth factor 14 (FGF14) gene underlies SCA27, an autosomal dominant spinocerebellar ataxia in humans. Mice with a targeted disruption of the Fgf14 locus (Fgf14(-/-)) develop ataxia resembling human SCA27. We tested the hypothesis that loss of FGF14 affects the firing properties of Purkinje neurons, which play an important role in motor control and coordination. Current clamp recordings from Purkinje neurons in cerebellar slices revealed attenuated spontaneous firing in Fgf14(-/-) neurons. Unlike in the wild type animals, more than 80% of Fgf14(-/-) Purkinje neurons were quiescent and failed to fire repetitively in response to depolarizing current injections. Immunohistochemical examination revealed reduced expression of Nav1.6 protein in Fgf14(-/-) Purkinje neurons. Together, these observations suggest that FGF14 is required for normal Nav1.6 expression in Purkinje neurons, and that the loss of FGF14 impairs spontaneous and repetitive firing in Purkinje neurons by altering the expression of Nav1.6 channels.

Functional evolutionary history of the mouse Fgf gene family

Fibroblast Growth Factors (FGFs) are polypeptides with diverse activities in development and physiology. The mammalian Fgf family can be divided into the intracellular Fgf11/12/13/14 subfamily (iFGFs), the hormone-like Fgf15/21/23 subfamily (hFGFs), and the canonical Fgf subfamilies, including Fgf1/2/5, Fgf3/4/6, Fgf7/10/22, Fgf8/17/18, and Fgf9/16/20. However, all Fgfs are evolutionarily related. We propose that an Fgf13-like gene is the ancestor of the iFgf subfamily and the most likely evolutionary ancestor of the entire Fgf family. Potential ancestors of the canonical and hFgf subfamilies, Fgf4-, Fgf5-, Fgf8-, Fgf9-, Fgf10-, and Fgf15-like, appear to have derived from an Fgf13-like ancestral gene. Canonical FGFs function in a paracrine manner, while hFGFs function in an endocrine manner. We conclude that the ancestral Fgfs for these subfamilies acquired this functional diversity before the evolution of vertebrates. During the evolution of early vertebrates, the Fgf subfamilies further expanded to contain three or four members in each subfamily.