Synchronization in STDP-driven memristive neural networks with time-varying topology.

Synchronization is a widespread phenomenon in the brain. Despite numerous studies, the specific parameter configurations of the synaptic network structure and learning rules needed to achieve robust and enduring synchronization in neurons driven by spike-timing-dependent plasticity (STDP) and temporal networks subject to homeostatic structural plasticity (HSP) rules remain unclear. Here, we bridge this gap by determining the configurations required to achieve high and stable degrees of complete synchronization (CS) and phase synchronization (PS) in time-varying small-world and random neural networks driven by STDP and HSP. In particular, we found that decreasing P (which enhances the strengthening effect of STDP on the average synaptic weight) and increasing F (which speeds up the swapping rate of synapses between neurons) always lead to higher and more stable degrees of CS and PS in small-world and random networks, provided that the network parameters such as the synaptic time delay [Formula: see text], the average degree [Formula: see text], and the rewiring probability [Formula: see text] have some appropriate values. When [Formula: see text], [Formula: see text], and [Formula: see text] are not fixed at these appropriate values, the degree and stability of CS and PS may increase or decrease when F increases, depending on the network topology. It is also found that the time delay [Formula: see text] can induce intermittent CS and PS whose occurrence is independent F. Our results could have applications in designing neuromorphic circuits for optimal information processing and transmission via synchronization phenomena.

Single-compartment model of a pyramidal neuron, fitted to recordings with current and conductance injection.

For single neuron models, reproducing characteristics of neuronal activity such as the firing rate, amplitude of spikes, and threshold potentials as functions of both synaptic current and conductance is a challenging task. In the present work, we measure these characteristics of regular spiking cortical neurons using the dynamic patch-clamp technique, compare the data with predictions from the standard Hodgkin-Huxley and Izhikevich models, and propose a relatively simple five-dimensional dynamical system model, based on threshold criteria. The model contains a single sodium channel with slow inactivation, fast activation and moderate deactivation, as well as, two fast repolarizing and slow shunting potassium channels. The model quantitatively reproduces characteristics of steady-state activity that are typical for a cortical pyramidal neuron, namely firing rate not exceeding 30 Hz; critical values of the stimulating current and conductance which induce the depolarization block not exceeding 80 mV and 3, respectively (both values are scaled by the resting input conductance); extremum of hyperpolarization close to the midpoint between spikes. The analysis of the model reveals that the spiking regime appears through a saddle-node-on-invariant-circle bifurcation, and the depolarization block is reached through a saddle-node bifurcation of cycles. The model can be used for realistic network simulations, and it can also be implemented within the so-called mean-field, refractory density framework.

Viruses - a major cause of amyloid deposition in the brain.

INTRODUCTION: Clinically, Alzheimer's disease (AD) is a syndrome with a spectrum of various cognitive disorders. There is a complete dissociation between the pathology and the clinical presentation. Therefore, we need a disruptive new approach to be able to prevent and treat AD. AREAS COVERED: In this review, the authors extensively discuss the evidence why the amyloid beta is not the pathological cause of AD which makes therefore the amyloid hypothesis not sustainable anymore. They review the experimental evidence underlying the role of microbes, especially that of viruses, as a trigger/cause for the production of amyloid beta leading to the establishment of a chronic neuroinflammation as the mediator manifesting decades later by AD as a clinical spectrum. In this context, the emergence and consequences of the infection/antimicrobial protection hypothesis are described. The epidemiological and clinical data supporting this hypothesis are also analyzed. EXPERT OPINION: For decades, we have known that viruses are involved in the pathogenesis of AD. This discovery was ignored and discarded for a long time. Now we should accept this fact, which is not a hypothesis anymore, and stimulate the research community to come up with new ideas, new treatments, and new concepts.

Idealized multiple-timescale model of cortical spreading depolarization initiation and pre-epileptic hyperexcitability caused by NaV1.1/SCN1A mutations.

