Skip to main content

Song exposure regulates known and novel microRNAs in the zebra finch auditory forebrain

Abstract

Background

In an important model for neuroscience, songbirds learn to discriminate songs they hear during tape-recorded playbacks, as demonstrated by song-specific habituation of both behavioral and neurogenomic responses in the auditory forebrain. We hypothesized that microRNAs (miRNAs or miRs) may participate in the changing pattern of gene expression induced by song exposure. To test this, we used massively parallel Illumina sequencing to analyse small RNAs from auditory forebrain of adult zebra finches exposed to tape-recorded birdsong or silence.

Results

In the auditory forebrain, we identified 121 known miRNAs conserved in other vertebrates. We also identified 34 novel miRNAs that do not align to human or chicken genomes. Five conserved miRNAs showed significant and consistent changes in copy number after song exposure across three biological replications of the song-silence comparison, with two increasing (tgu-miR-25, tgu-miR-192) and three decreasing (tgu-miR-92, tgu-miR-124, tgu-miR-129-5p). We also detected a locus on the Z sex chromosome that produces three different novel miRNAs, with supporting evidence from Northern blot and TaqMan qPCR assays for differential expression in males and females and in response to song playbacks. One of these, tgu-miR-2954-3p, is predicted (by TargetScan) to regulate eight song-responsive mRNAs that all have functions in cellular proliferation and neuronal differentiation.

Conclusions

The experience of hearing another bird singing alters the profile of miRNAs in the auditory forebrain of zebra finches. The response involves both known conserved miRNAs and novel miRNAs described so far only in the zebra finch, including a novel sex-linked, song-responsive miRNA. These results indicate that miRNAs are likely to contribute to the unique behavioural biology of learned song communication in songbirds.

Background

Songbirds are important models for exploring the neural and genomic mechanisms underlying vocal communication, social experience and learning (reviewed in [1]). Songbirds communicate using both innate calls and learned vocalizations (songs), and unique specializations of the brain evolved to support this behavior (reviewed in [2]). In the zebra finch, only the male produces songs, although both sexes process and discriminate specific songs [3–6]. The genome is actively engaged by song communication, as first shown in an early demonstration of how gene responses in the brain discriminate among different song stimuli [7]. The genomic response is not a simple correlate of neural activity and it can vary significantly according to the salience and behavioral context of the experience [8–13]. Recent studies using microarray technology have now shown that song exposure affects the expression of thousands of genes in the auditory forebrain [14, 15]. Repeated exposure to one song leads to an altered gene expression profile, correlated with habituation of both the behavioral and immediate genomic responses to that specific song. These observations suggest the involvement of large and dynamic transcriptional network in the recognition and memory of complex vocal signals [14].

MicroRNAs (miRNAs or miRs) are emerging as potential control points in transcriptional networks, and may be particularly important for the evolution of brain and behavior. Many miRNAs are expressed in the brain [16], often in different patterns in different species [17–19]. Brain miRNAs undergo dramatic changes in expression during development [20–22] and aging [23] and have been functionally implicated in neurological disease [24]. They may also function in the normal physiological operation of the nervous system as suggested by evidence for involvement of miR-132 and miR-219 in circadian clock regulation [25] and miR-134 in control of dendritic translation [26, 27].

Here we apply massively parallel Illumina sequencing to probe the involvement of miRNAs in the processing of song experience in the zebra finch auditory forebrain. We begin by identifying 155 different miRNA sequences and the genomic loci of their precursor sequences in the zebra finch genome, including 34 miRNA genes that have not been detected in the genomes of other species. We then ask whether the miRNA content changes after song exposure and find robust evidence of miRNA responses to song playbacks. We also assess correlations between expression changes of a novel miRNA and its predicted target mRNAs during song habituation. The results indicate an active role for miRNAs in the neural processing of a natural perceptual experience - hearing the sound of another bird singing.

Results

The miRNAs of the zebra finch auditory forebrain

We carried out Illumina small RNA sequencing (RNA-seq) on the small RNA (~18-30 nucleotides) fraction of total RNA isolated from adult zebra finch auditory forebrain. Ultimately, we performed 6 Illumina runs on 6 different RNA samples, to assess the effects of song exposure (next section). First we describe the overall small RNA profile obtained by combining the results of all the runs, representing 36 adult zebra finches (equal numbers of males and females). A total of 20 million reads were obtained (Table 1) and aligned to reference miRNA sequences from other species (miRBase version 13.0). Overall we identified 107 non-redundant miRNAs representing 52% of sequences that have been previously identified in chicken, rodent and human. The remaining sequences mapping to the piRNA database were denoted as piRNA reads (~30%) (Additional File 1, Table S1).

Table 1 Summary statistics for the read alignments

Reads that did not align to known RNAs were assessed for miRNA potential through a novel miRNA discovery pipeline described by Creighton et al.[28] which tests for properties that are characteristic of known miRNAs. These properties include the following: 1) The mature sequence must map to the stem region of the hairpin sequence of the putative precursor extracted from the zebra finch genome. 2) The mature miRNA sequence must map to the precursor such that it can be processed following the Drosha processing rules [29]. All novel miRNA candidates that map to the loop region and/or lack appropriate Drosha processing sites are failed. 3) Known miRNAs have stable 5'-ends that vary at the most by +/- 1 nucleotide. 4) By contrast the 3'-ends of miRNAs are highly heterogeneous in length due to imprecise Dicer processing [29, 30] and exhibit non-templated nucleotide sequence changes due to RNA editing [29–31]. 5) Consequently, the putative precursor must give a strong signal of sequence alignments in a tight area of 18-25 nucleotides. Small RNA sequences that are distributed fairly evenly along the entire length of the precursor are rejected since they likely represent degraded products of a large RNA. The candidates that also demonstrate the presence of the miRNA star sequence (miR*) mapping on the opposite side of the mature miRNA and occurring at a lower abundance in the deep sequencing data are considered to be confirmed novel miRNAs in zebra finch. Using this pipeline (Figure 1) we discovered 48 putative novel miRNAs that map on the zebra finch genome to a stem loop structure that folds with a minimum free energy of < -20 kcal/mol [32]. The complete analysis and mapping information for all the novel miRNA candidates is given in Additional File 1, Tables S2 and S3.

Figure 1
figure 1

Pipeline with yields for analysis of putative novel miRNAs. 52 small RNA sequences did not align to miRBase reference sequences and were assessed for miRNA potential. 48 sequences passed the minimum criteria and were categorized into three groups according to strength of evidence (sequences are color-coded in Additional File 1, Table S3, as indicated). Seven (7) are confirmed novel miRNAs since they had all the characteristics of known miRNAs and in addition also had a less abundant miR* sequence that maps on the opposite side of the stem from the putative novel miRNA. These are labelled green in Additional File 1, Table S3. Twenty-one (21) putative novel miRNAs are highly confident (labelled blue) since they also shared characteristics of known miRNAs but no sequence was found aligning to the miR* region. Given that the miR and miR* sequences for most known miRNAs have vastly different copy numbers such that the miR* sequence is sometimes not found, the highly confident candidates are also highly likely to be genuine novel miRNAs, Twenty (20) candidates (labelled grey) had a subset of the characteristics of known miRNAs but not all and therefore were deemed potential candidates that require more evidence.

All novel miRNA candidates were mapped to genomic loci in the zebra finch genome assembly [33], and also to human and chicken genomes using the BLAT function of the UCSC Genome Browser (Additional File 1, Table S3). In the zebra finch genome, the loci include both annotated exons and introns as well as unannotated intergenic regions. Thirty-four (34) novel microRNAs uncovered from zebra finch are not presently detected in the human or chicken genome assemblies. Eleven (11) map to genome positions in chicken, and six to positions in the human (with three of these found in human but not chicken assemblies). Tgu-mir-2976 maps to three loci in the finch and 14 in the chicken, indicating a probable expansion of this miRNA in the chicken lineage. This putative novel miRNA is not currently detected in the human assembly HG18. Tgu-mir-2985 is intriguing as it is located within two stem loops within the introns of two functionally related genes: the glutamate receptor subunits GRIA2 and GRIA4 in all three genomes.

