|
|
||||||||
1McGovern Institute and Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, Massachusetts; and 2Department of Physics and Department of Psychology, Pennsylvania State University, University Park, Pennsylvania
Submitted 6 September 2006; accepted in final form 16 December 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The extremely stereotyped vocal behavior of songbirds is implemented by a discrete and anatomically well-characterized neural circuit (for a review, see Brenowitz et al. 1997
). The song control system consists of two major pathways: the descending motor pathway and the anterior forebrain pathway (AFP) (Fig. 1A). The descending motor pathway consists of brain nuclei necessary for song production: HVC (used as a proper name), robust nucleus of the arcopallium (RA) and the brain stem motor nucleus innervating the muscles of the syrinx (Nottebohm et al. 1976
; Vicario 1991
; Vicario and Nottebohm 1988
). Nucleus HVC also projects to the anterior forebrain pathway (AFP). The AFP consists of basal ganglia (area X), the medial nucleus of the dorsolateral thalamus (DLM), and lateral magnocellular nucleus of the anterior nidopallium (LMAN), which projects to RA (Bottjer et al. 1989
; Luo and Perkel 1999a
,b
; Okuhata and Saito 1987
). Lesions of the AFP leave song production intact in adult birds but cause severe deficits of vocal learning in juvenile birds (Bottjer et al. 1984
; Scharff and Nottebohm 1991
; Sohrabji et al. 1990
). Neurons in HVC and the AFP exhibit auditory responses in anesthetized birds and have especially strong responses to the bird's own song (BOS) (Doupe and Konishi 1991
; Margoliash 1986
). Evidence available to date suggests that auditory responses of AFP neurons are driven by X-projecting HVC (HVC(X)) neurons (Doupe and Konishi 1991
; Vicario and Yohay 1993
; Williams and Nottebohm 1985
).
|
Because of the crucial importance of nucleus HVC for song production, and the possible role of the HVC projection to the AFP for song learning, HVC is one of the most-studied nuclei in the song control system. Neural types and synaptic connections within HVC have been characterized (Mooney 2000
; Mooney and Prather 2005
), and both multi- and single-unit recordings of singing-related activity in HVC have been obtained (Hahnloser et al. 2002
; McCasland 1987
; Yu and Margoliash 1996
). However, systematic characterization of the singing-related activity of identified HVC neurons has not been done. In this work, we identified three different classes of HVC neurons in zebra finches by antidromic stimulation and recorded their patterns of activity during singing. We characterized the patterns of activity of each class of neurons and analyzed the relation of the neural activity to the acoustic structure of the song.
Because HVC(X) neurons are implicated in providing auditory input to the AFP, we asked a question whether their singing-related activity has an auditory component. Auditory activity of HVC(X) neurons could constitute a vocal evaluation signal, which is used for vocal learning and error correction. To elucidate an auditory-related component in the activity of HVC(X) neurons, we recorded the firing patterns of HVC(X) neurons in juvenile birds during normal singing and during singing in the presence of distorted auditory feedback (DAF). By comparing firing patterns of HVC(X) neurons during singing with and without DAF, we tested the hypothesis that these neurons carry an auditory signal during singing.
| METHODS |
|---|
|
|
|---|
Subjects were adult (>90 days posthatch) and juvenile (7090 days posthatch) zebra finches. The care and experimental manipulation of the animals was carried out in accordance with guidelines of the National Institutes of Health and has been reviewed and approved by the local Institutional Animal Care and Use Committee. Data were collected from eight adult and five juvenile zebra finches. Before surgery, anesthesia was produced with 13% isofluorane in oxygen. Single-unit recordings were made with a three-channel miniature motorized microdrive using previously described techniques (Fee and Leonardo 2001
). Bipolar stimulating electrodes were implanted in RA and Area X for antidromic stimulation of HVC (Fig. 1A). Antidromic identification was done by stimulation RA or Area X just above the threshold for eliciting a reliable evoked spike. Classification of neurons as RA-projecting (HVC(RA)), X-projecting (HVC(X)), or putative interneurons (HVC(i), neurons not projecting either to RA or X) was done by measuring the latency variability of the first evoked spike and was verified using collision protocol. Detailed description of the antidromic identification protocol is presented elsewhere (Fee et al. 2004
; Hahnloser et al. 2006
).
Alignment of acoustic and neural patterns
If one makes a raster plot of the firing pattern of a neuron by aligning the spike trains at the beginning of the song motif, variability in the duration of the motif results in typically 1020 ms jitter in the spike trains at the end of the motif. To take into account this source of variability, spike trains were aligned using each syllable as a time reference. A single song motif was chosen as a template for time alignment. All other motifs were aligned to the template motif using syllable onsets and offsets as reference points and using piecewise-linear time warping (Glaze and Troyer 2006
; Leonardo 2002
). Using the same piecewise-linear time-warping, spike times were mapped onto the template motif.
This alignment procedure dramatically reduced the spike time variability due to the variations in the song motif duration. In the rest of the paper, unless noted otherwise, all the data analysis was done using time-warped spike trains aligned to the song template.
Discretized firing patterns, population activity of neurons and characterization of syllable song pattern
Discretized mean firing patterns ni(t) were computed for each neuron by counting the mean number of spikes occurring during nonoverlapping time windows
t = 2 ms long
![]() |
![]() |
Calculation of coherence
To characterize the relation between population activity of interneurons and the syllable structure of the song, we computed the coherence between the syllable pattern S(t) and the population average of interneurons Ipop(t). The advantage of working with the coherence versus time-domain correlation function (i.e., working in frequency domain vs. time domain) is the simplicity of the estimating confidence levels for the null hypothesis (Jarvis and Mitra 2001
).
Direct spectral estimators of S(t) and Ipop(t) were computed using the multi-taper method (Thompson 1982
) with time-bandwidth product NW = 7/2. Functions S(t) and Ipop(t) were convolved with taper windows and zero-padded to nfft = 512 points, and then fast Fourier transforms
i(
) and
pop(t) were computed. Analysis was done for five birds in which five or more interneurons per bird were recorded. Data from different birds were regarded as different samples. The coherency
(
) was computed by summing cross-spectral products over the samples and over the tapers and normalizing (Thompson 1982
). The coherence C(
) is defined as C(
) = |
(
)|2. The (p x 100)th percentile confidence intervals for the null hypothesis [i.e., that S(t) and Ipop(t) are independent] of coherence C(
) are given by
![]() |
Correlation coefficients and time-domain correlation functions
Pairwise correlations between the firing patterns of different neurons were characterized by computing Pearson correlation coefficients between the discretized mean firing patterns of neurons (defined in the preceding text). When computing correlation coefficients between the activity of an individual interneuron ni(t) and the population activity of interneurons, we did not include the interneuron itself in the estimator of the population activity to avoid artificial increase of the correlation
![]() |
We also analyzed the correlations between the population activities of HVC(RA) and HVC(X) neurons. We used two different methods for combining data from different birds: we computed correlation coefficients between the populations of HVC(X) and HVC(RA) neurons and averaged them over birds, and we computed correlation coefficients between the populations of HVC(X) and HVC(RA) neurons and averaged over birds weighting the correlation coefficients by the number of HVC(RA) neurons recorded in each birdthis way the birds in which more neurons were recorded received more weight. To assess the significance of the computed correlations, we randomly time-shifted the activity of individual neurons to construct shuffled population activities of HVC(X) and HVC(RA) neurons. Then we computed the correlation coefficients between shuffled population activities of HVC(RA) and HVC(X) neurons. The shuffling procedure was repeated n = 10,000 times, and the P = 0.05 confidence intervals for the correlations of the shuffled data were computed.
To study whether the activity of the projection neurons is modulated coherently with the syllable temporal pattern, we computed the mean correlation functions between the activity of a neuron and the syllable patterns of the song: CRA-S(
) =
corr[RAk(t), Sk(t
)]
k and CX-S(
) =
corr[Xk(t), Sk(t
)]
k (averaging is over all neurons in all birds).
To assess the significance of this correlation, we generated randomized data sets by introducing random shifts tkrand to the mean firing patterns of projection neurons and computing the mean correlation functions between the syllable pattern S(t) and the shifted neural firing patterns: CRA-Srand(
) =
corr[RAk(t + tkrand), Sk(t
)]
k and CX-Srand(
) =
corr[Xk(t + tkrand), Sk(t
)]
k. Because the relative timing between the neural activity and the syllable pattern is a priori unknown, we used a fairly conservative criterion for statistical significance: the mean correlation functions CRA-S(
) and CX-S(
) were considered significant if they exceeded the P = 0.01 confidence intervals for the distribution of maxima or if they were below the P = 0.01 confidence interval for the distribution of minima of the functions CRA-Srand(
) and CX-Srand(
) in the interval 200 ms < t < 200 ms.
Correlations between multiple bursts of X-projecting neurons and sound: burst-triggered similarity
We analyzed whether multiple bursts of HVC(X) neurons are related to the occurrence of similar sounds in the song. To analyze whether repeating bursts of individual HVC(X) neurons are preceded or followed by similar sounds, we computed the burst-triggered similarity function averaged over all pairs of bursts of individual HVC(X) neurons
![]() |
We used three different methods for constructing the similarity matrix between sounds in the song: syllable-syllable similarity, which only takes into account syllable pattern of the song and ignores fine-scale acoustic correlation; similarity of fine acoustic structure correlation using sound spectrograms; and similarity using acoustic features of the song (Tchernichovski et al. 2000
).
SYLLABLE-SYLLABLE SIMILARITY.
To take into account similarity associated with syllable structure of the song only, we used the following rule for constructing the syllable-syllable acoustic similarity matrix
![]() |
![]() |
SPECTRAL SIMILARITY.
Song spectrograms were computed over time window
t = 8 ms using the multi-taper method with NW = 3/2. If P(t,f) is time- and frequency-dependent acoustic power, the logarithmic spectrogram (called for shortness "spectrogram" in the following text) is B(t f) = log P(t,f). One can define the matrix of acoustic similarities by using the spectra at times t1 and t2: C(t1,t2) = corr[B(t1,f), B(t2,f)]. The acoustic similarity matrix C(t1,t2) defined this way has a lot of structure due to the syllable pattern of the song. The syllable-pattern structure in C(t1,t2) makes it difficult to analyze the similarities in fine acoustic structure (i.e., similar notes). To circumvent this problem, we calculated the mean temporal and spectral modes of the spectrogram B(t,f), normalized the spectrum at each time by the total acoustic energy and then subtracted the projection of the spectrum onto the mean spectral mode (Leonardo and Fee 2005
)
![]() |
(t) = log P(
P(t,f)
f) is the mean temporal mode,
(f) = log (
P(t,f)
t) is the mean spectral mode of the song and fmax = 8 kHz. Using the normalized spectrogram Bnorm(t,f), we calculated the spectral similarity matrix [this quantity was named correlation matrix in (Leonardo and Fee 2005
Matrix CS(t1,t2) has most of its structure due to similar notes (e.g.,harmonic stacks). The burst-triggered spectral similarity function averaged over multiple bursts of HVC(X) neurons is defined as
![]() |
FEATURE-BASED SIMILARITY.
We also computed acoustic features of the songs and constructed a feature-based similarity matrix
![]() |
are time-dependent scaled acoustic features (Tchernichovski et al. 2000
![]() |
![]() |
ASSESSMENT OF SIGNIFICANCE.
To assess the significance of the burst-triggered similarities, we generated pairs of bursts distributed uniformly during the song motif and computed the burst-triggered similarity functions using simulated bursts Crand(
).
For the burst-triggered syllable-syllable similarity, we repeated the procedure n = 3,000 times and computed the distribution of maxima and minima of Crand(
) in the interval 200 ms <
< 200 ms. Values of observed burst-triggered similarity functions were considered significant when they either exceeded the P = 0.01 confidence interval for the distribution of maxima of Crand(
) or were below the P = 0.01 confidence interval for the distribution of minima of Crand(
).
For the burst-triggered spectral similarity and the feature-based similarity, we computed the distribution of Crand(
= 0). The values of observed burst-triggered spectral and feature-based similarity functions were considered significant when they were outside the P = 0.01 confidence intervals for the distribution of Crand(
= 0). This significance criterion is less stringent than the one used for burst-triggered syllable-syllable similarity, which is justified because the observed burst-triggered spectral and feature similarities are not significant as assessed by confidence intervals for the distribution of Crand(
= 0).
Correlation times of vocal dynamics and neural dynamics in HVC
To compute acoustic correlation times we used the spectral similarity matrix Cspectr(t1, t2) defined in the preceding text. For computing the neural similarity matrix, we constructed a neural activity vector
(t) (the ith component of this vector is the mean firing pattern of the ith interneuron ni(t)) and computed Pearson correlation coefficient between the components of vector
(t) at different times: CN(t1,t2) = corr(
(t1),
(t2)).
The procedure for computing the correlation time for matrices CS(t1,t2) and CN(t1,t2) is as previously described (Leonardo and Fee 2005
). For each point (t0,t0) on the matrix diagonal, we determined the maximal size of the square around this point for which all matrix values inside the square are larger than the threshold value C = 0.5. The size of that square is the characteristic time
(t0). We then searched for the peaks of the function
(t0), starting from the largest peak, and for each peak occurring at time tp, we extended the peak value of the characteristic time
peak over the time interval
peak long so that the correlation time tcorr [t
(tp
peak/2,tp +
peak/2)] =
peak. Such procedure results in the correlation time being constant over the length of the interval where the values of the matrix are high (for example, sound correlation time tacoustic(t) is constant during a harmonic stack and equal to the duration of the harmonic stack).
Distorted auditory feedback
For the generation of DAF during neural recordings, loud bursts of broad-band noise were played back to the bird during 50% of song renditions (sound volume: 8090 dB SPL, the burst duration typically 70100 ms). Using the Golay pulse calibration (Golay 1961
), noise cancellation was performed and song signals were recovered (Leonardo 2002
).
Assessments of effects of DAF
Characterization of systematic changes in song structure as a result of DAF exposure was done by computing a spectral similarity score between each rendition of a syllable and a template syllable produced 34 days before the onset of DAF. The similarity score between two syllables was defined as a Pearson correlation coefficient between the multi-taper spectrograms (Thompson 1982
) of the syllables (time-bandwidth product: NW = 3/2). Spectra were averaged over three tapers and FFTs were computed over time windows 8-ms-long overlapping by 7 ms. Spectral components between 0.3 and 8 kHz were used to compute the similarity scores. This procedure for computing similarity scores typically yielded values of 0.80.9 for different renditions of the same syllable and values of 0.40.7 for different syllables (depending on how similar the different syllables are).
To assess whether DAF caused significant degradation of the acoustic structure of the song syllable, we computed similarity scores of each syllable to the template syllable and separated the scores in two groups: similarity scores on catch trials on the days when DAF was presented, starting on the second day of DAF presentation and through the first day after DAF presentation stopped, and similarity scores on the days when no DAF was present, both before and after the period of DAF presentation. The distributions of scores in these two groups were compared with each other, and the degradation of acoustic structure was considered significant if the similarity scores during the period of DAF presentation were significantly (P < 0.05, 1-sided KS test) lower than when no DAF was presented.
Comparison of firing patterns of HVC(X) neurons with and without DAF
To compare the firing patterns of HVC(X) neurons with and without DAF, for each burst of each HVC(X) neuron recorded, we computed several burst parameters: burst onset times during the song motif tonset, burst duration tdur, the number of spikes in a burst nburst, the mean firing rate during the burst and the interspike intervals in a burst (for the calculation of interspike intervals, we used nontime-warped spike times). We compared the distributions of each of these parameters during normal singing and during singing with DAF. We used the two-sided KS test and the one-way ANOVA to assess statistical differences in the distributions of the burst parameters.
| RESULTS |
|---|
|
|
|---|
We antidromically identified and obtained single-unit recordings from 28 HVC(RA) neurons (8 birds), 103 HVC(X) neurons (7 birds), and 52 putative HVC interneurons (8 birds) during singing. The three types of neurons have very distinct patterns of activity during singing (Fig. 2).
|
An additional eight units were identified by antidromic stimulation as HVC(RA) neurons, but these were not spontaneously active nor were they active during singing or during any calls that could be elicited. Because spikes were only observed from these neurons during antidromic stimulation, there is less certainty about the classification of these units
The burst duration of HVC(RA) neurons was 5.1 ± 1.8 ms (minimum = 3.0 ms, maximum = 9.1 ms), and the bursts contained 4.3 ± 1.3 spike (minimum = 2.0, maximum = 6.8; ranges are means ± 1 SD unless specified otherwise). The average firing rate during the bursts was 650 ± 150 Hz (minimum = 330 Hz, maximum = 890 Hz). The bursts were extremely reliable and were tightly locked to the vocalization the root mean square (RMS) jitter of the first spike is 0.73 ± 0.3 ms (minimum = 0.4 ms, maximum = 2.5 ms). The ratio of the SD of the number of spikes in a burst to the mean number of spikes in a burst was 0.10 ± 0.05.
HVC(X) neurons also fired sparse bursts during singing. Most HVC(X) neurons (72 cells) exhibited stereotyped high-frequency bursts during singing that were tightly locked to the vocalization. The bursts of these HVC(X) neurons were similar to the bursts of HVC(RA) neurons, but HVC(X) neurons fired up to 4 bursts per motif (1.3 ± 0.9 burst/motif). The burst duration was 6.3 ± 3.3 ms (minimum = 2.7 ms, maximum = 14 ms), the number of spikes in a burst was 3.7 ± 1.4 (minimum = 1.75, maximum = 7.0), and the average firing rate during the burst was 434 ± 130 Hz (minimum = 110 Hz, maximum = 750 Hz). The rms temporal jitter of the first spike in a burst was 1.0 ± 0.5 ms (minimum = 0.3 ms, maximum = 2.3 ms; 114 bursts in 72 neurons). The ratio of the SD of the number of spikes in a burst to the mean number of spikes in a burst is 0.16 ± 0.12. In awake nonsinging birds, HVC(X) neurons fired single spikes at an average rate of 1.5 ± 0.5 Hz.
About 10% of HVC(X) cells (n = 13) exhibited bursts with lower firing rate and higher variability in a number of parameters. These cells burst up to three times per motif and their bursts consisted of 3.2 ± 1.6 spikes. However, the bursts of these neurons were longer (burst duration: 10 ± 5 ms) and had a lower firing rate (150 ± 72 Hz), and the number of spikes in a burst was more variable: the ratio of the SD of the number of spikes to the mean number of spikes is 0.64 ± 0.32. The bursts of these HVC(X) cells were less tightly locked to the vocalization; the rms jitter of the first spike in the burst was 8.2 ± 7.8 ms. A scatter plot of the rms jitter of the first spike in a burst versus average firing rate during the burst for HVC(X) neurons is shown in Fig. 2B. Bursts with higher temporal jitter tended to have lower firing rates. Although the two groups of HVC(X) neurons had significantly different firing patterns, the sample size is too small to determine if this less stereotyped group of neurons form a distinct cluster in the space of burst characteristics.
In addition, 18 HVC(X) neurons were either completely silent during singing or produced sparse, irregular single spikes (<1 spike/s). The distribution of the number of bursts generated by HVC(X) neurons per song motif is presented in Fig. 2C.
Similarly to HVC neurons previously described (Yu and Margoliash 1996
) HVC interneurons (HVCi) were active throughout the vocalization and had tonic patterns of activity. Their firing patterns were much less stereotyped than the firing patterns of projection neurons, but there is a stable underlying time-dependent firing rate pattern (Fig. 2A). Interneurons generated 69 ± 31 spike/song motif (minimum = 20, maximum = 208). During singing they had average firing rates of 95 ± 40 Hz (minimum = 25 Hz, maximum = 262 Hz). The SD of the number of spikes of an interneuron per song motif is 6 ± 3 spike. In an awake nonsinging bird, HVC interneurons had an average spontaneous firing rate of 8 ± 6 Hz.
Population activity of different classes of HVC neurons
We asked whether the activity of different classes of HVC neurons is synchronized during singing. Because of the stereotypical patterns of activity, the firing patterns of different neurons are strongly correlated but this does not necessarily imply synchrony. Multiunit data suggest that there is a synchronized pattern of population activity in HVC (Vu et al. 1998
). However, only single-unit recordings from identified neurons can tell us whether neurons of different classes tend to fire in synchrony.
The average correlation coefficient between the activity of an interneuron and the population activity of other interneurons was 0.14 ± 0.11 (see METHODS), larger than the pairwise correlations between individual interneurons (1-way KS test, P < 2 x 105). These data show that the firing patterns of individual interneurons are weakly synchronized.
To address the question of whether neurons of different classes are synchronized with each other, we analyzed pairwise zero time-lag correlations between the firing patterns of HVC neurons. The mean and the SD of the pairwise correlation coefficients for different types of neurons are presented in Table 1. To test whether the synchrony in the firing patterns of individual neurons is statistically significant, we introduced random time shifts to the spike trains and calculated pairwise correlation coefficients for the shifted spike trains. We used KS test to test the hypothesis that the distributions of correlation coefficients for nonshifted and shifted spike trains are different. Statistically significant correlations were found between pairs of interneurons (KS test, P < 108) and between interneurons and HVC(X) neurons (P < 5 x 104). Thus activity of HVC(X) neurons is synchronized with the coherent modulation of interneuron activity in HVC.
|
Relation of activity of HVC neurons to syllable patterning
We analyzed the relation between the population activity of the three classes of HVC neurons and the syllable pattern of the song. At the simplest level, zebra finch song consists of song syllables separated by intervals of silence. We ask whether firing patterns of HVC neurons are related to the pattern of intervals of sound and silence in the song motif. For this question, we disregard the spectral content of the song and characterize the song with syllables represented by 1 and silences by 0. We refer to this pattern of sounds and silences as the syllable pattern of the song. The syllable pattern S(t), and the population averages of the three classes of neurons [RApop(t), Xpop(t) and Ipop(t), see METHODS] for one of the birds (bird 2) are shown in Fig. 3A.
|
0.25 (P < 0.01). The frequency band between 8 and 15 Hz corresponds to the time scale of the song syllables (
70130 ms), showing that the population activity of interneurons was modulated coherently with the syllable pattern of the song.
The correlation function between individual interneurons and S(t), averaged over all recorded interneurons from all birds, is shown in Fig. 3C. On average, the activity of an individual interneuron was weakly correlated with the syllable pattern of the song (correlation coefficient
0.08) and leads the song syllable pattern by 3050 ms. The peak of the correlation between Ipop(t) and S(t) (at a lag of 45 ms) was 0.2, averaged over birds. Thus the population activity of interneurons was modulated by
20% and was largest preceding song syllables by
45 ms.
To assess the syllable-related modulation of projection neurons, we computed the mean correlation between the firing pattern of a neuron of a given class and the syllable pattern and assessed its significance (see METHODS). The mean correlation functions of HVC(RA) and HVC(X) neurons with the syllable pattern of the song are presented in Fig. 3C. Although the mean correlation function for HVC(RA) neurons exhibits a peak at
40 ms, the value of the peak does not reach statistical significance. It is possible that this peak is real, but a larger dataset would be needed to establish significance.
In summary, the population activity of HVC interneurons exhibited syllable-related modulation, and lead the syllable pattern of the song by
45 ms. The population activity of HVC(RA) and HVC(X) neurons was not significantly modulated with the song syllable pattern in our dataset.
Relation of activity of HVC(X) neurons to acoustic structure of the song
Given the importance of the AFP for song learning and the proposed role of the HVC(X) neurons in the vocal error correction, we asked whether the activity of HVC(X) neurons is related to the acoustic structure of the vocalizations. Although we saw no effects of the song syllable pattern on the population activity of HVC(X) neurons, the occurrence of multiple bursts from HVC(X) neurons allows us to test the possibility that repeated bursts of a given HVC(X) neuron occur in some definite relation to repeated sounds in the song. For example, if the activity of HVC(X) neurons were auditory-related, one might expect that multiple bursts of a given neuron would be preceded by similar sounds. Alternatively, if activity of HVC(X) neurons has a premotor relation to a particular sound in the song, one might expect multiple bursts to precede similar sounds during the song motif (Fig. 4A).
|
![]() |
). The function Csyll(
) has a peak at t
40 ms significant at P < 0.01. This indicates that multiple bursts of an HVC(X) neuron tend to precede similar elements of the song syllable pattern (either sounds or silences). Thus their temporal relation to the syllable pattern is premotor-like.
We next turned to the question of whether multiple HVC(X) bursts are related to similar spectral structure in the song. We computed the burst-triggered spectral similarity function averaged over pairs of bursts of HVC(X) neurons, Cspectr(
), and assessed its significance (Fig. 4C, see METHODS). The spectral similarity did not have any pronounced structure as a function of
, i.e., multiple bursts of an individual HVC(X) neuron were not preceded or followed by similar spectral structure in the song. Analysis of feature-based similarity likewise shows that multiple bursts of HVC(X) neurons are neither preceded nor followed by similar acoustic features in the song.
In summary, multiple bursts of a given HVC(X) neuron are correlated with the syllable structure of the song; Bursts of a given HVC(X) neuron tend to precede similar syllable pattern elements (either sounds or silences) by
40 ms. We find no significant correlations of repeated bursts of individual HVC(X) neuron with repeated spectral structure in the song.
Time scales of vocal dynamics and neural activity in HVC
Vocal output during the bird's song can change on varying time scalesfrom <10 ms during rapid acoustic transitions to
100 ms during harmonic stacks. It has been shown that neural dynamics in RA evolve on fixed time scale (
8 ms) that is not related to the time scale of modulation in the vocalization (Leonardo and Fee 2005
). Consistent with this observation, and with experiments suggesting that RA bursts are driven from HVC (Hahnloser et al. 2002
), we expect that the time scale of HVC projection neuron dynamics will likewise not depend on the timescale of modulation in the vocalizations. Unfortunately, because of their sparse firing patterns, it is difficult to use HVC projection neurons to estimate the time scales of HVC dynamics at all times during the song motif.
In contrast, HVC interneurons are active continuously, and provide a means of estimating the time scale of neural dynamics in HVC. Correlation matrices of sound and interneuron firing patterns (see METHODS) for one of the birds are presented in Fig. 5, A and B. Figure 5C shows time dependence of the correlation time of the song
tacoustic(t) and the correlation time of activity of interneurons
tneural(t). The values of the correlation time
tneural(t) are consistent with the width (FWHM = 3.6 ms) of the central peak of the average autocorrelation function of the individual interneurons (Fig. 5C, inset). The linear regression of neural correlation time shifted forward by
= 20 ms
tneural(t +
) versus acoustic correlation time
tacoustic(t) (Fig. 5D, straight line) had essentially a zero slope, showing that premotor activity in HVC changed with a characteristic time uncorrelated with the characteristic timescales of the vocal output. We varied delay time
between 70 ms (which may correspond to auditory delay) to +70 ms (which exceeds delay time between premotor activity in HVC and vocal output). At all values of
between 70 and +70 ms, there was no significant positive correlation between
tneural(t +
) and
tacoustic(t); at all values 70 ms <
<70 ms the slope of the linear fit of
tneural(t +
) versus
tacoustic(t) was between 0.015 and 0.015.
|
Effects of distorted auditory feedback on the firing patterns of HVC(X) neurons
The most direct test of auditory sensitivity of HVC(X) neurons during singing is to alter what the bird hears while it sings and look for changes in the firing patterns of these neurons. We recorded the firing patterns of HVC(X) neurons in juvenile birds during normal singing and during singing in the presence of DAF (Leonardo 2002
; Leonardo and Konishi 1999
). By comparing the firing patterns of HVC(X) neurons during singing with and without DAF, we directly tested the hypothesis that these neurons carry auditory signal during singing.
Singing-related single-unit activity of a total of 30 HVC(X) neurons (3 juvenile birds) was recorded in the presence and in the absence of distorted auditory feedback. Data were taken when the birds were 7090 days posthatch. The firing patterns of HVC(X) neurons were quite similar with and without distorted feedback, and there were no apparent changes in the activity of these cells (Fig. 6C). The highly stereotypical and temporally precise firing patterns of HVC(X) neurons should make detectable even slight activity changes related to DAF, such as
1 extra spike in a burst, or change of the burst onset time by
1 ms, or changes in the interspike interval during the burst by
300 us. For each burst of each HVC(X) neuron studied, we computed several burst parameters: burst onset times during the song motif, burst duration, the number of spikes in a burst, and the mean firing rate during the burst and interspike intervals in a burst. A summary of the analysis of the burst onset times, burst durations and number of spikes in the bursts of HVC(X) neurons are presented in Fig. 7. Of the 43 bursts analyzed, only 2 bursts showed a statistically significant difference in burst onset time or burst duration (P < 0.05, 2-sided KS test). This would be expected by chanceif there is no change in the firing patterns, analysis of
5% of the bursts would yield statistical significance at P < 0.05 level. None of the 43 bursts analyzed showed statistical differences at P < 0.01 level. We conclude that the singing related firing patterns of HVC(X) neurons were not altered by on-line distorted auditory feedback.
|
|
In zebra finches, playback of loud sounds to the bird during singing does not immediately change the song but leads to gradual song degradation on a time scale from weeks to months (Leonardo and Konishi 1999
). Such degradation has been observed for playback of song syllables or bursts of broadband noise, both referred to as DAF. This suggests that the auditory information acquired during singing leads to plastic changes in the motor program.
We tested the effects of broadband noise DAF on the song of juvenile birds (n = 5 birds). The DAF was triggered by the birds vocalizations with high probability (P = 0.98) for several days, and changes in the songs were analyzed on catch trials.
After several days of distorted feedback exposure, the vocalizations of all five birds changed markedly (Fig. 6, A and B; supplementary Fig. 1).1 The changes observed included syllable stuttering (1 bird), abnormal repetition of an introductory note (2 birds), and changes in the acoustic structure of the song syllables (3 birds); one bird exhibited both introductory note stuttering and syllable structure degradation. Every bird studied exhibited some disruptive effect of DAF. In the three birds that showed degradation of the syllable structure, we found significant (P < 0.05) degradation of 4 syllables of total of 12 syllables in these three birds.
The observed changes in the song were reversibleboth the syllable sequence and the syllable structure recovered to near prefeedback levels after the distorted feedback had been removed. Observation of the chronic and reversible effects of the distorted feedback on the birds song indicates that the distorted feedback protocol used to test the auditory responses of HVC(X) neurons was effective in causing plastic changes in the song control system.
These changes in the vocalizations were similar to those reported previously for adult birds (Leonardo and Konishi 1999
) but occurred much fasterover the course of several days. The higher rate of song degradation is likely due to the fact that we used juvenile birds, consistent with previous observations that song degradation following deafening occurs faster in younger birds (Lombardino and Nottebohm 2000
).
| DISCUSSION |
|---|
|
|
|---|
The three classes of HVC neurons have distinct firing patterns during singing: whereas interneurons have tonic firing patterns, the firing patterns of projection neurons are phasic and sparse. Taking advantage of the stereotyped song motif to align the firing patterns of neurons, we addressed the questions of synchrony among neurons and the relation between neural activity in HVC and vocal output.
Motivated by the idea that HVC(X) neurons carry a vocal performance evaluation signal to the AFP during singing, we directly tested this possibility. We altered the acoustic environment of a juvenile bird in a way that causes gradual song deterioration and studied the effects of such distorted acoustic feedback on the firing patterns of HVC(X) neurons.
In the following text, we summarize and interpret our findings. Our findings present several lines of evidence suggesting that the activity of HVC(X) neurons during singing transmits information about timing of the premotor sequence, but not auditory information, to the AFP.
Singing-related spiking activity of identified HVC neurons
Projection HVC neurons generate sparse and stereotypical patterns of bursts during singing. Bursts of both HVC(RA) and HVC(X) neurons are extremely tightly locked to the vocalization and frequently have rms jitter of <1 ms. The sequence of bursts of HVC projection neurons is one of the most temporally precise neural sequences found in nature to date. In contrast, HVC interneurons have tonic and more variable patterns of activity. We find that there is significant, albeit weak, synchrony between individual interneurons. There is also statistically significant synchrony between the populations of interneurons and HVC(X) neurons.
Neural activity in HVC and the syllable pattern of the song
We find that the population activity of HVC interneurons is modulated coherently with the syllable structure of the song. The population activity of interneurons leads the syllable pattern by
45 ms. It has been noted previously that the onset of HVC multi-unit activity precedes the onset of the contact calls (which are learned vocalizations similar in acoustic structures to song syllables) by
50 ms (McCasland 1987
). This time delay between the onset of neural activity and the onset of sound has been interpreted to correspond to a neural delay between activity in HVC and vocal output (Troyer and Doupe 2000
). Estimation of the neural delay between the activity in HVC and the vocal output is important because this delay contributes to the total delay between premotor activity and auditory feedback, which, in turn, is regarded as a fundamental difficulty in motor learning (Dave and Margoliash 2000
; Troyer and Doupe 2000
). Paired recordings in HVC and RA suggest that the bursts in RA are driven from HVC and the delay between bursts of HVC(RA) neurons and RA neurons is
5 ms (Hahnloser et al. 2002
). Direct stimulation of RA during a harmonic stack results in the transient vocal output perturbation that starts
1015 ms after the stimulation pulse (Fee et al. 2004
). These results suggest that the delay between bursts of HVC(RA) neurons and the vocal output is
1520 ms, substantially smaller than 50 ms previously inferred from the onsets of multi-unit HVC activity. The functional significance of the syllable-related modulation of neural activity in HVC is unclear. However, the modulation may reflect activity in HVC involved in preparing the premotor circuitry or vocal organ for production of the next song syllable. For example, it may reflect external input to HVC involved in interhemispherical synchronization of HVC activity (Schmidt 2003
).
Relation of bursts of HVC(X) neurons to acoustic structure of the song
The idea that HVC transmits auditory information to the AFP suggests that bursts of HVC(X) neurons may be related to specific acoustic elements in the song. Our recordings of firing patterns of HVC(X) neurons allow for testing of this hypothesis: auditory-related activity of HVC(X) neurons implies that multiple bursts of an HVC(X) neuron would follow similar acoustic elements. This idea is not supported by our observations: we find that multiple bursts of a given HVC(X) neuron are not associated with similar notes in the song. However, multiple bursts of a given HVC(X) neuron tend to precede similar elements of the syllable pattern (either sound or silence) by
40 ms. The lack of relation of bursts of an individual HVC(X) neuron to similar spectral structure in the song does not support the hypothesis that HVC(X) neurons transmit auditory information to the AFP.
Neural activity and vocal output evolve on uncorrelated time scales
We use the population activity of interneurons to assess the time scales of neural dynamics and to study the relation of timescales of neural dynamics to the characteristic times of acoustic output. We find that neural dynamics in HVC occurs on a fixed fast (
5 ms) time scales, whereas vocal dynamics occurs on variable time scales (from <10 to
100 ms) uncorrelated with the time scales of neural dynamics. This observation is consistent with the idea that HVC codes for temporal order in the song rather than for sound: whereas vocal output can remain essentially constant for
100 ms (e.g., during a harmonic stack), the premotor sequence in HVC evolves on a substantially shorter time scale.
Absence of auditory responses of HVC(X) neurons during singing
Lesion studies, together with the auditory responses exhibited in the AFP under anesthesia, have led to the hypothesis that the AFP processes auditory information to produce an error signal used to guide vocal learning. In a simple implementation of this model, HVC might provide auditory information to the AFP, which then induces a gradual modification of the motor program. In this study, we used DAF to directly test the idea that HVC(X) neurons are responsive to auditory signals or errors during singing. Several previous experiments have addressed the issue of auditory feedback in the AFP (Hessler and Doupe 1999a
; Leonardo 2004
; McCasland and Konishi 1981
), but none of them were carried out in juvenile birds, and none have directly examined the role of HVC(X) neurons in the processing of singing-related auditory signals.
We recorded from HVC(X) neurons during singing with and without distorted feedback in juvenile birds and found no difference in the activity of HVC(X) neurons with and without DAF. Absence of auditory sensitivity of HVC(X) neurons during singing is in surprising contrast with their robust responses to playback of the bird's own song (Mooney 2000
), which are highly sensitive to distortions of the auditory input (Theunissen and Doupe 1998
). During singing, HVC(X) neurons generate extremely stereotyped patterns of spikes, thus even subtle changes in their firing patterns should be easily detectable.
Absence of auditory responses of HVC(X) neurons would be expected if our distorted auditory feedback were not effective in inducing song plasticity. To rule out this possibility, we carried out recordings in young (7090 days posthatch) birds the songs of which can undergo systematic changes over a time course of several hours (Deregnaucourt et al. 2003
; Tchernichovski et al. 2001
). We also verified that our DAF protocol was efficient in causing plastic changes in the motor program. Despite the effectiveness of DAF in song degradation and young age of birds, we observe no changes in the firing patterns of HVC(X) neurons in the presence of DAF.
Absence of auditory responses of HVC(X) neurons during singing raises questions about the possible role of HVC input to the AFP in vocal learning. It is conceivable that auditory responses of HVC(X) neurons during singing are too weak to be detected in the sample of spike trains which can be obtained with our technique (typically 1520 song motifs). However, it is not clear why vocal error correction would be carried out with such small signals, i.e., with low signal to noise. Given the exquisite vocal learning capabilities of the songbird, it seems advantageous for the auditory and error-correction circuits to have high signal to noise and for these circuits to effect gradual changes in the vocal output. Another possibility is that HVC transmits auditory information about vocal performance to the AFP off-line, i.e., not during singing. There is some evidence from behavioral experiments suggesting that song learning may be happening on-line as well as during sleep (Deregnaucourt et al. 2005
; Tchernichovski et al. 2001
). Regardless of when plasticity in the motor circuit occurs, auditory information about the bird's vocal performance must be collected during singing. Lack of auditory responses to DAF in HVC might indicate, for example, that song-related auditory information may be stored in auditory areas upstream from HVC (e.g., Field L, NCM, cmHV, NIf) and transmitted through HVC to the AFP some time after singing, for example, during sleep (Dave and Margoliash 2000
).
Activity of HVC(X) neurons may constitute precise timing signal
A fundamentally different possibility is that HVC may not be involved in auditory processing, but is fundamentally involved in generating and transmitting timing information about the ongoing song. We have previously proposed that activity of RA-projecting HVC neurons may constitute a representation of temporal order in the vocal motor sequence (Hahnloser et al. 2002
). Given that HVC(X) neurons also generate precisely timed sparse sequences of bursts, similar to those observed in HVC(RA) neurons, it is compelling to imagine that HVC(X) neurons perform a similar function of transmitting timing information to the basal ganglia. Further support of this idea is the lack of auditory responsiveness of HVC(X) neurons during singing, as well as the lack of correlation between the firing of these neurons and spectral structure of the song. Although the precise function of the AFP in song learning is unknown, it seems likely that information about timing of the ongoing motor sequence during learning will play an essential role in that function.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
1 The online version of this article contains supplemental data. ![]()
Address for reprint requests and other correspondence: M. S. Fee, MIT, 46-5133A, 77 Massachusetts Ave., Cambridge, MA 02139
| REFERENCES |
|---|
|
|
|---|
Bottjer SW, Miesner EA, Arnold AP. Forebrain lesions disrupt development but not maintenance of song in passerine birds. Science 224: 901903, 1984.
Brainard MS, Doupe AJ. What songbirds teach us about learning. Nature 417: 351358, 2002.[CrossRef][Medline]
Brenowitz EA, Margoliash D, Nordeen KW. An introduction to birdsong and the avian song system. J Neurobiol 33: 495500, 1997.[CrossRef][Web of Science][Medline]
Dave AS, Margoliash D. Song replay during sleep and computational rules for sensorimotor vocal learning. Science 290: 812816, 2000.
Deregnaucourt S, Maul K, Feher O, Mitra PP, Tchernichovski O. Circadian fluctuations in syllable structure during song development (Abstract). SFN 2003.
Deregnaucourt S, Mitra PP, Feher O, Pytte C, Tchernichovski O. How sleep affects the developmental learning of bird song. Nature 433: 710716, 2005.[CrossRef][Medline]
Doupe AJ, Konishi M. Song-selective auditory circuits in the vocal control system of the zebra finch. Proc Natl Acad Sci USA 88: 1133911343, 1991.
Fee MS, Leonardo A. Miniature motorized microdrive and commutator system for chronic neural recordings in small animals. J Neurosci Methods 112: 8394, 2001.[CrossRef][Web of Science][Medline]
Fee MS, Kozhevnikov AA, Hahnloser RHR. Neural mechanisms of sequence generation in the songbird. Ann NY Acad Sci 1016: 153170, 2004.[CrossRef][Web of Science][Medline]
Glaze CM, Troyer TW. Temporal structure in zebra finch song: implications for motor coding. J Neurosci 26: 9911005, 2006.
Golay MJE. Complementary series. IRE Trans Inform Theory II-7: 8287, 1961.[CrossRef]
Hahnloser RHR, Kozhevnikov AA, Fee MS. An ultra-sparse code underlies the generation of neural sequences in a songbird. Nature 419: 6570, 2002.[CrossRef][Medline]
Hahnloser RHR, Kozhevnikov AA, Fee MS. Sleep-related neural activity in a premotor and a basal-ganglia pathway of the songbird. J Neurophysiol 96: 794812, 2006.
Hessler N, Doupe A. Singing-related neural activity in a dorsal forebrain-basal ganglia circuit of adult zebra finches. J Neurosci 19: 1046110481, 1999.
Jarvis MR, Mitra PP. Sampling properties of the spectrum and coherency of sequences of action potentials. Neural Comput 13: 717749, 2001.[CrossRef][Web of Science][Medline]
Konishi M. The role of auditory feedback in the control of vocalizations in the white-crowned sparrow. Z Tierpsychol 22: 770783, 1965.[Medline]
Leonardo A. Neural dynamics underlying complex behavior in a songbird (PhD thesis). In: Computation and Neural Systems Program. Pasadena: California Institute of Technology. 2002.
Leonardo A. Experimental test of the birdsong error-correction model. Proc Nat Acad Sci USA 101: 1693516940, 2004.
Leonardo A, Fee MS. Ensemble coding of vocal control in birdsong. J Neurosci 25: 652661, 2005.
Leonardo A, Konishi M. Decrystallization of adult birdsong by perturbation of auditory feedback. Nature 399: 466470, 1999.[CrossRef][Medline]
Lombardino A, Nottebohm F. Age at deafening affects the stability of learned song in adult male zebra finches. J Neurosci 20: 50545064, 2000.[Medline]
Luo M, Perkel DJ. Long-range GABAergic projection in a circuit essential for vocal learning. J Comp Neurol 403: 6884, 1999a.[CrossRef][Web of Science][Medline]
Luo M, Perkel DJ. A GABAergic, strongly inhibitory projection to a thalamic nucleus in the zebra finch song system. J Neurosci 19: 67006711, 1999b.
Margoliash D. Preference for autogenous song by auditory neurons in a song system nucleus in the white-crowned sparrow. J Neurosci 6: 16431661, 1986.[Abstract]
McCasland JS. Neuronal control of bird song production. J Neurosci 7: 2339, 1987.[Abstract]
McCasland J, Konishi M. Interaction between auditory and motor activities in an avian song control nucleus. Proc Natl Acad Sci USA 78: 78157819, 1981.
Mooney R. Different subthreshold mechanisms underlie song selectivity in identified HVc neurons of the zebra finch. J Neurosci 201: 54205436, 2000.
Mooney R, Prather JF. The HVC microcircuit: the synaptic basis for interactions between song motor and vocal plasticity pathways. J Neurosci 25: 19521964, 2005.
Nottebohm F, Stokes TM, Leonard CM. Central control of song in the canary, Serinus canarius. J Comp Neurol 165: 457486, 1976.[CrossRef]
Okuhata S, Saito N. Synaptic connections of thalamo-cerebral vocal nuclei of the canary. Brain Res Bull 18: 3544, 1987.[CrossRef][Web of Science][Medline]
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 11: 28962913, 1991.[Abstract]
Schmidt M. Pattern of interhemispheric synchronization in HVC during singing correlates with key transitions in the song pattern. J Neurophysiol 90: 39313949, 2003.
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 53: 5163, 1990.[CrossRef][Web of Science][Medline]
Tchernichovski O, Mitra P, Lints T, Nottebohm F. Dynamics of the vocal imitation process: how a zebra finch learns its song. Science 291: 25642569, 2001.
Tchernichovski O, Nottebohm F, Ho CE, Pesaran B, Mitra PP. A procedure for an automated measurement of song similarity. Anim Behav 59: 11671176, 2000.[CrossRef][Web of Science][Medline]
Theunissen FE, Doupe AJ. Temporal and spectral sensitivity of complex auditory neurons in the nucleus HVc of male zebra finches. J Neurosci 20: 23152331, 1998.
Thompson DJ. Spectrum Estimation and Harmonic Analysis. Proc IEEE 70: 10551096, 1982.[CrossRef]
Troyer T, Doupe A. An associational model of birdsong sensorimotor learning. I. Efference copy and the learning of song syllables. J Neurophysiol 84: 12041223, 2000.
Vicario DS. Organization of the zebra finch song control system: II. Functional organization of outputs from nucleus Robustus archistriatalis. J Comp Neurol 309: 486494, 1991.[CrossRef][Web of Science][Medline]
Vicario DS, Nottebohm F. Organization of the zebra finch song control system. I. Representation of syringeal muscles in the hypoglossal nucleus. J Comp Neurol 271: 346354, 1988.[CrossRef][Web of Science][Medline]
Vicario DS, Yohay KH. Song-selective auditory input to a forebrain vocal control nucleus in the zebra finch. J Neurobiol 24: 488505, 1993.[CrossRef][Web of Science][Medline]
Vu ET, Schmidt MF, Mazurek ME. Interhemispheric coordination of premotor neural activity during singing in adult zebra finches. J Neurosci 18: 90889098, 1998.
Williams H, Nottebohm F. Auditory responses in avian vocal motor neurons: a motor theory for song perception in birds. Science 229: 279282, 1985.
Yu AC, Margoliash D. Temporal hierarchical control of singing in birds. Science 273: 18711875, 1996.
This article has been cited by other articles:
![]() |
J. Meitzen, A. L. Weaver, E. A. Brenowitz, and D. J. Perkel Plastic and Stable Electrophysiological Properties of Adult Avian Forebrain Song-Control Neurons across Changing Breeding Conditions J. Neurosci., May 20, 2009; 29(20): 6558 - 6567. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. H. Kao, B. D. Wright, and A. J. Doupe Neurons in a Forebrain Nucleus Required for Vocal Plasticity Rapidly Switch between Precise Firing and Variable Bursting Depending on Social Context J. Neurosci., December 3, 2008; 28(49): 13232 - 13247. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. R. Griff, M. Mafhouz, A. Perrut, and M. A. Chaput Comparison of Identified Mitral and Tufted Cells in Freely Breathing Rats: I. Conduction Velocity and Spontaneous Activity Chem Senses, November 1, 2008; 33(9): 779 - 792. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. T. Sakata and M. S. Brainard Online Contributions of Auditory Feedback to Neural Activity in Avian Song Control Circuitry J. Neurosci., October 29, 2008; 28(44): 11378 - 11390. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. F. Roberts, M. E. Klein, M. F. Kubke, J. M. Wild, and R. Mooney Telencephalic Neurons Monosynaptically Link Brainstem and Forebrain Premotor Networks Necessary for Song J. Neurosci., March 26, 2008; 28(13): 3479 - 3489. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. A. Thompson, W. Wu, R. Bertram, and F. Johnson Auditory-Dependent Vocal Recovery in Adult Male Zebra Finches Is Facilitated by Lesion of a Forebrain Pathway That Includes the Basal Ganglia J. Neurosci., November 7, 2007; 27(45): 12308 - 12320. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Meitzen, I. T. Moore, K. Lent, E. A. Brenowitz, and D. J. Perkel Steroid Hormones Act Transsynaptically within the Forebrain to Regulate Neuronal Phenotype and Song Stereotypy J. Neurosci., October 31, 2007; 27(44): 12045 - 12057. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. L. Person and D. J. Perkel Pallidal Neuron Activity Increases during Sensory Relay through Thalamus in a Songbird Circuit Essential for Learning J. Neurosci., August 8, 2007; 27(32): 8687 - 8698. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |