henhaobio
Sign in

bioRxivNeurosciencePreprint

Neural-vocal phase coupling reveals structured timing in birdsong production

Fiamma L. Leites, B. R. R. BoarettoiD, Cristina MasolleriD, Ana AmadoriD

PDFDOI

Abstract

Understanding how neural population activity is temporally coordinated with behavior remains a central challenge in neuroscience. Songbirds provide a powerful model system for addressing this question because learned vocal production requires precise coordination among neural dynamics, temporally structured motor output, and auditory feedback. However, quantifying neural-vocal interactions is challenging because both neural and acoustic signals are rhythmic, noisy, and highly nonstationary. Here, we investigate neural-vocal coordination during spontaneous canary singing using simultaneous recordings of neural population activity in a forebrain region of the song system and vocal behavior. Using a phase-resolved cross-correlation framework combined with surrogate-based statistical validation, we quantify neural-vocal interactions in short and highly variable song segments. Our analysis reveals that neural-vocal interactions are organized into distinct temporal regimes comprising positive, near-zero, and negative lags, consistent with neural activity preceding, accompanying, or following vocal output. The coexistence of these regimes is consistent with the integrative role of the recorded region, which receives auditory input, contributes to premotor control, and participates in the neural circuitry supporting song learning and the ongoing maintenance of adult song. We further find that correlated and anticorrelated interactions coexist throughout singing, with anticorrelated interactions consistently concentrated around near-zero lags. These anticorrelations identify periods in which decreases in neural population activity are closely aligned with sound production, revealing biologically relevant information that is obscured by analyses performed over complete song renditions. Together, these results uncover a robust temporal structure linking neural population activity to vocal behavior and provide a broadly applicable framework for extracting transient neural-behavioral interactions from complex biological signals. Author summaryModern neuroscience can simultaneously record the activity of large neural populations, yet extracting meaningful relationships between neural activity and behavior remains challenging because natural behaviors are highly variable. Songbirds provide a unique opportunity to study this problem because their learned vocalizations share key features with human speech while remaining experimentally accessible. Here, we analyzed simultaneous recordings of neural population activity and vocal behavior in freely singing canaries. Instead of averaging neural activity across entire songs, we examined brief time windows and combined local correlation analysis with statistical tests to identify reliable interactions between the brain and behavior. We found that these interactions switch among several preferred timing patterns: neural activity can precede vocal output, occur at nearly the same time, or follow acoustic events. This diversity is consistent with the combined roles of the recorded brain region. Our work provides an intuitive and broadly applicable framework for uncovering transient neural-behavioral interactions that remain hidden by traditional time-averaged approaches. Because it requires only simultaneous recordings of two time series and makes minimal assumptions about their dynamics, the framework can be applied to many biological systems involving complex temporal signals.

In-depth summary

Claude Opus 5.5 · xhigh effort · 8 Nov 2026 · from the full text of version 1 (22 Jul 2026) · CC BY 4.0

In three singing canaries, sliding-window cross-correlation with surrogate testing shows HVC multi-unit activity and song envelope couple at a few preferred lags (positive, near-zero, negative), correlated or anticorrelated.

Rhythmic, nonstationary neural and vocal signals often give whole-song cross-correlations with several comparable peaks, and the paper notes that Hilbert phase locking, coherence and Granger causality rely on narrowband or stationary signals. The framework assigns a lag only in short windows where a peak is prominent, beats surrogate thresholds and appears symmetrically in both cross-correlation directions, letting transient coupling regimes be tracked through a song.