NaV1.1 (SCN1A) is a voltage-gated sodium channel mainly expressed in GABAergic neurons. Loss of function mutations of NaV1.1 lead to epileptic disorders, while gain of function mutations cause a migraine in which cortical spreading depolarizations (CSDs) are involved. It is still debated how these opposite effects initiate two different manifestations of neuronal hyperactivity: epileptic seizures and CSD. To investigate this question, we previously built a conductance-based model of two neurons (GABAergic and pyramidal), with dynamic ion concentrations (Lemaire et al. in PLoS Comput Biol 17(7):e1009239, 2021. https://doi.org/10.1371/journal.pcbi.1009239 ). When implementing either NaV1.1 migraine or epileptogenic mutations, ion concentration modifications acted as slow processes driving the system to the corresponding pathological firing regime. However, the large dimensionality of the model complicated the exploitation of its implicit multi-timescale structure. Here, we substantially simplify our biophysical model to a minimal version more suitable for bifurcation analysis. The explicit timescale separation allows us to apply slow-fast theory, where slow variables are treated as parameters in the fast singular limit. In this setting, we reproduce both pathological transitions as dynamic bifurcations in the full system. In the epilepsy condition, we shift the spike-terminating bifurcation to lower inputs for the GABAergic neuron, to model an increased susceptibility to depolarization block. The resulting failure of synaptic inhibition triggers hyperactivity of the pyramidal neuron. In the migraine scenario, spiking-induced release of potassium leads to the abrupt increase of the extracellular potassium concentration. This causes a dynamic spike-terminating bifurcation of both neurons, which we interpret as CSD initiation.

Geometry of spiking patterns in early visual cortex: a topological data analytic approach.

In the brain, spiking patterns live in a high-dimensional space of neurons and time. Thus, determining the intrinsic structure of this space presents a theoretical and experimental challenge. To address this challenge, we introduce a new framework for applying topological data analysis (TDA) to spike train data and use it to determine the geometry of spiking patterns in the visual cortex. Key to our approach is a parametrized family of distances based on the timing of spikes that quantifies the dissimilarity between neuronal responses. We applied TDA to visually driven single-unit and multiple single-unit spiking activity in macaque V1 and V2. TDA across timescales reveals a common geometry for spiking patterns in V1 and V2 which, among simple models, is most similar to that of a low-dimensional space endowed with Euclidean or hyperbolic geometry with modest curvature. Remarkably, the inferred geometry depends on timescale and is clearest for the timescales that are important for encoding contrast, orientation and spatial correlations.

Interaction Mechanism Between the HSV-1 Glycoprotein B and the Antimicrobial Peptide Amyloid-ß.

Background: Unravelling the mystery of Alzheimer's disease (AD) requires urgent resolution given the worldwide increase of the aging population. There is a growing concern that the current leading AD hypothesis, the amyloid cascade hypothesis, does not stand up to validation with respect to emerging new data. Indeed, several paradoxes are being discussed in the literature, for instance, both the deposition of the amyloid-ß peptide (Aß) and the intracellular neurofibrillary tangles could occur within the brain without any cognitive pathology. Thus, these paradoxes suggest that something more fundamental is at play in the onset of the disease and other key and related pathomechanisms must be investigated. Objective: The present study follows our previous investigations on the infectious hypothesis, which posits that some pathogens are linked to late onset AD. Our studies also build upon the finding that Aß is a powerful antimicrobial agent, produced by neurons in response to viral infection, capable of inhibiting pathogens as observed in in vitro experiments. Herein, we ask what are the molecular mechanisms in play when Aß neutralizes infectious pathogens? Methods: To answer this question, we probed at nanoscale lengths with FRET (Förster Resonance Energy Transfer), the interaction between Aß peptides and glycoprotein B (responsible of virus-cell binding) within the HSV-1 virion. Results: The experiments show an energy transfer between Aß peptides and glycoprotein B when membrane is intact. No energy transfer occurs after membrane disruption or treatment with blocking antibody. Conclusion: We concluded that Aß insert into viral membrane, close to glycoprotein B, and participate in virus neutralization.

Cross-scale excitability in networks of quadratic integrate-and-fire neurons.

From the action potentials of neurons and cardiac cells to the amplification of calcium signals in oocytes, excitability is a hallmark of many biological signalling processes. In recent years, excitability in single cells has been related to multiple-timescale dynamics through canards, special solutions which determine the effective thresholds of the all-or-none responses. However, the emergence of excitability in large populations remains an open problem. Here, we show that the mechanism of excitability in large networks and mean-field descriptions of coupled quadratic integrate-and-fire (QIF) cells mirrors that of the individual components. We initially exploit the Ott-Antonsen ansatz to derive low-dimensional dynamics for the coupled network and use it to describe the structure of canards via slow periodic forcing. We demonstrate that the thresholds for onset and offset of population firing can be found in the same way as those of the single cell. We combine theoretical analysis and numerical computations to develop a novel and comprehensive framework for excitability in large populations, applicable not only to models amenable to Ott-Antonsen reduction, but also to networks without a closed-form mean-field limit, in particular sparse networks.

Immunosenescence and Altered Vaccine Efficiency in Older Subjects: A Myth Difficult to Change.

Organismal ageing is associated with many physiological changes, including differences in the immune system of most animals. These differences are often considered to be a key cause of age-associated diseases as well as decreased vaccine responses in humans. The most often cited vaccine failure is seasonal influenza, but, while it is usually the case that the efficiency of this vaccine is lower in older than younger adults, this is not always true, and the reasons for the differential responses are manifold. Undoubtedly, changes in the innate and adaptive immune response with ageing are associated with failure to respond to the influenza vaccine, but the cause is unclear. Moreover, recent advances in vaccine formulations and adjuvants, as well as in our understanding of immune changes with ageing, have contributed to the development of vaccines, such as those against herpes zoster and SARS-CoV-2, that can protect against serious disease in older adults just as well as in younger people. In the present article, we discuss the reasons why it is a myth that vaccines inevitably protect less well in older individuals, and that vaccines represent one of the most powerful means to protect the health and ensure the quality of life of older adults.

Classification of bursting patterns: A tale of two ducks.

Bursting is one of the fundamental rhythms that excitable cells can generate either in response to incoming stimuli or intrinsically. It has been a topic of intense research in computational biology for several decades. The classification of bursting oscillations in excitable systems has been the subject of active research since the early 1980s and is still ongoing. As a by-product, it establishes analytical and numerical foundations for studying complex temporal behaviors in multiple timescale models of cellular activity. In this review, we first present the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable systems. We recall a complementary mathematical classification approach by Bertram and colleagues, and then by Golubitsky and colleagues, which, together with the Rinzel-Izhikevich proposals, provide the state-of-the-art foundations to these classifications. Beyond classical approaches, we review a recent bursting example that falls outside the previous classification systems. Generalizing this example leads us to propose an extended classification, which requires the analysis of both fast and slow subsystems of an underlying slow-fast model and allows the dissection of a larger class of bursters. Namely, we provide a general framework for bursting systems with both subthreshold and superthreshold oscillations. A new class of bursters with at least 2 slow variables is then added, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. Key to this mechanism are so-called canard or duck orbits, organizing the underpinning excitability structure. We describe the 2 main families of folded-node bursters, depending upon the phase (active/spiking or silent/nonspiking) of the bursting cycle during which folded-node dynamics occurs. We classify both families and give examples of minimal systems displaying these novel bursting patterns. Finally, we provide a biophysical example by reinterpreting a generic conductance-based episodic burster as a folded-node burster, showing that the associated framework can explain its subthreshold oscillations over a larger parameter region than the fast subsystem approach.

Initiation of migraine-related cortical spreading depolarization by hyperactivity of GABAergic neurons and NaV1.1 channels.

Spreading depolarizations (SDs) are involved in migraine, epilepsy, stroke, traumatic brain injury, and subarachnoid hemorrhage. However, the cellular origin and specific differential mechanisms are not clear. Increased glutamatergic activity is thought to be the key factor for generating cortical spreading depression (CSD), a pathological mechanism of migraine. Here, we show that acute pharmacological activation of NaV1.1 (the main Na+ channel of interneurons) or optogenetic-induced hyperactivity of GABAergic interneurons is sufficient to ignite CSD in the neocortex by spiking-generated extracellular K+ build-up. Neither GABAergic nor glutamatergic synaptic transmission were required for CSD initiation. CSD was not generated in other brain areas, suggesting that this is a neocortex-specific mechanism of CSD initiation. Gain-of-function mutations of NaV1.1 (SCN1A) cause familial hemiplegic migraine type-3 (FHM3), a subtype of migraine with aura, of which CSD is the neurophysiological correlate. Our results provide the mechanism linking NaV1.1 gain of function to CSD generation in FHM3. Thus, we reveal the key role of hyperactivity of GABAergic interneurons in a mechanism of CSD initiation, which is relevant as a pathological mechanism of Nav1.1 FHM3 mutations, and possibly also for other types of migraine and diseases in which SDs are involved.

Modeling NaV1.1/SCN1A sodium channel mutations in a microcircuit with realistic ion concentration dynamics suggests differential GABAergic mechanisms leading to hyperexcitability in epilepsy and hemiplegic migraine.

Loss of function mutations of SCN1A, the gene coding for the voltage-gated sodium channel NaV1.1, cause different types of epilepsy, whereas gain of function mutations cause sporadic and familial hemiplegic migraine type 3 (FHM-3). However, it is not clear yet how these opposite effects can induce paroxysmal pathological activities involving neuronal networks' hyperexcitability that are specific of epilepsy (seizures) or migraine (cortical spreading depolarization, CSD). To better understand differential mechanisms leading to the initiation of these pathological activities, we used a two-neuron conductance-based model of interconnected GABAergic and pyramidal glutamatergic neurons, in which we incorporated ionic concentration dynamics in both neurons. We modeled FHM-3 mutations by increasing the persistent sodium current in the interneuron and epileptogenic mutations by decreasing the sodium conductance in the interneuron. Therefore, we studied both FHM-3 and epileptogenic mutations within the same framework, modifying only two parameters. In our model, the key effect of gain of function FHM-3 mutations is ion fluxes modification at each action potential (in particular the larger activation of voltage-gated potassium channels induced by the NaV1.1 gain of function), and the resulting CSD-triggering extracellular potassium accumulation, which is not caused only by modifications of firing frequency. Loss of function epileptogenic mutations, on the other hand, increase GABAergic neurons' susceptibility to depolarization block, without major modifications of firing frequency before it. Our modeling results connect qualitatively to experimental data: potassium accumulation in the case of FHM-3 mutations and facilitated depolarization block of the GABAergic neuron in the case of epileptogenic mutations. Both these effects can lead to pyramidal neuron hyperexcitability, inducing in the migraine condition depolarization block of both the GABAergic and the pyramidal neuron. Overall, our findings suggest different mechanisms of network hyperexcitability for migraine and epileptogenic NaV1.1 mutations, implying that the modifications of firing frequency may not be the only relevant pathological mechanism.

Hyperactivation of monocytes and macrophages in MCI patients contributes to the progression of Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) is the most common neurodegenerative disease ultimately manifesting as clinical dementia. Despite considerable effort and ample experimental data, the role of neuroinflammation related to systemic inflammation is still unsettled. While the implication of microglia is well recognized, the exact contribution of peripheral monocytes/macrophages is still largely unknown, especially concerning their role in the various stages of AD. OBJECTIVES: AD develops over decades and its clinical manifestation is preceded by subjective memory complaints (SMC) and mild cognitive impairment (MCI); thus, the question arises how the peripheral innate immune response changes with the progression of the disease. Therefore, to further investigate the roles of monocytes/macrophages in the progression of AD we assessed their phenotypes and functions in patients at SMC, MCI and AD stages and compared them with cognitively healthy controls. We also conceptualised an idealised mathematical model to explain the functionality of monocytes/macrophages along the progression of the disease. RESULTS: We show that there are distinct phenotypic and functional changes in monocyte and macrophage populations as the disease progresses. Higher free radical production upon stimulation could already be observed for the monocytes of SMC patients. The most striking results show that activation of peripheral monocytes (hyperactivation) is the strongest in the MCI group, at the prodromal stage of the disease. Monocytes exhibit significantly increased chemotaxis, free radical production, and cytokine production in response to TLR2 and TLR4 stimulation. CONCLUSION: Our data suggest that the peripheral innate immune system is activated during the progression from SMC through MCI to AD, with the highest levels of activation being in MCI subjects and the lowest in AD patients. Some of these parameters may be used as biomarkers, but more holistic immune studies are needed to find the best period of the disease for clinical intervention.

