Exploring the molecular basis of neuronal excitability in a vocal learner
BMC Genomics volume 20, Article number: 629 (2019)
Vocal learning, the ability to learn to produce vocalizations through imitation, relies on specialized brain circuitry known in songbirds as the song system. While the connectivity and various physiological properties of this system have been characterized, the molecular genetic basis of neuronal excitability in song nuclei remains understudied. We have focused our efforts on examining voltage-gated ion channels to gain insight into electrophysiological and functional features of vocal nuclei. A previous investigation of potassium channel genes in zebra finches (Taeniopygia guttata) revealed evolutionary modifications unique to songbirds, as well as transcriptional specializations in the song system [Lovell PV, Carleton JB, Mello CV. BMC Genomics 14:470 2013]. Here, we expand this approach to sodium, calcium, and chloride channels along with their modulatory subunits using comparative genomics and gene expression analysis encompassing microarrays and in situ hybridization.
We found 23 sodium, 38 calcium, and 33 chloride channel genes (HGNC-based classification) in the zebra finch genome, several of which were previously unannotated. We determined 15 genes are missing relative to mammals, including several genes (CLCAs, BEST2) linked to olfactory transduction. The majority of sodium and calcium but few chloride channels showed differential expression in the song system, among them SCN8A and CACNA1E in the direct motor pathway, and CACNG4 and RYR2 in the anterior forebrain pathway. In several cases, we noted a seemingly coordinated pattern across multiple nuclei (SCN1B, SCN3B, SCN4B, CACNB4) or sparse expression (SCN1A, CACNG5, CACNA1B).
The gene families examined are highly conserved between avian and mammalian lineages. Several cases of differential expression likely support high-frequency and burst firing in specific song nuclei, whereas cases of sparse patterns of expression may contribute to the unique electrophysiological signatures of distinct cell populations. These observations lay the groundwork for manipulations to determine how ion channels contribute to the neuronal excitability properties of vocal learning systems.
Motor learning, the process by which motor skills are acquired and perfected through practice, requires fine-tuning of sensorimotor circuit elements to ultimately produce precise motor output. One remarkable example is vocal learning, a trait that allows individuals to learn their vocalizations through auditorily guided vocal practice. Vocal learning is demonstrated only by some mammals (humans, cetaceans, bats, and possibly pinnipeds and elephants), and three groups of birds (parrots, hummingbirds, and songbirds) . Songbirds provide a particularly powerful model for studying the role of specific genes within vocal learning circuits. Decades of study in the zebra finch have revealed much about neuroanatomical substrates for vocal learning, including the connectivity and electrophysiological properties of a discrete set of vocal nuclei called the song system [2, 3]. Single-unit neural recordings and modeling have described the electrophysiological characteristics and underlying conductances of several song system neuron types [4,5,6,7,8,9,10,11,12,13]. Less established is how the regulated expression of genes - specifically ion channel genes - gives rise to the physiological and excitable properties of song system neurons.
The zebra finch song system is composed of the posterior direct motor pathway (DMP), necessary for song production [14, 15], and the anterior forebrain pathway (AFP) necessary for song learning and adult song variability [16,17,18,19]. Nucleus HVC (proper name) projects to both the DMP and AFP (Fig. 1), and receives inputs from nucleus interfacialis of the nidopallium (NIf)  and thalamic nucleus uvaeformis (Uva) . In the DMP, terminals from HVC synapse onto neurons of the robust nucleus of the arcopallium (RA), which are considered analogous to layer 5/6 motor neurons of mammalian laryngeal motor cortex . RA then projects to the midbrain’s dorsomedial (DM) nucleus of the intercollicular complex, to the tracheosyringeal subdivision of the XII cranial nerve nucleus (nXIIts), which innervates the syrinx (the avian vocal organ), and to medullary respiratory centers . In the AFP, striatal Area X receives projections from HVC and projects to the medial nucleus of dorsolateral thalamus (DLM), which projects to the lateral magnocellular nucleus of the anterior nidopallium (LMAN) . LMAN projects to Area X  but also to RA, converging with inputs from the DMP . Detailed knowledge of this circuitry facilitates the identification of the molecular genetic determinants of its electrophysiological properties.
Given their direct effects on membrane potential, cell excitability, and neuronal firing, ion channels genes are among the most important determinants of neuronal function [27,28,29]. Their regulation and modulation is important in establishing electrophysiological properties of specific circuits, and may also contribute to nonsynaptic mechanisms of learning and memory [30,31,32]. As such, much can be learned about circuit function and excitability properties by identifying regional differences in ion channel gene expression. We previously examined the genomics and brain expression of potassium channel genes in zebra finch  and uncovered evidence of differential expression for several key regulators of resting membrane and action potential repolarization. Here, we extend our analysis to sodium, calcium, and chloride channel families, all containing genes involved in determining neuronal excitability properties. Whereas sodium currents drive the rising phase of the action potential (AP) to ensure its initiation and propagation, calcium channels perform a wide array of functions such as regulating synaptic vesicle release, activating second messenger systems, and integrating signals over dendrites [27, 34, 35]. While less is known about chloride channels, some members of this family play a role in stabilizing resting membrane potential .
We identified 23 sodium, 38 calcium, and 33 chloride channel genes in the zebra finch genome - the majority with clear orthologs in mammals, indicating their broad conservation across vertebrates. We also established complements of ion channel genes that are differentially expressed within the major song system nuclei compared to the surround. These data represent the most comprehensive characterization of these ion channel gene families in a bird species to date. They also provide novel insights into how molecular specializations associated with neuronal excitability may have evolved to support a complex, learned vocal behavior.
Determining the full complement of sodium, calcium, and chloride channel genes in the zebra finch genome
We identified the sets of sodium, calcium, and chloride channel genes in zebra finch using a modified strategy of that previously used for potassium channels . We first compiled the full set of human orthologs for these genes by navigating gene family hierarchies in HGCN (Ion channels > Ion channels by channel type > Sodium channels / Calcium channels / Chloride channels) and retrieving all gene names within all subfamilies for each of these three ion channel families . Starting with human, arguably one of the most complete and best annotated genomes, allowed us to define a comprehensive set of genes to be examined. This effort identified a total of 108 ion channel genes, which included 23 sodium, 43 calcium, and 41 chloride channel and channel-related genes. This comprehensive set included genes encoding ion channels which have been linked to neuronal transmission, as well as the auxiliary subunits that modulate these channels. It also includes epithelial sodium channel (SCNN) genes, sperm-associated cation channel (CATSPER) genes, and volume-regulated anion channel subunits (LRRC8), even though these families have not been implicated in neuronal excitability and their expression in mammals is generally biased toward non-neuronal tissues [38,39,40].
We next searched for the orthologous genes in zebra finch using the following pipeline: (1) Obtain from NCBI (National Center for Biotechnology Information) the human records for the HGNC (HUGO Gene Nomenclature Committee) gene list; (2) Align the human models to the zebra finch genome to identify potential orthologous loci; (3) Verify the synteny of the top-scoring alignment; 4) Demonstrate that the zebra finch gene aligns preferentially to the corresponding human locus; and (5) Demonstrate that secondary alignments represent other family members. These steps were also applied to identify chicken orthologs (“Chicken locus” column of Additional file 1: Table S1) in the galGal5 (GCA_000002315.3 ) or updated galGal6 (GRCg6a; GCA_000002315.5) chicken genome. Each step is detailed next and graphically summarized in Additional file 2.
(1) We queried the NCBI Entrez Gene database with the names of the 108 genes in the starting HGNC gene list, to retrieve the corresponding human gene records. All 108 human genes were found in NCBI. (2) We retrieved the human RefSeq nucleotide sequences of the longest isoform for each annotated gene and used these as queries to search the zebra finch genome (taeGut1/WUGSC 3.2.4; GCA_000151805; Jul 2008 ) using the UCSC Browser BLAT , noting all significant hits, defined as BLAT scores above 50 (top BLAT scores for human-to-taeGut1 alignments are reported in Additional file 1: Table S1). (3) We compared the synteny of the top scoring alignment in finch with that of the human query. We manually examined the immediately flanking genes on both sides of the alignment and query, if necessary expanding our analysis to genes farther up- or downstream. (4) The gene prediction (Ensembl release 95 ) at the top scoring finch locus was aligned back to human (using BLAT and/or BLAST) and the top-scoring alignment in human examined to confirm it was to the expected locus (same as initial query). For zebra finch gene models with multiple transcripts in Ensembl, we selected the longest variant. (5) All significant secondary alignments found in the finch genome were verified as done for the top hits (synteny comparison and back-alignment to human). This step served both to further confirm the orthology of the top hit and to identify possible novel paralogs or paralogs missing from the human gene list (e.g. TPCN3) in the zebra finch genome.
In the majority of cases, the top-scoring alignments in both finch and humans were reciprocally aligning models with conserved synteny across species, and all significant secondary alignments were to other genes of the same family, thus confirming orthology to human of the primary alignment in finch. These zebra finch Ensembl models that passed our synteny and alignment criteria are reported in Table 1 (n = 61 genes, no symbols next to gene name and the “Zebra finch locus” column contains only Ensembl model(s)). We detected no cases of multiple top-scoring alignments with similar or different syntenies, thus we conclude that no novel paralogs and/or segmental duplications (e.g. of a cluster of syntenic genes) are present in the zebra finch genome for the gene families examined. In a few cases no Ensembl model was present but synteny of the top scoring finch locus was conserved with the human ortholog (n = 2; chromosome location in taeGut1 is indicated in the “Zebra finch locus” column of Table 1). Lastly, secondary high scoring BLAT alignments to chrUn (chromosome Unknown) were not examined further, as those cases represent allelic variants in the zebra finch genome (as detailed in ).
Several groups of genes required further verification efforts and/or variations of the main pipeline described above; these cases are discussed next and graphically summarized in Additional file 2.
Case 1: Gene clusters
Some sodium channel (SCNxA) genes are arranged in clusters; in such cases, BLAT alignments tend to span multiple family members, or identify another family member in a given cluster. For these genes, we compared the syntenies within and around the clusters as a group across species, and closely examined the number, orientation and alignments of the individual genes in each cluster before deciding orthology. Both clusters in this family were found conserved between finch and humans (n = 7 genes in 2 clusters: SCN5A, SCN10A, SCN11A, and SCN1A, SCN2A, SCN3A, SCN9A).
Case 2: Gene in taeGut1 but synteny diverges from human
Some genes met the reciprocal top-scoring alignment criteria but their synteny did not match that of human. To examine them further, we compared the synteny in zebra finch and human to that of orthologous loci in other birds, including Bengalese finch, medium ground finch, and tit (other passeriforms), falcon (outgroup to psitacopasseridae), chicken and/or quail (galliformes, basal neoaves), tinamou, ostrich, and emu (ratites, basal avian), lizard, alligator and turtle (non-avian sauropsids, outgroup to birds), Xenopus (amphibian, outgroup to amniotes), and platypus and opossum (basal mammals). We only used non-human orthologs that met two criteria: (1) conserved synteny and (2) aligned preferentially to the correct human locus. For 4 genes, this phylogenetic approach revealed an indirect relationship between the syntenies of the human query and the top scoring alignment in finch (indicated by a “$” in Table 1; details in “Synteny notes” column of Additional file 1: Table S1; examples in Figs. 2 and 3). For another 2 genes (ASIC2 and CLIC2), human and zebra finch syntenies could not be fully linked through phylogenetic relationships, indicating weaker evidence of orthology (indicated by a “#” in Table 1).
Case 3: Gene detection in taeGut1 required non-human ortholog queries
In other cases, there were no significant alignments of the human query to taeGut1. For these genes, we pursued further searches in taeGut1 using as queries orthologs from the same list of other species discussed above in Case 2. For all 6 genes in this category, the top-scoring alignment identified a zebra finch Ensembl model that itself aligned to the correct locus in human, and whose synteny was conserved with human (indicated by a “◆” in Table 1 and a chromosome location in the “Zebra finch locus” column). In the case of TPCN3, which lacks a human model, we used orthologs from other species (chicken, dog) as queries to confirm the locus (initially identified from secondary alignments of human TPCN1 and 2) in taeGut1, and confirmed that synteny in zebra finch was conserved with those species.
Case 4: Synteny verification required examination of PacBio assembly
In yet other cases, a locus with an Ensembl model in taeGut1 was identified through alignments of human or other species’ models and that Ensembl model aligned back to the correct locus in human but lacked syntenic context. This occurred when the top-scoring alignment in finch was located between gaps and/or on an unplaced scaffold or a chromosome “random” (e.g. chr13_random), or when the only significant alignment was to chrUn. For these genes, we next used the zebra finch Ensembl model to conduct searches of a newer zebra finch assembly (Tgut_diploid_1.0; GCA_002008985; Korlach et al. 2017) built exclusively from PacBio reads (Pacific Biosciences). This PacBio assembly has higher contiguity, fewer gaps (individual scaffolds have no gaps), and often more genes and intergenic regions than the Sanger, 454, or Illumina assemblies . Since this PacBio assembly currently does not have gene predictions, the PacBio scaffolds found to contain the top scoring alignments were examined further by BLASTing against RefSeq databases . Briefly, we BLAST searched the NCBI collections of avian and non-avian RefSeqs separately, using as query the identified PacBio scaffolds, broken down into 10–40 kb segments depending on gene size and gene density. The ion channel gene of interest was confirmed as present in zebra finch if the BLAST alignment output had high-scoring alignments of: a) the correct orthologs from multiple non-avian species at the same locus, and b) neighboring genes to that locus that confirmed a conserved synteny with human or other species as detailed in Case 2. All 8 genes in this category are indicated by the presence of both an Ensembl model and a PacBio location in the “Zebra finch locus” column of Table 1). We note that in 2 of these 8 cases (ASIC3 and CATSPERE), the top scoring alignment back to human was not to the correct locus, but both Ensembl models are very partial relative to their full-length human ortholog. In such cases the query seems to be targeting a conserved region or motif rather than the entire gene, thus the reciprocal alignment is not an accurate criterion.
Case 5: Gene located only in PacBio
In several other cases, no locus could be found through alignments of human or other species’ orthologs to taeGut1. For these genes, we searched the PacBio assembly with the human model and/or other species’ orthologs and confirmed synteny using BLAST to explore the surrounding scaffold sequence (see example in Additional file 3), as described in Cases 3 and 4 above. As there are no models in the PacBio assembly, we did not perform alignments of the zebra finch locus to human. In all cases in this category, we found a correct ortholog (n = 10; genes have a PacBio location instead of an Ensembl model in the “Zebra finch locus” column of Table 1). In 5 of these 10 cases (indicated by “▽” in Table 1 and Additional file 1: Table S1), we found only short segments of the ion channel gene between its syntenic genes, indicating severe truncation. Because these PacBio scaffolds have no internal missing sequence due to gaps, such truncations point to likely examples of pseudogenization.
Case 6: Gene could not be found in either zebra finch assembly
When no significant hits could be recovered from alignments of orthologs in either taeGut1 or the zebra finch PacBio assembly, we followed a pipeline for candidate missing genes, specifically: (1) We BLAST searched the finch PacBio assembly using models of syntenic genes from comparative species as queries. When the syntenic genes were found on the same scaffold, we conducted additional BLAST searches of short segments (1–2 kb) of their intergenic region at a time. We interpret cases in which syntenic genes were adjacent but no traces of the ion channel gene could be found intergenically as high likelihood cases of gene loss, noting again that finch PacBio scaffolds are gapless. (2) To better understand the origin of such losses, we systematically searched for the gene of interest and syntenic genes in all birds with an assembled genome deposited in NCBI, initially based on gene name searches in Entrez Gene. In cases where other avian species’ scaffolds contained syntenic genes but lacked a model for the gene of interest, we conducted BLAST searches of the intergenic region as described in Case 4. In cases where we could not find evidence for the gene in any bird species, we applied the same general strategy above to search for the gene of interest and its syntenic context in representative outgroups (e.g. alligator, lizard), to trace the most likely pattern of gene loss in groups ancestral to birds (for example, see Fig. 3). We note, however, that there are often gaps in the intergenic regions of these Illumina-based assemblies, which limits the certainty of conclusions regarding gene loss in species currently lacking a PacBio assembly.
In all cases where the Ensembl model at the zebra finch locus was annotated, that annotation was correct. In cases where an Ensembl model was present but annotated as a “novel gene”, we annotated it accordingly (n = 34 models; indicated by a ‘^’ next to the Ensembl model in Table 1).
Assessing gene model completeness and expanding gene models with additional sequence
While inspecting taeGut1 loci to verify orthology, we found many Ensembl predictions that seemed very partial. As a first pass to estimate gene model quality and completeness, we calculated a length ratio by dividing the length in bases of each zebra finch model by that of the orthologous human model (values in Additional file 1: Table S1; frequency histogram in Additional file 4A). As with the orthology analysis, we selected the longest transcript variant for zebra finch and human. If there were multiple non-overlapping zebra finch Ensembl models at a locus (i.e. split models), their lengths were summed.
We also used an alignment-based method to evaluate the completeness of zebra finch models. Because chicken (galGal6) had a higher coverage (82x) and contig N50 (17,655,422) than taeGut1 (5.5x; 38,639) and thus had more reliable and complete gene predictions, we used chicken models from Ensembl (release 95). RYR1 was excluded from this analysis as it is > 100 kb long, fragmented into at least 8 models across 8 different scaffolds that lack syntenic context in galGal5, and the only verified model in taeGut1 is highly partial and in a gappy region of chrUn. For all other genes with orthology-verified models in both zebra finch and chicken (n = 79 genes, models in both “Zebra finch locus” and “Chicken locus” columns of Additional file 1: Table S1), we first obtained the transcript sequences of the orthologous chicken Ensembl models. As with zebra finch, we selected the longest transcript variant and concatenated split models. We then aligned the chicken and corresponding zebra finch models to taeGut1 using the UCSC Browser BLAT and exported the results as psl files. Using a custom Python (version 3.6.7) script , we compared the top-scoring alignments of chicken and zebra finch models for each locus to find blocks of aligned chicken model sequence that did not overlap with any blocks of aligned zebra finch model sequence. These blocks of aligned chicken sequence missing from the zebra finch Ensembl gene predictions highlight the incompleteness of the zebra finch models and define sequences that can be used to expand these incomplete models (Additional file 5). We created a taeGut1 BED track that displays these additional sequence blocks in the UCSC genome browser (Additional file 6). For each gene in this analysis, we calculated the number of bases recovered, as well as a “percent of post-recovery length” (original model length / [original model length + number of bases recovered]; values in Additional file 1: Table S1; frequency histogram in Additional file 4B). Genes were considered problematic if the original Ensembl model length was less than 90% of the percent post-recovery length (indicated by a “*” in the “Zebra finch locus” column of Table 1, and values reported as % post-recovery length < 90% in Additional file 1: Table S1).
Animal subjects and brain tissue
Adult male zebra finches were obtained from a commercial supplier (Magnolia Bird Farm, Pasadena, CA) and acclimated in our animal facility in same-sex group cages for at least 2 weeks prior to the onset of experiments. The work described in this study was approved by the OHSU IACUC and is in accordance with NIH guidelines. To minimize the confound of song-induced gene expression, all birds were placed overnight in acoustic isolation chambers and monitored for singing for a period of 2 h after lights-on. Verified non-singing birds were sacrificed by decapitation and brains were rapidly dissected, blocked in TissueTec OCT (Sakura Finetek USA, Inc.; Torrance, CA), then flash-frozen in a slurry of isopropanol and crushed dry ice. Frozen brains were sectioned in the sagittal plane on a Leica CM1850 cryostat at 10 μm and melted onto charged microscope slides (Colorfrost Plus; Thermo Fisher Scientific; Waltham, MA). Slides were stored at − 80 °C until use.
cDNA clone selection and riboprobe synthesis
All probes were derived from clones selected from the ESTIMA zebra finch cDNA library . All genes present in taeGut1 were examined in the UCSC Genome Browser for the presence of ESTIMA clones. For genes found only in the PacBio assembly, we BLASTed each corresponding scaffold region against all NCBI zebra finch EST collections to evaluate potential evidence of expression and identify suitable cDNA clones (high alignment scores and alignments that mirrored the exon structure of aligned orthologs). Preference was given to clones representing the 3’UTR region, to minimize cross-alignment of conserved coding sequences with other gene family members. Candidate ESTs were aligned to the zebra finch genome using BLAT (UCSC Genome Browser) to confirm mapped unambiguous mapping to the target gene without significant hits to other loci. Clone IDs are provided for all genes for which an in situ was run (“EST evidence” column of Additional file 1: Table S1).
We followed the protocol by Carleton et al.  to generate non-radioactive, digoxigenin (DIG)-labeled riboprobes. In brief, isolated plasmid DNA was restriction digested using BSSHII (New England Biolabs; Ipswich, MA), purified using Invitrogen’s PureLink PCR kit (Invitrogen; Carlsbad, CA), and run on an agarose gel to verify digestion, DNA content, and correct template size. Antisense probes were synthesized by incubating template with T3 RNA polymerase (Promega; Madison, WI) and DIG labeling mix (Roche Applied Science) for 2 h at 37 °C, purified using Sephadex G-50 columns and stored at − 80 °C until use.
In situ hybridization and gene expression analysis
Non-radioactive in situ hybridizations to zebra finch sagittal brain sections were performed as in , using VectaMount Permanent Mounting solution for coverslipping. Expression for each gene was evaluated in sections from 2 to 3 different birds from a total of 10 birds. Brightfield microscopy was used to evaluate slides, and high-resolution images were acquired and uploaded to the Zebra Finch Brain Expression Atlas, ZEBrA . Genes that showed greater or lesser mRNA signal within a song nucleus relative to their surround were designated as higher or lower expression, respectively. Incubation with no probes or antisense riboprobes for genes of known expression pattern (e.g. GAD2) were routinely included in the hybridizations as negative and positive controls, respectively. The hybridization conditions used have been previously shown not to generate significant signal to various sense-strand riboprobes.
We also evaluated differential brain gene expression of ion channels based on data from four previous microarray experiments where the song nucleus of interest and a nearby region were microdissected using laser capture microscopy. The HVC dataset was generated using spotted glass ESTIMA:Song collection “20 k” cDNA microarrays as detailed in  and HVC was contrasted with the immediately ventral nidopallial shelf. Details on sample preparation, mRNA isolation, probe synthesis, and hybridization can be found in . The RA, Area X, and nXIIts datasets were generated using Agilent microarrays spotted with oligonucleotides. RA and nXIIts samples were hybridized to the Duke University Taeniopygia guttata 45 K oligo array and Area X samples were hybridized to the 20 K Agilent-019785 Custom Zebra Finch Microarray. RA was contrasted to ventrolateral arcopallium, Area X was contrasted to ventral striato-pallidum, and nXIIts was contrasted with the supra-spinal medullary nucleus (SSP). Detailed methods of sample preparation, mRNA isolation, array hybridization, and first-pass analyses are available for RA , Area X , and nXIIts .
To determine the expression of ion channel and related genes in the zebra finch song system, these microarray datasets were reanalyzed as detailed in [54, 55]. Briefly, this involved mapping the oligonucleotide sequences onto taeGut1, visualizing them as a BED track uploaded to the UCSC Genome Browser, manually annotating oligo and cDNA entries, and using in situ hybridization data from the ZEBrA database for establishing validated significance criteria of differential expression. Here we evaluated the differential microarray data for each ion channel gene present in taeGut1 by navigating to the locus to confirm oligo and cDNA annotation and verify the microarray data outcome for that gene as higher, lower, or non-differential expression in each song nucleus analyzed compared to its contrasting region. This effort led to annotation of previously unannotated oligos or correction of mis-annotated ones, contributing to the manual curation efforts outlined in [54, 55]. Differential expression tables constructed in Excel are presented in Figs. 4 5, 6. When the microarray outcome was in disagreement with the in situ hybridization data, the latter designation was used as the final outcome (indicated by “$” in Figs. 4 5, 6). Some of these discrepancies were likely due to the song nuclei and comparison regions being in the same sagittal plane for in situs but not for microarrays (in the latter, RA was compared to ventrolateral arcopallium, and Area X was compared to ventromedial striato-pallidum). In some other cases, probes were not available for in situs, so expression data reflects microarray data only (indicated by asterisks in Figs. 4 5, 6).
Sodium, calcium, and chloride channel genes in the zebra finch genome
Genes that encode or modulate sodium, calcium, and chloride ion channels play important roles in determining the electrophysiological properties of brain circuits. In order to understand the potential contributions of these channels to functional circuits underlying vocal learning and production in songbirds, we started with a comprehensive analysis to define the full complement of these ion channel families in the zebra finch genome. Starting with a list of 108 ion channel genes previously identified in human plus one additional gene that lacks a human model (TPCN3), we used BLAT alignments and synteny analysis to identify orthologous loci in zebra finch and chicken. We identified 94 (23 sodium, 38 calcium, 33 chloride) of these genes in zebra finch (Table 1) and 98 in chicken (Additional file 1: Table S1), the latter encompassing four genes (BEST2, BEST4, CLCA1, and CLCA2) not found in zebra finch (Table 2 and Additional file 1: Table S2). Of the identified zebra finch genes, 54 were predicted and correctly annotated by Ensembl (taeGut3.2.4 release 95), 27 genes comprising 34 models were annotated as “novel genes” in Ensembl (indicated by a “^” in the “Zebra finch locus” column of Table 1), and 14 had no Ensembl predictions in taeGut1. Of these 14, four genes were found in taeGut1 through BLAT alignments of human or other species’ orthologs (indicated by a taeGut1 locus only in the “Zebra finch locus” column of Table 1.) The remaining 10 genes were not in taeGut1, but were identified in the recently released zebra finch PacBio assembly (indicated by a PacBio locus only in the “Zebra finch locus” column of Table 1) through alignments of orthologous models. For five of these genes only very partial alignments were found, suggesting gene truncation (indicated by a “▽” in Table 1.) In most cases we established orthology by demonstrating reciprocal top-scoring alignments and similar syntenies in birds and humans. For 8 genes, we linked the zebra finch synteny with human through phylogenetic relationships (indicated by a “$” in Table 1), but for 2 genes we could not fully link the synteny in zebra finch with that of human (indicated by a “#” in Table 1, detail provided in “Synteny notes” column of Additional file 1: Table S1). For a total of 14 genes in this study, we used orthologs from other species to find the locus in zebra finch (indicated by a “◆” in Table 1.) While no novel paralogs were found in zebra finch, we confirmed the presence of TPCN3 in birds and reptiles, a gene that is severely truncated in humans and chimp and completely missing from mouse and rat .
For 15 mammalian genes we did not find orthologs in zebra finch, but we established a possible origin of the loss by exploring outgroup species (Table 2; details in Additional file 1: Table S2). These genes are discussed further in the context of their respective subfamilies. We note these conclusions about gene losses are tentative given that many avian genomes are incomplete and phylogenetic representation of sequenced genomes is somewhat limited.
Sodium channel genes (Table 1A)
Sodium voltage-gated channel alpha (α) subunits (SCNxA):
The subunits encoded by these genes form the pore of voltage-gated sodium channels , which play important roles in the initiation and rising phase of action potential (AP) generation, and are critical for AP propagation . This family originated from a singular ancestral form that has undergone multiple duplications and diversified broadly [57, 58]. We found 9 of the 10 mammalian SCNxA genes in zebra finch and chicken. The missing gene was SNC7A, which arose in mammals and presumably functions not as a voltage-gated channel but rather as a sodium sensor .
Sodium voltage-gated channel beta subunits (SCNxB):
Four genes encoding non-pore-forming, auxiliary β subunits have been described in mammals . They are thought to modulate sodium currents by associating covalently (SCN1B and 3B) or non-covalently (SCN2B and 4B) with α subunits to form heteromeric complexes. We found all four β subunit genes in zebra finch and chicken, and note that in the latter, it was previously identified in a PacBio assembly  . SCN1B was missing from taeGut1 but limited evidence suggested its presence in songbirds . We have now identified the complete gene and syntenic context in the PacBio assembly through ortholog alignments and verification of conserved synteny (Fig. 2 and Additional file 3). The synteny of SCN1B throughout avian phylogeny was difficult to resolve as many avian models were either on short scaffolds with little syntenic context or on scaffolds with many gaps, but we identified a consistent group of syntenic genes across lineages.
Acid-sensing ion channels (ASIC):
ASICs are proton-activated, nonselective cation channels that pass sodium, lithium, and potassium ions with the highest affinity for sodium . All 5 mammalian ASIC genes were found in both zebra finch and chicken genomes.
Sodium leak channel (NALCN):
NALCN is a single gene that encodes a widely-expressed, nonselective, voltage-insensitive cation channel. This channel carries a persistent sodium leak current that sets excitability thresholds in virtually all neurons . It is present in zebra finch and chicken genomes.
Epithelial sodium channels (SCNN):
Channel subunits in the SCNN family are encoded by four genes and drive active sodium reabsorption from extracellular fluid in epithelial tissues . We found that all four genes are present in zebra finch and chicken genomes, confirming Hanukoglu et al.  and strengthening the evidence for orthology by confirming synteny.
Calcium channel genes (Table 1B)
Calcium voltage-gated channel alpha (α) -1 subunits (CACNA1x):
Genes in this family encode the pore-forming α1 subunits that carry calcium across the plasma membrane . In mammals, α1 subunits belong to three major subfamilies: Cav1 (L-type) channels (4 genes), Cav2 (P/Q-, N-, and R- type) channels (3 genes), and Cav3 (T-type) channels (3 genes). All genes were found in zebra finch and chicken, with CACNA1A and CACNA1F in the PacBio assembly.
Calcium voltage-gated channel alpha-2-delta (α2/δ) subunits (CACNA2Dx):
Calcium channels are modified by non-conducting subunits, including α2 and δ. α2/δ genes are unique in that each gene codes for both the α2 and δ subunits: α2 is an extracellular glycoprotein that forms a disulfide linkage to the δ subunit, and the δ subunit keeps the complex tethered to the plasma membrane . All four human α2/δ genes were found in zebra finch and chicken.
Calcium voltage-gated channel auxiliary beta (β) subunits (CACNBx):
These modulatory subunits are entirely intracellular and attached to the calcium channel complex through binding sites on the α1 subunit . All four human genes were found in zebra finch and chicken. CACNB1 in zebra finch was found only in the PacBio assembly.
Calcium channel auxiliary gamma (γ) subunits (CACNGx):
Gamma subunits are glycoproteins with four transmembrane domains that associate with calcium channel complexes in the membrane . Of the eight human genes, three (CACNG6, CACNG7, and CACNG8) occur in a syntenic block within a region of human Chr19 that is largely missing in birds, possibly due to chromosomal rearrangements in ancestral archosaurs and/or dinosaurs . Subsequent assessments have identified CACNG7 and CACNG8 but not CACNG6 in the galGal5 chicken assembly . Here we have confirmed with high confidence that CACNG6 is missing in zebra finch and other avian lineages based on BLAST searches and synteny verification of bird genomes including chicken and zebra finch PacBio assemblies (Fig. 3), suggesting an ancestral avian loss. We also found evidence that CACNG7 is present and likely complete in zebra finch, while CACNG8 is present but appears truncated in both zebra finch and chicken.
Sperm associated cation channels (CATPSER):
This family encodes four α subunits (CATSPER1–4) and five auxiliary subunits (CATSPER B, D, E, G, & Z) essential to sperm activation and motility [70, 71]. In zebra finch, we found fragments of CATSPER3, B, D, and E, but no trace of CATSPER1, 2, 4, G and Z. Chicken shares the same pattern, but retains fragments of CATSPER2. Because CATSPER1, − 4, and -G are present in ratites, their absence in zebra finch and chicken appears to reflect loss in Neognaths. Our conclusions are partially consistent with Cai et al. , Chung et al. , and Warren et al. , but due to the higher quality PacBio scaffolds and examination of additional birds, we provide stronger evidence with regards to truncation and phylogenetic patterns of gene loss. We note that the chicken gene annotated as CATSPER2 in NCBI (Gene ID: 395985) is misannotated – the correct annotation is CACNA1S.
Other calcium channel genes (ITPR, RYR, TPCN):
Other calcium channels may impact neuronal function by regulating internal calcium stores, including inositol 1,4,5-trisphosphate receptors (ITPR), ryanodine receptors (RYR), and two-pore channels (TPCN) [56, 73, 74]. Each family has 3 genes, all of which were identified in zebra finch and chicken. TPCN3 appears to have been lost in rodents and truncated in primates, including humans . Interestingly, the truncated gene appears to have been duplicated in human, as can be seen in the UCSC Browser (GRCh38/hg38 chr2:109947776–109,974,086 and chr2:110397405–110,423,734).
Chloride Channel genes (Table 1C)
Voltage-dependent chloride channels (CLCN):
This family is composed of channels that open and conduct chloride ions in response to changes in membrane potential . Mammals have 10 CLCN genes, of which 9 were identified in zebra finch and chicken. Two adjacent human genes - CLCNKA and CLCNKB - have high sequence identity and appear to be the result of a duplication specific to mammals. Both chicken and finch have a single gene at this locus that mirrors the ancestral condition. Because of their high similarity, sequence alignments could not resolve which human paralog is the ancestral form present in birds or reptiles. As the nomenclature in ancestral species is inconsistent, we suggest the name “CLCNK” for the single gene found at this locus in birds and reptiles.
Chloride intracellular channels (CLIC):
These channels are integrated into the membranes of organelles and participate in processes like membrane trafficking and cytoskeleton dynamics . They can also exist in a soluble state in the cytoplasm, the function of which is poorly understood. Of the 7 functional genes in mammals, we found all but CLIC1 in chicken and zebra finch. We could not conclusively establish a CLIC1 loss because the syntenic genes are at the ends of different scaffolds in both zebra finch and chicken PacBio assemblies. It thus remains possible that CLIC1 is not in these assemblies due to incompleteness of the scaffolds. Consistent with this possibility, CLIC1 and syntenic genes are present in alligator and in kiwi (a ratite). Alternatively, the gene was lost in neognaths, as kiwi is currently the only bird with a detectable CLIC1.
Calcium-activated chloride channels (ANO, BEST):
Some anoctamin (ANO1 & ANO2) and bestrohpin (BEST2) genes encode channels that conduct chloride currents in response to surges in intracellular calcium, but the function of most other family members remains speculative [76, 77]. Mammals have 10 ANO genes, of which we found 9 in chicken and zebra finch, and 4 BEST genes, of which we found 4 in chicken and 2 in zebra finch.
ANO7 lacks models in zebra finch and chicken, but traces were detected by aligning orthologs to the zebra finch and chicken PacBio assemblies. The synteny in these and other birds (kiwi, goose, golden eagle, and rifleman) is conserved with mammals and alligators. Given that only a few 3′ exons from orthologous models of this gene aligned to the chicken and zebra finch PacBio assemblies, combined with the lack of EST-based expression evidence, we conclude that ANO7 has likely become a pseudogene in some avian lineages, including oscines and galliformes.
BEST4 is present and complete in galliforms (chicken) and psittaciformes (budgie), but appears pseudogenized in passerines, as only short segments of the 5′ and 3′ UTR regions were detected in zebra finch (PacBio), Bengalese finch, great tit, and starling. While BEST2 seems complete in chicken and songbirds like Tibetan tit and starling, we did not detect it at the corresponding locus in the zebra finch PacBio assembly. The upstream genes are in zebra finch similarly as in other birds, reptiles, and mammals, while the genes downstream of this locus reflect a synteny shared by Tibetan tit and starling that is presumably unique to songbirds. Because there are no BEST2 models and syntenic genes are isolated on short scaffolds in Bengalese and medium ground finch, we cannot distinguish between a loss in the finch lineage and a zebra finch-specific loss. Similar to CACNG6/7/8, BEST2 is located in a region of human chromosome 19 associated with extensive rearrangements and gene losses in birds .
Chloride channel accessory subunits (CLCA):
Mammals have several non-pore-forming accessory subunits with a broad functional repertoire that includes modulating chloride channels, cell adhesion, and tumor suppression . In humans, three functional CLCAs are in a syntenic cluster (CLCA1, CLCA2, and CLCA4) along with a pseudogene (CLCA3P). We found at least one CLCA gene with conserved synteny in psittaciformes (budgie), accipitriformes (eagles), and more basal birds (chicken, guinea fowl, and tinamou). In zebra finch and other passerines (e.g. golden-collared manakin), the immediately syntenic genes were adjacent to one another, with no trace of CLCA and no gaps in the intergenic region. From this, we conclude that CLCA genes were likely lost in Passeriformes.
Volume-regulated chloride channels (LRRC8):
LRRC8 channels are activated by cell swelling and may influence extracellular signaling by transporting neurotransmitters [79, 80]. Mammals have five LRRC8 subunits, and all but one (LRRC8E) were found in zebra finch and chicken. This gene appears to have been lost in birds .
Other chloride channel genes:
Cystic fibrosis transmembrane conductance regulator (CFTR) is a chloride channel gated by ATP . Chloride nucleotide-sensitive channel 1A (CLNS1A) is expressed at relatively high levels in human  and mouse  brain tissue, and strongly implicated in cell volume regulation . Both genes are present in zebra finch and chicken.
Evaluating zebra finch ion channel gene models
While examining zebra finch gene orthology, we noticed that many human models were noticeably longer than the corresponding zebra finch Ensembl models, suggesting that the latter may be incomplete predictions. To record this difference, we calculated a ratio (Additional file 1: Table S1 and 3A) of zebra finch to human transcript model length, and found that 28 of the 80 existing zebra finch Ensembl models were less than 90% of the human length, and 19 were less than 75%, even after summing the lengths of partial models, when present. To further evaluate model completeness, we compared the alignments of orthologous chicken and zebra finch Ensembl models at each gene locus in taeGut1 (n = 79 genes with an Ensembl model in both species). Because Ensembl chicken models are from a higher quality genome and are generally more complete than zebra finch models, this method revealed sequence blocks where the chicken model aligned but that are not present in the zebra finch Ensembl models (example in Additional file 5). We recovered these additional sequences for 66 of the 79 genes analyzed (Additional file 1: Table S1, “Additional bases recovered” > 0), and created a taeGut1 BED track for their visualization in the UCSC Browser (Additional file 6). We also calculated a percent of post-recovery length (Additional file 1: Table S1; Additional file 4B) and found 36 genes for which the original model length was less than 90% of the post-recovery length (indicated by a “*” in the “Zebra finch locus” column of Table 1), and 16 that were less than 75%. This effort demonstrates that many zebra finch Ensembl models are incomplete, but additional genomic sequence related to those genes exists and should be taken into account in future studies.
Differential expression of sodium, calcium, and chloride channel genes in the song system
To assess the expression of ion channel genes in the song system of adult male zebra finches, we analyzed previously generated microarray data and conducted in situ hybridizations on brain sections containing the 4 major telencephalic song nuclei (HVC, RA, LMAN, and Area X), and a brainstem motor nucleus, nXIIts. Because our focus was on vocal control circuitry, not all brain areas of relevance for vocal learning (such as auditory forebrain) were assessed. This analysis encompassed 55 genes (14 Na, 25 Ca, 16 Cl) that were part of the microarrays and/or for which cDNA clones were available in zebra finch. Summaries of the data are presented in Figs. 4, 5, 6, and robust examples of differential in situ patterns are presented in Figs. 7 8, 9, 10, and 11 Additional file 7. In situ images can also be accessed through ZEBrA . For a few genes, in situs also revealed differential expression in DLM and Uva. For all genes analyzed, the terms higher, lower and non-differential below refer to expression compared to the respective surrounds (nidopallium for HVC, LMAN and NIf, arcopallium for RA, striatum for Area X, dorsal thamalus for DLM and Uva, rostral medulla for nXIIts).
The four alpha subunits with the highest brain expression in mammals (SCN1A/2A/3A/8A) showed varied patterns within the song system (Fig. 4). SCN1A was higher in RA (Fig. 7b), HVC (Fig. 8b), LMAN (Fig. 9j), Area X, and DLM (Fig. 10b). SCN2A was lower in RA (Additional file 7), Area X, LMAN, and DLM (Fig. 10c), but non-differential in HVC (Fig. 8c). SCN3A was lower in RA (Fig. 7d), HVC (Additional file 7), and LMAN. SCN8A was higher in RA (Additional file 7) and HVC (Additional file 7). In contrast, SCN5A was non-differential in all nuclei, and SCN9A was non-differential except in nXIIts where it was higher and the only alpha subunit to show differential expression in this nucleus.
SCN1B showed higher expression in all telencephalic nuclei, including NIf, as well as in DLM (Fig. 11). In contrast, SCN2B was non-differential across song nuclei (e.g. RA in Fig. 7c) except for higher expression in HVC (Fig. 8c). SCN3B was an exquisite negative marker of all telencephalic song nuclei (Figs. 7d, 8d, 9b, and 10c), but higher in nXIIts. Interestingly, SCN3B was highly expressed in a sparse cell population in Area X (Additional file 7). SCN4B was higher in RA (Fig. 7e), HVC, and LMAN (Fig. 9k). Area X showed lower expression of ASIC1 and ASIC4, and ASIC2 was lower in nXIIts. NALCN, the only sodium channel that is a leak channel, showed lower expression in Area X.
Among Cav1s, CACNA1C was higher in Area X and nXIIts, while CACNA1D was lower in RA, Area X, and LMAN. Among Cav2s, CACNA1E was lower in RA (Fig. 7f) and LMAN, while CACNA1B was higher in HVC and Area X (Fig. 9c) where it was particularly strong in a sparse population of Area X cells (Additional file 7). Differential expression of Cav3’s was biased toward lower expression, except for higher expression of CACNA1I (Cav3.3) in HVC (Additional file 7), the only calcium channel α subunit to exhibit opposing expression in two different song nuclei. All three Cav3 channels were lower in RA, the only nucleus to show differential expression in the same direction of every gene within a channel subfamily.
CACNA2D1 and CACNA2D2 were complementary in their differential expression: CACNA2D1 was lower in HVC (Fig. 8f) and LMAN, and CACNA2D2 higher in HVC and LMAN (Fig. 9d). In Area X, CACNA2D2 was non-differential but strongly expressed in a sparse cell population, with a similar pattern throughout the striatum (Fig. 9d). CACNA2D2 was also higher in RA, and CACNA2D3 higher only in Area X.
The only calcium channel β subunit to show lower expression was CACNB2, which was turned down in RA, HVC (Additional file 7), and LMAN (Fig. 9l), but higher in DLM (Additional file 7). CACNB2 was also generally lower throughout much of the anterior striatum, but a sparse cell population showed enhanced expression (Additional file 7). CACNB3 was higher in HVC and LMAN and CACNB4 higher in HVC (Fig. 8e) and Area X (Fig. 9e).
Differential expression of γ subunits in the song system was predominantly in the AFP. CACNG5 was lower in LMAN (Fig. 9g) and higher in Area X, and CACNG3 was lower in Area X and LMAN. CACNG5 also showed enhanced expression in a sparse cell population of Area X (Additional file 7), and in Uva (Fig. 10d). CACNG4 was lower in HVC, LMAN and Area X (Fig. 9f), but higher in nXIIts. All γ subunit genes examined were non-differential in RA (e.g. Additional file 7), except for CACNG7, which was higher in RA as well as in HVC and Area X.
Differential expression for TPCN and RYR genes was always higher, including RYR2 in HVC (Additional file 7), Area X and LMAN (Fig. 9h), and RYR3, TPCN − 1 and − 2 in HVC. Only nXIIts showed differential expression of ITPRs, which were lower.
The majority of identified chloride channels lacked zebra finch ESTs. Among assessed genes, most were non-differential, except for three genes with higher expression - ANO5, CLIC4, and CLNS1A - which were all higher in HVC, with ANO5 and CLIC4 also higher in RA, and one gene with lower expression, LRRC8, in RA and Area X. Thus, differential expression of chloride channel genes was biased toward nuclei of the DMP and higher expression.
In spite of the wealth of information on song system connectivity and neurophysiology, an understanding of the molecular basis for the intrinsic electrophysiological properties of its constituent neurons is lacking. Here we report on the genomics and song system expression of sodium, calcium, and chloride ion channel genes in the zebra finch, including pore-forming and modulatory subunits. This effort annotated 34 “novel gene” Ensembl models, detected additional sequences for several partial gene models in taeGut1, and established what we believe is the full complement of these gene families in zebra finch and chicken, with no novel paralogs detected. PacBio assemblies were particularly helpful to: (1) identify previously undetected genes; (2) confirm with high confidence gene truncation; (3) define orthology based on synteny; (4) suggest a likely phylogenetic path to previously reported avian gene losses or truncations [41, 61]. A large number of zebra finch gene models are partial relative to human and chicken models, but we detected additional sequences in taeGut1 that could expand those models. Incomplete gene predictions (indicated by “*” in Table 1) cast doubt on “off-the-shelf” dN/dS values from automated algorithms. Thus, conclusions regarding positive selection or rapid gene evolution in these families (e.g. [42, 84, 85]) should be treated as provisional until more complete gene models are available.
Each song nucleus examined differentially expresses a unique complement of ion channels, suggesting that diverse expression of these genes is fundamental to shaping excitability within the song system. With few exceptions, differential expression of individual genes was in the same direction across nuclei. Area X and HVC showed enhanced expression of more calcium channel genes than other song nuclei, suggesting a more prominent role of calcium currents in governing electrophysiological properties of the AFP. In contrast, chloride channel genes were mostly non-differential. Compared to other nuclei, nXIIts differentially expressed the fewest number of genes. We note that most genes that lacked brain expression evidence in zebra finch (e.g. SCN4A, CACNG1, and CLCN1) are primarily expressed in non-neural tissues (e.g. skeletal muscle) in other organisms [38, 39].
Some genes were strongly expressed in sparse cells and might contribute to defining the unique neuronal populations reported in Area X [12, 86, 87], HVC [8, 13], and RA . We acknowledge, however, that several of the genes in this study participate in processes other than excitability (e.g. cell signaling and metabolism), and expression in specific cell types or non-neuronal cells (e.g. glia) was not assessed. Furthermore, protein-based assays and coexpression analysis are needed to determine the subcellular localization of the channels and their possible interactions. We did not evaluate splice variants, but note that in mammals, genes like T-type calcium channels give rise to variants with distinct gating properties and spatiotemporal regulation . We also note that while differential expression provides clues as to regulatory mechanisms, non-differentially expressed ion channels may still contribute to excitability. Below we discuss the main implications of the brain expression findings with regards to gene families.
Based on its known role in in mammals, the higher expression of SCN1A (Nav1.1) in all major song nuclei suggests an essential function of somatodendritic signal integration and action potential (AP) thresholding in the song system . SCN1A also confers fast spiking in a diversity of GABAergic interneurons [90,91,92]. Because in most nuclei SCN1A was enhanced in small, sparse cells reminiscent of interneurons, the enhanced expression of SCN1A may contribute to the fast-spiking behavior of interneurons that has been described in multiple song nuclei [8,9,10,11, 86]. In mammals, SCN8A (Nav1.6) facilitates AP initiation and propagation [89, 93,94,95]. Its higher expression in RA and HVC could thus reinforce rapid and reliable axonal propagation in the DMP. SCN3A (Nav1.3) brain expression is high early in mammalian development but is progressively replaced by SCN1A, SCN2A, and SCN8A [96,97,98]. In adult zebra finch, SCN3A is widely expressed throughout the brain but lower in RA, HVC and LMAN - a species discrepancy that highlights the diversity in sodium channel expression across phylogenies.
β subunits alter cell surface expression, gating properties, and voltage sensitivities of α subunits , and high levels of β subunits are likely necessary to shuttle and modulate a high volume of voltage-gated sodium channels. Thus, it is not surprising that some α and β subunits were highly expressed within the same song nuclei. In mammals, SCN1B increases the surface expression and current densities of α subunits [60, 99, 100]. It also supports resurgent sodium currents that enable high-frequency firing , and reinforces AP initiation by targeting SCN8A to the axon initial segment . Additionally, SCN1B can modulate voltage-gated potassium (Kv) family channels [102,103,104], several of which are higher in multiple song nuclei, including Kv1.1, Kv4.2, and Kv4.3 . Collectively these features suggest that SCN1B may be an important regulator of neuronal excitability in song nuclei by supporting fast, reliable spiking. HVC was the only nucleus to show higher expression of both SCN1B and SCN2B, consistent with the higher expression of α subunits in HVC, and the observation that the sodium current-boosting effects of SCN2B require SCN1B coexpression [60, 99].
SCN4B induces resurgent sodium current, an important adaptation for neurons that exhibit high-frequency firing [95, 105,106,107,108,109]. Its high expression suggests that resurgent currents might support the high-frequency and bursty firing characteristic of neurons in RA [110, 111], HVC [112, 113], and LMAN . SCN4B also supports persistent sodium currents , which might account for the non-inactivating sodium conductances of HVC neurons . While the effects of SCN3B on excitability are inconsistent [115,116,117,118,119,120], its high brain expression in adult zebra finches contrasts with rodents, where expression is high during development and lower in adulthood . Much like SCN3A, this stark divergence exemplifies how gene expression patterns can differ in birds versus mammals. As for ASICs, little is known about the role of acid signaling in neurotransmission, although there is evidence of their role in dendritic spine remodeling . Lastly, leak currents carried by NALCN play roles in rhythmic and spontaneous firing [63, 122], thus NALCN differential expression might contribute to the heterogeneity in spontaneous firing rates across Area X cell types .
Calcium channels influence input processing in dendrites and signal transmission in synaptic terminals, which links them to plasticity mechanisms that serve learning and memory . Voltage-gated calcium channels (Cav’s) are vital to dendritic processing, burst firing, and neurotransmitter release. Among pore-forming alpha subunits, L-type (Cav1) channels conduct large, long-lasting currents in response to strong depolarization . The higher expression of the L-type channel CACNA1C suggests that integration of synaptic input is important in Area X and nXIIts. CACNA1C can also associate with BK channel KCa1.1 (KCNMA1) to modulate repolarization and spiking frequency [124, 125]. While not differential, KCNMA1 is expressed in Area X  and may interact with CACNA1C to shape repolarization and repetitive firing. Compared to CACNA1C, CACNA1D activates at lower membrane potentials and is slower to inactivate . Thus, its lower expression may refine the timing and restrict the duration of depolarization-induced calcium influx in RA, Area X, and LMAN.
T-type (Cav3) channels create low-threshold calcium potentials that are integral to driving rebound and burst firing in neurons [88, 123, 127]. The lower expression of all T-type channels in RA is intriguing given that low expression of T-type channels is associated with single or tonic AP firing , yet RA neurons exhibit bursty, high-frequency firing [110, 111, 128]. This suggests that high-frequency firing in RA might instead be driven by voltage-gated sodium channel beta subunits (SCN1B, SCN4B) and/or potassium channel subunits like Kv3.1 [27, 129], which show high expression in RA [33, 50], or by hyperpolarization-activated (HCN; Ih) channels as suggested by electrophysiology studies of HVC . HVC was the only nucleus to highly express a T-type channel (CACNA1I / Cav3.3), providing a candidate for the low-threshold calcium conductances recorded in HVC’s Area X-projecting neurons and interneurons [4, 130]. Relative to other T-type channels, CACNA1I channels are the most resistant to attenuation during high-frequency firing, exhibit current amplitude facilitation, and have more positive thresholds of activation and inactivation [88, 131, 132]. These properties might help to regulate excitability in HVC, which shows temporally-precise bursting during production of mature song . CACNA1I is also an attractive candidate for development studies given that T-type calcium channels may underlie changes in intrinsic electrophysiological properties of HVC neurons during song learning .
CACNA1B, an N-type calcium channel, has a large single-channel conductance and is closely associated with the calcium-sensitive presynaptic vesicle release machinery [133,134,135,136]. The distinctive higher expression of CACNA1B in HVC and Area X suggests that these song nuclei are specialized for rapid and efficient neurotransmission. CACNA1E, an R-type calcium channel, functions in both dendritic calcium signaling and neurotransmitter release [35, 137]. In mammalian hippocampal neurons, these channels activate potassium channels within dendritic spines to limit depolarization in response to excitatory post-synaptic potentials [138, 139]. Additionally, CACNA1E−/− mice are more resistant to seizure , possibly due to dampened excitability . Depending on CACNA1E subcellular localization, lower expression in RA may translate to larger calcium transients in dendritic spines and/or reduce the potential for run-away excitability.
The modulatory α2δ (CACNA2D) and β (CACNB) subunits increase calcium currents as well as vesicle release probability [66, 142,143,144,145]. Their widespread high expression, particularly in the AFP, might thus relate to more efficient neurotransmission. The nearly identical expression patterns of α2δ and β in HVC and LMAN contrast with the divergent patterns of α1 subunits in these nuclei, suggesting that differences in calcium conductance might arise from unique combinations of α1 subunits rather than different auxiliary subunits. In turn, the higher expression of both modulatory α2δ subunits and pore-forming α1 subunits in Area X and HVC suggests that high levels of both these subunits are required. Finally, evidence that CACNA2D3 participates in neurexin-mediated retrograde signaling  suggests Area X might utilize a α2δ-dependent, inhibitory synaptic feedback mechanism that is unique among song nuclei.
Most γ (CACNG) subunits are transmembrane AMPA receptor (AMPAR) regulatory proteins, or TARPs, which modulate gating and synaptic targeting of AMPARs [147, 148]. A sparse population of Area X cells that resemble the large DLM-projecting neurons in Area X  showed enhanced expression of CACNG5, a TARP unique in its ability to increase AMPAR desensitization and deactivation rates . Intriguingly, a population of sparse, large cells in Area X also expresses GluR4, the AMPAR subunit with the fastest kinetics , further suggesting that temporally precise signaling through AMPARs is particularly important for Area X projection neurons. In contrast, CACNG4 and CACNG8 are the most effective TARPs in slowing the desensitization, deactivation, and mEPSC decay rates of AMPARs . The lower expression of CACNG4 in HVC, Area X, and LMAN and lack of a functional CACNG8 gene suggest that these nuclei rely on alternate mechanisms to elicit slow AMPA kinetics, such as differential splicing of AMPA subunits into their “flip” or “flop” forms . CACNG7 selectively shuttles to the plasma membrane and boosts currents of calcium-permeable AMPARs while inhibiting membrane expression of calcium-impermeable ones [151, 152]. Thus, biasing AMPAR pools toward calcium permeability and enhancing calcium flow through AMPARs might be a shared specialization in RA, HVC, and Area X. In addition, CACNG7 can dampen CACNA1B (Cav2.2) expression in heterologous expression systems, but not in sympathetic neurons . Accordingly, CACNA1B was not lower in any nuclei that showed higher CACNG7 expression. Finally, the non-TARP CACNG1 and CACNG6 decrease calcium currents by interacting with voltage-gated calcium channels. Birds lack CACNG6, but CACNG1 could potentially compensate for this gene loss.
The high expression of multiple intracellular calcium channels suggests that mobilization of internal calcium stores is important for HVC. While the role of TPCNs in neural processing is unknown, inositol 1,2,5-trisphosphate (IP3) receptors (ITPRs) and ryanodine receptors (RYRs) are implicated in plasticity, learning, and memory . The higher expression of RYR2 and CACNA1C in Area X is also intriguing, given that RYR2 interacts with CACNA1C in muscle cells to coordinate excitation-contraction coupling . Area X cells might thus funnel calcium ions through plasma membrane-bound CACNA1C to ER-bound RYR channels as a way to amplify calcium signaling. Lastly, ITPRs and RYRs do not overlap in location or direction of expression, suggesting the recruitment of distinct intracellular calcium stores and signaling pathways across song nuclei .
Of the gene families studied here, chloride channels are the least characterized and had a low incidence of differential expression. Considering their role in stabilizing the resting potential [36, 157], the higher expression of chloride channels in HVC and RA might help the membrane potential to reset despite large voltage fluctuations.
Passeriformes appear to lack CLCA genes. In mammals, these calcium-activated proteases are expressed predominantly in intestine and airway epithelium, and only a few are expressed in brain . The absence of CLCAs suggests that songbirds rely on other mechanisms to regulate chloride transport in respiratory, digestive, and excretory tissues. Some CLCAs are expressed in olfactory sensory neurons (OSNs) in mammals, where they modulate chloride currents involved in olfactory signaling [158, 159]. Contrary to previous notions, birds use olfactory cues [160,161,162] and have large repertoires of olfactory receptor genes [42, 163,164,165], but downstream signal transduction mechanisms are unknown. In mammals, ~ 90% of the transduction current in OSNs is carried by the calcium-sensitive chloride channel ANO2 (TMEM16B) [166,167,168]. Because ANO2 is present in zebra finch and other avian species, birds appear to have the necessary machinery for olfactory transduction. Moreover, BEST2 - once the strongest candidate for transduction in OSNs  - is missing from zebra finch, consistent with the conclusion that BEST2 is not required for olfaction . In sum, while the main carrier of olfactory transduction current (ANO2) may be conserved between mammals and birds, the lack of CLCAs suggests that regulation of this current may be markedly different in Passeriformes.
We report here on the genomic identity and song system expression of three major ion channel families. The majority of these genes are conserved between birds and mammals, and the expression data paint a picture of abundant differential patterns, especially for sodium and calcium channel genes, in circuitry that supports the acquisition and production of learned vocalizations. We note that while changes in synaptic and intrinsic properties during the song learning period have been reported [11, 110, 111, 171], further studies of developmental regulation of ion channel genes in song nuclei are needed [172, 173]. Finally, this investigation establishes a foundation and framework for future studies to test the contribution of specific ion channel genes to distinct excitability properties of vocal nuclei in songbirds.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and Additional files 1, 2, 3, 4, 5, 6, 7. A GitHub repository for the sequence recovery analysis is available from . Complete microarray datasets are available from  and details regarding curation methods can be found in . High-resolution images of most of the in situ data presented here are publicly available on ZEBrA .
Anterior forebrain pathway
Medial nucleus of the dorsolateral thalamus
Dorsomedial nucleus of the intercollicular complex
Direct motor pathway
Lateral magnocellular nucleus of the anterior nidopallium
Nucleus interfacialis of the nidopallium
Tracheosyringeal portion of the hypoglossal nerve nucleus
Robust nucleus of the arcopallium
Thalamic nucleus uvaeformis
Petkov CI, Jarvis ED. Birds, primates, and spoken language origins: behavioral phenotypes and neurobiological substrates. Front Evol Neurosci. 2012;4:12.
Zeigler HP, Marler P. Neuroscience of birdsong: Cambridge University Press; 2008.
Zeigler HP, Marler P. Behavioral neurobiology of bird song: New York Academy of Sciences; 2004.
Daou A, Ross MT, Johnson F, Hyson RL, Bertram R. Electrophysiological characterization and computational models of HVC neurons in the zebra finch. J Neurophysiol. 2013;110(5):1227–45.
Kadakia N, Armstrong E, Breen D, Morone U, Daou A, Margoliash D, et al. Nonlinear statistical data assimilation for HVC[formula: see text] neurons in the avian song system. Biol Cybern. 2016;110(6):417–34.
Meliza CD, Kostuk M, Huang H, Nogaret A, Margoliash D, Abarbanel HD. Estimating parameters and predicting membrane voltages with conductance-based neuron models. Biol Cybern. 2014;108(4):495–516.
Nogaret A, Meliza CD, Margoliash D, Abarbanel HD. Automatic construction of predictive neuron models through large scale assimilation of electrophysiological data. Sci Rep. 2016;6:32749.
Dutar P, Vu HM, Perkel DJ. Multiple cell types distinguished by physiological, pharmacological, and anatomic properties in nucleus HVc of the adult zebra finch. J Neurophysiol. 1998;80(4):1828–38.
Spiro JE, Dalva MB, Mooney R. Long-range inhibition within the zebra finch song nucleus RA can coordinate the firing of multiple projection neurons. J Neurophysiol. 1999;81(6):3007–20.
Livingston FS, Mooney R. Development of intrinsic and synaptic properties in a forebrain nucleus essential to avian song learning. J Neurosci. 1997;17(23):8997–9009.
Goldberg JH, Fee MS. Singing-related neural activity distinguishes four classes of putative striatal neurons in the songbird basal ganglia. J Neurophysiol. 2010;103(4):2002–14.
Luo M, Perkel DJ. A GABAergic, strongly inhibitory projection to a thalamic nucleus in the zebra finch song system. J Neurosci. 1999;19(15):6700–11.
Kubota M, Taniguchi I. Electrophysiological characteristics of classes of neuron in the HVc of the zebra finch. J Neurophysiol. 1998;80(2):914–23.
Nottebohm F, Kelley DB, Paton JA. Connections of vocal control nuclei in the canary telencephalon. J Comp Neurol. 1982;207(4):344–57.
Nottebohm F, Stokes TM, Leonard CM. Central control of song in the canary, Serinus canarius. J Comp Neurol. 1976;165(4):457–86.
Bottjer SW, Miesner EA, Arnold AP. Forebrain lesions disrupt development but not maintenance of song in passerine birds. Science. 1984;224(4651):901–3.
Brainard MS, Doupe AJ. Interruption of a basal ganglia-forebrain circuit prevents plasticity of learned vocalizations. Nature. 2000;404(6779):762–6.
Sohrabji F, Nordeen EJ, Nordeen KW. Selective impairment of song learning following lesions of a forebrain nucleus in the juvenile zebra finch. Behav Neural Biol. 1990;53(1):51–63.
Scharff C, Nottebohm F. A comparative study of the behavioral deficits following lesions of various parts of the zebra finch song system: implications for vocal learning. J Neurosci. 1991;11(9):2896–913.
Cardin JA, Raksin JN, Schmidt MF. Sensorimotor nucleus NIf is necessary for auditory processing but not vocal motor output in the avian song system. J Neurophysiol. 2005;93(4):2157–66.
Williams H, Vicario DS. Temporal patterning of song production: participation of nucleus uvaeformis of the thalamus. J Neurobiol. 1993;24(7):903–12.
Pfenning AR, Hara E, Whitney O, Rivas MV, Wang R, Roulhac PL, et al. Convergent transcriptional specializations in the brains of humans and song-learning birds. Science. 2014;346(6215):1256846.
Vicario DS. Organization of the zebra finch song control system: II. Functional organization of outputs from nucleus Robustus archistriatalis. J Comp Neurol. 1991;309(4):486–94.
Bottjer SW, Halsema KA, Brown SA, Miesner EA. Axonal connections of a forebrain nucleus involved with vocal learning in zebra finches. J Comp Neurol. 1989;279(2):312–26.
Vates GE, Nottebohm F. Feedback circuitry within a song-learning pathway. Proc Natl Acad Sci U S A. 1995;92(11):5139–43.
Johnson F, Sablan MM, Bottjer SW. Topographic organization of a forebrain pathway involved with vocal learning in zebra finches. J Comp Neurol. 1995;358(2):260–78.
Bean BP. The action potential in mammalian central neurons. Nat Rev Neurosci. 2007;8(6):451–65.
Hille B. Ion channels of excitable membranes. Sunderland: Sinauer; 2001.
Zheng J, Trudeau MC. Handbook of ion channels: CRC Press; 2015.
Mozzachiodi R, Byrne JH. More than synaptic plasticity: role of nonsynaptic plasticity in learning and memory. Trends Neurosci. 2010;33(1):17–26.
Zhang W, Linden DJ. The other side of the engram: experience-driven changes in neuronal intrinsic excitability. Nat Rev Neurosci. 2003;4(11):885–900.
Ross MT, Flores D, Bertram R, Johnson F, Wu W, Hyson RL. Experience-dependent intrinsic plasticity during auditory learning. J Neurosci. 2019;39(7):1206–21.
Lovell PV, Carleton JB, Mello CV. Genomics analysis of potassium channel genes in songbirds reveals molecular specializations of brain circuits for the maintenance and production of learned vocalizations. BMC Genomics. 2013;14:470.
Yu FH, Catterall WA. Overview of the voltage-gated sodium channel family. Genome Biol. 2003;4(3):207.
Catterall WA. Voltage-gated calcium channels. Cold Spring Harb Perspect Biol. 2011;3(8):a003947.
Jentsch TJ, Stein V, Weinreich F, Zdebik AA. Molecular structure and physiological function of chloride channels. Physiol Rev. 2002;82(2):503–68.
Yates B, Braschi B, Gray KA, Seal RL, Tweedie S, Bruford EA. Genenames.org: the HGNC and VGNC resources in 2017. Nucleic Acids Res 2017;45(D1):D619-DD25.
Fagerberg L, Hallstrom BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics. 2014;13(2):397–406.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Yue F, Cheng Y, Breschi A, Vierstra J, Wu W, Ryba T, et al. A comparative encyclopedia of DNA elements in the mouse genome. Nature. 2014;515(7527):355–64.
Warren WC, Hillier LW, Tomlinson C, Minx P, Kremitzki M, Graves T, et al. A new chicken genome assembly provides insight into avian genome structure. G3 (Bethesda). 2017;7(1):109–17.
Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, et al. The genome of a songbird. Nature. 2010;464(7289):757–62.
Kent WJ. BLAT--the BLAST-like alignment tool. Genome Res. 2002;12(4):656–64.
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46(D1):D754–D61.
Korlach J, Gedman G, Kingan SB, Chin CS, Howard JT, Audet JN, et al. De novo PacBio long-read and phased avian genome assemblies correct and add to reference genes generated with intermediate and short reads. Gigascience. 2017;6(10):1–16.
Boratyn GM, Camacho C, Cooper PS, Coulouris G, Fong A, Ma N, et al. BLAST: a more efficient report with usability improvements. Nucleic Acids Res. 2013;41(Web Server issue:W29–33.
Friedrich SR. Seq-recovery: initial release (version v1.0.0): Zenodo. 2019. https://doi.org/10.5281/zenodo.2578221
Replogle K, Arnold AP, Ball GF, Band M, Bensch S, Brenowitz EA, et al. The songbird Neurogenomics (SoNG) initiative: community-based tools and strategies for study of brain gene function and evolution. BMC Genomics. 2008;9:131.
Carleton JB, Lovell PV, McHugh A, Marzulla T, Horback KL, Mello CV. An optimized protocol for high-throughput in situ hybridization of zebra finch brain. Cold Spring Harb Protoc. 2014;2014(12):1249–58.
ZEBrA: Zebra finch brain expression atlas. Available from: http://www.zebrafinchatlas.org. Accessed 23 Feb 2019.
Lovell PV, Clayton DF, Replogle KL, Mello CV. Birdsong “transcriptomics”: neurochemical specializations of the oscine song system. PLoS One. 2008;3(10):e3440.
Hilliard AT, Miller JE, Horvath S, White SA. Distinct neurogenomic states in basal ganglia subregions relate differently to singing behavior in songbirds. PLoS Comput Biol. 2012;8(11):e1002773.
Hara E, Rivas MV, Ward JM, Okanoya K, Jarvis ED. Convergent differential regulation of parvalbumin in the brains of vocal learners. PLoS One. 2012;7(1):e29457.
Lovell PV, Huizinga NA, Friedrich SR, Wirthlin M, Mello CV. The constitutive differential transcriptome of a brain circuit for vocal learning. BMC Genomics. 2018;19(1):231.
Lovell PV, Huizinga NA, Getachew A, Mees B, Friedrich SR, Wirthlin M, et al. Curation of microarray oligonucleotides and corresponding ESTs/cDNAs used for gene expression analysis in zebra finches. BMC Res Notes. 2018;11(1):309.
Zhu MX, Evans AM, Ma J, Parrington J, Galione A. Two-pore channels for integrative ca signaling. Commun Integr Biol. 2010;3(1):12–7.
Zakon HH. Adaptive evolution of voltage-gated sodium channels: the first 800 million years. Proc Natl Acad Sci U S A. 2012;109(Suppl 1):10619–25.
Goldin AL. Evolution of voltage-gated Na(+) channels. J Exp Biol. 2002;205(Pt 5:575–84.
Widmark J, Sundstrom G, Ocampo Daza D, Larhammar D. Differential evolution of voltage-gated sodium channels in tetrapods and teleost fishes. Mol Biol Evol. 2011;28(1):859–71.
Calhoun JD, Isom LL. The role of non-pore-forming β subunits in physiology and pathophysiology of voltage-gated sodium channels. Voltage Gated Sodium Channels 221. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. p. 51–89.
Lovell PV, Wirthlin M, Wilhelm L, Minx P, Lazar NH, Carbone L, et al. Conserved syntenic clusters of protein coding genes are missing in birds. Genome Biol. 2014;15(12):565.
Zha XM. Acid-sensing ion channels: trafficking and synaptic function. Mol Brain. 2013;6(1):1-13.
Ren D. Sodium leak channels in neuronal excitability and rhythmic behaviors. Neuron. 2011;72(6):899–911.
Garty H, Palmer LG. Epithelial sodium channels: function, structure, and regulation. Physiol Rev. 1997;77(2):359–96.
Hanukoglu I, Hanukoglu A. Epithelial sodium channel (ENaC) family: phylogeny, structure-function, tissue distribution, and associated inherited diseases. Gene. 2016;579(2):95–132.
Zamponi GW, Striessnig J, Koschak A, Dolphin AC. The physiology, pathology, and pharmacology of voltage-gated calcium channels and their future therapeutic potential. Pharmacol Rev. 2015;67(4):821–70.
Dolphin AC. The alpha2delta subunits of voltage-gated calcium channels. Biochim Biophys Acta. 2013;1828(7):1541–9.
Buraei Z, Yang J. The β subunit of voltage-gated Ca2+ channels. Physiol Rev. 2010;90(4):1461–506.
Chen R-S, Deng T-C, Garcia T, Sellers ZM, Best PM. Calcium channel γ subunits: a functionally diverse protein family. Cell Biochem Biophys. 2007;47(2):178–86.
Chung JJ, Miki K, Kim D, Shim SH, Shi HF, Hwang JY, et al. CatSperzeta regulates the structural continuity of sperm ca(2+) signaling domains and is required for normal fertility. Elife. 2017;6:1-25.
Ren D, Xia J. Calcium signaling through CatSper channels in mammalian fertilization. Physiology (Bethesda). 2010;25(3):165–75.
Cai X, Clapham DE. Evolutionary genomics reveals lineage-specific gene loss and rapid evolution of a sperm-specific ion channel complex: CatSpers and CatSperbeta. PLoS One. 2008;3(10):e3569.
Berridge MJ. The inositol trisphosphate/calcium signaling pathway in health and disease. Physiol Rev. 2016;96(4):1261–96.
Berridge MJ. Inositol trisphosphate and calcium signalling. Nature. 1993;361(6410):315–25.
Singh H. Two decades with dimorphic chloride intracellular channels (CLICs). FEBS Lett. 2010;584(10):2112–21.
Kunzelmann K. TMEM16, LRRC8A, bestrophin: chloride channels controlled by ca(2+) and cell volume. Trends Biochem Sci. 2015;40(9):535–43.
Pedemonte N, Galietta LJ. Structure and function of TMEM16 proteins (anoctamins). Physiol Rev. 2014;94(2):419–59.
Loewen ME, Forsyth GW. Structure and function of CLCA proteins. Physiol Rev. 2005;85(3):1061–92.
Voss FK, Ullrich F, Munch J, Lazarow K, Lutter D, Mah N, et al. Identification of LRRC8 heteromers as an essential component of the volume-regulated anion channel VRAC. Science. 2014;344(6184):634–8.
Lutter D, Ullrich F, Lueck JC, Kempa S, Jentsch TJ. Selective transport of neurotransmitters and modulators by distinct volume-regulated LRRC8 anion channels. J Cell Sci. 2017;130(6):1122–33.
Abascal F, Zardoya R. LRRC8 proteins share a common ancestor with pannexins, and may form hexameric channels involved in cell-cell communication. Bioessays. 2012;34(7):551–60.
Gadsby DC, Vergani P, Csanady L. The ABC protein turned chloride channel whose failure causes cystic fibrosis. Nature. 2006;440(7083):477–83.
Hubert MD, Levitan I, Hoffman MM, Zraggen M, Hofreiter ME, Garber SS. Modulation of volume regulated anion current by I(Cln). Biochim Biophys Acta. 2000;1466(1–2):105–14.
Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, et al. Comparative genomics reveals insights into avian genome evolution and adaptation. Science. 2014;346(6215):1311–20.
Nam K, Mugal C, Nabholz B, Schielzeth H, Wolf JB, Backstrom N, et al. Molecular evolution of genes in avian genomes. Genome Biol. 2010;11(6):R68.
Farries MA, Perkel DJ. A telencephalic nucleus essential for song learning contains neurons with physiological characteristics of both striatum and globus pallidus. J Neurosci. 2002;22(9):3776–87.
Goldberg JH, Adler A, Bergman H, Fee MS. Singing-related neural activity distinguishes two putative pallidal cell types in the songbird basal ganglia: comparison to the primate internal and external pallidal segments. J Neurosci. 2010;30(20):7088–98.
Cain SM, Snutch TP. Contributions of T-type calcium channel isoforms to neuronal firing. Channels (Austin). 2010;4(6):475–82.
Vacher H, Mohapatra DP, Trimmer JS. Localization and targeting of voltage-dependent ion channels in mammalian central neurons. Physiol Rev. 2008;88(4):1407–47.
Kalume F, Yu FH, Westenbroek RE, Scheuer T, Catterall WA. Reduced sodium current in Purkinje neurons from Nav1.1 mutant mice: implications for ataxia in severe myoclonic epilepsy in infancy. J Neurosci. 2007;27(41):11065–74.
Ogiwara I, Miyamoto H, Morita N, Atapour N, Mazaki E, Inoue I, et al. Nav1.1 localizes to axons of parvalbumin-positive inhibitory interneurons: a circuit basis for epileptic seizures in mice carrying an Scn1a gene mutation. J Neurosci. 2007;27(22):5903–14.
Yu FH, Mantegazza M, Westenbroek RE, Robbins CA, Kalume F, Burton KA, et al. Reduced sodium current in GABAergic interneurons in a mouse model of severe myoclonic epilepsy in infancy. Nat Neurosci. 2006;9(9):1142–9.
Hu W, Tian C, Li T, Yang M, Hou H, Shu Y. Distinct contributions of Na(v)1.6 and Na(v)1.2 in action potential initiation and backpropagation. Nat Neurosci. 2009;12(8):996–1002.
Boiko T, Rasband MN, Levinson SR, Caldwell JH, Mandel G, Trimmer JS, et al. Compact myelin dictates the differential targeting of two sodium channel isoforms in the same axon. Neuron. 2001;30(1):91–104.
Leao RN, Naves MM, Leao KE, Walmsley B. Altered sodium currents in auditory neurons of congenitally deaf mice. Eur J Neurosci. 2006;24(4):1137–46.
Felts PA, Yokoyama S, Dib-Hajj S, Black JA, Waxman SG. Sodium channel alpha-subunit mRNAs I, II, III, NaG, Na6 and hNE (PN1): different expression patterns in developing rat nervous system. Brain Res Mol Brain Res. 1997;45(1):71–82.
Beckh S, Noda M, Lubbert H, Numa S. Differential regulation of three sodium channel messenger RNAs in the rat central nervous system during development. EMBO J. 1989;8(12):3611–6.
Cheah CS, Westenbroek RE, Roden WH, Kalume F, Oakley JC, Jansen LA, et al. Correlations in timing of sodium channel expression, epilepsy, and sudden death in Dravet syndrome. Channels (Austin). 2013;7(6):468–72.
Kazarinova-Noyes K, Malhotra JD, McEwen DP, Mattei LN, Berglund EO, Ranscht B, et al. Contactin associates with Na+ channels and increases their functional expression. J Neurosci. 2001;21(19):7517–25.
Kim DY, Carey BW, Wang H, Ingano LA, Binshtok AM, Wertz MH, et al. BACE1 regulates voltage-gated sodium channels and neuronal activity. Nat Cell Biol. 2007;9(7):755–64.
Brackenbury WJ, Calhoun JD, Chen C, Miyazaki H, Nukina N, Oyama F, et al. Functional reciprocity between Na+ channel Nav1.6 and beta1 subunits in the coordinated regulation of excitability and neurite outgrowth. Proc Natl Acad Sci U S A. 2010;107(5):2283–8.
Marionneau C, Carrasquillo Y, Norris AJ, Townsend RR, Isom LL, Link AJ, et al. The sodium channel accessory subunit Navbeta1 regulates neuronal excitability through modulation of repolarizing voltage-gated K(+) channels. J Neurosci. 2012;32(17):5716–27.
Nguyen HM, Miyazaki H, Hoshi N, Smith BJ, Nukina N, Goldin AL, et al. Modulation of voltage-gated K+ channels by the sodium channel beta1 subunit. Proc Natl Acad Sci U S A. 2012;109(45):18577–82.
Deschenes I, Tomaselli GF. Modulation of Kv4.3 current by accessory subunits. FEBS Lett. 2002;528(1–3):183–8.
Grieco TM, Malhotra JD, Chen C, Isom LL, Raman IM. Open-channel block by the cytoplasmic tail of sodium channel beta4 as a mechanism for resurgent sodium current. Neuron. 2005;45(2):233–44.
Bant JS, Raman IM. Control of transient, resurgent, and persistent current by open-channel block by Na channel beta4 in cultured cerebellar granule neurons. Proc Natl Acad Sci U S A. 2010;107(27):12357–62.
Kim JH, Kushmerick C, von Gersdorff H. Presynaptic resurgent Na+ currents sculpt the action potential waveform and increase firing reliability at a CNS nerve terminal. J Neurosci. 2010;30(46):15479–90.
Barbosa C, Tan Z-Y, Wang R, Xie W, Strong JA, Patel RR, et al. Navβ4 regulates fast resurgent sodium currents and excitability in sensory neurons. Mol Pain. 2015;11(1):60.
Raman IM, Bean BP. Resurgent sodium current and action potential formation in dissociated cerebellar Purkinje neurons. J Neurosci. 1997;17(12):4517–26.
Adret P, Margoliash D. Metabolic and neural activity in the song system nucleus robustus archistriatalis: effect of age and gender. J Comp Neurol. 2002;454(4):409–23.
Olveczky BP, Otchy TM, Goldberg JH, Aronov D, Fee MS. Changes in the neural control of a complex motor sequence during learning. J Neurophysiol. 2011;106(1):386–97.
Aronov D, Andalman AS, Fee MS. A specialized forebrain circuit for vocal babbling in the juvenile songbird. Science. 2008;320(5876):630–4.
Hahnloser RHR, Kozhevnikov AA, Fee MS. An ultra-sparse code underlies the generation of neural sequences in a songbird. Nature. 2002;419(6902):65–70.
Hessler NA, Doupe AJ. Singing-related neural activity in a dorsal forebrain-basal ganglia circuit of adult zebra finches. J Neurosci. 1999;19(23):10461–81.
Morgan K, Stevens EB, Shah B, Cox PJ, Dixon AK, Lee K, et al. Beta 3: an additional auxiliary subunit of the voltage-sensitive sodium channel that modulates channel gating with distinct kinetics. Proc Natl Acad Sci U S A. 2000;97(5):2308–13.
Cummins TR, Aglieco F, Renganathan M, Herzog RI, Dib-Hajj SD, Waxman SG. Nav1.3 sodium channels: rapid repriming and slow closed-state inactivation display quantitative differences after expression in a mammalian cell line and in spinal sensory neurons. J Neurosci. 2001;21(16):5952–61.
Cusdin FS, Nietlispach D, Maman J, Dale TJ, Powell AJ, Clare JJ, et al. The sodium channel beta3-subunit induces multiphasic gating in NaV1.3 and affects fast inactivation via distinct intracellular regions. J Biol Chem. 2010;285(43):33404–12.
Merrick EC, Kalmar CL, Snyder SL, Cusdin FS, Yu EJ, Sando JJ, et al. The importance of serine 161 in the sodium channel beta3 subunit for modulation of Na(V)1.2 gating. Pflugers Arch. 2010;460(4):743–53.
Ko SH, Lenkowski PW, Lee HC, Mounsey JP, Patel MK. Modulation of Na(v)1.5 by beta1-- and beta3-subunit co-expression in mammalian cells. Pflugers Arch. 2005;449(4):403–12.
Meadows LS, Chen YH, Powell AJ, Clare JJ, Ragsdale DS. Functional modulation of human brain Nav1.3 sodium channels, expressed in mammalian cells, by auxiliary beta 1, beta 2 and beta 3 subunits. Neuroscience. 2002;114(3):745–53.
Shah BS, Stevens EB, Pinnock RD, Dixon AK, Lee K. Developmental expression of the novel voltage-gated sodium channel auxiliary subunit beta3, in rat CNS. J Physiol. 2001;534(Pt 3:763–76.
Lu B, Su Y, Das S, Liu J, Xia J, Ren D. The neuronal channel NALCN contributes resting sodium permeability and is required for normal respiratory rhythm. Cell. 2007;129(2):371–83.
Tsien RW, Lipscombe D, Madison DV, Bley KR, Fox AP. Multiple types of neuronal calcium channels and their selective modulation. Trends Neurosci. 1988;11(10):431–8.
Berkefeld H, Fakler B. Repolarizing responses of BKCa-Cav complexes are distinctly shaped by their Cav subunits. J Neurosci. 2008;28(33):8238–45.
Berkefeld H, Sailer CA, Bildl W, Rohde V, Thumfart JO, Eble S, et al. BKCa-Cav channel complexes mediate rapid and localized Ca2+−activated K+ signaling. Science. 2006;314(5799):615–20.
Xu W, Lipscombe D. Neuronal ca(V)1.3alpha(1) L-type channels activate at relatively hyperpolarized membrane potentials and are incompletely inhibited by dihydropyridines. J Neurosci. 2001;21(16):5944–51.
Cain SM, Snutch TP. T-type calcium channels in burst-firing, network synchrony, and epilepsy. Biochim Biophys Acta. 2013;1828(7):1572–8.
Yu AC, Margoliash D. Temporal hierarchical control of singing in birds. Science. 1996;273(5283):1871–5.
Rudy B, McBain CJ. Kv3 channels: voltage-gated K+ channels designed for high-frequency repetitive firing. Trends Neurosci. 2001;24(9):517–26.
Kubota M, Saito N. Sodium- and calcium-dependent conductances of neurones in the zebra finch hyperstriatum ventrale pars caudale in vitro. J Physiol. 1991;440:131–42.
Chemin J, Monteil A, Perez-Reyes E, Bourinet E, Nargeot J, Lory P. Specific contribution of human T-type calcium channel isotypes (alpha(1G), alpha(1H) and alpha(1I)) to neuronal excitability. J Physiol. 2002;540(Pt 1:3–14.
Kozlov AS, McKenna F, Lee JH, Cribbs LL, Perez-Reyes E, Feltz A, et al. Distinct kinetics of cloned T-type Ca2 + channels lead to differential Ca2 + entry and frequency-dependence during mock action potentials. Eur J Neurosci. 1999;11(12):4149–58.
Wheeler DB, Randall A, Tsien RW. Roles of N-type and Q-type Ca2+ channels in supporting hippocampal synaptic transmission. Science. 1994;264(5155):107–11.
Sheng ZH, Rettig J, Takahashi M, Catterall WA. Identification of a syntaxin-binding site on N-type calcium channels. Neuron. 1994;13(6):1303–13.
Stanley EF. Single calcium channels and acetylcholine release at a presynaptic nerve terminal. Neuron. 1993;11(6):1007–11.
Weber AM, Wong FK, Tufford AR, Schlichter LC, Matveev V, Stanley EF. N-type Ca2+ channels carry the largest current: implications for nanodomains and transmitter release. Nat Neurosci. 2010;13(11):1348–50.
Simms BA, Zamponi GW. Neuronal voltage-gated calcium channels: structure, function, and dysfunction. Neuron. 2014;82(1):24–45.
Bloodgood BL, Sabatini BL. Nonlinear regulation of unitary synaptic signals by CaV2.3 voltage-sensitive calcium channels located in dendritic spines. Neuron. 2007;53(2):249–60.
Wang K, Kelley MH, Wu WW, Adelman JP, Maylie J. Apamin boosting of synaptic potentials in CaV2.3 R-type Ca2+ channel null mice. PLoS One. 2015;10(9):e0139332.
Weiergraber M, Henry M, Krieger A, Kamp M, Radhakrishnan K, Hescheler J, et al. Altered seizure susceptibility in mice lacking the ca(v)2.3 E-type Ca2+ channel. Epilepsia. 2006;47(5):839–50.
Zaman T, Lee K, Park C, Paydar A, Choi JH, Cheong E, et al. Cav2.3 channels are critical for oscillatory burst discharges in the reticular thalamus and absence epilepsy. Neuron. 2011;70(1):95–108.
Jarvis SE, Zamponi GW. Trafficking and regulation of neuronal voltage-gated calcium channels. Curr Opin Cell Biol. 2007;19(4):474–82.
Dolphin AC. Voltage-gated calcium channels and their auxiliary subunits: physiology and pathophysiology and pharmacology. J Physiol. 2016;594(19):5369–90.
Hoppa MB, Lana B, Margas W, Dolphin AC, Ryan TA. alpha2delta expression sets presynaptic calcium channel abundance and release probability. Nature. 2012;486(7401):122–5.
Taylor CP, Garrido R. Immunostaining of rat brain, spinal cord, sensory neurons and skeletal muscle for calcium channel alpha2-delta (alpha2-delta) type 1 protein. Neuroscience. 2008;155(2):510–21.
Tong XJ, Lopez-Soto EJ, Li L, Liu H, Nedelcu D, Lipscombe D, et al. Retrograde synaptic inhibition is mediated by alpha-Neurexin binding to the alpha2delta subunits of N-type calcium channels. Neuron. 2017;95(2):326–40 e5.
Jackson AC, Nicoll RA. The expanding social network of ionotropic glutamate receptors: TARPs and other transmembrane auxiliary subunits. Neuron. 2011;70(2):178–99.
Kato AS, Gill MB, Yu H, Nisenbaum ES, Bredt DS. TARPs differentially decorate AMPA receptors to specify neuropharmacology. Trends Neurosci. 2010;33(5):241–8.
Wada K, Sakaguchi H, Jarvis ED, Hagiwara M. Differential expression of glutamate receptors in avian neural pathways for learned vocalization. J Comp Neurol. 2004;476(1):44–64.
Cho CH, St-Gelais F, Zhang W, Tomita S, Howe JR. Two families of TARP isoforms that have distinct effects on the kinetic properties of AMPA receptors and synaptic currents. Neuron. 2007;55(6):890–904.
Yamazaki M, Fukaya M, Hashimoto K, Yamasaki M, Tsujita M, Itakura M, et al. TARPs gamma-2 and gamma-7 are essential for AMPA receptor expression in the cerebellum. Eur J Neurosci. 2010;31(12):2204–20.
Studniarczyk D, Coombs I, Cull-Candy SG, Farrant M. TARP gamma-7 selectively enhances synaptic expression of calcium-permeable AMPARs. Nat Neurosci. 2013;16(9):1266–74.
Moss FJ, Viard P, Davies A, Bertaso F, Page KM, Graham A, et al. The novel product of a five-exon stargazin-related gene abolishes ca(V)2.2 calcium channel expression. EMBO J. 2002;21(7):1514–23.
Baker KD, Edwards TM, Rickard NS. The role of intracellular calcium stores in synaptic plasticity and memory consolidation. Neurosci Biobehav Rev. 2013;37(7):1211–39.
Mouton J, Ronjat M, Jona I, Villaz M, Feltz A, Maulet Y. Skeletal and cardiac ryanodine receptors bind to the ca(2+)-sensor region of dihydropyridine receptor alpha(1C) subunit. FEBS Lett. 2001;505(3):441–4.
Walton PD, Airey JA, Sutko JL, Beck CF, Mignery GA, Sudhof TC, et al. Ryanodine and inositol trisphosphate receptors coexist in avian cerebellar Purkinje neurons. J Cell Biol. 1991;113(5):1145–57.
Nilius B, Droogmans G. Amazing chloride channels: an overview. Acta Physiol Scand. 2003;177(2):119–47.
Mura CV, Delgado R, Delgado MG, Restrepo D, Bacigalupo J. A CLCA regulatory protein present in the chemosensory cilia of olfactory sensory neurons induces a ca(2+)-activated cl(−) current when transfected into HEK293. BMC Neurosci. 2017;18(1):61.
Gonzalez-Silva C, Vera J, Bono MR, Gonzalez-Billault C, Baxter B, Hansen A, et al. Ca2+−activated cl- channels of the ClCa family express in the cilia of a subset of rat olfactory sensory neurons. PLoS One. 2013;8(7):e69295.
Krause ET, Bischof H-J, Engel K, Golüke S, Maraci Ö, Mayer U, et al. Chapter two - olfaction in the Zebra finch (Taeniopygia guttata): what is known and further perspectives. In: Naguib M, Barrett L, Healy SD, Podos J, Simmons LW, Zuk M, editors. Advances in the study of behavior. 50: Academic Press; 2018. p. 37–85.
Amo L, Galván I, Tomás G, Sanz JJ. Predator odour recognition and avoidance in a songbird. Funct Ecol. 2008;22(2):289–93.
Petit C, Hossaert-McKey M, Perret P, Blondel J, Lambrechts MM. Blue tits use selected plants and olfaction to maintain an aromatic environment for nestlings. Ecol Lett. 2002;5(4):585–9.
Khan I, Yang Z, Maldonado E, Li C, Zhang G, Gilbert MT, et al. Olfactory receptor subgenomes linked with broad ecological adaptations in Sauropsida. Mol Biol Evol. 2015;32(11):2832–43.
Steiger SS, Fidler AE, Valcu M, Kempenaers B. Avian olfactory receptor gene repertoires: evidence for a well-developed sense of smell in birds? Proc Biol Sci. 2008;275(1649):2309–17.
Steiger SS, Kuryshev VY, Stensmyr MC, Kempenaers B, Mueller JC. A comparison of reptilian and avian olfactory receptor gene repertoires: species-specific expansion of group gamma genes in birds. BMC Genomics. 2009;10:446.
Stephan AB, Shum EY, Hirsh S, Cygnar KD, Reisert J, Zhao H. ANO2 is the cilial calcium-activated chloride channel that may mediate olfactory amplification. Proc Natl Acad Sci U S A. 2009;106(28):11776–81.
Dibattista M, Pifferi S, Boccaccio A, Menini A, Reisert J. The long tale of the calcium activated cl(−) channels in olfactory transduction. Channels (Austin). 2017;11(5):399–414.
Pifferi S, Cenedese V, Menini A. Anoctamin 2/TMEM16B: a calcium-activated chloride channel in olfactory transduction. Exp Physiol. 2012;97(2):193–9.
Pifferi S, Pascarella G, Boccaccio A, Mazzatenta A, Gustincich S, Menini A, et al. Bestrophin-2 is a candidate calcium-activated chloride channel involved in olfactory transduction. Proc Natl Acad Sci U S A. 2006;103(34):12929–34.
Pifferi S, Dibattista M, Sagheddu C, Boccaccio A, Al Qteishat A, Ghirardi F, et al. Calcium-activated chloride currents in olfactory sensory neurons from mice lacking bestrophin-2. J Physiol. 2009;587(Pt 17:4265–79.
Ross MT, Flores D, Bertram R, Johnson F, Hyson RL. Neuronal intrinsic physiology changes during development of a learned behavior. eNeuro. 2017;4(5):1-16.
Olson CR, Hodges LK, Mello CV. Dynamic gene expression in the song system of zebra finches during the song learning period. Dev Neurobiol. 2015;75(12):1315–38.
Mori C, Wada K. Audition-independent vocal crystallization associated with intrinsic developmental gene expression dynamics. J Neurosci. 2015;35(3):878–89.
We wish to thank Henrique Von Gersdorff and Benjamin Zemel for proofreading the manuscript and providing feedback. We also wish to thank Alexa Buckner for preparing and uploading high-resolution scans of in situ data to ZEBrA, and Alison Thomas for her help with gene model evaluation. Finally, we would like to thank Wes Warren and Erich Jarvis for providing us the chicken and zebra finch PacBio genomes before they became available on NCBI.
This research was supported by the National Institutes of Health (R24_GM120464 and R21_DC014432) and by the National Science Foundation (NSF-143602 and Graduate Research Fellowship DGE-1448072). These funding bodies played no role in the design of the study, collection, analysis, or interpretation of data, nor in writing or editing the manuscript.
All procedures involving animals conformed to NIH guidelines for the use and care of animals in research and were approved by OHSU’s Institutional Animal Care Use Committee (IACUC).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Documents synteny analysis, chicken orthologs, model completeness, sequence recovery data, and EST evidence for ion channel genes in the zebra finch genome. Table S2. details evidence supporting gene loss conclusions. (XLSX 29 kb)
Overview of ortholog identification pipeline. Each box represents a step toward identifying an ortholog in zebra finch. Arrows connecting boxes indicate the most common workflows. Details of main pipeline (yellow) and all variations (Cases 1–6), including the use of zebra finch PacBio (pink) and other species (yellow), can be found in methods section. Dotted line indicates Case 4, where there is limited synteny information in taeGut1 and PacBio is required for verification (PDF 607 kb)
Confirming orthology in the PacBio assembly using BLAST. Graphic summary of BLAST results demonstrating that SCN1B is present on gap-less PacBio scaffold MUGN01000920.1 with conserved exon structure and the conserved syntenic genes BCL3 and a fragment of RMB42 (see Fig. 2). Numbers indicate location (bases) along the scaffold. All alignments shown are non-avian RefSeqs. Double dashes indicate additional contracted sequence between selected regions shown. (PDF 234 kb)
Quantitative analysis of zebra finch model completeness. A) Frequency histogram of genes by zebra finch/human model length ratio. B) Frequency histogram of genes by percent of post-recovery length. (PDF 183 kb)
Visualization of sequence recovery for a gene with split models. The zebra finch gene, CACNA2D2, is displayed using the UCSC Genome Browser. A zoomed-in region from within the green rectangle is shown in the bottom panel to highlight the detailed structure of alignments. This gene has three partial zebra finch Ensembl models (dark red track) in a region with numerous gaps (black Gap tack). Note the dips in sequence quality scores (light blue track) surrounding the gaps. Alignment of a more complete chicken model (magenta arrowheads) reveals additional sequence blocks that are missing from the zebra finch model, displayed in the SeqRecovery BED track highlighted in yellow. Non-zebra finch RefSeqs (dark blue tracks) provide further support for blocks of additional sequence recovered through alignments of the chicken model. (PDF 336 kb)
Sequence recovery BED track. A taeGut1 BED track of the sequence recovered from aligning chicken models. (TXT 19 kb)
Expression of select ion channel genes in song nuclei. Representative in situ hybridization photomicrographs of select ion channel genes in song nuclei RA, HVC, Area X, and DLM. All genes are differentially expressed except for CACNA2D1 and CACNG5 in RA, which appear to label sparse populations of cells. Area X panels show select genes with enhanced expression in sparse populations of cells. Camera lucida drawings indicating the location of these nuclei can be found in Fig. 7a for RA, Fig. 8a for HVC, Fig. 9a for Area X, and Fig. 10a for DLM. Gene abbreviations are given in Table 1. All scale bars = 500 μm. (PDF 5215 kb)
About this article
Cite this article
Friedrich, S.R., Lovell, P.V., Kaser, T.M. et al. Exploring the molecular basis of neuronal excitability in a vocal learner. BMC Genomics 20, 629 (2019). https://doi.org/10.1186/s12864-019-5871-2
- Vocal learning
- Zebra finch
- Neuronal excitability
- Ion channels
- Gene expression
- Action potential
- Comparative genomics