A generalized linear model for estimating spectrotemporal receptive fields from responses to natural sounds

In the auditory system, the stimulus-response properties of single neurons are often described in terms of the spectrotemporal receptive field (STRF), a linear kernel relating the spectrogram of the sound stimulus to the instantaneous firing rate of the neuron. Several algorithms have been used to estimate STRFs from responses to natural stimuli; these algorithms differ in their functional models, cost functions, and regularization methods. Here, we characterize the stimulus-response function of auditory neurons using a generalized linear model (GLM). In this model, each cell's input is described by: 1) a stimulus filter (STRF); and 2) a post-spike filter, which captures dependencies on the neuron's spiking history. The output of the model is given by a series of spike trains rather than instantaneous firing rate, allowing the prediction of spike train responses to novel stimuli. We fit the model by maximum penalized likelihood to the spiking activity of zebra finch auditory midbrain neurons in response to conspecific vocalizations (songs) and modulation limited (ml) noise. We compare this model to normalized reverse correlation (NRC), the traditional method for STRF estimation, in terms of predictive power and the basic tuning properties of the estimated STRFs. We find that a GLM with a sparse prior predicts novel responses to both stimulus classes significantly better than NRC. Importantly, we find that STRFs from the two models derived from the same responses can differ substantially and that GLM STRFs are more consistent between stimulus classes than NRC STRFs. These results suggest that a GLM with a sparse prior provides a more accurate characterization of spectrotemporal tuning than does the NRC method when responses to complex sounds are studied in these neurons.     

Species differences in auditory processing dynamics in songbird auditory telencephalon

The caudomedial nidopallium (NCM) of songbirds is a telencephalic area involved in the auditory processing and memorization of complex vocal communication signals. We used pure tone stimuli and multiunit electrophysiological recordings in awake birds to investigate whether the basic properties of song-responding circuits in NCM differ between canaries and zebra finches, two species whose songs are markedly different in their spectral and temporal organization. We found that the responses in zebra finch NCM are characterized by broad tuning and sustained responses that may facilitate the integration of zebra finch song syllables and call notes that are of long duration and have a broad harmonic structure. In contrast, we found that the responses in canary NCM show narrower tuning and less sustained responses over the time periods analyzed. These characteristics may contribute to enhanced processing of the narrow-band whistles, rapid trills, and steep frequency modulations that are prominent features of canary song. These species differences are much less pronounced in field L2, the direct thalamorecipient region that represents a preceding station in the central avian auditory pathway. NCM responses did not differ across sexes of either species, but field L2 did show wider tuning in zebra finch females relative to males. In sum, species differences in the response properties of NCM likely reflect selectivity for the acoustic elements of each species' vocal repertoire.     

Developmental experience alters information coding in auditory midbrain and forebrain neurons

In songbirds, species identity and developmental experience shape vocal behavior and behavioral responses to vocalizations. The interaction of species identity and developmental experience may also shape the coding properties of sensory neurons. We tested whether responses of auditory midbrain and forebrain neurons to songs differed between species and between groups of conspecific birds with different developmental exposure to song. We also compared responses of individual neurons to conspecific and heterospecific songs. Zebra and Bengalese finches that were raised and tutored by conspecific birds, and zebra finches that were cross‐tutored by Bengalese finches were studied. Single‐unit responses to zebra and Bengalese finch songs were recorded and analyzed by calculating mutual information (MI), response reliability, mean spike rate, fluctuations in time‐varying spike rate, distributions of time‐varying spike rates, and neural discrimination of individual songs. MI quantifies a response's capacity to encode information about a stimulus. In midbrain and forebrain neurons, MI was significantly higher in normal zebra finch neurons than in Bengalese finch and cross‐tutored zebra finch neurons, but not between Bengalese finch and cross‐tutored zebra finch neurons. Information rate differences were largely due to spike rate differences. MI did not differ between responses to conspecific and heterospecific songs. Therefore, neurons from normal zebra finches encoded more information about songs than did neurons from other birds, but conspecific and heterospecific songs were encoded equally. Neural discrimination of songs and MI were highly correlated. Results demonstrate that developmental exposure to vocalizations shapes the information coding properties of songbird auditory neurons. © 2009 Wiley Periodicals, Inc. Develop Neurobiol 70: 235–252, 2010.     