miRNA responses to song exposure

When zebra finches are exposed to playback of a song they have not heard recently, changes occur in the expression of many different mRNAs as detected 30 min after stimulus onset [14]. To determine whether specific miRNAs also change in expression, we counted the Illumina reads in samples of RNA pooled from the auditory forebrain of birds either 30 min after onset of song playback (Song group) or from matched controls (Silence group). In our first such experiment, the birds in both groups were all males (n = 6 each). The read count for each miRNA in each sample was normalized to the total number of usable reads mapped in that sample. We then calculated the ratio of the normalized count in the Song-stimulated condition compared to the Silence condition and performed a Fisher's exact test (with correction for multiple testing) to evaluate whether the ratio differed significantly from the range of expected values at a 95% confidence interval. In the initial experiment with males, 49 of the known conserved miRNAs showed a significant difference, with 28 decreasing and 21 increasing in the group exposed to song (Additional File 1, Table S4).

To address the biological reproducibility of the miRNA responses to song more broadly, we then repeated the small RNA-seq comparison two additional times using new groups of birds. In the second experiment, we used only females, and in the third we used an equal mix of males and females. In total, therefore, we performed three independent "song-silence" pairwise comparisons by small RNA-seq, with an overall sex balance but different sex ratios in each individual comparison. These second and third experiments were done six months after the first and Illumina technology had improved by this time so that we obtained twice as many read counts (Table 1) - but again we normalized to the total mapped read number in each individual sample for our statistic analyses. As in the first experiment, we again observed differential read counts for roughly a third of the miRNAs, but the identities of the miRNAs affected were somewhat different in each comparison. This is summarized graphically as a Venn diagram (Additional File 2, Figure S1), and comprehensive read count data are presented in Additional File 1, Table S4. Across all three experiments, five conserved miRNAs showed changes that were both significant and in same direction in all comparisons (Table 2). For a number of other miRNAs, including let-7f, an apparent effect of song exposure was measured in all three experiments but the direction of change was not consistent (Additional File 1, Table S4).

Table 2 Conserved miRNAs with consistent responses to song exposure

We performed TaqMan assays on RNA from additional birds, probing for eleven of the "significantly affected" miRNAs, and obtained fluorescent signals in PCR for ten. In nine out of ten cases, we observed the same direction of song response by TaqMan as in the small RNA-seq experiment, although the P-value by TaqMan was below 0.05 in only five cases (tgu-miR-124, tgu-miR-29a, tgu-miR-92, tgu-129-5p, and tgu-miR-2954-3p, Additional File 1, Table S4). The lack of statistical significance in the TaqMan assay for the others could reflect differences in the sensitivity and resolution of Illumina vs. TaqMan assays, or the operation of other uncontrolled factors in our experiments that lead to variability in the expression of some miRNAs.

The transcriptional response in the auditory forebrain of zenk and other mRNAs is specific to song relative to non-song auditory stimuli [6, 7, 34, 35]. To test for song-specificity of the miRNA response, we conducted a further TaqMan experiment assessing the levels of six miRNAs (tgu-miR-124, tgu-miR-92, tgu-miR-129-5p, and three miRNAs derived from the tgu-miR-2954 locus, next section), in birds who had heard either a normal song or a carefully matched non-song acoustic stimulus, "song enveloped noise" (SEN). SEN has the same amplitude envelope as the song from which it is derived but spectral content has been randomized so it does not sound like a song [34]. By TaqMan PCR, we confirmed that normal song induced a larger increase in zenk mRNA in these birds than did SEN (Additional File 2, Figure S3 panel D). In these same animals, normal song, but not SEN, triggered a significant decrease in the levels of tgu-miR-124, tgu-mir-129-5p, tgu-miR-92 and tgu-miR-2954-3p (Additional File 2, Figure S3 panels A-C, H). Thus we conclude that there is indeed a unique miRNA response in the auditory forebrain that is selective for song over non-song acoustic stimuli.

A complex sex-linked miRNA locus in zebra finch and other birds

The novel miRNA, tgu-mir-2954, that was detected most frequently in our Illumina assays maps to the sense strand of an intron in the XPA gene, on the Z chromosome (Figure 2A). The precursor hairpin contains reads from both arms, thus meeting our bioinformatic criteria for a confirmed miRNA (Figure 2B). By contrast to most known miRNAs, the numbers of reads from both 5' and 3' arms were found at similar copy numbers, suggesting that both arms may make functional mature miRNAs. BLAST analysis of the mir-2954 hairpin precursor sequence against the NCBI nr database identified a putative mature miRNA in chicken (gi|145279910|emb|AM691163.1|), and BLAT analysis of a collection of transcripts from crocodile and 11 other bird species [36] detected mir-2954 transcripts in 2 non-passerine species (two hummingbirds) and 3 passerine species (the American crow, the pied flycatcher, and the golden collared manakin) (Additional File 2, Figure S2). There was no BLAT hit in the crocodile, the remaining 3 non-passerine birds (Emu, budgerigar, and ringneck dove), and 3 passerine species (collared flycatcher, blue tit and Eastern phoebe). The lack of a hit does not necessarily mean absence of the gene as these datasets represent incomplete transcriptomes derived by 454 sequencing [36]. These results clarify that the sequence is not unique to the zebra finch or passerines, but may nevertheless have a restricted distribution within birds.

Figure 2
figure 2

The genome locus for tgu-mir-2954 produces three different miRNAs. A. Alignments via the UCSC Genome Browser of the three detected miRNAs to the intron of the zebra finch XPA gene. B. Hairpin precursors for the three miRNAs. C. Northern blot analysis using an RNA probe complementary to novel miRNA tgu-miR-2954-3p.

To validate the existence of these two miRNAs in zebra finch, we performed TaqMan analyses for both, using their reverse complements as controls. Interestingly, we got significant expression values not only for the predicted miRNAs but also for one of the reverse-complement miRNAs (tgu-miR-2954R-5p) although no significant song regulation for miR-2954R-5p was found (Additional File 2, Figure S3 panels I-J). With respect to the XPA gene within which this locus is embedded (Figure 2A), these data suggest that precursor-miRNA-stem loops are produced from both the sense (same orientation as XPA) and antisense strands. The stem loop precursor processed by Drosha from the sense RNA (tgu-mir-2954) generates two active miRNAs from its both arms (tgu-miR-2954-3p and tgu-miR-2954-5p). The stem loop precursor processed by Drosha from the antisense RNA (tgu-mir-2954R) generates at least one active miRNA (tgu-miR-2954R-5p) from its 5' end sequence.

We carried out Northern analysis on tgu-miR-2954-3p, which is the miRNA that has the highest number of read counts detected in our Illumina assays among the three miRNAs from the tgu-mir-2954 locus. A robust signal at ~22 nucleotides is evident in mixed-sex pools of RNA from birds hearing either song or silence, and a weaker signal is also detectable in two female-only pools of RNA (Figure 2C). Greater expression in males is consistent with the ZZ genotype of males and the lack of efficient sex chromosome dosage compensation in the zebra finch [37, 38].

By TaqMan as well as by Illumina, we observed an apparent sex difference in the direction of the response of tgu-miR-2954-3p to song - up in males and down in females (Figure 3 and Additional File 1, Table S4). This suggests this locus may be under complex regulation, integrating information about sex, auditory or social experience and perhaps also other factors related to XPA gene expression.