Targeting Impaired Antimicrobial Immunity in the Brain for the Treatment of Alzheimer's Disease.

Alzheimer's disease (AD) is the most common form of dementia and aging is the most common risk factor for developing the disease. The etiology of AD is not known but AD may be considered as a clinical syndrome with multiple causal pathways contributing to it. The amyloid cascade hypothesis, claiming that excess production or reduced clearance of amyloid-beta (Aß) and its aggregation into amyloid plaques, was accepted for a long time as the main cause of AD. However, many studies showed that Aß is a frequent consequence of many challenges/pathologic processes occurring in the brain for decades. A key factor, sustained by experimental data, is that low-grade infection leading to production and deposition of Aß, which has antimicrobial activity, precedes the development of clinically apparent AD. This infection is chronic, low grade, largely clinically silent for decades because of a nearly efficient antimicrobial immune response in the brain. A chronic inflammatory state is induced that results in neurodegeneration. Interventions that appear to prevent, retard or mitigate the development of AD also appear to modify the disease. In this review, we conceptualize further that the changes in the brain antimicrobial immune response during aging and especially in AD sufferers serve as a foundation that could lead to improved treatment strategies for preventing or decreasing the progression of AD in a disease-modifying treatment.

Why we should use topological data analysis in ageing: Towards defining the "topological shape of ageing".