Key results

  1. SuggestedWithin each bird, significant neural-vocal lags (τ*) are not spread continuously. They fall into one to three preferred lag classes per bird and coupling sign, and where several classes occur they are separated by low-density gaps.

    Gaussian KDE (bandwidth factor 0.65; classes kept only if they hold at least 10% of points). Class medians: Bird 1 correlations ≈ 23 ms, anticorrelations ≈ −65, ≈ 0 and ≈ 65 ms. Bird 2 correlations ≈ −8, ≈ 42 and ≈ 114 ms, anticorrelations ≈ 15 ms. Bird 3 correlations ≈ −2 and ≈ 54 ms, anticorrelations ≈ 7.5 ms. Medians and IQRs are in Table S1. · Neural-vocal interactions cluster into discrete…; Class detection using kernel density estimation; Supporting information

  2. SuggestedSignificant coupling between HVC multi-unit activity and the song envelope appears only in limited stretches of a song. Sign and lag switch between stretches, but similar C* and τ* values persist over consecutive windows, forming temporally localized regimes.

    Sliding 1200 ms windows (N = 1200 at 1 kS/s), τmax = 200 ms. A lag was assigned only for peaks with prominence above 0.1 and |C*| > 0.5 that exceeded thresholds from 100 IAAFT surrogates and 1000 random window pairs, with CXY/CYX peak lags within 10 ms of each other. Time courses shown for Bird 3 (sample 837), Bird 2 (sample 240) and Bird 1 (sample 356). · Phase-resolved analyses reveal temporal…; Cross-correlation analysis; Statistical significance analysis; Figure legends

  3. ShownCross-correlation computed over a whole song frequently gives several peaks of similar height at positive and negative lags, so a unique neural-vocal lag is often ambiguous and a single whole-song lag cannot adequately describe the relationship.

    In the three representative whole-song CXY(τ) and CYX(τ) examples (τmax = 1000 ms; Birds 1, 2 and 3, Fig. 2), the functions show several peaks of comparable amplitude rather than one dominant maximum. · Neural population activity exhibits structured…; Cross-correlation analysis; Figure legends

  4. SuggestedAnticorrelated coupling (C* < 0) recurs as a near-zero-lag class in all three birds. Example segments cover every combination of coupling sign (positive or negative) and lag direction (positive, near-zero, negative).

    The authors place an anticorrelation class near τ* ≈ 0 in every bird: medians ≈ 0 ms (Bird 1), ≈ 15 ms (Bird 2) and ≈ 7.5 ms (Bird 3), overlapping across birds. Fig. 6 shows example segments with C* > 0.5 and with C* < −0.5, each at τ* > 0, τ* ≈ 0 and τ* < 0. · Neural-vocal interactions cluster into discrete…; Distinct phase regimes are associated with…; Discussion; Supporting information

  5. SuggestedTwo birds have correlated classes at roughly 40-60 ms with HVC activity leading the envelope. The authors see this as compatible with HVC-to-vocal-output delays through downstream motor pathways. Infrequent lags beyond 100 ms occur during slowly repeated or long syllables.

    Correlation class medians are ≈ 42 ms (Bird 2) and ≈ 54 ms (Bird 3). The lag distributions have tails with |τ*| > 100 ms, seen in Fig. 5 and Fig. S2. The motor-pathway reading comes from the Discussion; no causal test was done. · Neural-vocal interactions cluster into discrete…; Discussion; Supporting information

Methods, data and reproducibility

Tetrode recordings from HVC of singing canaries were reduced to multi-unit activity (MUA) and cross-correlated with the root-mean-square song envelope (ENV) in sliding 1200 ms windows. Only prominent, surrogate-significant peaks consistent between CXY and CYX were kept. Lags were then grouped into classes by kernel density estimation and tested with linear mixed-effects models.

  • Simultaneous extracellular HVC recordings (tetrode MUA) and song audio during spontaneous singing; whole songs analysed, including baseline silence before and after · canary (Serinus canaria) · n = three adult male canaries; four MUA traces per song
Methods & tools
MATLAB, IAAFT surrogates (Iterated Amplitude Adjusted Fourier Transform; unnamed Python library), Intan RHD2000 USB interface, RHD2132 16-channel amplifier, custom tungsten tetrode arrays (California Fine Wire), manual microdrive, MAX4466 amplifier board with electret microphone, third-order Butterworth high-pass filter, Gaussian kernel density estimation, linear mixed-effects models with likelihood ratio tests
Reagents & environment
Preprocessing in MATLAB; IAAFT surrogates generated with a Python library that is not named; no software versions given.
Reproducibility gaps
The text has no code or data availability statement and no links. The number of songs analysed per bird is not stated.
Limitations (4)
  • The analysis is correlational and does not establish causal relationships between HVC activity and vocal output. · stated by the authors
  • Window size, maximum lag, prominence and correlation thresholds depend on the dataset and must be tuned to each signal's temporal structure, noise and sampling; no single setting is optimal across datasets. · stated by the authors
  • Only three birds, and the text does not say how many songs or windows went into each lag distribution. Lag classes come from KDE peak picking (bandwidth factor 0.65), with no test against a continuous or unimodal null. · noted in this summary
  • Several settings (MUA smoothing, the 10 ms symmetry tolerance) were chosen by inspection. The linear mixed-effects comparisons of correlated vs anticorrelated lags are described in Methods, but their results are not reported in the text. · noted in this summary

Read this if you want to detect transient, lag-specific coupling between a neural population signal and a rhythmic, nonstationary behavioural signal using sliding-window cross-correlation with surrogate tests.

Citations
3

Reproducibility ledger

Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.

Sign in to file a reproducibility report.