Figure 3
figure 3

Analysis of miRNAs produced at the tgu-mir-2954 locus. TaqMan and Illumina RNA-seq data generated from independent sets of birds (n = 6 in each data set) for expression from the tgu-mir-2954 locus. A) TaqMan results, where the relative gene expression of each individual bird (open circle) was obtained by using the 2^-ddCt method [98]; the relative gene expression of either Silence (white bar) or Song (gray bar) group was the mean of six individuals; the P value was calculated by paired t test since each song stimulated animal was explicitly paired with a silence control animal collected simultaneously. B) Read counts from the Illumina RNA-seq for miR-2954-3p and miR-2954-5p (also shown in the Additional File 1, Table S4).

To gain insight into the potential functional role of tgu-miR-2954-3p in the response to song, we used a conservative strategy to predict gene targets that are both conserved in birds and responsive to song exposure in the zebra finch. Potential targets of miRNAs are described as mRNAs that have sequences that can undergo Watson-Crick base pairing with the 5'-seed (nucleotide 2-7) of the miRNA [39]. For target prediction we applied the TargetScan (5.1) algorithm using the chicken genome as an initial reference, and then confirmed presence of the target sequence in the zebra finch. For evidence of song responsiveness, we used the data set of Dong et al. [14]. Eight genes met all these criteria (Table 3) and are thus both song-responsive and also subject to regulation by tgu-miR-2954-3p. These genes all have functions in control of cell proliferation or neurite outgrowth (see below).

Table 3 Song-regulated targets of tgu-miR-2954-3p

Discussion

Here we show that a natural perceptual experience, hearing the sound of another bird singing, alters the profile of miRNAs in parts of the songbird brain responsible for auditory perception, integration and memory. The song-regulated population includes both known (conserved) and novel miRNAs. We highlight one sex-linked song-responsive miRNA and identify mRNAs that are potential targets of its action during song exposure. Thus miRNAs may have roles in the information processing functions of the brain, in addition to their roles in brain development and evolution.

To demonstrate this, we first catalogued the miRNAs expressed in the adult zebra finch auditory forebrain. We used massively parallel Illumina sequencing of small RNAs to perform this cataloguing efficiently. In addition to known conserved miRNAs, our analysis identified 48 small RNA sequences that meet the structural criteria for miRNAs but had not been described in miRBase in any organism at the time of our analysis. Fourteen of these are detected in the chicken or human genome assemblies and may give rise to miRNAs that have not yet been described elsewhere due to low copy number, restricted tissue distribution or other factors. The remaining novel miRNAs, 34 in number, may be unique to the zebra finch or the songbird lineage. Few studies have attempted de novo identification of miRNAs from the brain [18] and ours is the first to report direct sequencing of songbird brain miRNAs. A previous study did identify precursor sequences for five conserved miRNAs in the developing zebra finch brain [40]. Also, in parallel with our own Illumina analysis, Li and her colleagues used 454 sequencing to identify miRNAs in the brain and liver of adult zebra finches. These different sets of annotations are compared and collated in a supplement to the analysis of the zebra finch genome assembly [33].

By comparing birds hearing novel song playbacks or silence, we found evidence for experience-dependent fluctuations in large numbers of miRNAs in the auditory forebrain. We performed three separate pairwise comparisons by Illumina, where all aspects of the experimental conditions were carefully counterbalanced between the two groups in each comparison. The three comparisons were not direct replications of each other, as each had a different sex ratio. Our reasons for varying the sex ratio were partly pragmatic (limited numbers of birds of the same sex that could be removed from our aviary) and partly analytical (males and females have different behavioral responses to songs). Some of the differences between the three sets of results may reflect real biological differences in the responses of males and females. Indeed, our Northern analysis of the tgu-miR-2954-3p confirms a sex difference in expression of this Z-linked miRNA gene. This is especially intriguing because we also obtained TaqMan evidence for both sense and antisense transcripts of this miRNA. One can imagine scenarios where different ratios of sense and antisense transcription occur in males (two copies of the gene) and females (one copy of the gene) with different consequences on the transcriptional networks affected by song exposure in the two sexes.

Ignoring the potential effects of sex, we identified five miRNAs that showed significant and consistent changes in response to song across all three Illumina comparisons. Three miRNAs consistently decreased after song (tgu-miR-92, tgu-miR-124, tgu-miR-129-5p) and two increased (tgu-miR-25, tgu-miR-192). The down-regulated miRNAs are at much higher abundance (> 1000 reads in each run) and perhaps for this reason we were more successful at detecting them and replicating their song regulation by TaqMan assay in subsequent experiments with additional groups of birds. The most abundant miRNA in our regulated set, tgu-miR-124, consistently met the statistical test for significant down-regulation by song, in each of six separate experiments (three Illumina comparisons, two TaqMan analyses in Additional File 1, Table S4, and the TaqMan comparison of song vs. SEN in Additional File 2, Figure 3).

In studies in other species, miR-124 has been linked to brain plasticity and development in several contexts. Chronic cocaine administration results in down-regulation of miR-124 in the rodent mesolimbic dopaminergic system [41]. In the developing chick neural tube, miR-124a is a component of a regulatory network that controls the transition between neural progenitors and post-mitotic neurons [42]. miR-124 also regulates adult neurogenesis, and its overexpression promotes neuronal differentiation [42, 43] and neurite outgrowth [44]. Intriguingly, in songbirds neurogenesis continues in the forebrain throughout adulthood, from a population of precursor cells that line the walls of the lateral ventricles and have the characteristics of neural stem cells [45–47]. The net rate of neuronal addition and loss in the adult songbird has been shown to depend on social and environmental influences [48–51]. Perhaps tgu-miR-124 is a regulatory link between experience and neurogenesis - further study of this fascinating possibility is clearly warranted.

Although miRNAs can have diverse functions, they often act by altering the concentrations of specific mRNAs they target via complementary base pairing. We used the TargetScan algorithm [52] to predict binding sites of tgu-miR-2954-3p in chicken genes, and then we confirmed the presence of the same conserved target sequence in the zebra finch genome assembly. We found eight targets that met these criteria and were also regulated by song in the Dong et al. microarray data [14]. These eight genes have a provocative coherence in their function, as they are all implicated in control of cell proliferation and neuronal differentiation. Six operate by affecting gene expression and chromatin remodeling as we briefly review here. ELAVL2 is a member of a protein family that binds AU-rich regions in the 3'UTR of genes such as c-fos and promotes the shift from cell proliferation into cellular differentiation [53–57]. TLK2 is a kinase tightly associated with DNA replication during cell division [58]. At least one of its targets, the histone chaperone Asf1, controls chromatin assembly, thus TLK2 activity can regulate transcription and elongation [59–61]. BTG1 is also regulated during the cell cycle [62]. It acts as a cofactor for Hoxb9, a transcription factor that controls cell proliferation and differentiation, and BTG1 reduces rates of cell proliferation [62–64]. CHD2 can potentially affect transcription of many genes by remodeling chromatin [65, 66]; disruption of CHD2 has profound consequences for development and is implicated in many human diseases [67–69]. HMGB1 is another DNA binding protein that facilitates transcription by altering chromatin structure to ease promoter binding [70–73]. Some of the genes regulated by HMGB1 may play a role in cell proliferation and migration [74, 75]. Neuronal migration and neurite outgrowth are affected by CRKL, a transcriptional activator that is a component of the reelin pathway [76–79]. Unlike the other six genes, NEGR1 and LINGO2 do not seem to alter transcription but they do have established roles in neuronal differentiation. NEGR1 affects cell-cell adhesion to modulate neurite outgrowth and synapse formation [80–82]. LINGO2 is one member of a family of transmembrane proteins that are involved in neural and axonal regeneration [83, 84]. The function of LINGO2 is untested, but expression of a related protein, LINGO1, is attenuated in cortical areas deprived of sensory input and is a partner in a signaling pathway that correlates with neuronal activity during a learning paradigm [85, 86].

Conclusions

In conclusion, these data reveal a network of miRNAs in the zebra finch's auditory forebrain, responsive to the experience of hearing another bird sing. The network includes well-characterized conserved miRNA known to have roles in neuronal differentiation (miR-124), and novel miRNAs that can target genes that control neuronal differentiation (tgu-miR-2954-3p). Our data suggest this miRNA network may influence the fundamental shift we have observed in the transcriptional and metabolic state of the auditory forebrain during the process of song-specific habituation [14, 87]. Further study of song responses in the zebra finch may reveal general insights into the neurogenomic mechanisms that underlie learning, memory and the ongoing adaptation to experience.

Methods

Song stimulation and brain dissections

Zebra finches were obtained from aviaries maintained at the University of Illinois. All procedures involving animals were conducted with the approval of the University of Illinois Institutional Animal Care and Use Committee. The birds were raised in a standard breeding aviary and were tutored under normal social conditions (i.e., by their parents or other adult birds in the breeding colony). All birds used in this study were adults (older than 90 days after hatching). The song playback procedures and brain dissections were performed exactly as in previous microarray analyses, using the same equipment [14, 88]. Briefly, each bird was put individually into a sound isolation chamber for 18 hours on the first day, and on the second day those in the song group heard 30 minutes of a song not heard previously ("novel song"). Matched controls collected in parallel heard no song playback ("silence"). Birds were sacrificed in song-silence pairs, so that 5 minutes before the end of the song playback to one bird, a bird in the silence group was sacrificed and its auditory forebrain was dissected and frozen in dry ice. Then the auditory forebrain of the song-stimulated bird was dissected and frozen in dry ice. The auditory forebrain dissection (also referred to as auditory lobule) is described in [89] and collects NCM (caudomedial nidopallium), CMM (caudomedial mesopallium) and the enclosed Field L subregions. At the end of the song stimulation procedure, all auditory forebrains were transferred and stored at -80C until RNA isolation. For the comparison of responses after overnight isolation to song versus SEN (Additional File 2, Figure S3), we used two matched stimuli derived from bird "C7" as previously described [34].

RNA Samples

For Illumina analyses: Total RNA was extracted using the mirVana miRNA Isolation Kit (Ambion) from three pairs of pooled auditory forebrain samples. 1) Males (samples S7 and S8): 6 birds per pool, collected in November 2008. 2) Females (samples S1 and S2): 6 birds per pool, collected in May 2009. 3) Mixed (samples S3 and S4), 3 males and 3 females each pool, collected in May 2009. Samples with odd numbers were from birds hearing song, and even number hearing silence.

For Northern analysis: Auditory forebrains of 22 birds (12 females and 10 males) were collected in April 2009, and total RNA was extracted by Tri-Reagent (Ambion). Male and female samples were pooled after extraction.

For TaqMan analysis: Analyses were performed on total RNA extracted either by mirVana or Tri-Reagent (Ambion), from the auditory forebrains of individual males or females, collected in April-August 2009, March 2010 or December 2010.

Illumina small RNA sequencing and novel miRNA discovery