Living systems are subject to the arrow of time; from birth, they undergo complex transformations (self-organization) in a constant battle for survival, but inevitably ageing and disease trap them to death. Can ageing be understood and eventually reversed? What tools can be employed to further our understanding of ageing? The present article is an invitation for biologists and clinicians to consider key conceptual ideas and computational tools (known to mathematicians and physicists), which potentially may help dissect some of the underlying processes of ageing and disease. Specifically, we first discuss how to classify and analyse complex systems, as well as highlight critical theoretical difficulties that make complex systems hard to study. Subsequently, we introduce Topological Data Analysis - a novel Big Data tool - which may help in the study of complex systems since it extracts knowledge from data in a holistic approach via topological considerations. These conceptual ideas and tools are discussed in a relatively informal way to pave future discussions and collaborations between mathematicians and biologists studying ageing.

Targeting Infectious Agents as a Therapeutic Strategy in Alzheimer's Disease.

Alzheimer's disease (AD) is the most prevalent dementia in the world. Its cause(s) are presently largely unknown. The most common explanation for AD, now, is the amyloid cascade hypothesis, which states that the cause of AD is senile plaque formation by the amyloid ß peptide, and the formation of neurofibrillary tangles by hyperphosphorylated tau. A second, burgeoning theory by which to explain AD is based on the infection hypothesis. Much experimental and epidemiological data support the involvement of infections in the development of dementia. According to this mechanism, the infection either directly or via microbial virulence factors precedes the formation of amyloid ß plaques. The amyloid ß peptide, possessing antimicrobial properties, may be beneficial at an early stage of AD, but becomes detrimental with the progression of the disease, concomitantly with alterations to the innate immune system at both the peripheral and central levels. Infection results in neuroinflammation, leading to, and sustained by, systemic inflammation, causing eventual neurodegeneration, and the senescence of the immune cells. The sources of AD-involved microbes are various body microbiome communities from the gut, mouth, nose, and skin. The infection hypothesis of AD opens a vista to new therapeutic approaches, either by treating the infection itself or modulating the immune system, its senescence, or the body's metabolism, either separately, in parallel, or in a multi-step way.

Canard-induced complex oscillations in an excitatory network.

In Neuroscience, mathematical modelling involving multiple spatial and temporal scales can unveil complex oscillatory activity such as excitable responses to an input current, subthreshold oscillations, spiking or bursting. While the number of slow and fast variables and the geometry of the system determine the type of the complex oscillations, canard structures define boundaries between them. In this study, we use geometric singular perturbation theory to identify and characterise boundaries between different dynamical regimes in multiple-timescale firing rate models of the developing spinal cord. These rate models are either three or four dimensional with state variables chosen within an overall group of two slow and two fast variables. The fast subsystem corresponds to a recurrent excitatory network with fast activity-dependent synaptic depression, and the slow variables represent the cell firing threshold and slow activity-dependent synaptic depression, respectively. We start by demonstrating canard-induced bursting and mixed-mode oscillations in two different three-dimensional rate models. Then, in the full four-dimensional model we show that a canard-mediated slow passage creates dynamics that combine these complex oscillations and give rise to mixed-mode bursting oscillations (MMBOs). We unveil complicated isolas along which MMBOs exist in parameter space. The profile of solutions along each isola undergoes canard-mediated transitions between the sub-threshold regime and the bursting regime; these explosive transitions change the number of oscillations in each regime. Finally, we relate the MMBO dynamics to experimental recordings and discuss their effects on the silent phases of bursting patterns as well as their potential role in creating subthreshold fluctuations that are often interpreted as noise. The mathematical framework used in this paper is relevant for modelling multiple timescale dynamics in excitable systems.

Excitable dynamics in neural and cardiac systems.

