skip to main content

Title: Cortical Localization of the Auditory Temporal Response Function from MEG via Non-convex Optimization
The magnetoencephalography (MEG) response to continuous auditory stimuli, such as speech, is commonly described using a linear filter, the auditory temporal response function (TRF). Though components of the sensor level TRFs have been well characterized, the underlying neural sources responsible for these components are not well understood. In this work, we provide a unified framework for determining the TRFs of neural sources directly from the MEG data, by integrating the TRF and distributed forward source models into one, and casting the joint estimation task as a Bayesian optimization problem. Though the resulting problem emerges as non-convex, we propose efficient solutions that leverage recent advances in evidence maximization. We demonstrate the effectiveness of the resulting algorithm in both simulated and experimentally recorded MEG data from humans.
; ; ;
Award ID(s):
Publication Date:
Journal Name:
2018 52nd Asilomar Conference on Signals, Systems, and Computers
Page Range or eLocation-ID:
373 - 378
Sponsoring Org:
National Science Foundation
More Like this
  1. Objective: The Temporal Response Function (TRF) is a linear model of neural activity time-locked to continuous stimuli, including continuous speech. TRFs based on speech envelopes typically have distinct components that have provided remarkable insights into the cortical processing of speech. However, current methods may lead to less than reliable estimates of single-subject TRF components. Here, we compare two established methods, in TRF component estimation, and also propose novel algorithms that utilize prior knowledge of these components, bypassing the full TRF estimation. Methods: We compared two established algorithms, ridge and boosting, and two novel algorithms based on Subspace Pursuit (SP) andmore »Expectation Maximization (EM), which directly estimate TRF components given plausible assumptions regarding component characteristics. Single-channel, multi-channel, and source-localized TRFs were fit on simulations and real magnetoencephalographic data. Performance metrics included model fit and component estimation accuracy. Results: Boosting and ridge have comparable performance in component estimation. The novel algorithms outperformed the others in simulations, but not on real data, possibly due to the plausible assumptions not actually being met. Ridge had slightly better model fits on real data compared to boosting, but also more spurious TRF activity. Conclusion: Results indicate that both smooth (ridge) and sparse (boosting) algorithms perform comparably at TRF component estimation. The SP and EM algorithms may be accurate, but rely on assumptions of component characteristics. Significance: This systematic comparison establishes the suitability of widely used and novel algorithms for estimating robust TRF components, which is essential for improved subject-specific investigations into the cortical processing of speech.« less
  2. Identifying the directed connectivity that underlie networked activity between different cortical areas is critical for understanding the neural mechanisms behind sensory processing. Granger causality (GC) is widely used for this purpose in functional magnetic resonance imaging analysis, but there the temporal resolution is low, making it difficult to capture the millisecond-scale interactions underlying sensory processing. Magne- toencephalography (MEG) has millisecond resolution, but only provides low-dimensional sensor-level linear mixtures of neural sources, which makes GC inference challenging. Conventional methods proceed in two stages: First, cortical sources are estimated from MEG using a source localization technique, followed by GC inference among themore »estimated sources. However, the spatiotemporal biases in estimating sources propagate into the subsequent GC analysis stage, may result in both false alarms and missing true GC links. Here, we introduce the Network Localized Granger Causality (NLGC) inference paradigm, which models the source dynamics as latent sparse multivariate autoregressive processes and estimates their parameters directly from the MEG measurements, integrated with source localization, and employs the resulting parameter estimates to produce a precise statistical characterization of the detected GC links. We offer several theoretical and algorithmic innovations within NLGC and further examine its utility via comprehensive simulations and application to MEG data from an auditory task involving tone processing from both younger and older participants. Our simulation studies reveal that NLGC is markedly robust with respect to model mismatch, network size, and low signal-to-noise ratio, whereas the conventional two-stage methods result in high false alarms and mis-detections. We also demonstrate the advantages of NLGC in revealing the cortical network- level characterization of neural activity during tone processing and resting state by delineating task- and age-related connectivity changes.« less
  3. Abstract Harmful algal blooms produce potent neurotoxins that accumulate in seafood and are hazardous to human health. Developmental exposure to the harmful algal bloom toxin, domoic acid (DomA), has behavioral consequences well into adulthood, but the cellular and molecular mechanisms of DomA developmental neurotoxicity are largely unknown. To assess these, we exposed zebrafish embryos to DomA during the previously identified window of susceptibility and used the well-known startle response circuit as a tool to identify specific neuronal components that are targeted by exposure to DomA. Exposure to DomA reduced startle responsiveness to both auditory/vibrational and electrical stimuli, and even atmore »the highest stimulus intensities tested, led to a dramatic reduction of one type of startle (short latency c-starts). Furthermore, DomA-exposed larvae had altered kinematics for both types of startle responses tested, exhibiting shallower bend angles and slower maximal angular velocities. Using vital dye staining, immunolabelling, and live imaging of transgenic lines, we determined that while the sensory inputs were intact, the reticulospinal neurons required for short latency c-starts were absent in most DomA-exposed larvae. Furthermore, axon tracing revealed that DomA-treated larvae also showed significantly reduced primary motor neuron axon collaterals. Overall, these results show that developmental exposure to DomA targets large reticulospinal neurons and motor neuron axon collaterals, resulting in measurable deficits in startle behavior. They further provide a framework for using the startle response circuit to identify specific neural populations disrupted by toxins or toxicants and to link these disruptions to functional consequences for neural circuit function and behavior.« less
  4. Abstract

    We combined magnetoencephalography (MEG), 7 T proton magnetic resonance spectroscopy (MRS), and 7 T fMRI during performance of a task in a group of 23 first episode psychosis (FEP) patients and 26 matched healthy controls (HC). We recorded both the auditory evoked response to 40 Hz tone clicks and the resting state in MEG. Neurometabolite levels were obtained from the anterior cingulate cortex (ACC). The fMRI BOLD response was obtained during the Stroop inhibitory control task. FEP showed a significant increase in resting state low frequency theta activity (p < 0.05; Cohend= 0.69), but no significant difference in the 40 Hz auditory evoked response compared tomore »HC. An across-groups whole brain analysis of the fMRI BOLD response identified eight regions that were significantly activated during task performance (p < 0.01, FDR-corrected); the mean signal extracted from those regions was significantly different between the groups (p = 0.0006;d = 1.19). In the combined FEP and HC group, there was a significant correlation between the BOLD signal during task performance and MEG resting state low frequency activity (p < 0.05). In FEP, we report significant alteration in resting state low frequency MEG activity, but no alterations in auditory evoked gamma band response, suggesting that the former is a more robust biomarker of early psychosis. There were no correlations between gamma oscillations and GABA levels in either HC or FEP. Finally, in this study, each of the three imaging modalities differentiated FEP from HC; fMRI with good and MEG and MRS with moderate effect size.

    « less
  5. Abstract Transfer RNA-derived fragments (tRFs) are a new class of small non-coding RNAs and play important roles in biological and physiological processes. Prediction of tRF target genes and binding sites is crucial in understanding the biological functions of tRFs in the molecular mechanisms of human diseases. We developed a publicly accessible web-based database, tRFtarget (, for tRF target prediction. It contains the computationally predicted interactions between tRFs and mRNA transcripts using the two state-of-the-art prediction tools RNAhybrid and IntaRNA, including location of the binding sites on the target, the binding region, and free energy of the binding stability with graphicmore »illustration. tRFtarget covers 936 tRFs and 135 thousand predicted targets in eight species. It allows researchers to search either target genes by tRF IDs or tRFs by gene symbols/transcript names. We also integrated the manually curated experimental evidence of the predicted interactions into the database. Furthermore, we provided a convenient link to the DAVID® web server to perform downstream functional pathway analysis and gene ontology annotation on the predicted target genes. This database provides useful information for the scientific community to experimentally validate tRF target genes and facilitate the investigation of the molecular functions and mechanisms of tRFs.« less