Thalamocortical Model Research Articles

Rapid context inference in a thalamocortical model using recurrent neural networks

Cognitive flexibility is a fundamental ability that enables humans and animals to exhibit appropriate behaviors in various contexts. The thalamocortical interactions between the prefrontal cortex (PFC) and the mediodorsal thalamus (MD) have been identified as crucial for inferring temporal context, a critical component of cognitive flexibility. However, the neural mechanism responsible for context inference remains unknown. To address this issue, we propose a PFC-MD neural circuit model that utilizes a Hebbian plasticity rule to support rapid, online context inference. Specifically, the model MD thalamus can infer temporal contexts from prefrontal inputs within a few trials. This is achieved through the use of PFC-to-MD synaptic plasticity with pre-synaptic traces and adaptive thresholding, along with winner-take-all normalization in the MD. Furthermore, our model thalamus gates context-irrelevant neurons in the PFC, thus facilitating continual learning. We evaluate our model performance by having it sequentially learn various cognitive tasks. Incorporating an MD-like component alleviates catastrophic forgetting of previously learned contexts and demonstrates the transfer of knowledge to future contexts. Our work provides insight into how biological properties of thalamocortical circuits can be leveraged to achieve rapid context inference and continual learning.

Thalamocortical contributions to hierarchical cognitive control.

Cognitive flexibility relies on hierarchically structured task representations that organize task contexts, relevant environmental features, and subordinate decisions. Despite ongoing interest in the human thalamus, its role in cognitive control has been understudied. This study explored thalamic representation and thalamocortical interactions that contribute to hierarchical cognitive control in humans. We found that several thalamic nuclei, including the anterior, mediodorsal, ventrolateral, and pulvinar nuclei, exhibited stronger evoked responses when subjects switch between task contexts. Decoding analysis revealed that thalamic activity encodes task contexts within the hierarchical task representations. To determine how thalamocortical interactions contribute to task representations, we developed a thalamocortical functional interaction model to predict task-related cortical representation. This data-driven model outperformed comparison models, particularly in predicting activity patterns in cortical regions that encode context representations. Collectively, our findings highlight the significant contribution of thalamic activity and thalamocortical interactions for contextually guided hierarchical cognitive control.

Emergent effects of synaptic connectivity on the dynamics of global and local slow waves in a large-scale thalamocortical network model of the human brain.

Slow-wave sleep (SWS), characterized by slow oscillations (SOs, <1Hz) of alternating active and silent states in the thalamocortical network, is a primary brain state during Non-Rapid Eye Movement (NREM) sleep. In the last two decades, the traditional view of SWS as a global and uniform whole-brain state has been challenged by a growing body of evidence indicating that SO can be local and can coexist with wake-like activity. However, the mechanisms by which global and local SOs arise from micro-scale neuronal dynamics and network connectivity remain poorly understood. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and SWS, and we investigated the role of connectivity in the spatio-temporal dynamics of sleep SO. We found that the overall strength and a relative balance between long and short-range synaptic connections determined the network state. Importantly, for a range of synaptic strengths, the model demonstrated complex mixed SO states, where periods of synchronized global slow-wave activity were intermittent with the periods of asynchronous local slow-waves. An increase in the overall synaptic strength led to synchronized global SO, while a decrease in synaptic connectivity produced only local slow-waves that would not propagate beyond local areas. These results were compared to human data to validate probable models of biophysically realistic SO. The model producing mixed states provided the best match to the spatial coherence profile and the functional connectivity estimated from human subjects. These findings shed light on how the spatio-temporal properties of SO emerge from local and global cortical connectivity and provide a framework for further exploring the mechanisms and functions of SWS in health and disease.

Possible mechanisms to improve sleep spindles via closed loop stimulation during slow wave sleep: A computational study.

Sleep spindles are one of the prominent EEG oscillatory rhythms of non-rapid eye movement sleep. In the memory consolidation, these oscillations have an important role in the processes of long-term potentiation and synaptic plasticity. Moreover, the activity (spindle density and/or sigma power) of spindles has a linear association with learning performance in different paradigms. According to the experimental observations, the sleep spindle activity can be improved by closed loop acoustic stimulations (CLAS) which eventually improve memory performance. To examine the effects of CLAS on spindles, we propose a biophysical thalamocortical model for slow oscillations (SOs) and sleep spindles. In addition, closed loop stimulation protocols are applied on a thalamic network. Our model results show that the power of spindles is increased when stimulation cues are applied at the commencing of an SO Down-to-Up-state transition, but that activity gradually decreases when cues are applied with an increased time delay from this SO phase. Conversely, stimulation is not effective when cues are applied during the transition of an Up-to-Down-state. Furthermore, our model suggests that a strong inhibitory input from the reticular (RE) layer to the thalamocortical (TC) layer in the thalamic network shifts leads to an emergence of spindle activity at the Up-to-Down-state transition (rather than at Down-to-Up-state transition), and the spindle frequency is also reduced (8-11 Hz) by thalamic inhibition.

Emergent effects of synaptic connectivity on the dynamics of global and local slow waves in a large-scale thalamocortical network model of the human brain.

Slow-wave sleep (SWS), characterized by slow oscillations (SO, <1Hz) of alternating active and silent states in the thalamocortical network, is a primary brain state during Non-Rapid Eye Movement (NREM) sleep. In the last two decades, the traditional view of SWS as a global and uniform whole-brain state has been challenged by a growing body of evidence indicating that SO can be local and can coexist with wake-like activity. However, the understanding of how global and local SO emerges from micro-scale neuron dynamics and network connectivity remains unclear. We developed a multi-scale, biophysically realistic human whole-brain thalamocortical network model capable of transitioning between the awake state and slow-wave sleep, and we investigated the role of connectivity in the spatio-temporal dynamics of sleep SO. We found that the overall strength and a relative balance between long and short-range synaptic connections determined the network state. Importantly, for a range of synaptic strengths, the model demonstrated complex mixed SO states, where periods of synchronized global slow-wave activity were intermittent with the periods of asynchronous local slow-waves. Increase of the overall synaptic strength led to synchronized global SO, while decrease of synaptic connectivity produced only local slow-waves that would not propagate beyond local area. These results were compared to human data to validate probable models of biophysically realistic SO. The model producing mixed states provided the best match to the spatial coherence profile and the functional connectivity estimated from human subjects. These findings shed light on how the spatio-temporal properties of SO emerge from local and global cortical connectivity and provide a framework for further exploring the mechanisms and functions of SWS in health and disease.

Research on the pathogenesis of Alzheimer's disease based on thalamocortical computational model.

Alpha rhythm slowing is an important electroencephalogram(EEG) feature associated with (AD). This study aims to understand the correlation between alpha band deceleration and molecular changes from the perspective of neural computing. Considering the effect of Aβ amyloid deposition on the inhibitory changes in the thalamic, a thalamic cortical model coupled with Aβ amyloid is established. The results show that Aβ amyloid deposition may induce neurotoxicity in thalamic reticular nucleus neurons, which results in inhibitory changes in the thalamus and slows the alpha rhythm of EEG output from the thalamus. In order to understand the pathogenesis more intuitively, some numerical simulations are provided to illustrate the obtained theories. This research is helpful to understand the pathogenesis of AD, so as to provide theoretical basis for the intervention and control of the disease.

Neurophysiological evidence that frontoparietal connectivity and GABA-A receptor changes underpin the antidepressant response to ketamine

Revealing the acute cortical pharmacodynamics of an antidepressant dose of ketamine in humans with depression is key to determining the specific mechanism(s) of action for alleviating symptoms. While the downstream effects are characterised by increases in plasticity and reductions in depressive symptoms—it is the acute response in the brain that triggers this cascade of events. Computational modelling of cortical interlaminar and cortico-cortical connectivity and receptor dynamics provide the opportunity to interrogate this question using human electroencephalography (EEG) data recorded during a ketamine infusion. Here, resting-state EEG was recorded in a group of 30 patients with major depressive disorder (MDD) at baseline and during a 0.44 mg/kg ketamine dose comprising a bolus and infusion. Fronto-parietal connectivity was assessed using dynamic causal modelling to fit a thalamocortical model to hierarchically connected nodes in the medial prefrontal cortex and superior parietal lobule. We found a significant increase in parietal-to-frontal AMPA-mediated connectivity and a significant decrease in the frontal GABA time constant. Both parameter changes were correlated across participants with the antidepressant response to ketamine. Changes to the NMDA receptor time constant and inhibitory intraneuronal input into superficial pyramidal cells did not survive correction for multiple comparisons and were not correlated with the antidepressant response. These results provide evidence that the antidepressant effects of ketamine may be mediated by acute fronto-parietal connectivity and GABA receptor dynamics. Furthermore, it supports the large body of literature suggesting the acute mechanism underlying ketamine’s antidepressant properties is related to GABA-A and AMPA receptors rather than NMDA receptor antagonism.

Regulatory role of excitatory interneurons by combining electrical stimulation for absence seizures in the coupled thalamocortical model

&lt;abstract&gt;&lt;p&gt;The role of excitatory interneurons (EINs) in the cortical has received increasing attention in the discussion of absence seizures. Numerous physiological experiments have confirmed the correlation between EIN and absence seizures. However, the dynamic mechanisms underlying this relationship are not well understood, and there are some challenges in selecting appropriate stimulation strategies for pyramidal clusters. In this study, we incorporated EIN into the previous Taylor model and developed an improved thalamocortical coupled model consisting of ten neuronal populations. Initially, we investigated the excitatory induction effect of EIN to pyramidal clusters and the external input of EIN. Then, four different targeted treatment approaches (deep brain stimulation (DBS), current balanced biphasic pulse (CBBP), 1:0 coordinated resetting stimulation (1:0 CRS), and 3:2 CRS) were applied to the pyramidal clusters. Moreover, we established two quantitative indices to evaluate the stimulation effects. The results showed that modifying the external input of EIN and the coupling strength projected onto the pyramidal clusters can effectively transition the system from an absence seizure state to other normal states. Additionally, inputs from the left compartment were found to reduce the generation of abnormal discharge regions in the right compartment. Furthermore, considering the treatment effects and current consumption, the 3:2 CRS stimulation strategy appeared to be the most suitable treatment approach for the pyramidal clusters. This work introduces a novel coupled model containing EIN, which contributes new theoretical foundations and insights for the future treatment of absence seizures.&lt;/p&gt;&lt;/abstract&gt;

Diencephalic organoids - A key to unraveling development, connectivity, and pathology of the human diencephalon.

The diencephalon, an integral component of the forebrain, governs a spectrum of crucial functions, ranging from sensory processing to emotional regulation. Yet, unraveling its unique development, intricate connectivity, and its role in neurodevelopmental disorders has long been hampered by the scarcity of human brain tissue and ethical constraints. Recent advancements in stem cell technology, particularly the emergence of brain organoids, have heralded a new era in neuroscience research. Although most brain organoid methodologies have hitherto concentrated on directing stem cells toward telencephalic fates, novel techniques now permit the generation of region-specific brain organoids that faithfully replicate precise diencephalic identities. These models mirror the complexity of the human diencephalon, providing unprecedented opportunities for investigating diencephalic development, functionality, connectivity, and pathophysiology in vitro. This review summarizes the development, function, and connectivity of diencephalic structures and touches upon developmental brain disorders linked to diencephalic abnormalities. Furthermore, it presents current diencephalic organoid models and their applications in unraveling the intricacies of diencephalic development, function, and pathology in humans. Lastly, it highlights thalamocortical assembloid models, adept at capturing human-specific aspects of thalamocortical connections, along with their relevance in neurodevelopmental disorders.

Modeling Synaptic Integration of Bursty and β Oscillatory Inputs in Ventromedial Motor Thalamic Neurons in Normal and Parkinsonian States.

The ventromedial motor thalamus (VM) is implicated in multiple motor functions and occupies a central position in the cortico-basal ganglia-thalamocortical loop. It integrates glutamatergic inputs from motor cortex (MC) and motor-related subcortical areas, and it is a major recipient of inhibition from basal ganglia. Previous in vitro experiments performed in mice showed that dopamine depletion enhances the excitability of thalamocortical (TC) neurons in VM due to reduced M-type potassium currents. To understand how these excitability changes impact synaptic integration in vivo, we constructed biophysically detailed mouse VM TC model neurons fit to normal and dopamine-depleted conditions, using the NEURON simulator. These models allowed us to assess the influence of excitability changes with dopamine depletion on the integration of synaptic inputs expected in vivo We found that VM neuron models in the dopamine-depleted state showed increased firing rates with the same synaptic inputs. Synchronous bursting in inhibitory input from the substantia nigra pars reticulata (SNR), as observed in parkinsonian conditions, evoked a postinhibitory firing rate increase with a longer duration in dopamine-depleted than control conditions, due to different M-type potassium channel densities. With β oscillations in the inhibitory inputs from SNR and the excitatory inputs from cortex, we observed spike-phase locking in the activity of the models in normal and dopamine-depleted states, which relayed and amplified the oscillations of the inputs, suggesting that the increased β oscillations observed in VM of parkinsonian animals are predominantly a consequence of changes in the presynaptic activity rather than changes in intrinsic properties.

Regulating epileptiform discharges by heterogeneous interneurons in thalamocortical model.

Inhibitory interneurons in the cortex are abundant and have diverse roles, classified as parvalbumin (PV), somatostatin (SOM), and vasoactive intestinal polypeptide (VIP) according to chemically defined categories. Currently, their involvement with seizures has been partially uncovered in physiological terms. Here, we propose a corticothalamic model containing heterogeneous interneurons to study the effects of various interneurons on absence seizure dynamics by means of optogenetic stimulation. First, the important role of feedforward inhibition caused by SRN→PV→PN projections on seizures is verified. Then, we demonstrate that light activation targeting either PV or SOM INs can control seizures. Finally, with different inhibition contributions from PV INs and SOM INs, the possible disinhibitory effect of blue light acting on VIP INs is mainly discussed. The results suggest that depending on the inhibition degree of both types, the disinhibition brought about by the VIP INs will trigger seizures, will control seizures, and will not work or cause the PNs to tend toward a high saturation state with high excitability. The circuit mechanism and the related bifurcation characteristics in various cases are emphatically revealed. In the model presented, in addition to Hopf and saddle-node bifurcations, the system may also undergo period-doubling and torus bifurcations under stimulus action, with more complex dynamics. Our work may provide a theoretical basis for understanding and further exploring the role of heterogeneous interneurons, in particular, the VIP INs, a novel target, in absence seizures.

The role of age-related sleep EEG changes in memory decline: experiments and computational modeling.

The slow oscillation (SO) observed during deep sleep is known to facilitate memory consolidation. However, the impact of age-related changes in sleep electroencephalography (EEG) oscillations and memory remains unknown. In this study, we aimed to investigate the contribution of age-related changes in sleep SO and its role in memory decline by combining EEG recordings and computational modeling. Based on the detected SO events, we found that older adults exhibit lower SO density, lower SO frequency, and longer Up and Down state durations during N3 sleep compared to young and middle-aged groups. Using a biophysically detailed thalamocortical network model, we simulated the "aged" brain as a partial loss of synaptic connections between neurons in the cortex. Our simulations showed that the changes in sleep SO properties in the "aged" brain, similar to those observed in older adults, resulting in impaired memory consolidation. Overall, this study provides mechanistic insights into how age-related changes modulate sleep SOs and memory decline.Clinical Relevance- This study contributes towards finding feasible biomarkers and target mechanism for designing therapy in older adults with memory deficits, such as Alzheimer's disease patients.

Propagation Effect of Epileptic Seizures in a Coupled Thalamocortical Network

Epileptic seizures have spatial features related to the propagation of seizure waves. As the main characteristic of absence seizures, 2–4[Formula: see text]Hz spike-wave discharges (SWDs) originate from the cortices and are maintained by the thalamus. In this study, we explore the onset and propagation effect of absence seizures based on a thalamocortical model. First, we develop a two-compartment model and consider the autapse of the thalamic reticular nucleus as a crucial parameter to investigate transition behaviors. Moreover, we present dynamical mechanisms through bifurcation analysis. Simulation results show that the absence seizures can be induced and advanced as the coupling strength increases. Second, we investigate excitatory and inhibitory coupling functions in a three-compartment model. Our research indicates that the excitatory coupling function can lead to SWDs when all the compartments are initially saturated. In the process of propagation, excitatory coupling also gives rise to SWDs in normal compartments, whereas inhibitory coupling plays a limited role. Finally, we reproduce the above results in a 10-compartment model and verify the robustness against the variation of the number of modules. This work may shed new light on the field of seizure propagation and provide potential dynamical mechanisms.

Dynamic mechanism of epileptic seizures induced by excitatory pyramidal neuronal population

&lt;abstract&gt;&lt;p&gt;The pyramidal neuronal population (PY) in the cerebral cortex is closely related to epilepsy, while the excitability of PY is directly affected by the excitatory interneurons (EIN), the inhibitory interneurons (IN), and the thalamic relay nucleus (TC). Based on this, we use the thalamocortical neural field model to explore the dynamic mechanism of system transition by taking the synaptic connection strengths of the above three nuclei on PY as the main factor affecting seizures. The results show that the excitatory effects of EIN on PY induce transitions from 1-spike and wave discharges (SWDs) to 2-spike and wave discharges (2-SWDs), the inhibitory effects of IN on PY induce transitions from saturated state to tonic oscillation state, and the excitatory effects of TC on PY induce transitions from clonic oscillation state to saturated state. According to the single-parameter bifurcation analysis, it is found that Hopf and fold limit cycle bifurcations are the key factors leading to the state transition. In addition, the state analysis of the three pathways is carried out in pairs. The results show that the system produces more types of epileptic seizures with the combined action of EIN and TC on PY. According to the two-parameter bifurcation curve, we obtain the stable parameter areas of tonic-clonic oscillations, SWDs, 2-SWDs and saturated discharges, and clearly find the reasonable transition path between tonic-clonic seizures and absence seizures. This may provide some theoretical guidance for the transmission and evolution of seizures.&lt;/p&gt;&lt;/abstract&gt;

Thalamocortical contributions to cognitive task activity.

Thalamocortical interaction is a ubiquitous functional motif in the mammalian brain. Previously (Hwang et al., 2021), we reported that lesions to network hubs in the human thalamus are associated with multi-domain behavioral impairments in language, memory, and executive functions. Here, we show how task-evoked thalamic activity is organized to support these broad cognitive abilities. We analyzed functional magnetic resonance imaging (MRI) data from human subjects that performed 127 tasks encompassing a broad range of cognitive representations. We first investigated the spatial organization of task-evoked activity and found a basis set of activity patterns evoked to support processing needs of each task. Specifically, the anterior, medial, and posterior-medial thalamus exhibit hub-like activity profiles that are suggestive of broad functional participation. These thalamic task hubs overlapped with network hubs interlinking cortical systems. To further determine the cognitive relevance of thalamic activity and thalamocortical functional connectivity, we built a data-driven thalamocortical model to test whether thalamic activity can be used to predict cortical task activity. The thalamocortical model predicted task-specific cortical activity patterns, and outperformed comparison models built on cortical, hippocampal, and striatal regions. Simulated lesions to low-dimensional, multi-task thalamic hub regions impaired task activity prediction. This simulation result was further supported by profiles of neuropsychological impairments in human patients with focal thalamic lesions. In summary, our results suggest a general organizational principle of how the human thalamocortical system supports cognitive task activity.

Differential thalamocortical interactions in slow and fast spindle generation: A computational model.

Cortical slow oscillations (SOs) and thalamocortical sleep spindles are two prominent EEG rhythms of slow wave sleep. These EEG rhythms play an essential role in memory consolidation. In humans, sleep spindles are categorized into slow spindles (8-12 Hz) and fast spindles (12-16 Hz), with different properties. Slow spindles that couple with the up-to-down phase of the SO require more experimental and computational investigation to disclose their origin, functional relevance and most importantly their relation with SOs regarding memory consolidation. To examine slow spindles, we propose a biophysical thalamocortical model with two independent thalamic networks (one for slow and the other for fast spindles). Our modeling results show that fast spindles lead to faster cortical cell firing, and subsequently increase the amplitude of the cortical local field potential (LFP) during the SO down-to-up phase. Slow spindles also facilitate cortical cell firing, but the response is slower, thereby increasing the cortical LFP amplitude later, at the SO up-to-down phase of the SO cycle. Neither the SO rhythm nor the duration of the SO down state is affected by slow spindle activity. Furthermore, at a more hyperpolarized membrane potential level of fast thalamic subnetwork cells, the activity of fast spindles decreases, while the slow spindles activity increases. Together, our model results suggest that slow spindles may facilitate the initiation of the following SO cycle, without however affecting expression of the SO Up and Down states.

Thalamocortical contributions to working memory processes during the n-back task