Competing sound sources reveal spatial effects in cortical processing

Why is spatial tuning in auditory cortex weak, even though location is important to object recognition in natural settings? This question continues to vex neuroscientists focused on linking physiological results to auditory perception. Here we show that the spatial locations of simultaneous, competing sound sources dramatically influence how well neural spike trains recorded from the zebra finch field L (an analog of mammalian primary auditory cortex) encode source identity. We find that the location of a birdsong played in quiet has little effect on the fidelity of the neural encoding of the song. However, when the song is presented along with a masker, spatial effects are pronounced. For each spatial configuration, a subset of neurons encodes song identity more robustly than others. As a result, competing sources from different locations dominate responses of different neural subpopulations, helping to separate neural responses into independent representations. These results help elucidate how cortical processing exploits spatial information to provide a substrate for selective spatial auditory attention.     

Mechanisms for temporal tuning and filtering by postsynaptic signaling pathways

Networks of signaling pathways perform complex temporal decoding functions in diverse biological systems, including the synapse, development, and bacterial chemotaxis. This paper examines temporal filtering and tuning properties of synaptic signaling pathways as a possible substrate for emergent temporal decoding. A mass action kinetic model of 16 synaptic signaling pathways was used to dissect out the contribution of these pathways in linear cascades and when coupled to form a network. The model predicts two primary mechanisms of temporal tuning of pathways: a weighted summation of responses of pathways with different timings and the presence of biochemical feedback loop(s) with emergent dynamics. Regulatory inputs act differently on these two tuning mechanisms. In the first case, regulators act like a gain-control on pathways with different intrinsic tuning. In the case of feedback loops, the temporal properties of the loop itself are changed. These basic tuning mechanisms may underlie specialized temporal tuning functions in more complex signaling systems in biology.     

Studying spike trains using a van Rossum metric with a synapse-like filter

Spike trains are unreliable. For example, in the primary sensory areas, spike patterns and precise spike times will vary between responses to the same stimulus. Nonetheless, information about sensory inputs is communicated in the form of spike trains. A challenge in understanding spike trains is to assess the significance of individual spikes in encoding information. One approach is to define a spike train metric, allowing a distance to be calculated between pairs of spike trains. In a good metric, this distance will depend on the information the spike trains encode. This method has been used previously to calculate the timescale over which the precision of spike times is significant. Here, a new metric is constructed based on a simple model of synaptic conductances which includes binding site depletion. Including binding site depletion in the metric means that a given individual spike has a smaller effect on the distance if it occurs soon after other spikes. The metric proves effective at classifying neuronal responses by stimuli in the sample data set of electro-physiological recordings from the primary auditory area of the zebra finch fore-brain. This shows that this is an effective metric for these spike trains suggesting that in these spike trains the significance of a spike is modulated by its proximity to previous spikes. This modulation is a putative information-coding property of spike trains.     

Response properties of the auditory telencephalon in songbirds change with recent experience and season

The caudomedial nidopallium (NCM) is a telencephalic auditory area that is selectively activated by conspecific vocalizations in zebra finches and canaries. We recently demonstrated that temporal and spectral dynamics of auditory tuning in NCM differ between these species [1]. In order to determine whether these differences reflect recent experience, we exposed separate groups of each species and sex to different housing conditions. Adult birds were housed either in an aviary with conspecifics (NORM), with heterospecifics (canary subjects in a zebra finch aviary, and vice versa: (CROSS)), or in isolation (ISO) for 9 days prior to testing. We then recorded extracellular multi-unit electrophysiological responses to simple pure tone stimuli (250-5000 Hz) in awake birds from each group and analyzed auditory tuning width using methods from our earlier studies. Relative to NORM birds, tuning was narrower in CROSS birds, and wider in ISO birds. The trend was greater in canaries, especially females. The date of recording was also included as a covariate in ANCOVAs that analyzed a larger set of the canary data, including data from birds tested outside of the breeding season, and treated housing condition and sex as independent variables. These tests show that tuning width was narrower early in the year and broader later. This effect was most pronounced in CROSS males. The degree of the short-term neural plasticity described here differs across sexes and species, and may reflect differences in NCM's anatomical and functional organization related to species differences in song characteristics, adult plasticity and/or social factors. More generally, NCM tuning is labile and may be modulated by recent experience to reflect the auditory processing required for behavioral adaptation to the current acoustic, social or seasonal context.     

Disentangling sub-millisecond processes within an auditory transduction chain

Every sensation begins with the conversion of a sensory stimulus into the response of a receptor neuron. Typically, this involves a sequence of multiple biophysical processes that cannot all be monitored directly. In this work, we present an approach that is based on analyzing different stimuli that cause the same final output, here defined as the probability of the receptor neuron to fire a single action potential. Comparing such iso-response stimuli within the framework of nonlinear cascade models allows us to extract the characteristics of individual signal-processing steps with a temporal resolution much finer than the trial-to-trial variability of the measured output spike times. Applied to insect auditory receptor cells, the technique reveals the sub-millisecond dynamics of the eardrum vibration and of the electrical potential and yields a quantitative four-step cascade model. The model accounts for the tuning properties of this class of neurons and explains their high temporal resolution under natural stimulation. Owing to its simplicity and generality, the presented method is readily applicable to other nonlinear cascades and a large variety of signal-processing systems.     

Spectrotemporal Response Properties of Core Auditory Cortex Neurons in Awake Monkey

So far, most studies of core auditory cortex (AC) have characterized the spectral and temporal tuning properties of cells in non-awake, anesthetized preparations. As experiments in awake animals are scarce, we here used dynamic spectral-temporal broadband ripples to study the properties of the spectrotemporal receptive fields (STRFs) of AC cells in awake monkeys. We show that AC neurons were typically most sensitive to low ripple densities (spectral) and low velocities (temporal), and that most cells were not selective for a particular spectrotemporal sweep direction. A substantial proportion of neurons preferred amplitude-modulated sounds (at zero ripple density) to dynamic ripples (at non-zero densities). The vast majority (>93%) of modulation transfer functions were separable with respect to spectral and temporal modulations, indicating that time and spectrum are independently processed in AC neurons. We also analyzed the linear predictability of AC responses to natural vocalizations on the basis of the STRF. We discuss our findings in the light of results obtained from the monkey midbrain inferior colliculus by comparing the spectrotemporal tuning properties and linear predictability of these two important auditory stages.     

M-channels modulate the intrinsic excitability and synaptic responses of layer 2/3 pyramidal neurons in auditory cortex

Neurons in the auditory cortex are believed to utilize temporal patterns of neural activity to accurately process auditory information but the intrinsic neuronal mechanism underlying the control of auditory neural activity is not known. The slowly activating, persistent K⁺ channel, also called M-channel that belongs to the Kv7 family, is already known to be important in regulating subthreshold neural excitability and synaptic summation in neocortical and hippocampal pyramidal neurons. However, its functional role in the primary auditory cortex (A1) has never been characterized. In this study, we investigated the roles of M-channels on neuronal excitability, short-term plasticity, and synaptic summation of A1 layer 2/3 regular spiking pyramidal neurons with whole-cell current-clamp recordings in vitro. We found that blocking M-channels with a selective M-channel blocker, XE991, significantly increased neural excitability of A1 layer 2/3 pyramidal neurons. Furthermore, M-channels controled synaptic responses of intralaminar-evoked excitatory postsynaptic potentials (EPSPs); XE991 significantly increased EPSP amplitude, decreased the rate of short-term depression, and increased the synaptic summation. These results suggest that M-channels are involved in controlling spike output patterns and synaptic responses of A1 layer 2/3 pyramidal neurons, which would have important implications in auditory information processing.     

