Unveiling serotonergic dysfunction of obsessive-compulsive disorder on prefrontal network dynamics: a computational perspective.

Serotonin (5-HT) regulates working memory within the prefrontal cortex network, which is crucial for understanding obsessive-compulsive disorder. However, the mechanisms how network dynamics and serotonin interact in obsessive-compulsive disorder remain elusive. Here, we incorporate 5-HT receptors (5-HT1A, 5-HT2A) and dopamine receptors into a multistable prefrontal cortex network model, replicating the experimentally observed inverted U-curve phenomenon. We show how the two 5-HT receptors antagonize neuronal activity and modulate network multistability. Reduced binding of 5-HT1A receptors increases global firing, while reduced binding of 5-HT2A receptors deepens attractors. The obtained results suggest reward-dependent synaptic plasticity mechanisms may attenuate 5-HT related network impairments. Integrating serotonin-mediated dopamine release into circuit, we observe that decreased serotonin concentration triggers the network into a deep attractor state, expanding the domain of attraction of stable nodes with high firing rate, potentially causing aberrant reverse learning. This suggests a hypothesis wherein elevated dopamine concentrations in obsessive-compulsive disorder might result from primary deficits in serotonin levels. Findings of this work underscore the pivotal role of serotonergic dysregulation in modulating synaptic plasticity through dopamine pathways, potentially contributing to learned obsessions. Interestingly, serotonin reuptake inhibitors and antidopaminergic potentiators can counteract the over-stable state of high-firing stable points, providing new insights into obsessive-compulsive disorder treatment.

A large-scale neuronal network modelling study: Stimulus size modulates gamma oscillations in the primary visual cortex by long-range connections.

Stimulus size modulation of neuronal firing activity is a fundamental property of the primary visual cortex. Numerous biological experiments have shown that stimulus size modulation is affected by multiple factors at different spatiotemporal scales, but the exact pathways and mechanisms remain incompletely understood. In this paper, we establish a large-scale neuronal network model of primary visual cortex with layer 2/3 to study how gamma oscillation properties are modulated by stimulus size and especially how long-range connections affect the modulation as realistic neuronal properties and spatial distributions of synaptic connections are considered. It is shown that long-range horizontal synaptic connections are sufficient to produce dimensional modulation of firing rates and gamma oscillations. In particular, with increasing grating stimulus size, the firing rate increases and then decreases, the peak frequency of gamma oscillations decreases and the spectral power increases. These are consistent with biological experimental observations. Furthermore, we explain in detail how the number and spatial distribution of long-range connections affect the size modulation of gamma oscillations by using the analysis of neuronal firing activity and synaptic current fluctuations. Our results provide a mechanism explanation for size modulation of gamma oscillations in the primary visual cortex and reveal the important and unique role played by long-range connections, which contributes to a deeper understanding of the cognitive function of gamma oscillations in visual cortex.

Retrospective cohort evaluation of non-HIV Castleman disease from a single academic center in Beijing, China.

The purpose is to ascertain the clinical impact of Castleman disease (CD) by reassessment of the real-world data from Peking University First Hospital (PKUFH). The results will contribute to the standardization of diagnosis and treatment on CDs. Based on the last 15-year retrospective real-world data from Peking University First Hospital (PKUFH), we reclassified and re-evaluated the clinical and pathological information of patients with pathologically suspected diagnosis of CD. A total of 203 patients were included in our study, in which the diagnosis of CD was confirmed in 189 cases, including 118 patients with unicentric CD (UCD, n = 118, 62.4%) and 71 patients with multicentric CD (MCD, n = 71, 37.6%). A total of 44.1% (n = 52) of UCDs in our cohort were complicated with paraneoplastic pemphigus (PNP). The treatment of UCD is primarily surgical, with a 5-year overall survival (OS) of 88.1%. Patients with PNP had a poorer prognosis than those without PNP (82.9% (95% CI 123-178) vs 92.8% (95% CI 168-196), log-rank p = 0.041). The rate of concurrent systemic symptoms was 74.6% (n = 53), and renal involvement occurred in 49.3% (n = 35) MCD patients. The MCD treatments were mainly chemotherapy regimens, with a 5-year OS of 77.6% (95% CI, 143-213). Patients with UCD demonstrate a better overall prognosis than patients with MCD. But the prognosis of those complicated with PNP was poor. The differential diagnosis of MCD is extensive. MCD treatment in China is heterogeneous. The inaccessibility of anti-IL-6-targeted drugs in China may contribute to the poor prognosis for patients with MCD.A preprint has previously been published (Guo et al. 34).

Theoretical Study on the O-H Fracture of Methanol on PtnCu4-n (n = 1, 2, 3) Catalysts with Different Coverages.

Direct methanol fuel cells (DMFCs) have attracted increasing attention as a very promising and important energy source. In this paper, density functional theory (DFT) is used to study the structure and O-H fracture mechanism of methanol adsorption on PtnCu4-n (111) (n = 1, 2, 3) binary metal catalyst surfaces under different coverages. By comparing the adsorption energy and dehydrogenation energy barriers of methanol, it is found that the adsorption strength and dehydrogenation energy barriers of methanol on Pt and Cu sites decreased with increasing coverage. At the same Pt and Cu ratio, methanol is more easily adsorbed on Cu sites. When Pt/Cu = 3:1 and 1:3, the PtCu binary catalyst has a significant impact on the energy barrier of breaking the O-H bond in methanol with the increase of coverage. Especially when Pt/Cu = 1:3 and the coverage is 1/4 ML, the energy barriers of O-H bond breaking in methanol on Pt and Cu sites are 0.63 and 0.61 eV, respectively, which are lower than that on pure Pt. It means that the Cu sites played a very important role in reducing the O-H fracture energy barrier of methanol. When Pt/Cu = 1:1, the change in the dehydrogenation energy barrier of methanol on Pt sites and Cu sites is not significant, indicating that the coverage has little effect on it.

Application of fecal collection device in Intensive Care Unit Patients with fecal incontinence Receiving Extracorporeal Membrane Oxygenation: A Comparison Cohort Study.

Objective: The objective of this study was to assess the effectiveness of fecal collection devices in preventing incontinence-associated dermatitis (IAD) and reducing skin care time in ICU patients with fecal incontinence undergoing Extracorporeal Membrane Oxygenation (ECMO). Methods: A nonrandomized comparison cohort (quasi-experimental) study with pre-post comparison was carried out in a general intensive care unit. 85 bedridden patients receiving ECMO with fecal incontinence (FI) in a general intensive care unit between June 2017 and May 2022 participated in the study and separated into two groups according to the fecal collection device they received. 40 were assigned to the Control group (structured IAD preventive care protocol alone) and 45 to the Intervention group (structured IAD preventive care protocol plus application of fecal collection device). The status of IAD was assessed using the Incontinence Associated Dermatitis Intervention Tool (IAD-IT). Fecal consistency was evaluated via the Bristol Stool Scale. Outcome measures included the nursing time for skin care and the incidence of IAD, and bleeding complications between the two groups during the period. Results: Participants in the Intervention group had fewer IAD occurrences than participants in the Control group (13.33% vs. 52.50%, P < .05). The patients in the Intervention group significantly reduced skincare time (63.30±14.09 min in the Control group versus 28.44±2.04 min in the Intervention group, P < .01). There was 3 turning complications for bleeding in the Intervention group and 11 in the Control group and had a significant reduction in urning complications(3 vs.11, P = .022). Conclusions: Applying a fecal collection device may reduce skincare time and reduce occurrences of IAD and bleeding related to turning position for skin care in ICU patients with FI associated with diarrhea receiving ECMO Support.This study offers a more efficient way to use the fecal collection device in ECMO patients.Future research needs to focus on the perianal skin in ECMO patients regarding fecal collection devices connected to continuous low-negative-pressure suction devices.

Functional identification of anthocyanin glucosyltransferase genes: a Ps3GT catalyzes pelargonidin to pelargonidin 3-O-glucoside painting the vivid red flower color of Paeonia.

MAIN CONCLUSION: Glycosylation from an anthocyanidin 3-O-glucosyltransferase Ps3GT (PsUGT78A27) facilitates the accumulation of pelargonidin 3-O-glucoside, which defines the vivid red flower color and occurs only in specific peony tree cultivars. Although tree peony cultivars of Chinese and Japanese both originated from China, vivid red color is only found in flowers of Japanese cultivars but not of Chinese cultivar groups. In this study, a Japanese tree peony cultivar 'Taiyoh' with vivid red petals and a Chinese tree peony cultivar 'Hu Hong' with reddish pink petals were chosen as the experimental materials. Flavonoids profiling indicated that pelargonidin 3-O-glucoside (Pg3G) detected only in Japanese cultivar contributed to vivid red color of tree peony petals, while pelargonidin 3,5-di-O-glucoside (Pg3G5G) found in both of Japanese and Chinese cultivars was responsible for pink flower color. Through the integration of full-length transcriptome sequencing and in vitro enzymatic activity analysis, two anthocyanin glucosyltransferase genes PsUGT78A27 and PsUGT75L45 were isolated from the petals of tree peony, and their encoding products exhibited enzymatic activities of pelargonidin 3-O-glucosyltransferase and anthocyanin 5-O-glucosyltransferase, respectively. Further quantitative real-time PCR revealed that PsUGT78A27 displayed high expression in petals of both cultivars and PsUGT75L45 was expressed at high levels in cultivar 'Hu Hong' only. Using a gene gun technique, the GFP fusion proteins of PsUGT78A27 and PsUGT75L45 were visualized to be cytoplasmic and nuclear localization in the epidermal cells of tree peony petals, and the glucosylation function of PsUGT78A27 and PsUGT75L45 to alter petal color of tree peony and herbaceous peony had been directly validated in vivo. These results demonstrated that PsUGT78A27 and PsUGT75L45 are key players for the presence or absence of vivid red flower color in tree peony cultivars. Our findings further elucidated the chemical and molecular mechanism of petal pigmentation of Paeonia and could help breed the Paeonia cultivars possessing novel flower colors.

Real-world outcomes and prognostic factors among patients with acute myeloid leukemia treated with allogeneic hematopoietic stem cell transplantation.

Acute myeloid leukemia (AML) is a malignant lymphohematopoietic tumor that ranks among the most frequent indications for allogeneic hematopoietic stem cell transplantation (allo-HSCT). This article aims to provide a comprehensive analysis of the application of allo-HSCT for AML and identify prognostic factors to enhance future treatment effect. This retrospective study collected data from 323 patients diagnosed with AML at Peking University First Hospital who underwent allo-HSCT between September 2003 and July 2022. The annual number of transplantations has steadily increased. Our center has observed a rise in the proportion of cytogenetic high-risk and measurable residual disease (MRD) positive patients since 2013, as well as an increase in the number of haploidentical transplantations. The overall leukocyte engraftment time has decreased over the past 20 years. Furthermore, both overall survival (OS) and disease-free survival (DFS) have significantly improved, while non-relapse mortality (NRM) has significantly decreased since 2013. Multivariate analysis identified transplantation before 2013, patients in complete remission (CR) 2 or non-CR, and recipients older than 50 years as risk factors for NRM, while patients in non-CR and patients with positive MRD are risk factors for recurrence. These findings offer insights into AML treatment outcomes in China, highlighting changes in transplantation practices and the need to reduce post-transplant relapse. Effective interventions, such as MRD monitoring and risk stratification schemes, are crucial for further enhancing transplant outcomes.

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.

Correlation between Flow Temperature and Average Molar Ionic Potential of Ash during Gasification of Coal and Phosphorus-Rich Biomass.

The co-gasification of biomass and coal is helpful for achieving the clean and efficient utilization of phosphorus-rich biomass. A large number of alkali and alkaline earth metals (AAEMs) present in the ash system of coal (or biomass) cause varying degrees of ash, slagging, and corrosion problems in the entrained flow gasifier. Meanwhile, phosphorus is present in the slag in the form of PO43-, which has a strong affinity for AAEMs (especially for Ca2+) to produce minerals dominated by calcium phosphates or alkaline Ca-phosphate, effectively mitigating the aforementioned problems. To investigate the changing behavior of the slag flow temperature (FT) under different CaO/P2O5 ratios, 72 synthetic ashes with varying CaO/P2O5 ratios at different Si/Al contents and compositions were prepared, and their ash fusion temperatures were tested. The effects of different CaO/P2O5 ratios on the FT were analyzed using FactSage thermodynamic simulation. A model for predicting slag FT at different CaO/P2O5 ratios was constructed on the basis of the average molar ionic potential (Ia) method and used to predict data reported from 19 mixed ashes in the literature. The results showed that Ia and FT gradually increased with a decreasing CaO/P2O5 ratio, and the main mineral types shifted from anorthite → mullite → berlinite, which reasonably explained the decrease in ash fusion temperatures in the mixed ash. The established model showed good adaptability to the prediction of 19 actual coal ash FTs in the literature; the deviation of the prediction was in the range of 40 °C. The model proposed between FT and Ia based on the different CaO/P2O5 ratios can be used to predict the low-rank coal and phosphorus-rich biomass and their mixed ashes.

NHC-alcohol adduct-mediated photocatalytic deoxygenation for the synthesis of 6-phenanthridine derivatives.

NHC-alcohol adduct-mediated deoxygenation of alcohols under photocatalytic conditions is described. This process provides various alkyl radicals, which can react with 2-isocyanobiaryls to afford 6-substituted phenanthridine derivatives in moderate to good yields. This method offered the first example of directly using alcohols as radical sources for 6-phenanthridine synthesis.

Impact of Double-Machine Replacement Protocol at Start of Continuous Renal Replacement Therapy in Vasopressor-Dependent Patients: A Retrospective Cohort Study.

INTRODUCTION: When starting continuous renal replacement therapy (CRRT), vasopressor-dependent patients are at risk of hemodynamic instability. Thus far, only a few studies have analyzed the impact of CRRT circuit replacement for vasopressor-dependent patients. Hence, we compared the effect of double-machine replacement protocol (DMRP) with single-machine replacement protocol (SMRP) for CRRT circuit replacement in vasopressor-dependent patients. METHODS: The medical records of 96 vasopressor-dependent patients treated with CRRT in the general intensive care unit of the Shunde Hospital, Southern Medical University, between January 2017 and April 2018 were retrospectively analyzed. The major measures of the SMRP included returning the blood to the patient and sealing access catheter with heparin and starting a new CRRT machine with a slow blood pump, while DMRP involved simultaneous drawing and return of blood with two machines using a slow blood pump for circuit replacement. The primary outcome measures were changes in vasopressor dose and hemodynamic parameters, and the secondary outcome measure was the pause time difference between the two groups during the period. RESULTS: A total of 53 patients were treated with SMRP and 43 patients with DMRP. Heart rate was higher in the SMRP group as compared to the DMRP group (p < 0.05). There were no significant changes in central venous pressure, mean arterial pressure, and vasopressor dose in either group (p > 0.05). The patients in the DMRP group had a significant reduction in CRRT pause time (5.62 ± 0.69 min in DMRP group vs. 37.01 ± 8.72 min in SMRP group, p < 0.01). The DMRP group needed a lower volume of circuit purging and priming fluid related to CRRT circuit replacement (0 mL in DMRP group vs. 463 mL in SMRP group). CONCLUSIONS: Implementation of the DMRP for CRRT circuit replacement had a slight hemodynamic effect on vasopressor-dependent patients. It also reduced the pause time and volume of circuit purging and priming fluid related to CRRT circuit replacement compared with SMRP.