Working memory allows individuals to temporally maintain and manipulate information that is no longer accessible from the sensorium. Whereas prior studies have detailed frontoparietal contributions to working memory processes, less emphasis has been placed on subcortical regions, in particular the human thalamus. The thalamus has a complex anatomy that consists of several distinct nuclei, many of which have dense anatomical connectivity with frontoparietal regions, and thus might play an important yet underspecified role for working memory. The goal of our study is to characterize the detailed functional neuroanatomy of the human thalamus and thalamocortical interactions during the n-back task. To that end, we analyzed an n-back fMRI dataset consisting of 395 subjects from the Human Connectome Project (HCP). We found that thalamic nuclei in the anterior, medial, ventral lateral, and posterior medial thalamus showed stronger evoked responses in response to higher working memory load. Activity in most thalamic nuclei were only modulated by working memory load, but not by categorical membership of the memorized stimuli, suggesting that thalamic function supports domain-general processing for working memory. To determine whether thalamocortical interactions contribute to cortical activity for working memory, we employed an activity flow mapping analysis to test whether thalamocortical interactions can predict cortical task activity patterns. In support, this data-driven thalamocortical interaction model explained a significant amount of variance in the observed cortical activity patterns modulated by working memory load. Our results suggest that the anterior, medial, and posterior medial thalamus, and their associated thalamocortical interactions, contribute to the modulations of distributed cortical activity during working memory.

The thalamocortical inhibitory network controls human conscious perception

Although conscious perception is a fundamental cognitive function, its neural correlates remain unclear. It remains debatable whether thalamocortical interactions play a decisive role in conscious perception. To clarify this, we used functional magnetic resonance imaging (fMRI) where flickering red and green visual cues could be perceived either as a non-fused colour or fused colour. Here we show significantly differentiated fMRI neurodynamics only in higher-order thalamocortical regions, compared with first-order thalamocortical regions. Anticorrelated neurodynamic behaviours were observed between the visual stream network and default-mode network. Its dynamic causal modelling consistently provided compelling evidence for the involvement of higher-order thalamocortical iterative integration during conscious perception of fused colour, while inhibitory control was revealed during the non-fusion condition. Taken together with our recent magnetoencephalography study, our fMRI findings corroborate a thalamocortical inhibitory model for consciousness, where both thalamic inhibitory regulation and integrative signal iterations across higher-order thalamocortical regions are essential for conscious perception.

Thalamic bursting and the role of timing and synchrony in thalamocortical signaling in the awake mouse

The thalamus controls transmission of sensory signals from periphery to cortex, ultimately shaping perception. Despite this significant role, dynamic thalamic gating and the consequences for downstream cortical sensory representations have not been well studied in the awake brain. We optogenetically modulated the ventro-posterior-medial thalamus in the vibrissa pathway of the awake mouse and measured spiking activity in the thalamus and activity in primary somatosensory cortex (S1) using extracellular electrophysiology and genetically encoded voltage imaging. Thalamic hyperpolarization significantly enhanced thalamic sensory-evoked bursting; however, surprisingly, the S1 cortical response was not amplified, but instead, timing precisionwas significantly increased, spatial activation more focused, and there was an increased synchronization of cortical inhibitory neurons. A thalamocortical network model implicates the modulation of precise timing of feedforward thalamic population spiking, presenting a highly sensitive, timing-based gating of sensory signaling to the cortex.

Role of Sleep in Formation of Relational Associative Memory.

Relational memory, the ability to make and remember associations between objects, is an essential component of mammalian reasoning. In relational memory tasks, it has been shown that periods of offline processing, such as sleep, are critical to making indirect associations. To understand biophysical mechanisms behind the role of sleep in improving relational memory, we developed a model of the thalamocortical network to test how slow-wave sleep affects performance on an unordered relational memory task. First, the model was trained in the awake state on a paired associate inference task, in which the model learned to recall direct associations. After a period of subsequent slow-wave sleep, the model developed the ability to recall indirect associations. We found that replay, during sleep, of memory patterns learned in awake increased synaptic connectivity between neurons representing the item that was overlapping between tasks and neurons representing the unlinked items of the different tasks; this forms an attractor that enables indirect memory recall. Our study predicts that overlapping items between indirectly associated tasks are essential for relational memory, and sleep can reactivate pathways to and from overlapping items to the unlinked objects to strengthen these pathways and form new relational memories.SIGNIFICANCE STATEMENT Experimental studies have shown that some types of associative memory, such as transitive inference and relational memory, can improve after sleep. Still, it remains unknown what specific mechanisms are responsible for these sleep-related changes. In this new work, we addressed this problem by building a thalamocortical network model that can learn relational memory tasks and that can be simulated in awake or sleep states. We found that memory traces learned in awake were replayed during slow waves of NREM sleep and revealed that replay increased connections to and from overlapping memory items to form new relational memories. Our work discovered specific mechanisms behind the role of sleep in associative memory and made testable predictions about how sleep augments associative learning.