Nucleotide Variation, Linkage Disequilibrium and Founder-Facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (K_ST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.     

The selectivity of canary HVC neurons for the Bird’s Own Song: Rate coding, temporal coding, or both?

The neuronal selectivity observed in the avian song system for the Bird’s Own Song progressively emerged as an extraordinary fruitful model to investigate the neural code underlying the recognition of complex stimuli and the occurrence of learned behaviors. In adult zebra finch, neurons from the HVC (used as a proper name) show very selective auditory responses, firing more to presentation of the Bird’s Own Song (BOS) than to reverse BOS or other conspecific songs. However, as adult zebra finches always produce the same stereotyped song, the presence of such highly selective neurons in birds with larger repertoire still remains an open question. Data presented here show that neurons selective for the BOS can be found in adult canary, a seasonal breeding bird which display a large repertoire. More precisely, we found that a large proportion of neurons (29/36) exhibits higher responses to presentation of the forward than to the reverse BOS, and that 22% of the cells were identified as selective on the basis of the d′ value. For a cell that was extensively studied, we evaluated to what extent temporal stimulus-related structure predicts the acoustic stimulus using linear or non-linear artificial neural networks (ANN). These analyses indicated that the temporal structure contained in spike trains characterizes more accurately the stimulus than the firing rate. The limitations of applying ANN analyses to electrophysiological data are discussed and potential solutions to increase the confidence in these analysis are proposed.     

The selectivity of canary HVC neurons for the Bird's Own Song: rate coding, temporal coding, or both?

The neuronal selectivity observed in the avian song system for the Bird's Own Song progressively emerged as an extraordinary fruitful model to investigate the neural code underlying the recognition of complex stimuli and the occurrence of learned behaviors. In adult zebra finch, neurons from the HVC (used as a proper name) show very selective auditory responses, firing more to presentation of the Bird's Own Song (BOS) than to reverse BOS or other conspecific songs. However, as adult zebra finches always produce the same stereotyped song, the presence of such highly selective neurons in birds with larger repertoire still remains an open question. Data presented here show that neurons selective for the BOS can be found in adult canary, a seasonal breeding bird which display a large repertoire. More precisely, we found that a large proportion of neurons (29/36) exhibits higher responses to presentation of the forward than to the reverse BOS, and that 22% of the cells were identified as selective on the basis of the d' value. For a cell that was extensively studied, we evaluated to what extent temporal stimulus-related structure predicts the acoustic stimulus using linear or non-linear artificial neural networks (ANN). These analyses indicated that the temporal structure contained in spike trains characterizes more accurately the stimulus than the firing rate. The limitations of applying ANN analyses to electrophysiological data are discussed and potential solutions to increase the confidence in these analysis are proposed.     

Intrinsic and synaptic properties of neurons in the vocal-control nucleus lMAN from in vitro slice preparations of juvenile and adult zebra finches

A common theme of diverse neural systems is that circuits that are important for initial acquisition of learning do not necessarily serve as a substrate for the long-term storage of that memory. The neural basis of vocal learning in songbirds provides an example of this phenomenon, since a circuit that is necessary for vocal production during initial stages of vocal development apparently plays no subsequent role in controlling learned vocalizations. This striking functional change suggests the possibility of marked physiological changes in synaptic transmission within this circuit. We therefore examined intrinsic and synaptic properties of neurons in the cortical nucleus lMAN (lateral magnocellular nucleus of the anterior neostriatum), which forms part of this developmentally regulated circuit, in an in vitro preparation of the zebra finch forebrain. Although both functional and morphological characteristics of these neurons change substantially during vocal development, we did not observe widespread, substantive changes in the electrophysiological characteristics of juvenile versus adult lMAN neurons examined in vitro. Overall, both the intrinsic properties and synaptic responses of lMAN neurons were similar in slices from juvenile birds (at ages when lesions of lMAN disrupt vocal production) and in slices from adult birds (when lMAN lesions have no effect on song production). However, one intrinsic property that did vary between juvenile and adult cells was spike duration, which was longer in juvenile cells, suggesting the potential for activation of second-messenger cascades and/or enhanced synaptic transmission onto target cells of lMAN neurons. The pattern of synaptic response observed in both juvenile and adult cells suggests that lMAN projection neurons receive direct excitatory afferent inputs, as well as disynaptic inhibitory inputs from interneurons within lMAN. Activation of inhibitory interneurons rapidly curtails the excitatory response seen in projection neurons. This inhibition was abolished by bicuculline, indicating that the inhibitory interneurons normally exert their postsynaptic response via GABA_A receptors on projection neurons. The inhibitory response could also be blocked by CNQX (6-cyano-7-nitroquinoxaline-2,3-dione), suggesting that the activation of inhibitory interneurons within lMAN may be governed primarily by AMPA receptors. © 1998 John Wiley & Sons, Inc. J Neurobiol 37: 642–658, 1998     

Robustness of neuronal tuning to binaural sound localization cues against age-related loss of inhibitory synaptic inputs

Sound localization relies on minute differences in the timing and intensity of sound arriving at both ears. Neurons of the lateral superior olive (LSO) in the brainstem process these interaural disparities by precisely detecting excitatory and inhibitory synaptic inputs. Aging generally induces selective loss of inhibitory synaptic transmission along the entire auditory pathways, including the reduction of inhibitory afferents to LSO. Electrophysiological recordings in animals, however, reported only minor functional changes in aged LSO. The perplexing discrepancy between anatomical and physiological observations suggests a role for activity-dependent plasticity that would help neurons retain their binaural tuning function despite loss of inhibitory inputs. To explore this hypothesis, we use a computational model of LSO to investigate mechanisms underlying the observed functional robustness against age-related loss of inhibitory inputs. The LSO model is an integrate-and-fire type enhanced with a small amount of low-voltage activated potassium conductance and driven with (in)homogeneous Poissonian inputs. Without synaptic input loss, model spike rates varied smoothly with interaural time and level differences, replicating empirical tuning properties of LSO. By reducing the number of inhibitory afferents to mimic age-related loss of inhibition, overall spike rates increased, which negatively impacted binaural tuning performance, measured as modulation depth and neuronal discriminability. To simulate a recovery process compensating for the loss of inhibitory fibers, the strength of remaining inhibitory inputs was increased. By this modification, effects of inhibition loss on binaural tuning were considerably weakened, leading to an improvement of functional performance. These neuron-level observations were further confirmed by population modeling, in which binaural tuning properties of multiple LSO neurons were varied according to empirical measurements. These results demonstrate the plausibility that homeostatic plasticity could effectively counteract known age-dependent loss of inhibitory fibers in LSO and suggest that behavioral degradation of sound localization might originate from changes occurring more centrally.     

Lateral sharpening of cortical frequency tuning by approximately balanced inhibition

Cortical inhibition plays an important role in shaping neuronal processing. The underlying synaptic mechanisms remain controversial. Here, in vivo whole-cell recordings from neurons in the rat primary auditory cortex revealed that the frequency tuning curve of inhibitory input was broader than that of excitatory input. This results in relatively stronger inhibition in frequency domains flanking the preferred frequencies of the cell and a significant sharpening of the frequency tuning of membrane responses. The less selective inhibition can be attributed to a broader bandwidth and lower threshold of spike tonal receptive field of fast-spike inhibitory neurons than nearby excitatory neurons, although both types of neurons receive similar ranges of excitatory input and are organized into the same tonotopic map. Thus, the balance between excitation and inhibition is only approximate, and intracortical inhibition with high sensitivity and low selectivity can laterally sharpen the frequency tuning of neurons, ensuring their highly selective representation.     

Comparative study of long-term reduction of high-frequency song syllables in deafened zebra and Bengalese finches: implications for the role of auditory feedback in song maintenance

Auditory feedback is necessary for adult song maintenance in both oscines and psittacines. Although belonging to phylogenically separated orders, deafened adult oscine Bengalese finches and psittacine budgerigars exhibit similarities in certain aspects of song changes. An interesting question is whether these birds share common mechanisms for song maintenance. Therefore, it is important to compare the effects of deafening on adult song patterns among and within orders. Although zebra and Bengalese finches are closely related oscine species, few studies have performed direct, long-term, quantitative comparisons of their songs after deafening because suitable song characteristics have not been identified. Based on our previous findings for Bengalese finch songs, we analyzed zebra finch songs over 9 months after deafening, focusing on changes in the number of syllables categorized according to fundamental frequencies. Deafened zebra finches demonstrated a gradual but significant decrease in high-frequency syllables and a tendency to increase low-frequency syllables, similar to deafened Bengalese finches. Although this change took longer in zebra finches, the altered proportion of syllables eventually stabilized. Results indicated that adult songs show similar aspects after auditory deprivation, and that neural mechanisms involved in the maintenance of high-frequency song syllables, using auditory feedback, may be present in both finches despite species differences.     

Variability and coding efficiency of noisy neural spike encoders

Encoding synaptic inputs as a train of action potentials is a fundamental function of nerve cells. Although spike trains recorded in vivo have been shown to be highly variable, it is unclear whether variability in spike timing represents faithful encoding of temporally varying synaptic inputs or noise inherent in the spike encoding mechanism. It has been reported that spike timing variability is more pronounced for constant, unvarying inputs than for inputs with rich temporal structure. This could have significant implications for the nature of neural coding, particularly if precise timing of spikes and temporal synchrony between neurons is used to represent information in the nervous system. To study the potential functional role of spike timing variability, we estimate the fraction of spike timing variability which conveys information about the input for two types of noisy spike encoders--an integrate and fire model with randomly chosen thresholds and a model of a patch of neuronal membrane containing stochastic Na(+) and K(+) channels obeying Hodgkin-Huxley kinetics. The quality of signal encoding is assessed by reconstructing the input stimuli from the output spike trains using optimal linear mean square estimation. A comparison of the estimation performance of noisy neuronal models of spike generation enables us to assess the impact of neuronal noise on the efficacy of neural coding. The results for both models suggest that spike timing variability reduces the ability of spike trains to encode rapid time-varying stimuli. Moreover, contrary to expectations based on earlier studies, we find that the noisy spike encoding models encode slowly varying stimuli more effectively than rapidly varying ones.     

Intrinsic bursting enhances the robustness of a neural network model of sequence generation by avian brain area HVC

Avian brain area HVC is known to be important for the production of birdsong. In zebra finches, each RA-projecting neuron in HVC emits a single burst of spikes during a song motif. The population of neurons is activated in a precisely timed, stereotyped sequence. We propose a model of these burst sequences that relies on two hypotheses. First, we hypothesize that the sequential order of bursting is reflected in the excitatory synaptic connections between neurons. Second, we propose that the neurons are intrinsically bursting, so that burst duration is set by cellular properties. Our model generates burst sequences similar to those observed in HVC. If intrinsic bursting is removed from the model, burst sequences can also be produced. However, they require more fine-tuning of synaptic strengths, and are therefore less robust. In our model, intrinsic bursting is caused by dendritic calcium spikes, and strong spike frequency adaptation in the soma contributes to burst termination.     

The Essential Complexity of Auditory Receptive Fields

Encoding properties of sensory neurons are commonly modeled using linear finite impulse response (FIR) filters. For the auditory system, the FIR filter is instantiated in the spectro-temporal receptive field (STRF), often in the framework of the generalized linear model. Despite widespread use of the FIR STRF, numerous formulations for linear filters are possible that require many fewer parameters, potentially permitting more efficient and accurate model estimates. To explore these alternative STRF architectures, we recorded single-unit neural activity from auditory cortex of awake ferrets during presentation of natural sound stimuli. We compared performance of > 1000 linear STRF architectures, evaluating their ability to predict neural responses to a novel natural stimulus. Many were able to outperform the FIR filter. Two basic constraints on the architecture lead to the improved performance: (1) factorization of the STRF matrix into a small number of spectral and temporal filters and (2) low-dimensional parameterization of the factorized filters. The best parameterized model was able to outperform the full FIR filter in both primary and secondary auditory cortex, despite requiring fewer than 30 parameters, about 10% of the number required by the FIR filter. After accounting for noise from finite data sampling, these STRFs were able to explain an average of 40% of A1 response variance. The simpler models permitted more straightforward interpretation of sensory tuning properties. They also showed greater benefit from incorporating nonlinear terms, such as short term plasticity, that provide theoretical advances over the linear model. Architectures that minimize parameter count while maintaining maximum predictive power provide insight into the essential degrees of freedom governing auditory cortical function. They also maximize statistical power available for characterizing additional nonlinear properties that limit current auditory models.