Extracting the transition network of epileptic seizure onset.

In presurgical monitoring, focal seizure onset is visually assessed from intracranial electroencephalogram (EEG), typically based on the selection of channels that show the strongest changes in amplitude and frequency. As epileptic seizure dynamics is increasingly considered to reflect changes in potentially distributed neural networks, it becomes important to also assess the interrelationships between channels. We propose a workflow to quantitatively extract the nodes and edges contributing to the seizure onset using an across-seizure scoring. We propose a quantification of the consistency of EEG channel contributions to seizure onset within a patient. The workflow is exemplified using recordings from patients with different degrees of seizure-onset consistency.

Spontaneous transitions to focal-onset epileptic seizures: A dynamical study.

Given the complex temporal evolution of epileptic seizures, understanding their dynamic nature might be beneficial for clinical diagnosis and treatment. Yet, the mechanisms behind, for instance, the onset of seizures are still unknown. According to an existing classification, two basic types of dynamic onset patterns plus a number of more complex onset waveforms can be distinguished. Here, we introduce a basic three-variable model with two time scales to study potential mechanisms of spontaneous seizure onset. We expand the model to demonstrate how coupling of oscillators leads to more complex seizure onset waveforms. Finally, we test the response to pulse perturbation as a potential biomarker of interictal changes.

Addressing Noise and Skewness in Interpretable Health-Condition Assessment by Learning Model Confidence.

Assessing the health condition has a wide range of applications in healthcare, military, aerospace, and industrial fields. Nevertheless, traditional feature-engineered techniques involve manual feature extraction, which are too cumbersome to adapt to the changes caused by the development of sensor network technology. Recently, deep-learning-based methods have achieved initial success in health-condition assessment research, but insufficient considerations for problems such as class skewness, noisy segments, and result interpretability make it difficult to apply them to real-world applications. In this paper, we propose a K-margin-based Interpretable Learning approach for health-condition assessment. In detail, a skewness-aware RCR-Net model is employed to handle problems of class skewness. Furthermore, we present a diagnosis model based on K-margin to automatically handle noisy segments by naturally exploiting expected consistency among the segments associated with each record. Additionally, a knowledge-directed interpretation method is presented to learn domain knowledge-level features automatically without the help of human experts which can be used as an interpretable decision-making basis. Finally, through experimental validation in the field of both medical and aerospace, the proposed method has a better generality and high efficiency with 0.7974 and 0.8005 F1 scores, which outperform all state-of-the-art deep learning methods for health-condition assessment task by 3.30% and 2.99%, respectively.

Influence of aspartic acid and lysine on the uptake of gold nanoparticles in rice.

The interactions between plants and nanomaterials (NMs) can shed light on the environmental consequences of nanotechnology. We used the major crop plant rice (Oryza sativa L.) to investigate the uptake of gold nanoparticles (GNPs) coated with either negatively or positively charged ligands, over a 5-day period, in the absence or presence of one of two amino acids, aspartic acid (Asp) or lysine (Lys), acting as components of rice root exudates. The presence of Asp or Lys influenced the uptake and distribution of GNPs in rice, which depended on the electrical interaction between the coated GNPs and each amino acid. When the electrical charge of the amino acid was the same as that of the surface ligand coated onto the GNPs, the GNPs could disperse well in nutrient solution, resulting in increased uptake of GNPs into rice tissue. The opposite was true where the charge on the surface ligand was different from that on the amino acid, resulting in agglomeration and reduced Au uptake into rice tissue. The behavior of GNPs in the hydroponic nutrient solution was monitored in terms of agglomeration, particle size distribution, and surface charge in the presence and absence of Asp or Lys, which depended strongly on the electrostatic interaction. Results from this study indicated that the species of root exudates must be taken into account in assessing the bioavailability of nanomaterials to plants.

Synchronous high-frequency oscillations in inhibitory-dominant network motifs consisting of three dentate gyrus-CA3 systems.

Studies on the structural-functional connectomes of the human brain have demonstrated the existence of synchronous firings in a specific brain network motif. In particular, synchronization of high-frequency oscillations (HFOs) has been observed in the experimental data sets of temporal lobe epilepsy (TLE). In addition, both clinical and experimental evidences have accumulated to demonstrate the effect of electrical stimulation on TLE, which, however, remains largely unexplored. In this work, we first employ our previously proposed dentate gyrus (DG)-CA3 network model to investigate the influence of an external electrical stimulus on the HFO transitions. The results indicate that the reinforcing stimulus can induce the HFO transitions of the DG-CA3 system from the gamma band to the fast ripples band. Along with that, the consistent oscillations of neurons within DG-CA3 can also be enhanced with the increasing of stimulus. Then, we expand into a simple motif of three coupled DG-CA3 systems in both the feedforward inhibition and feedback inhibition connections, to investigate the synchronous evolutions of HFOs by regulating both the stimulation strength and inhibitory function. It is shown that the comprehensive effects, which lead to band transition, are independent of the motif configurations. The enhanced external electrical stimulus weakens the synchronism and correlation of connected motifs. In contrast, we demonstrate that the increased inhibitory coupling could facilitate correlation to some extent. Overall, our work highlights the possible origin of synchronous HFOs of hippocampal motifs governed by external inputs and inhibitory connection, which might contribute to a better understanding of the interplay between synchronization dynamics and epileptic structure in the human brain.

Effects of brain-derived neurotrophic factor and noise on transitions in temporal lobe epilepsy in a hippocampal network.

Brain-derived neurotrophic factor (BDNF) has recently been implicated in the modulation of receptor activation leading to dynamic state transitions in temporal lobe epilepsy (TLE). In addition, the crucial role of neuronal noise in these transitions has been studied in electrophysiological experiments. However, the precise role of these factors during seizure generation in TLE is not known. Building on a previously proposed model of an epileptogenic hippocampal network, we included the actions of BDNF-regulated receptors and intrinsic noise. We found that the effects of both BDNF and noise can increase the activation of N-methyl-D-aspartate receptors leading to excessive C a 2 + flux, which induces abnormal fast spiking and bursting. Our results indicate that the combined effects have a strong influence on the seizure-generating network, resulting in higher firing frequency and amplitude. As correlations between firing increase, the synchronization of the entire network increases, a marker of the ictogenic transitions from normal to seizures-like dynamics. Our work on the effects of BDNF dynamics in a noisy environment might lead to an improved model-based understanding of the pathological mechanisms in TLE.

Predicting seizure by modeling synaptic plasticity based on EEG signals - a case study of inherited epilepsy.

This paper explores the internal dynamical mechanisms of epileptic seizures through quantitative modeling based on full brain electroencephalogram (EEG) signals. Our goal is to provide seizure prediction and facilitate treatment for epileptic patients. Motivated by an earlier mathematical model with incorporated synaptic plasticity, we studied the nonlinear dynamics of inherited seizures through a differential equation model. First, driven by a set of clinical inherited electroencephalogram data recorded from a patient with diagnosed Glucose Transporter Deficiency, we developed a dynamic seizure model on a system of ordinary differential equations. The model was reduced in complexity after considering and removing redundancy of each EEG channel. Then we verified that the proposed model produces qualitatively relevant behavior which matches the basic experimental observations of inherited seizure, including synchronization index and frequency. Meanwhile, the rationality of the connectivity structure hypothesis in the modeling process was verified. Further, through varying the threshold condition and excitation strength of synaptic plasticity, we elucidated the effect of synaptic plasticity to our seizure model. Results suggest that synaptic plasticity has great effect on the duration of seizure activities, which support the plausibility of therapeutic interventions for seizure control.

Stimulus-induced transitions between spike-wave discharges and spindles with the modulation of thalamic reticular nucleus.

It is believed that thalamic reticular nucleus (TRN) controls spindles and spike-wave discharges (SWD) in seizure or sleeping processes. The dynamical mechanisms of spatiotemporal evolutions between these two types of activity, however, are not well understood. In light of this, we first use a single-compartment thalamocortical neural field model to investigate the effects of TRN on occurrence of SWD and its transition. Results show that the increasing inhibition from TRN to specific relay nuclei (SRN) can lead to the transition of system from SWD to slow-wave oscillation. Specially, it is shown that stimulations applied in the cortical neuronal populations can also initiate the SWD and slow-wave oscillation from the resting states under the typical inhibitory intensity from TRN to SRN. Then, we expand into a 3-compartment coupled thalamocortical model network in linear and circular structures, respectively, to explore the spatiotemporal evolutions of wave states in different compartments. The main results are: (i) for the open-ended model network, SWD induced by stimulus in the first compartment can be transformed into sleep-like slow UP-DOWN and spindle states as it propagates into the downstream compartments; (ii) for the close-ended model network, weak stimulations performed in the first compartment can result in the consistent experimentally observed spindle oscillations in all three compartments; in contrast, stronger periodic single-pulse stimulations applied in the first compartment can induce periodic transitions between SWD and spindle oscillations. Detailed investigations reveal that multi-attractor coexistence mechanism composed of SWD, spindles and background state underlies these state evolutions. What's more, in order to demonstrate the state evolution stability with respect to the topological structures of neural network, we further expand the 3-compartment coupled network into 10-compartment coupled one, with linear and circular structures, and nearest-neighbor (NN) coupled network as well as its realization of small-world (SW) topology via random rewiring, respectively. Interestingly, for the cases of linear and circular connetivities, qualitatively similar results were obtained in addition to the more irregularity of firings. However, SWD can be eventually transformed into the consistent low-amplitude oscillations for both NN and SW networks. In particular, SWD evolves into the slow spindling oscillations and background tonic oscillations within the NN and SW network, respectively. Our modeling and simulation studies highlight the effect of network topology in the evolutions of SWD and spindling oscillations, which provides new insights into the mechanisms of cortical seizures development.

The pacemaker role of thalamic reticular nucleus in controlling spike-wave discharges and spindles.

Absence epilepsy, characterized by 2-4 Hz spike-wave discharges (SWDs), can be caused by pathological interactions within the thalamocortical system. Cortical spindling oscillations are also demonstrated to involve the oscillatory thalamocortical rhythms generated by the synaptic circuitry of the thalamus and cortex. This implies that SWDs and spindling oscillations can share the common thalamocortical mechanism. Additionally, the thalamic reticular nucleus (RE) is hypothesized to regulate the onsets and propagations of both the epileptic SWDs and sleep spindles. Based on the proposed single-compartment thalamocortical neural field model, we firstly investigate the stimulation effect of RE on the initiations, terminations, and transitions of SWDs. It is shown that the activations and deactivations of RE triggered by single-pulse stimuli can drive the cortical subsystem to behave as the experimentally observed onsets and self-abatements of SWDs, as well as the transitions from 2-spike and wave discharges (2-SWDs) to SWDs. In particular, with increasing inhibition from RE to the specific relay nucleus (TC), rich transition behaviors in cortex can be obtained through the upstream projection path, RE→TC→Cortex. Although some of the complex dynamical patterns can be expected from the earlier single compartment thalamocortical model, the effect of brain network topology on the emergence of SWDs and spindles, as well as the transitions between them, has not been fully investigated. We thereby develop a spatially extended 3-compartment coupled network model with open-/closed-end connective configurations, to investigate the spatiotemporal effect of RE on the SWDs and spindles. Results show that the degrees of activations of RE1 can induce the rich spatiotemporal evolution properties including the propagations from SWDs to spindles within different compartments and the transitions between them, through the RE1→TC1→Cortex1 and Cortex1→Cortex2→Cortex3 projecting paths, respectively. Overall, those results imply that RE possesses the pacemaker function in controlling SWDs and spindling oscillations, which computationally provide causal support for the involvement of RE in absence seizures and sleep spindles.