Computational modeling of threat learning reveals links with anxiety and neuroanatomy in humans
Abstract
Influential theories implicate variations in the mechanisms supporting threat learning in the severity of anxiety symptoms. We use computational models of associative learning in conjunction with structural imaging to explicate links among the mechanisms underlying threat learning, their neuroanatomical substrates, and anxiety severity in humans. We recorded skin-conductance data during a threat-learning task from individuals with and without anxiety disorders (N=251; 8-50 years; 116 females). Reinforcement-learning model variants quantified processes hypothesized to relate to anxiety: threat conditioning, threat generalization, safety learning, and threat extinction. We identified the best-fitting models for these processes and tested associations among latent learning parameters, whole-brain anatomy, and anxiety severity. Results indicate that greater anxiety severity related specifically to slower safety learning and slower extinction of response to safe stimuli. Nucleus accumbens gray-matter volume moderated learning-anxiety associations. Using a modeling approach, we identify computational mechanisms linking threat learning and anxiety severity and their neuroanatomical substrates.
Data availability
We cannot share the full dataset due to the NIH IRB requirements, which require participants to explicitly consent to their data being shared publicly. An important element in that is to protect patients who agree to participate in studies that relate to their psychopathology. Such consent was not acquired from most participants; as such, we cannot upload our complete dataset in its raw or deidentified form, or derivatives of the data, since we will be violating IRB protocols. Still, a subset of participants did consent to data sharing and we have uploaded their data as noted in the revised manuscript (https://github.com/rany-abend/threat_learning_eLife). Researchers interested in potentially acquiring access to the data could contact Dr. Daniel Pine (pined@mail.nih.gov), Chief of the Emotion and Development Branch at NIH, with a research proposal; as per IRB rules, the IRB may approve adding such researchers as Associate Investigators if a formal collaboration is initiated. No commercial use of the data is a lowed. The modeling and imaging analyses have now been uploaded in full as source code files.
Article and author information
Author details
Funding
National Institutes of Health (ZIAMH002781-15)
- Daniel S Pine
National Institutes of Health (R00MH091183)
- Jennifer C Britton
Brain and Behavior Research Foundation (28239)
- Rany Abend
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Human subjects: Human Subjects: Yes Ethics Statement: Written informed consent was obtained from adult (greater than or equal to 18 years) participants as we l as parents, and written assent was obtained from youth. Procedures were approved by the NIMH Institutional Review Board (protocol 01-M-0192).
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 2,182
- views
-
- 385
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Neuroscience
We determined the intersubject association between the rhythmic entrainment abilities of human subjects during a synchronization-continuation tapping task (SCT) and the macro- and microstructural properties of their superficial (SWM) and deep (DWM) white matter. Diffusion-weighted images were obtained from 32 subjects who performed the SCT with auditory or visual metronomes and five tempos ranging from 550 to 950 ms. We developed a method to determine the density of short-range fibers that run underneath the cortical mantle, interconnecting nearby cortical regions (U-fibers). Notably, individual differences in the density of U-fibers in the right audiomotor system were correlated with the degree of phase accuracy between the stimuli and taps across subjects. These correlations were specific to the synchronization epoch with auditory metronomes and tempos around 1.5 Hz. In addition, a significant association was found between phase accuracy and the density and bundle diameter of the corpus callosum (CC), forming an interval-selective map where short and long intervals were behaviorally correlated with the anterior and posterior portions of the CC. These findings suggest that the structural properties of the SWM and DWM in the audiomotor system support the tapping synchronization abilities of subjects, as cortical U-fiber density is linked to the preferred tapping tempo and the bundle properties of the CC define an interval-selective topography.
-
- Computational and Systems Biology
- Neuroscience
In this study, we develop new reverse engineering (RE) techniques to identify the organization of the synaptic inputs generating firing patterns of populations of neurons. We tested these techniques in silico to allow rigorous evaluation of their effectiveness, using remarkably extensive parameter searches enabled by massively-parallel computation on supercomputers. We chose spinal motoneurons as our target neural system, since motoneurons process all motor commands and have well-established input-output properties. One set of simulated motoneurons was driven by 300,000+ simulated combinations of excitatory, inhibitory, and neuromodulatory inputs. Our goal was to determine if these firing patterns had sufficient information to allow RE identification of the input combinations. Like other neural systems, the motoneuron input-output system is likely non-unique. This non-uniqueness could potentially limit this RE approach, as many input combinations can produce similar outputs. However, our simulations revealed that firing patterns contained sufficient information to sharply restrict the solution space. Thus, our RE approach successfully generated estimates of the actual simulated patterns of excitation, inhibition, and neuromodulation, with variances accounted for ranging from 75–90%. It was striking that nonlinearities induced in firing patterns by the neuromodulation inputs did not impede RE, but instead generated distinctive features in firing patterns that aided RE. These simulations demonstrate the potential of this form of RE analysis. It is likely that the ever-increasing capacity of supercomputers will allow increasingly accurate RE of neuron inputs from their firing patterns from many neural systems.