Fifteen micrograms of total RNA from auditory forebrain of song bird samples described above were gel-fractionated to isolate 18-30 nt small RNAs. 3' and 5' adapters were ligated to the small RNAs and constructs amplified following RT-PCR following the conditions specified in the small RNA kit (FC-102-1009, Illumina) protocol. The small RNA library was sequenced using a Solexa/Illumina GA-1 Genome analyzer. Small RNA sequences were analyzed through a high-throughput computational pipeline described by [28, 29, 90, 91]. To identify zebra finch miRNAs that are also conserved in chicken, human and mouse, we performed a local Smith-Waterman alignment of each unique sequence read against each of the mature miRNAs in miRBase version 13.0 for each of these species. We allowed for a 3 base overhang on the 5' end and a 6 base overhang on the 3' end. In the case of redundantly aligning reads, mature miRNA sequences were equally apportioned among each of the hairpins. For each sample, all sequence reads were aligned to a reference set of precursor miRNAs from miRBase version 13.0. The reads that did not align to any known miRNA were passed to our novel miRNA discovery platform as previously described [28]. Briefly, each sequence is first mapped to the reference genome sequence (WUGSC 3.2.4) and 200 bases of flanking sequence are extracted to further define the putative hairpin. This extracted sequence is then folded using the Vienna RNA folding package [92] and those sequences that form a plausible hairpin are selected as potential novel miRNA hairpins. These candidates are filtered through a set of three Ambros criteria: 1) the mature putative miRNA sequence must rest on one side of a single hairpin; 2) the putative miRNA sequence must bind relatively tightly within the hairpin stem containing no large or energetically unfavorable loops; and 3) the putative hairpin must have a miRNA-appropriate energy (free energy below -20 kcal/mol). All sequences that passed were then carefully curated to determine if Drosha and Dicer processing could yield the resulting mature sequence from the predicted hairpin. These candidates are then divided into four different categories: "not likely", "potential", "high confidence", and "confirmed" (as in red, gray, blue and green colors in Additional File 1, Tables S2 and S3). Candidates that are flagged red as "not likely" either failed to map in a pile of sequences in a very tight space of 15-25 nt of the predicted hairpin (e.g. were scattered evenly across the full length of the hairpin), mapped within the loop of the hairpin, or mapped to known tRNAs or rRNAs. Candidates that passed all of the above criteria, and also mapped within a hairpin with predicted Drosha and Dicer cut sites were categorized as "high confidence" (blue annotation in Additional File 1, Tables S2 and S3). All high confidence candidates for which we detected both the mature sequence and the putative star sequence from the same hairpin we categorized as "confirmed" (green annotation in Additional File 1, Table S3). In addition to miRNA precursors, the reads were also mapped to the reference zebra finch genome using the Pash software package [93, 94], and uploaded to the Genboree platform (http://www.genboree.com) to identify potential mappings to piRNAs, snoRNAs and other annotations in addition to miRNAs (data shown in Additional File 1, Table S1). PiRNAs (i.e., Piwi-interacting RNAs) have a central role in the maintenance of the integrity of genomes through the silencing of transposable elements [95]. SnoRNAs (small nucleolar RNAs) function in site-specific ribosomal RNA modification, rRNA processing and more recently have been found to guide alternate splicing and RNA editing of mRNA transcripts [96].

TaqMan qPCR

To measure the mature miRNA, the TaqMan MicroRNA Assay Kit (Applied Biosystems) was used according to the manufacturer's instructions. Probe sequences used for each target miRNA are given in Table 4.

Table 4 Probes used for Taqman analysis of specific miRNA sequences

Northern Blot Analysis

Northern blotting to confirm novel miRNA tgu-miR-2954-3p was performed by modifying the protocol of [97]. 2 μg of total RNA was heated at 65°C for 5 min with 2X loading dye (Ambion), quenched on ice, and loaded on a 15% TBE Urea gel (Invitrogen). Total RNA was separated by electrophoresis at 200V for 50 min. The gel was stained with with EtBr in 1x TBE (4 μL of 10 mg/ml EtBr per 100 ml of 1x TBE) for 3 minutes with gentle shaking and transferred to nylon membrane for 90 min at 200V using 1X TBE buffer at room temperature. The membrane was cross-linked at 1200 kJ for 45 seconds. RNA probes were synthesized for tgu-miR-2954-3p probe 5' - UGCUAGGAGUGGAAUGGGGAU G - 3' by Integrated DNA Technologies. Radio labeling was carried out in a reaction of 12.0ul dH2O + 2.0ul PNK buffer + 1.0ul (100ng/ul) probe + 1.0ul PNK polymerase (Promega) + 4.0ul P32-gamma-ATP (10mCi/ml) (PerkinElmer). The reaction was incubated at 37°C for 1 hour and inactivated at 65°C for 10 min. The probe was purified using Nick columns from GE following manufacturer's instructions. The membranes were pre-hybridized for 30 min with 20 ml of pre-hybridization buffer (5X SSC + 20 mM NaPO4 + 7X SDS + 2X Denhardt (pre warmed) at 60° C) in a rotating hybridization oven. Hybridization was carried out at 50°C in a rotating incubator for 24h. The membranes were washed for 10 min at 50°C with 20-30mL of wash buffer (2X SSC + 0.5% SDS). When background was ~0.5 cpm, the membranes were wrapped in saran wrap and exposed at -80°C for ~72h.

References

  1. Clayton DF, Balakrishnan CN, London SE: Integrating genomes, brain and behavior in the study of songbirds. Curr Biol. 2009, 19 (18): R865-873. 10.1016/j.cub.2009.07.006.

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Jarvis ED: Learned birdsong and the neurobiology of human language. Ann N Y Acad Sci. 2004, 1016: 749-777. 10.1196/annals.1298.038.

    PubMed  PubMed Central  Google Scholar 

  3. Miller DB: Acoustic Basis of Mate Recognition by Female Zebra Finches (Taeniopygia-Guttata). Animal Behaviour. 1979, 27 (May): 376-380.

    Google Scholar 

  4. Miller DB: Long-Term Recognition of Fathers Song by Female Zebra Finches. Nature. 1979, 280 (5721): 389-391. 10.1038/280389a0.

    Google Scholar 

  5. Clayton NS: Song Discrimination-Learning in Zebra Finches. Animal Behaviour. 1988, 36: 1016-1024. 10.1016/S0003-3472(88)80061-7.

    Google Scholar 

  6. Stripling R, Kruse AA, Clayton DF: Development of song responses in the zebra finch caudomedial neostriatum: Role of genomic and electrophysiological activities. Journal of Neurobiology. 2001, 48 (3): 163-180. 10.1002/neu.1049.

    CAS  PubMed  Google Scholar 

  7. Mello CV, Vicario DS, Clayton DF: Song presentation induces gene expression in the songbird forebrain. Proc Natl Acad Sci USA. 1992, 89 (15): 6818-6822. 10.1073/pnas.89.15.6818.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Mello C, Nottebohm F, Clayton D: Repeated exposure to one song leads to a rapid and persistent decline in an immediate early gene's response to that song in zebra finch telencephalon. J Neurosci. 1995, 15 (10): 6919-6925.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Jarvis ED, Scharff C, Grossman MR, Ramos JA, Nottebohm F: For whom the bird sings: context-dependent gene expression. Neuron. 1998, 21 (4): 775-788. 10.1016/S0896-6273(00)80594-2.

    CAS  PubMed  Google Scholar 

  10. Clayton DF: The genomic action potential. Neurobiol Learn Mem. 2000, 74 (3): 185-216. 10.1006/nlme.2000.3967.

    CAS  PubMed  Google Scholar 

  11. Kruse AA, Stripling R, Clayton DF: Context-specific habituation of the zenk gene response to song in adult zebra finches. Neurobiol Learn Mem. 2004, 82 (2): 99-108. 10.1016/j.nlm.2004.05.001.

    CAS  PubMed  Google Scholar 

  12. Vignal C, Andru J, Mathevon N: Social context modulates behavioural and brain immediate early gene responses to sound in male songbird. Eur J Neurosci. 2005, 22 (4): 949-955. 10.1111/j.1460-9568.2005.04254.x.

    PubMed  Google Scholar 

  13. Woolley SC, Doupe AJ: Social context-induced song variation affects female behavior and gene expression. PLoS Biol. 2008, 6 (3): e62-10.1371/journal.pbio.0060062.

    PubMed  PubMed Central  Google Scholar 

  14. Dong S, Replogle KL, Hasadsri L, Imai BS, Yau PM, Rodriguez-Zas S, Southey BR, Sweedler JV, Clayton DF: Discrete molecular states in the brain accompany changing responses to a vocal signal. Proc Natl Acad Sci USA. 2009, 106 (27): 11364-11369. 10.1073/pnas.0812998106.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. London SE, Dong S, Replogle K, Clayton DF: Developmental shifts in gene expression in the auditory forebrain during the sensitive period for song learning. Dev Neurobiol. 2009, 69 (7): 437-450. 10.1002/dneu.20719.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. Cao X, Yeo G, Muotri AR, Kuwabara T, Gage FH: Noncoding RNAs in the mammalian central nervous system. Annu Rev Neurosci. 2006, 29: 77-103. 10.1146/annurev.neuro.29.051605.112839.

    CAS  PubMed  Google Scholar 

  17. Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH: Differences in vertebrate microRNA expression. Proc Natl Acad Sci USA. 2006, 103 (39): 14385-14389. 10.1073/pnas.0603529103.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Berezikov E, Thuemmler F, van Laake LW, Kondova I, Bontrop R, Cuppen E, Plasterk RH: Diversity of microRNAs in human and chimpanzee brain. Nat Genet. 2006, 38 (12): 1375-1377. 10.1038/ng1914.

    CAS  PubMed  Google Scholar 

  19. Bak M, Silahtaroglu A, Moller M, Christensen M, Rath MF, Skryabin B, Tommerup N, Kauppinen S: MicroRNA expression in the adult mouse central nervous system. RNA. 2008, 14 (3): 432-444. 10.1261/rna.783108.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Krichevsky AM, King KS, Donahue CP, Khrapko K, Kosik KS: A microRNA array reveals extensive regulation of microRNAs during brain development. RNA. 2003, 9 (10): 1274-1281. 10.1261/rna.5980303.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Miska EA, Alvarez-Saavedra E, Townsend M, Yoshii A, Sestan N, Rakic P, Constantine-Paton M, Horvitz HR: Microarray analysis of microRNA expression in the developing mammalian brain. Genome Biol. 2004, 5 (9): R68-10.1186/gb-2004-5-9-r68.

    PubMed  PubMed Central  Google Scholar 

  22. Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V: Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol. 2004, 5 (3): R13-10.1186/gb-2004-5-3-r13.

    PubMed  PubMed Central  Google Scholar 

  23. Li N, Bates DJ, An J, Terry DA, Wang E: Up-regulation of key microRNAs, and inverse down-regulation of their predicted oxidative phosphorylation target genes, during aging in mouse brain. Neurobiol Aging. 2009

    Google Scholar 

  24. Schratt G: Fine-tuning neural gene expression with microRNAs. Curr Opin Neurobiol. 2009, 19 (2): 213-219. 10.1016/j.conb.2009.05.015.

    CAS  PubMed  Google Scholar 

  25. Cheng HY, Papp JW, Varlamova O, Dziema H, Russell B, Curfman JP, Nakazawa T, Shimizu K, Okamura H, Impey S, et al: microRNA modulation of circadian-clock period and entrainment. Neuron. 2007, 54 (5): 813-829. 10.1016/j.neuron.2007.05.017.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Schratt GM, Tuebing F, Nigh EA, Kane CG, Sabatini ME, Kiebler M, Greenberg ME: A brain-specific microRNA regulates dendritic spine development. Nature. 2006, 439 (7074): 283-289. 10.1038/nature04367.

    CAS  PubMed  Google Scholar 

  27. Fiore R, Khudayberdiev S, Christensen M, Siegel G, Flavell SW, Kim TK, Greenberg ME, Schratt G: Mef2-mediated transcription of the miR379-410 cluster regulates activity-dependent dendritogenesis by fine-tuning Pumilio2 protein levels. EMBO J. 2009, 28 (6): 697-710. 10.1038/emboj.2009.10.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Creighton CJ, Reid JG, Gunaratne PH: Expression profiling of microRNAs by deep sequencing. Brief Bioinform. 2009, 10 (5): 490-497. 10.1093/bib/bbp019.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Reid JG, Nagaraja AK, Lynn FC, Drabek RB, Muzny DM, Shaw CA, Weiss MK, Naghavi AO, Khan M, Zhu H, et al: Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5'-seed/cleavage/anchor regions and stabilize predicted mmu-let-7a:mRNA duplexes. Genome Res. 2008, 18 (10): 1571-1581. 10.1101/gr.078246.108.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, et al: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008, 18 (4): 610-621. 10.1101/gr.7179508.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al: A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007, 129 (7): 1401-1414. 10.1016/j.cell.2007.04.040.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Lee RC, Ambros V: An extensive class of small RNAs in Caenorhabditis elegans. Science. 2001, 294 (5543): 862-864. 10.1126/science.1065329.

    CAS  PubMed  Google Scholar 

  33. Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, Searle S, White S, Vilella AJ, Fairley S, et al: The genome of a songbird. Nature. 2010, 464 (7289): 757-762. 10.1038/nature08819.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Park KH, Clayton DF: Influence of restraint and acute isolation on the selectivity of the adult zebra finch zenk gene response to acoustic stimuli. Behav Brain Res. 2002, 136 (1): 185-191. 10.1016/S0166-4328(02)00129-8.

    PubMed  Google Scholar 

  35. Bailey D, Wade J: Differential expression of the immediate early genes FOS and ZENK following auditory stimulation in the juvenile male and female zebra finch. Brain Res Mol Brain Res. 2003, 116 (1-2): 147-154.

    CAS  PubMed  Google Scholar 

  36. Kánstner A, Wolf JBW, Backström N, Whitney O, Balakrishnan CN, Day L, Edwards SV, Janes DE, Schlinger BA, Wilson RK, et al: Comparative genomics based on massive parallel transcriptome sequencing reveals patterns of substitution and selection across 10 bird species. Molecular Ecology. 2010, 19 (SUPPL. 1): 266-276.

    Google Scholar 

  37. Itoh Y, Melamed E, Yang X, Kampf K, Wang S, Yehya N, Van Nas A, Replogle K, Band MR, Clayton DF, et al: Dosage compensation is less effective in birds than in mammals. J Biol. 2007, 6 (1): 2-10.1186/jbiol53.

    PubMed  PubMed Central  Google Scholar 

  38. Itoh Y, Replogle K, Kim YH, Wade J, Clayton DF, Arnold AP: Sex bias and dosage compensation in the zebra finch versus chicken genomes: General and specialized patterns among birds. Genome Research. 2010, 20 (4): 512-518. 10.1101/gr.102343.109.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Li X, Wang XJ, Tannenhauser J, Podell S, Mukherjee P, Hertel M, Biane J, Masuda S, Nottebohm F, Gaasterland T: Genomic resources for songbird research and their use in characterizing gene expression during brain development. Proc Natl Acad Sci USA. 2007, 104 (16): 6834-6839. 10.1073/pnas.0701619104.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Chandrasekar V, Dreyer JL: microRNAs miR-124, let-7d and miR-181a regulate cocaine-induced plasticity. Mol Cell Neurosci. 2009, 42 (4): 350-362. 10.1016/j.mcn.2009.08.009.

    CAS  PubMed  Google Scholar 

  42. Visvanathan J, Lee S, Lee B, Lee JW, Lee SK: The microRNA miR-124 antagonizes the anti-neural REST/SCP1 pathway during embryonic CNS development. Genes Dev. 2007, 21 (7): 744-749. 10.1101/gad.1519107.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. Cheng LC, Pastrana E, Tavazoie M, Doetsch F: miR-124 regulates adult neurogenesis in the subventricular zone stem cell niche. Nat Neurosci. 2009, 12 (4): 399-408. 10.1038/nn.2294.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Yu JY, Chung KH, Deo M, Thompson RC, Turner DL: MicroRNA miR-124 regulates neurite outgrowth during neuronal differentiation. Exp Cell Res. 2008, 314 (14): 2618-2633. 10.1016/j.yexcr.2008.06.002.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. Goldman SA, Nottebohm F: Neuronal production, migration and differentiation in a vocal control nucleus of the adult female canary brain. Proceedings of the National Academy of Sciences of the United States of America. 1983, 80: 2390-2394. 10.1073/pnas.80.8.2390.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. Alvarez-Buylla A, Theelen M, Nottebohm F: Proliferation "hot spots" in adult avian ventricular zone reveal radial cell division. Neuron. 1990, 5 (1): 101-109. 10.1016/0896-6273(90)90038-H.

    CAS  PubMed  Google Scholar 

  47. Alvarez-Buylla A, Kirn JR: Birth, migration, incorporation, and death of vocal control neurons in adult songbirds. Journal of Neurobiology. 1997, 33: 585-601. 10.1002/(SICI)1097-4695(19971105)33:5<585::AID-NEU7>3.0.CO;2-0.

    CAS  PubMed  Google Scholar 

  48. Barnea A: Interactions between environmental changes and brain plasticity in birds. General and Comparative Endocrinology. 2009, 163 (1-2): 128-134. 10.1016/j.ygcen.2009.03.031.

    CAS  PubMed  Google Scholar 

  49. Kirn JR: The relationship of neurogenesis and growth of brain regions to song learning. Brain and Language. 2010, 115 (1): 29-44. 10.1016/j.bandl.2009.09.006.

    PubMed  Google Scholar 

  50. Nottebohm F, O'Loughlin B, Gould K, Yohay K, Alvarez-Buylla A: The life span of new neurons in a song control nucleus of the adult canary brain depends on time of year when these cells are born. Proceedings of the National Academy of Sciences of the United States of America. 1994, 91 (17): 7849-7853. 10.1073/pnas.91.17.7849.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. Wilbrecht L, Crionas A, Nottebohm F: Experience affects recruitment of new neurons but not adult neuron number. Journal of Neuroscience. 2002, 22 (3): 825-831.

    CAS  PubMed  Google Scholar 

  52. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.

    CAS  PubMed  Google Scholar 

  53. Levine TD, Gao F, King PH, Andrews LG, Keene JD: Hel-N1: an autoimmune RNA-binding protein with specificity for 3' uridylate-rich untranslated regions of growth factor mRNAs. Mol Cell Biol. 1993, 13 (6): 3494-3504.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Abe R, Yamamoto K, Sakamoto H: Target specificity of neuronal RNA-binding protein, Mel-N1: direct binding to the 3' untranslated region of its own mRNA. Nucleic Acids Res. 1996, 24 (11): 2011-2016. 10.1093/nar/24.11.2011.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Ma WJ, Chung S, Furneaux H: The Elav-like proteins bind to AU-rich elements and to the poly(A) tail of mRNA. Nucleic Acids Res. 1997, 25 (18): 3564-3569. 10.1093/nar/25.18.3564.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Akamatsu W, Okano HJ, Osumi N, Inoue T, Nakamura S, Sakakibara S, Miura M, Matsuo N, Darnell RB, Okano H: Mammalian ELAV-like neuronal RNA-binding proteins HuB and HuC promote neuronal development in both the central and the peripheral nervous systems. Proc Natl Acad Sci USA. 1999, 96 (17): 9885-9890. 10.1073/pnas.96.17.9885.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. Hambardzumyan D, Sergent-Tanguy S, Thinard R, Bonnamain V, Masip M, Fabre A, Boudin H, Neveu I, Naveilhan P: AUF1 and Hu proteins in the developing rat brain: implication in the proliferation and differentiation of neural progenitors. J Neurosci Res. 2009, 87 (6): 1296-1309. 10.1002/jnr.21957.

    CAS  PubMed  Google Scholar 

  58. Sillje HH, Takahashi K, Tanaka K, Van Houwe G, Nigg EA: Mammalian homologues of the plant Tousled gene code for cell-cycle-regulated kinases with maximal activities linked to ongoing DNA replication. EMBO J. 1999, 18 (20): 5691-5702. 10.1093/emboj/18.20.5691.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Sillje HH, Nigg EA: Identification of human Asf1 chromatin assembly factors as substrates of Tousled-like kinases. Curr Biol. 2001, 11 (13): 1068-1073. 10.1016/S0960-9822(01)00298-6.

    CAS  PubMed  Google Scholar 

  60. Blackwell TK, Walker AK: Transcription elongation: TLKing to chromatin?. Curr Biol. 2003, 13 (23): R915-916. 10.1016/j.cub.2003.11.013.

    CAS  PubMed  Google Scholar 

  61. Carrera P, Moshkin YM, Gronke S, Sillje HH, Nigg EA, Jackle H, Karch F: Tousled-like kinase functions with the chromatin assembly pathway regulating nuclear divisions. Genes Dev. 2003, 17 (20): 2578-2590. 10.1101/gad.276703.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. Rouault JP, Puisieux A, Samarut C, Guehenneux F, Berthet C, Rimokh R, Falette N, Magaud JP: Involvement of the BTG genes family in the control of cell cycle and DNA repair. Experimental Hematology. 1997, 25 (8): 229-229.

    Google Scholar 

  63. Corjay MH, Kearney MA, Munzer DA, Diamond SM, Stoltenborg JK: Antiproliferative gene BTG1 is highly expressed in apoptotic cells in macrophage-rich areas of advanced lesions in Watanabe heritable hyperlipidemic rabbit and human. Laboratory Investigation. 1998, 78 (7): 847-858.

    CAS  PubMed  Google Scholar 

  64. Li F, Liu J, Park ES, Jo M, Curry TE: The B cell translocation gene (BTG) family in the rat ovary: hormonal induction, regulation, and impact on cell cycle kinetics. Endocrinology. 2009, 150 (8): 3894-3902. 10.1210/en.2008-1650.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. Hall JA, Georgel PT: CHD proteins: a diverse family with strong ties. Biochem Cell Biol. 2007, 85 (4): 463-476. 10.1139/O07-063.

    CAS  PubMed  Google Scholar 

  66. Marfella CG, Imbalzano AN: The Chd family of chromatin remodelers. Mutat Res. 2007, 618 (1-2): 30-40.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Bandres E, Malumbres R, Cubedo E, Honorato B, Zarate R, Labarga A, Gabisu U, Sola JJ, Garcia-Foncillas J: A gene signature of 8 genes could identify the risk of recurrence and progression in Dukes' B colon cancer patients. Oncology Reports. 2007, 17 (5): 1089-1094.

    CAS  PubMed  Google Scholar 

  68. Kulkarni S, Nagarajan P, Wall J, Donovan DJ, Donell RL, Ligon AH, Venkatachalam S, Quade BJ: Disruption of chromodomain helicase DNA binding protein 2 (CHD2) causes scoliosis. Am J Med Genet A. 2008, 146A (9): 1117-1127. 10.1002/ajmg.a.32178.

    CAS  PubMed  PubMed Central  Google Scholar 

  69. Nagarajan P, Onami TM, Rajagopalan S, Kania S, Donnell R, Venkatachalam S: Role of chromodomain helicase DNA-binding protein 2 in DNA damage response signaling and tumorigenesis. Oncogene. 2009, 28 (8): 1053-1062. 10.1038/onc.2008.440.

    CAS  PubMed  PubMed Central  Google Scholar 

  70. Bustin M, Reeves R: High-mobility-group chromosomal proteins: Architectural components that facilitate chromatin function. Progress in Nucleic Acid Research and Molecular Biology. 1996, 54: 35-100. Vol 54

    CAS  PubMed  Google Scholar 

  71. Grasser KD: HMG1 and HU proteins: architectural elements in plant chromatin. Trends in Plant Science. 1998, 3 (7): 260-265. 10.1016/S1360-1385(98)01259-X.

    Google Scholar 

  72. Hall J, Thomas KL, Everitt BJ: Cellular imaging of zif268 expression in the hippocampus and amygdala during contextual and cued fear memory retrieval: Selective activation of hippocampal CA1 neurons during the recall of contextual memories. Journal of Neuroscience. 2001, 21 (6): 2186-2193.

    CAS  PubMed  Google Scholar 

  73. Bustin M: At the crossroads of necrosis and apoptosis: signaling to multiple cellular targets by HMGB1. Sci STKE. 2002, 2002 (151): pe39-

    PubMed  Google Scholar 

  74. Guazzi S, Strangio A, Franzi AT, Bianchi ME: HMGB1, an architectural chromatin protein and extracellular signalling factor, has a spatially and temporally restricted expression pattern in mouse brain. Gene Expression Patterns. 2003, 3 (1): 29-33. 10.1016/S1567-133X(02)00093-5.

    CAS  PubMed  Google Scholar 

  75. Bassi R, Giussani P, Anelli V, Colleoni T, Pedrazzi M, Patrone M, Viani P, Sparatore B, Melloni E, Riboni L: HMGB1 as an autocrine stimulus in human T98G glioblastoma cells: role in cell growth and migration. Journal of Neuro-Oncology. 2008, 87 (1): 23-33. 10.1007/s11060-007-9488-y.

    CAS  PubMed  Google Scholar 

  76. Ballif BA, Arnaud L, Arthur WT, Guris D, Imamoto A, Cooper JA: Activation of a Dab1/CrkL/C3G/Rap1 pathway in Reelin-stimulated neurons. Curr Biol. 2004, 14 (7): 606-610. 10.1016/j.cub.2004.03.038.

    CAS  PubMed  Google Scholar 

  77. Yip YP, Kronstadt-O'Brien P, Capriotti C, Cooper JA, Yip JW: Migration of sympathetic preganglionic neurons in the spinal cord is regulated by reelin-dependent Dab1 tyrosine phosphorylation and CrkL. Journal of Comparative Neurology. 2007, 502 (4): 635-643. 10.1002/cne.21318.

    CAS  Google Scholar 

  78. Matsuki T, Pramatarova A, Howell BW: Reduction of Crk and CrkL expression blocks reelin-induced dendritogenesis. J Cell Sci. 2008, 121 (Pt 11): 1869-1875.

    CAS  PubMed  Google Scholar 

  79. Hubbard TJ, Aken BL, Ayling S, Ballester B, Beal K, Bragin E, Brent S, Chen Y, Clapham P, Clarke L, et al: Ensembl 2009. Nucleic Acids Res. 2009, 37 (Database): D690-697. 10.1093/nar/gkn828.

    CAS  PubMed  Google Scholar 

  80. Brennan PA, Schellinck HM, Keverne EB: Patterns of expression of the immediate-early gene egr-1 in the accessory olfactory bulb of female mice exposed to pheromonal constituents of male urine. Neuroscience. 1999, 90 (4): 1463-1470. 10.1016/S0306-4522(98)00556-9.

    CAS  PubMed  Google Scholar 

  81. Schafer M, Brauer AU, Savaskan NE, Rathjen FG, Brummendorf T: Neurotractin/kilon promotes neurite outgrowth and is expressed on reactive astrocytes after entorhinal cortex lesion. Molecular and Cellular Neuroscience. 2005, 29 (4): 580-590. 10.1016/j.mcn.2005.04.010.

    PubMed  Google Scholar 

  82. Hashimoto T, Yamada M, Maekawa S, Nakashima T, Miyata S: IgLON cell adhesion molecule Kilon is a crucial modulator for synapse number in hippocampal neurons. Brain Research. 2008, 1224: 1-11.

    CAS  PubMed  Google Scholar 

  83. Ishii N, Wanaka A, Tohyama M: Increased expression of NLRR-3 mRNA after cortical brain injury in mouse. Brain Res Mol Brain Res. 1996, 40 (1): 148-152.

    CAS  PubMed  Google Scholar 

  84. Bormann P, Roth LWA, Andel D, Ackermann M, Reinhard E: zfNLRR, a novel leucine-rich repeat protein is preferentially expressed during regeneration in zebrafish. Molecular and Cellular Neuroscience. 1999, 13 (3): 167-179. 10.1006/mcne.1999.0742.

    CAS  PubMed  Google Scholar 

  85. Josephson A, Trifunovski A, Scheele C, Widenfalk J, Wahlestedt C, Brene S, Olson L, Spenger C: Activity-induced and developmental downregulation of the Nogo receptor. Cell Tissue Res. 2003, 311 (3): 333-342.

    CAS  PubMed  Google Scholar 

  86. Endo T, Spenger C, Tominaga T, Brene S, Olson L: Cortical sensory map rearrangement after spinal cord injury: fMRI responses linked to Nogo signalling. Brain. 2007, 130 (Pt 11): 2951-2961.

    PubMed  Google Scholar 

  87. Dong S, Clayton DF: Habituation in songbirds. Neurobiol Learn Mem. 2009, 92 (2): 183-188. 10.1016/j.nlm.2008.09.009.

    PubMed  Google Scholar 

  88. Replogle K, Arnold AP, Ball GF, Band M, Bensch S, Brenowitz EA, Dong S, Drnevich J, Ferris M, George JM, 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-10.1186/1471-2164-9-131.

    PubMed  PubMed Central  Google Scholar 

  89. Cheng HY, Clayton DF: Activation and habituation of extracellular signal-regulated kinase phosphorylation in zebra finch auditory forebrain during song presentation. Journal of Neuroscience. 2004, 24 (34): 7503-7513. 10.1523/JNEUROSCI.1405-04.2004.

    CAS  PubMed  Google Scholar 

  90. Nagaraja AK, Andreu-Vieyra C, Franco HL, Ma L, Chen R, Han DY, Zhu H, Agno JE, Gunaratne PH, DeMayo FJ, et al: Deletion of Dicer in somatic cells of the female reproductive tract causes sterility. Mol Endocrinol. 2008, 22 (10): 2336-2352. 10.1210/me.2008-0142.

    CAS  PubMed  PubMed Central  Google Scholar 

  91. Ma L, Buchold GM, Greenbaum MP, Roy A, Burns KH, Zhu H, Han DY, Harris RA, Coarfa C, Gunaratne PH: Correction: GASZ Is Essential for Male Meiosis and Suppression of Retrotransposon Expression in the Male Germline. PLoS Genet. 2009, 5 (12):

  92. Schuster P, Fontana W, Stadler PF, Hofacker IL: From sequences to shapes and back: a case study in RNA secondary structures. Proc Biol Sci. 1994, 255 (1344): 279-284. 10.1098/rspb.1994.0040.

    CAS  PubMed  Google Scholar 

  93. Kalafus KJ, Jackson AR, Milosavljevic A: Pash: efficient genome-scale sequence anchoring by Positional Hashing. Genome Res. 2004, 14 (4): 672-678. 10.1101/gr.1963804.

    CAS  PubMed  PubMed Central  Google Scholar 

  94. Coarfa C, Milosavljevic A: Pash 2.0: scaleable sequence anchoring for next-generation sequencing technologies. Pac Symp Biocomput. 2008, 102-113.

    Google Scholar 

  95. Thomson T, Lin H: The biogenesis and function of PIWI proteins and piRNAs: progress and prospect. Annu Rev Cell Dev Biol. 2009, 25: 355-376. 10.1146/annurev.cellbio.24.110707.175327.

    CAS  PubMed  PubMed Central  Google Scholar 

  96. Royo H, Cavaille J: Non-coding RNAs in imprinted gene clusters. Biol Cell. 2008, 100 (3): 149-166. 10.1042/BC20070126.

    CAS  PubMed  Google Scholar 

  97. Gu P, Reid JG, Gao X, Shaw CA, Creighton C, Tran PL, Zhou X, Drabek RB, Steffen DL, Hoang DM, et al: Novel microRNA candidates and miRNA-mRNA pairs in embryonic stem (ES) cells. PLoS One. 2008, 3 (7): e2548-10.1371/journal.pone.0002548.

    PubMed  PubMed Central  Google Scholar 

  98. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Sarah London for useful discussions and contributions to the text. Supported by NIH RO1 NS045264 and RO1 NS051820 (to D.F.C).

Note Added in Proof:

The novel miRNA referred to here as "miR-2954-3p" is now identified in miRBase as "miR-2954". The novel miRNA referred to here as "miR-2954-5p" is now identified in miRBase as "miR-2954*".

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to David F Clayton.

Additional information

Authors' contributions

PHG coordinated the work of Illumina RNA-seq and prepared the manuscript. YL conducted the song exposure experiments and subsequent dissections and RNA extractions, performed TaqMan qPCR, analyzed differentially expressed miRNAs and participated in drafting the manuscript. ALB performed Illumina RNA-seq and Northern blot. JD helped analyze expression data of Illumina RNA-seq and TaqMan qPCR. CC, JBT, CJC, JHK, and AM participated in mapping and analyzing Illumina RNA-seq data. MW and SGJ helped with miRNA sequence annotation. DFC designed and coordinated the study and drafted the manuscript. All authors read and approved the manuscript.

Preethi H Gunaratne, Ya-Chi Lin contributed equally to this work.

Electronic supplementary material

12864_2010_10144_MOESM1_ESM.XLS

Additional file 1: Supplemental tables.xls. This one file contains all four Supplemental Tables, each as a separate worksheet. Table S1 ("1 overview") is a summary of Illumina sequence read alignments for six pools of RNA from zebra finch auditory forebrain responding to song versus silence, and shows the distribution of sequence reads in relation to multiple genomes and multiple annotations in the current genomic databases. Table S2 ("2 novel hairpins") gives detailed alignments of putative pre-miRNAs and read sequences. Table S3 ("3 novel genes") shows annotations of all novel miRNA loci mapped in genome assemblies of zebra finch, chicken or human. Table S4 ("4 all read counts") gives read counts and current annotation in miRBase of all conserved and novel miRNAs, with statistics. (XLS 2 MB)

12864_2010_10144_MOESM2_ESM.DOC

Additional file 2: Supplemental figures.doc. This one file contains all three supplemental figures. Figure S1 is a Venn diagram of numbers of miRNAs with significant differential expression in response to novel song in three Illumina experiments. Figure S2 shows a comparative mapping in other avian transcriptomes of tgu-mir-2954. Figure S3 demonstrates the song-specificity of the miRNA response, using TaqMan to compare the levels of specific miRNAs in animals from groups that heard song, matching song-enveloped noise, or silence. (DOC 698 KB)

Authors’ original submitted files for images

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Gunaratne, P.H., Lin, YC., Benham, A.L. et al. Song exposure regulates known and novel microRNAs in the zebra finch auditory forebrain. BMC Genomics 12, 277 (2011). https://doi.org/10.1186/1471-2164-12-277

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-12-277

Keywords