Electrophysiology publications

Sleep changes were studied in mice (n = 59) from early adolescence to adulthood (postnatal days P19–111). REM sleep declined steeply in early adolescence, while total sleep remained constant and NREM sleep increased slightly. Four hours of sleep deprivation starting at light onset were performed from ages P26 through adulthood (>P60). Following this acute sleep deprivation all mice slept longer and with more consolidated sleep bouts, while NREM slow wave activity (SWA) showed high interindividual variability in the younger groups, and increased consistently only after P42. Three parameters together explained up to 67% of the variance in SWA rebound in frontal cortex, including weight-adjusted age and increase in alpha power during sleep deprivation, both of which positively correlated with the SWA response. The third, and strongest predictor was the SWA decline during the light phase in baseline: mice with high peak SWA at light onset, resulting in a large SWA decline, were more likely to show no SWA rebound after sleep deprivation, a result that was also confirmed in parietal cortex. During baseline, however, SWA showed the same homeostatic changes in adolescents and adults, declining in the course of sleep and increasing across periods of spontaneous wake. Thus, we hypothesize that, in young adolescent mice, a ceiling effect and not the immaturity of the cellular mechanisms underlying sleep homeostasis may prevent the SWA rebound when wake is extended beyond its physiological duration.

The relationship between structure and function is an invaluable context with which to explore biological mechanisms of normal and dysfunctional hearing. The systematic and topographic representation of frequency originates at the cochlea, and is retained throughout much of the central auditory system. The cochlear nucleus (CN), which initiates all ascending auditory pathways, represents an essential link for understanding frequency organization. A model of the CN that maps frequency representation in 3D would facilitate investigations of possible frequency specializations and pathologic changes that disturb frequency organization. Toward this goal, we reconstructed in 3D the trajectories of labeled auditory nerve (AN) fibers following multiunit recordings and dye injections in the anteroventral CN of the CBA/J mouse. We observed that each injection produced a continuous sheet of labeled AN fibers. Individual cases were normalized to a template using 3D alignment procedures that revealed a systematic and tonotopic arrangement of AN fibers in each subdivision with a clear indication of isofrequency laminae. The combined dataset was used to mathematically derive a 3D quantitative map of frequency organization throughout the entire volume of the CN. This model, available online ( http://3D.ryugolab.com/ ), can serve as a tool for quantitatively testing hypotheses concerning frequency and location in the CN. J. Comp. Neurol. 521:1510–1532, 2013. © 2012 Wiley Periodicals, Inc.

Objective: Functional electrical stimulation (FES) approaches often utilize an open-loop controller to drive state transitions. The addition of sensory feedback may allow for closed-loop control that can respond effectively to perturbations and muscle fatigue. Approach: We evaluated the use of natural sensory nerve signals obtained with penetrating microelectrode arrays in lumbar dorsal root ganglia (DRG) as real-time feedback for closed-loop control of FES-generated hind limb stepping in anesthetized cats. Main results: Leg position feedback was obtained in near real-time at 50 ms intervals by decoding the firing rates of more than 120 DRG neurons recorded simultaneously. Over 5 m of effective linear distance was traversed during closed-loop stepping trials in each of two cats. The controller compensated effectively for perturbations in the stepping path when DRG sensory feedback was provided. The presence of stimulation artifacts and the quality of DRG unit sorting did not significantly affect the accuracy of leg position feedback obtained from the linear decoding model as long as at least 20 DRG units were included in the model. Significance: This work demonstrates the feasibility and utility of closed-loop FES control based on natural neural sensors. Further work is needed to improve the controller and electrode technologies and to evaluate long-term viability.

The thalamus is crucial in determining the sensory information conveyed to cortex. In the visual system, the thalamic lateral geniculate nucleus (LGN) is generally thought to encode simple center-surround receptive fields, which are combined into more sophisticated features in cortex, such as orientation and direction selectivity. However, recent evidence suggests that a more diverse set of retinal ganglion cells projects to the LGN. We therefore used multisite extracellular recordings to define the repertoire of visual features represented in the LGN of mouse, an emerging model for visual processing. In addition to center-surround cells, we discovered a substantial population with more selective coding properties, including direction and orientation selectivity, as well as neurons that signal absence of contrast in a visual scene. The direction and orientation selective neurons were enriched in regions that match the termination zones of direction selective ganglion cells from the retina, suggesting a source for their tuning. Together, these data demonstrate that the mouse LGN contains a far more elaborate representation of the visual scene than current models posit. These findings should therefore have a significant impact on our understanding of the computations performed in mouse visual cortex.

Given the extraordinary ability of humans and animals to recognize communication signals over a background of noise, describing noise invariant neural responses is critical not only to pinpoint the brain regions that are mediating our robust perceptions but also to understand the neural computations that are performing these tasks and the underlying circuitry. Although invariant neural responses, such as rotation-invariant face cells, are well described in the visual system, high-level auditory neurons that can represent the same behaviorally relevant signal in a range of listening conditions have yet to be discovered. Here we found neurons in a secondary area of the avian auditory cortex that exhibit noise-invariant responses in the sense that they responded with similar spike patterns to song stimuli presented in silence and over a background of naturalistic noise. By characterizing the neurons' tuning in terms of their responses to modulations in the temporal and spectral envelope of the sound, we then show that noise invariance is partly achieved by selectively responding to long sounds with sharp spectral structure. Finally, to demonstrate that such computations could explain noise invariance, we designed a biologically inspired noise-filtering algorithm that can be used to separate song or speech from noise. This novel noise-filtering method performs as well as other state-of-the-art de-noising algorithms and could be used in clinical or consumer oriented applications. Our biologically inspired model also shows how high-level noise-invariant responses could be created from neural responses typically found in primary auditory cortex.

Inner hair cells (IHCs) of the cochlea use ribbon synapses to transmit auditory information faithfully to spiral ganglion neurons (SGNs). In the present study, we used genetic disruption of the presynaptic scaffold protein bassoon in mice to manipulate the morphology and function of the IHC synapse. Although partial-deletion mutants lacking functional bassoon (Bsn ΔEx4/5 ) had a near-complete loss of ribbons from the synapses (up to 88% ribbonless synapses), gene-trap mutants (Bsn gt ) showed weak residual expression of bassoon and 56% ribbonless synapses, whereas the remaining 44% had a loosely anchored ribbon. Patch-clamp recordings and synaptic Ca V 1.3 immunolabeling indicated a larger number of Ca 2+ channels for Bsn gt IHCs compared with Bsn ΔEx4/5 IHCs and for Bsn gt ribbon-occupied versus Bsn gt ribbonless synapses. An intermediate phenotype of Bsn gt IHCs was also found by membrane capacitance measurements for sustained exocytosis, but not for the size of the readily releasable vesicle pool. The frequency and amplitude of EPSCs were reduced in Bsn ΔEx4/5 mouse SGNs, whereas their postsynaptic AMPA receptor clusters were largely unaltered. Sound coding in SGN, assessed by recordings of single auditory nerve fibers and their population responses in vivo , was similarly affected in Bsn gt and Bsn ΔEx4/5 mice. Both genotypes showed impaired sound onset coding and reduced evoked and spontaneous spike rates. In summary, reduced bassoon expression or complete lack of full-length bassoon impaired sound encoding to a similar extent, which is consistent with the comparable reduction of the readily releasable vesicle pool. This suggests that the remaining loosely anchored ribbons in Bsn gt IHCs were functionally inadequate or that ribbon independent mechanisms dominated the coding deficit.

Objectives: Presbyacusis, one of the most common ailments of the elderly, is often treated with hearing aids, which serve to reintroduce some or all of those sounds lost to peripheral hearing loss. However, little is known about the underlying changes to the ear and brain as a result of such experience with sound late in life. The present study attempts to model this process by rearing aged CBA mice in an augmented acoustic environment (AAE). Design: Aged (22-23 months) male (n = 12) and female (n = 9) CBA/CaJ mice were reared in either 6 weeks of low-level (70 dB SPL) broadband noise stimulation (AAE) or normal vivarium conditions. Changes as a function of the treatment were measured for behavior, auditory brainstem response thresholds, hair cell cochleograms, and gamma aminobutyric acid neurochemistry in the key central auditory structures of the inferior colliculus and primary auditory cortex. Results: The AAE-exposed group was associated with sex-specific changes in cochlear pathology, auditory brainstem response thresholds, and gamma aminobutyric acid neurochemistry. Males exhibited significantly better thresholds and reduced hair cell loss (relative to controls) whereas females exhibited the opposite effect. AAE was associated with increased glutamic acid decarboxylase (GAD67) levels in the inferior colliculus of both male and female mice. However, in primary auditory cortex AAE exposure was associated with increased GAD67 labeling in females and decreased GAD67 in males. Conclusions: These findings suggest that exposing aged mice to a low-level AAE alters both peripheral and central properties of the auditory system and these changes partially interact with sex or the degree of hearing loss before AAE. Although direct application of these findings to hearing aid use or auditory training in aged humans would be premature, the results do begin to provide direct evidence for the underlying changes that might be occurring as a result of hearing aid use late in life. These results suggest the aged brain retains significantly anatomical, electrophysiological, and neurochemical plasticity.

Auditory-cued behavioral training can alter neural circuits in primary auditory cortex (A1), but the mechanisms and consequences of experience-dependent cortical plasticity are not fully understood. To address this issue, we trained adult rats to detect a 5 kHz target in order to receive a food reward. After 14 days training we identified three locations within A1: (i) the region representing the characteristic frequency (CF) 5 kHz, (ii) a nearby region with CF ∼10 kHz, and (iii) a more distant region with CF ∼20 kHz. In order to compare functional connectivity in A1 near to, vs. far from, the representation of the target frequency, we placed a 16-channel multiprobe in middle- (∼10 kHz) and high- (∼20 kHz) CF regions and obtained current-source density (CSD) profiles evoked by a range of tone stimuli (CF±1-3 octaves in quarter-octave steps). Our aim was to construct "CSD receptive fields" (CSD RFs) in order to determine the laminar and spectral profile of tone-evoked current sinks, and infer changes to thalamocortical and intracortical inputs. Behavioral training altered CSD RFs at the 10 kHz, but not 20 kHz, site relative to CSD RFs in untrained control animals. At the 10 kHz site, current sinks evoked by the target frequency were enhanced in layer 2/3, but the initial current sink in layer 4 was not altered. The results imply training-induced plasticity along intracortical pathways connecting the target representation with nearby cortical regions. Finally, we related behavioral performance (sensitivity index, d') to CSD responses in individual animals, and found a significant correlation between the development of d' over training and the amplitude of the target-evoked current sink in layer 2/3. The results suggest that plasticity along intracortical pathways is important for auditory learning.

The basolateral amygdala (BLA) modulates memory, particularly for arousing or emotional events, during post-training periods of consolidation. It strengthens memories whose substrates in part or whole are stored remotely, in structures such as the hippocampus, striatum and cerebral cortex. However, the mechanisms by which the BLA influences distant memory traces are unknown, largely because of the need for identifiable target mnemonic representations. Associative tuning plasticity in the primary auditory cortex (A1) constitutes a well-characterized candidate specific memory substrate that is ubiquitous across species, tasks and motivational states. When tone predicts reinforcement, the tuning of cells in A1 shifts toward or to the signal frequency within its tonotopic map, producing an over-representation of behaviorally important sounds. Tuning shifts have the cardinal attributes of forms of memory, including associativity, specificity, rapid induction, consolidation and long-term retention and are therefore likely memory representations. We hypothesized that the BLA strengthens memories by increasing their cortical representations. We recorded multiple unit activity from A1 of rats that received a single discrimination training session in which two tones (2.0 s) separated by 1.25 octaves were either paired with brief electrical stimulation (400 ms) of the BLA (CS+) or not (CS-). Frequency response areas generated by presenting a matrix of test tones (0.5-53.82 kHz, 0-70 dB) were obtained before training and daily for 3 weeks post-training. Tuning both at threshold and above threshold shifted predominantly toward the CS+ beginning on day 1. Tuning shifts were maintained for the entire 3 weeks. Absolute threshold and bandwidth decreased, producing less enduring increases in sensitivity and selectivity. BLA-induced tuning shifts were associative, highly specific and long-lasting. We propose that the BLA strengthens memory for important experiences by increasing the number of neurons that come to best represent that event. Traumatic, intrusive memories might reflect abnormally extensive representational networks due to hyper-activity of the BLA consequent to the release of excessive amounts of stress hormones.

Evidence regarding the functional subdivisions of human auditory cortex has been slow to converge on a definite model. In part, this reflects inadequacies of current understanding of how the cortex represents temporal information in acoustic signals. To address this, we investigated spatiotemporal properties of auditory responses in human posterolateral superior temporal (PLST) gyrus to acoustic click-train stimuli using intracranial recordings from neurosurgical patients. Subjects were patients undergoing chronic invasive monitoring for refractory epilepsy. The subjects listened passively to acoustic click-train stimuli of varying durations (160 or 1,000 ms) and rates (4–200 Hz), delivered diotically via insert earphones. Multicontact subdural grids placed over the perisylvian cortex recorded intracranial electrocorticographic responses from PLST and surrounding areas. Analyses focused on averaged evoked potentials (AEPs) and high gamma (70–150 Hz) event-related band power (ERBP). Responses to click trains featured prominent AEP waveforms and increases in ERBP. The magnitude of AEPs and ERBP typically increased with click rate. Superimposed on the AEPs were frequency-following responses (FFRs), most prominent at 50-Hz click rates but still detectable at stimulus rates up to 200 Hz. Loci with the largest high gamma responses on PLST were often different from those sites that exhibited the strongest FFRs. The data indicate that responses of non-core auditory cortex of PLST represent temporal stimulus features in multiple ways. These include an isomorphic representation of periodicity (as measured by the FFR), a representation based on increases in non-phase-locked activity (as measured by high gamma ERBP), and spatially distributed patterns of activity.

The potassium channel antagonist 4-aminopyridine (4-AP) improves a variety of motor abnormalities associated with disorders of the cerebellum. The most rigorous quantitative data relate to 4-AP's ability to improve eye movement deficits in humans referable to dysfunction of the cerebellar flocculus. Largely based on work in the ataxic mouse mutant tottering (which carries a mutation of the Cacna1a gene of the P/Q voltage-activated calcium channel), 4-AP is hypothesized to function by enhancing excitability or rhythmicity of floccular Purkinje cells. We tested this hypothesis by determining whether systemic or intrafloccular administration of 4-AP would ameliorate the eye movement deficits in tottering that are attributable to flocculus dysfunction, including the reductions in amplitude of the yaw-axis vestibulo-ocular reflex (VOR) and vision-enhanced vestibulo-ocular reflex (VVOR), and the optokinetic reflex (OKR) about yaw and roll axes. Because tottering's deficits increase with age, both young and elderly mutants were tested to detect any age-dependent 4-AP effects. 4-AP failed to improve VOR, VVOR, and OKR gains during sinusoidal stimuli, although it may have reduced the tendency of the mutants' responses to VOR and VVOR to decline over the course of a one-hour recording session. For constant-velocity optokinetic stimuli, 4-AP generated some enhancement of yaw OKR and upward-directed roll OKR, but the effects were also seen in normal C57BL/6 controls, and thus do not represent a specific reversal of the electrophysiological consequences of the tottering mutation. Data support a possible extra-floccular locus for the effects of 4-AP on habituation and roll OKR. Unilateral intrafloccular 4-AP injections did not affect ocular motility, except to generate mild eye elevations, consistent with reduced floccular output. Because 4-AP did not produce the effects expected if it normalized outputs of floccular Purkinje cells, there is a need for further studies to elucidate the drug's mechanism of action on cerebellar motor dysfunction.

The efficacy of spinal cord stimulators is dependent on the ability of the device to functionally activate targeted structures within the spinal cord, while avoiding activation of near-by non-targeted structures. In theory, these objectives can best be achieved by delivering electrical stimuli directly to the surface of the spinal cord. The current experiments were performed to study the influence of different stimulating electrode positions on patterns of spinal cord electrophysiological activation. A custom-designed spinal cord neurostimulator was used to investigate the effects of lead position and stimulus amplitude on cortical electrophysiological responses to spinal cord stimulation. Brain recordings were obtained from subdural grids placed in four adult sheep. We systematically varied the position of the stimulating lead relative to the spinal cord and the voltage delivered by the device at each position, and then examined how these variables influenced cortical responses. A clear relationship was observed between voltage and electrode position, and the magnitude of high gamma-band oscillations. Direct stimulation of the dorsal column contralateral to the grid required the lowest voltage to evoke brain responses to spinal cord stimulation. Given the lower voltage thresholds associated with direct stimulation of the dorsal column, and its possible impact on the therapeutic window, this intradural modality may have particular clinical advantages over standard epidural techniques now in routine use.

The central thalamus plays an important role in the regulation of arousal and allocation of attentional resources in the performance of even simple tasks. To assess the contribution of central thalamic neurons to short-term adjustments of attentional effort, we analyzed 166 microelectrode recordings obtained from two rhesus monkeys performing a visuomotor simple reaction time task with a variable foreperiod. Multiunit responses showed maintained firing rate elevations during the variable delay period of the task in ∼24% of recording sites. Simultaneously recorded local field potentials demonstrated significant decreases in power at ∼10–20 Hz and increases in power at 30–100 Hz during the delay period when compared against precue baselines. Comparison of the spectral power of local field potentials during the delay period of correct and incorrect trials showed that, during incorrect trials, similar, but reduced, shifts of spectral power occurred within the same frequency bands. Sustained performance of even simple tasks requires regulation of arousal and attention that combine in the concept of “attentional effort”. Our findings suggest that central thalamic neurons regulate task performance through brief changes in firing rates and spectral power changes during task-relevant short-term shifts of attentional effort. Increases in attentional effort may be reflected in changes within the central thalamic local populations, where correct task performance associates with more robust maintenance of firing rates during the delay period. Such ongoing fluctuations of central thalamic activity likely reflect a mix of influences, including variations in moment-to-moment levels of motivation, arousal, and availability of cognitive resources.

Much is known about the neuronal cell types and circuitry of the mammalian respiratory brainstem and its role in normal, quiet breathing. Our understanding of the role of respiration in the context of vocal production, however, is very limited. Songbirds contain a well-defined neural circuit, known as the song system, which is necessary for song production and is strongly coupled to the respiratory system. A major target of this system is nucleus parambigualis (PAm) in the ventrolateral medulla, a structure that controls inspiration by way of its bulbospinal projections but is also an integral part of the song-pattern generation circuit by way of its “thalamocortical” projections to song-control nuclei in the telencephalon. We have mapped out PAm to characterize the cell types and its functional organization. Extracellular single units were obtained in anesthetized adult male zebra finches while measuring air sac pressure to monitor respiration. Single units were characterized by their discharge patterns and the phase of the activity in the respiratory cycle. Several classes of neurons were identified and were analogous to those reported for mammalian medullary respiratory neurons. The majority of the neurons in PAm was classified as inspiratory augmenting or preinspiratory, although other basic discharge patterns were observed as well. The well-characterized connectivity of PAm within the vocal motor circuit and the similarity of its neural firing patterns to the rostral ventral respiratory group and pre-Bötzinger complex of mammals make it an ideal system for investigating the integration of breathing and vocalization.

One in 15 school age children have dyslexia, which is characterized by phoneme-processing problems and difficulty learning to read. Dyslexia is associated with mutations in the gene KIAA0319. It is not known whether reduced expression of KIAA0319 can degrade the brain's ability to process phonemes. In the current study, we used RNA interference (RNAi) to reduce expression of Kiaa0319 (the rat homolog of the human gene KIAA0319) and evaluate the effect in a rat model of phoneme discrimination. Speech discrimination thresholds in normal rats are nearly identical to human thresholds. We recorded multiunit neural responses to isolated speech sounds in primary auditory cortex (A1) of rats that received in utero RNAi of Kiaa0319. Reduced expression of Kiaa0319 increased the trial-by-trial variability of speech responses and reduced the neural discrimination ability of speech sounds. Intracellular recordings from affected neurons revealed that reduced expression of Kiaa0319 increased neural excitability and input resistance. These results provide the first evidence that decreased expression of the dyslexia-associated gene Kiaa0319 can alter cortical responses and impair phoneme processing in auditory cortex.

The convergence of eye-specific thalamic inputs to visual cortical neurons forms the basis of binocular vision. Inputs from the same eye that signal light increment (On) and decrement (Off) are spatially segregated into subregions, giving rise to cortical receptive fields (RFs) that are selective for stimulus orientation. Here we map RFs of binocular neurons in the mouse primary visual cortex using spike-triggered average. We find that subregions of the same sign (On-On and Off-Off) preferentially overlap between the 2 monocular RFs, leading to binocularly matched orientation tuning. We further demonstrate that such subregion correspondence and the consequent matching of RF orientation are disrupted in mice reared in darkness during development. Surprisingly, despite the lack of all postnatal visual experience, a substantial degree of subregion correspondence still remains. In addition, dark-reared mice show normal monocular RF structures and binocular overlap. These results thus reveal the specific roles of experience-dependent and -independent processes in binocular convergence and refinement of On and Off inputs onto single cortical neurons.

In their natural environment, animals face complex and highly dynamic olfactory input. Thus vertebrates as well as invertebrates require fast and reliable processing of olfactory information. Parallel processing has been shown to improve processing speed and power in other sensory systems and is characterized by extraction of different stimulus parameters along parallel sensory information streams. Honeybees possess an elaborate olfactory system with unique neuronal architecture: a dual olfactory pathway comprising a medial projection-neuron (PN) antennal lobe (AL) protocerebral output tract (m-APT) and a lateral PN AL output tract (l-APT) connecting the olfactory lobes with higher-order brain centers. We asked whether this neuronal architecture serves parallel processing and employed a novel technique for simultaneous multiunit recordings from both tracts. The results revealed response profiles from a high number of PNs of both tracts to floral, pheromonal, and biologically relevant odor mixtures tested over multiple trials. PNs from both tracts responded to all tested odors, but with different characteristics indicating parallel processing of similar odors. Both PN tracts were activated by widely overlapping response profiles, which is a requirement for parallel processing. The l-APT PNs had broad response profiles suggesting generalized coding properties, whereas the responses of m-APT PNs were comparatively weaker and less frequent, indicating higher odor specificity. Comparison of response latencies within and across tracts revealed odor-dependent latencies. We suggest that parallel processing via the honeybee dual olfactory pathway provides enhanced odor processing capabilities serving sophisticated odor perception and olfactory demands associated with a complex olfactory world of this social insect.

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

Whole-cell patch neurophysiology and pharmacological manipulations have provided unprecedented insight into the functions of central neurons, but their combined use has been largely restricted to in vitro preparations. We describe a method for performing whole-cell patch recording and focal application of pharmacological agents in vivo. A key feature of this technique involves iontophoresis of glutamate to establish proximity of drug and recording pipettes. We show data from iontophoresis of glutamate during extracellular and whole-cell recordings made in vivo from auditory neurons in the midbrain of the leopard frog, Rana pipiens, and the effects of blocking GABA(A) receptors while making a whole-cell recording. This methodology should accelerate our understanding of the roles of particular neurotransmitter systems in normal and pathological conditions, and facilitate investigation of the in vivo effects of drugs and the mechanisms underlying computations.

Photoacoustic tomography (PAT) offers three-dimensional (3D) structural and functional imaging of living biological tissue with label-free, optical absorption contrast. These attributes lend PAT imaging to a wide variety of applications in clinical medicine and preclinical research. Despite advances in live animal imaging with PAT, there is still a need for 3D imaging at centimeter depths in real-time. We report the development of four dimensional (4D) PAT, which integrates time resolutions with 3D spatial resolution, obtained using spherical arrays of ultrasonic detectors. The 4D PAT technique generates motion pictures of imaged tissue, enabling real time tracking of dynamic physiological and pathological processes at hundred micrometer-millisecond resolutions. The 4D PAT technique is used here to image needle-based drug delivery and pharmacokinetics. We also use this technique to monitor 1) fast hemodynamic changes during inter-ictal epileptic seizures and 2) temperature variations during tumor thermal therapy.