The application of mathematics, physics and engineering to medical research is continuously growing; interactions among these disciplines have become increasingly important and have contributed to an improved understanding of clinical and biological phenomena, with implications for disease prevention, diagnosis and treatment. This special issue presents examples of this synergy, with a particular focus on the investigation of cardiac and neural excitability. This issue includes 24 original research papers and covers a broad range of topics related to the physiological and pathophysiological function of the brain and the heart. Studies span scales from isolated neurons and small networks of neurons to whole-organ dynamics for the brain and from cardiac subcellular domains and cardiomyocytes to one-dimensional tissues for the heart. This preface is part of the Special Issue on "Excitable Dynamics in Neural and Cardiac Systems".

Modeling cortical spreading depression induced by the hyperactivity of interneurons.

Cortical spreading depression (CSD) is a wave of transient intense neuronal firing leading to a long lasting depolarizing block of neuronal activity. It is a proposed pathological mechanism of migraine with aura. Some forms of migraine are associated with a genetic mutation of the Nav1.1 channel, resulting in its gain of function and implying hyperexcitability of interneurons. This leads to the counterintuitive hypothesis that intense firing of interneurons can cause CSD ignition. To test this hypothesis in silico, we developed a computational model of an E-I pair (a pyramidal cell and an interneuron), in which the coupling between the cells in not just synaptic, but takes into account also the effects of the accumulation of extracellular potassium caused by the activity of the neurons and of the synapses. In the context of this model, we show that the intense firing of the interneuron can lead to CSD. We have investigated the effect of various biophysical parameters on the transition to CSD, including the levels of glutamate or GABA, frequency of the interneuron firing and the efficacy of the KCC2 co-transporter. The key element for CSD ignition in our model was the frequency of interneuron firing and the related accumulation of extracellular potassium, which induced a depolarizing block of the pyramidal cell. This constitutes a new mechanism of CSD ignition.

Metastable Resting State Brain Dynamics.

Metastability refers to the fact that the state of a dynamical system spends a large amount of time in a restricted region of its available phase space before a transition takes place, bringing the system into another state from where it might recur into the previous one. beim Graben and Hutt (2013) suggested to use the recurrence plot (RP) technique introduced by Eckmann et al. (1987) for the segmentation of system's trajectories into metastable states using recurrence grammars. Here, we apply this recurrence structure analysis (RSA) for the first time to resting-state brain dynamics obtained from functional magnetic resonance imaging (fMRI). Brain regions are defined according to the brain hierarchical atlas (BHA) developed by Diez et al. (2015), and as a consequence, regions present high-connectivity in both structure (obtained from diffusion tensor imaging) and function (from the blood-level dependent-oxygenation-BOLD-signal). Remarkably, regions observed by Diez et al. were completely time-invariant. Here, in order to compare this static picture with the metastable systems dynamics obtained from the RSA segmentation, we determine the number of metastable states as a measure of complexity for all subjects and for region numbers varying from 3 to 100. We find RSA convergence toward an optimal segmentation of 40 metastable states for normalized BOLD signals, averaged over BHA modules. Next, we build a bistable dynamics at population level by pooling 30 subjects after Hausdorff clustering. In link with this finding, we reflect on the different modeling frameworks that can allow for such scenarios: heteroclinic dynamics, dynamics with riddled basins of attraction, multiple-timescale dynamics. Finally, we characterize the metastable states both functionally and structurally, using templates for resting state networks (RSNs) and the automated anatomical labeling (AAL) atlas, respectively.

Conductance-Based Refractory Density Approach for a Population of Bursting Neurons.

The conductance-based refractory density (CBRD) approach is a parsimonious mathematical-computational framework for modelling interacting populations of regular spiking neurons, which, however, has not been yet extended for a population of bursting neurons. The canonical CBRD method allows to describe the firing activity of a statistical ensemble of uncoupled Hodgkin-Huxley-like neurons (differentiated by noise) and has demonstrated its validity against experimental data. The present manuscript generalises the CBRD for a population of bursting neurons; however, in this pilot computational study, we consider the simplest setting in which each individual neuron is governed by a piecewise linear bursting dynamics. The resulting population model makes use of slow-fast analysis, which leads to a novel methodology that combines CBRD with the theory of multiple timescale dynamics. The main prospect is that it opens novel avenues for mathematical explorations, as well as, the derivation of more sophisticated population activity from Hodgkin-Huxley-like bursting neurons, which will allow to capture the activity of synchronised bursting activity in hyper-excitable brain states (e.g. onset of epilepsy).