NeuroCa: integrated framework for systematic analysis of spatiotemporal neuronal activity patterns from large-scale optical recording data

Harnessing the potential of human induced pluripotent stem cells, functional assays and machine learning for neurodevelopmental disorders

Neurodevelopmental disorders (NDDs) affect 4.7% of the global population and are associated with delays in brain development and a spectrum of impairments that can lead to lifelong disability and even mortality. Identification of biomarkers for accurate diagnosis and medications for effective treatment are lacking, in part due to the historical use of preclinical model systems that do not translate well to the clinic for neurological disorders, such as rodents and heterologous cell lines. Human-induced pluripotent stem cells (hiPSCs) are a promising in vitro system for modeling NDDs, providing opportunities to understand mechanisms driving NDDs in human neurons. Functional assays, including patch clamping, multielectrode array, and imaging-based assays, are popular tools employed with hiPSC disease models for disease investigation. Recent progress in machine learning (ML) algorithms also presents unprecedented opportunities to advance the NDD research process. In this review, we compare two-dimensional and three-dimensional hiPSC formats for disease modeling, discuss the applications of functional assays, and offer insights on incorporating ML into hiPSC-based NDD research and drug screening.

A 3D neuronal network read-out interface with high recording performance using a neuronal cluster patterning on a microelectrode array

In recent years, in vitro three-dimensional (3D) neuronal network models utilizing extracellular matrices have been advancing. To understand the network activity from these models, attempts have been made to measure activity in multiple regions simultaneously using a microelectrode array (MEA). Although there hve been many attempts to measure the activity of 3D networks using 2-dimensional (2D) MEAs, the physical coupling between the 3D network and the microelectrodes was not stable and needed to be improved. In this study, we proposed a neuronal cluster interface that improves the active channel ratio of commercial 2D MEAs, enabling reliable measurement of 3D network activity. To achieve this, neuronal clusters, which consist of a small number of neurons, were patterned on microelectrodes and used as mediators to transmit the signal between the 3D network and the microelectrodes. We confirmed that the patterned neuronal clusters enhanced the active channel ratio and SNR(signal-to-noise-ratio) about 3D network recording and stimulation for a month. Our interface was able to functionally connect with 3D networks and measure the 3D network activity without significant alternation of activity characteristics. Finally, we demonstrated that our interface can be used to analyze the differences in the dynamics of 3D and 2D networks and to construct the 3D clustered network. This method is expected to be useful for studying the functional activity of various 3D neuronal network models, offering broad applications for the use of these models.

ALDH5A1-deficient iPSC-derived excitatory and inhibitory neurons display cell type specific alterations

Succinic semialdehyde dehydrogenase deficiency (SSADHD) is a neurometabolic disorder caused by ALDH5A1 mutations presenting with autism and epilepsy. SSADHD leads to impaired GABA metabolism and results in accumulation of GABA and γ-hydroxybutyrate (GHB), which alter neurotransmission and are thought to lead to neurobehavioral symptoms. However, why increased inhibitory neurotransmitters lead to seizures remains unclear. We used induced pluripotent stem cells from SSADHD patients (one female and two male) and differentiated them into GABAergic and glutamatergic neurons. SSADHD iGABA neurons show altered GABA metabolism and concomitant changes in expression of genes associated with inhibitory neurotransmission. In contrast, glutamatergic neurons display increased spontaneous activity and upregulation of mitochondrial genes. CRISPR correction of the pathogenic variants or SSADHD mRNA expression rescue various metabolic and functional abnormalities in human neurons. Our findings uncover a previously unknown role for SSADHD in excitatory human neurons and provide unique insights into the cellular and molecular basis of SSADHD and potential therapeutic interventions.

High temporal resolution transparent thermoelectric temperature sensors for photothermal effect sensing

We propose inkjet-printed high-speed and transparent temperature sensors based on the thermoelectric effect for direct monitoring of the photothermal effect. They consist of highly transparent organic thermoelectric materials that allow excellent biocompatibility and sub-ms temporal resolution, simultaneously. Our transparent thermoelectric temperature sensors can be used to advance various photothermal biomedical applications.

NRSF/REST lies at the intersection between epigenetic regulation, miRNA-mediated gene control and neurodevelopmental pathways associated with Intellectual disability (ID) and Schizophrenia

Genetic evidence indicates disrupted epigenetic regulation as a major risk factor for psychiatric disorders, but the molecular mechanisms that drive this association remain to be determined. EHMT1 is an epigenetic repressor that is causal for Kleefstra Syndrome (KS), a genetic disorder linked with neurodevelopmental disorders and associated with schizophrenia. Here, we show that reduced EHMT1 activity decreases NRSF/REST protein leading to abnormal neuronal gene expression and progression of neurodevelopment in human iPSC. We further show that EHMT1 regulates NRSF/REST indirectly via repression of miRNA and leads to aberrant neuronal gene regulation and neurodevelopment timing. Expression of a NRSF/REST mRNA that lacks the miRNA-binding sites restores neuronal gene regulation to EHMT1 deficient cells. Significantly, the EHMT1-regulated miRNA gene set not only controls NRSF/REST but is enriched for association for Intellectual Disability (ID) and schizophrenia. This reveals a broad molecular interaction between H3K9 demethylation, NSRF/REST regulation and risk for ID and Schizophrenia.

Cortical neuronal hyperexcitability and synaptic changes in SGCE mutation-positive myoclonus dystonia

Myoclonus Dystonia is a childhood-onset hyperkinetic movement disorder with a combined motor and psychiatric phenotype. It represents one of the few autosomal dominant inherited dystonic disorders and is caused by mutations in the ε-sarcoglycan (SGCE) gene. Work to date suggests that dystonia is caused by disruption of neuronal networks, principally basal ganglia-cerebello-thalamo-cortical circuits. Investigation of cortical involvement has primarily focused on disruption to interneuron inhibitory activity, rather than the excitatory activity of cortical pyramidal neurons. Here, we have sought to examine excitatory cortical glutamatergic activity using two approaches; the CRISPR/Cas9 editing of a human embryonic cell line, generating an SGCE compound heterozygous mutation, and three patient-derived induced pluripotent stem cell lines (iPSC) each gene edited to generate matched wild-type SGCE control lines. Differentiation towards a cortical neuronal phenotype demonstrated no significant differences in neither early- (PAX6, FOXG1) nor late-stage (CTIP2, TBR1) neurodevelopmental markers. However, functional characterisation using Ca2+ imaging and MEA approaches identified an increase in network activity, while single-cell patch clamp studies found a greater propensity towards action potential generation with larger amplitudes and shorter half-widths associated with SGCE-mutations. Bulk-RNA-seq analysis identified gene ontological enrichment for neuron projection development, synaptic signalling, and synaptic transmission. Examination of dendritic morphology found SGCE-mutations to be associated with a significantly higher number of branches and longer branch lengths, together with longer ion-channel dense axon initial segments, particularly towards the latter stages of differentiation (D80 and D100). Gene expression and protein quantification of key synaptic proteins (synaptophysin, synapsin and PSD95), AMPA and NMDA receptor subunits found no significant differences between the SGCE-mutation and matched wild-type lines. By contrast, significant changes to synaptic adhesion molecule expression were identified, namely higher pre-synaptic neurexin-1 and lower post-synaptic neuroligin-4 levels in the SGCE mutation carrying lines. Our study demonstrates an increased intrinsic excitability of cortical glutamatergic neuronal cells in the context of SGCE mutations, coupled with a more complex neurite morphology and disruption to synaptic adhesion molecules. These changes potentially represent key components to the development of the hyperkinetic clinical phenotype observed in Myoclonus Dystonia, as well a central feature to the wider spectrum of dystonic disorders, potentially providing targets for future therapeutic development.

Analysis of co-isogenic prion protein deficient mice reveals behavioral deficits, learning impairment, and enhanced hippocampal excitability

Background

Cellular prion protein (PrP^C) is a cell surface GPI-anchored protein, usually known for its role in the pathogenesis of human and animal prionopathies. However, increasing knowledge about the participation of PrP^C in prion pathogenesis contrasts with puzzling data regarding its natural physiological role. PrP^C is expressed in a number of tissues, including at high levels in the nervous system, especially in neurons and glial cells, and while previous studies have established a neuroprotective role, conflicting evidence for a synaptic function has revealed both reduced and enhanced long-term potentiation, and variable observations on memory, learning, and behavior. Such evidence has been confounded by the absence of an appropriate knock-out mouse model to dissect the biological relevance of PrP^C, with some functions recently shown to be misattributed to PrP^C due to the presence of genetic artifacts in mouse models. Here we elucidate the role of PrP^C in the hippocampal circuitry and its related functions, such as learning and memory, using a recently available strictly co-isogenic Prnp^0/0 mouse model (Prnp^ZH3/ZH3).

Results

We performed behavioral and operant conditioning tests to evaluate memory and learning capabilities, with results showing decreased motility, impaired operant conditioning learning, and anxiety-related behavior in Prnp^ZH3/ZH3 animals. We also carried in vivo electrophysiological recordings on CA3-CA1 synapses in living behaving mice and monitored spontaneous neuronal firing and network formation in primary neuronal cultures of Prnp^ZH3/ZH3 vs wildtype mice. PrP^C absence enhanced susceptibility to high-intensity stimulations and kainate-induced seizures. However, long-term potentiation (LTP) was not enhanced in the Prnp^ZH3/ZH3 hippocampus. In addition, we observed a delay in neuronal maturation and network formation in Prnp^ZH3/ZH3 cultures.

Conclusion

Our results demonstrate that PrP^C promotes neuronal network formation and connectivity. PrP^C mediates synaptic function and protects the synapse from excitotoxic insults. Its deletion may underlie an epileptogenic-susceptible brain that fails to perform highly cognitive-demanding tasks such as associative learning and anxiety-like behaviors.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01203-0.

Acoustically triggered mechanotherapy using genetically encoded gas vesicles

Recent advances in molecular engineering and synthetic biology provide biomolecular and cell-based therapies with a high degree of molecular specificity, but limited spatiotemporal control. Here we show that biomolecules and cells can be engineered to deliver potent mechanical effects at specific locations inside the body through ultrasound-induced inertial cavitation. This capability is enabled by gas vesicles, a unique class of genetically encodable air-filled protein nanostructures. We show that low-frequency ultrasound can convert these biomolecules into micrometre-scale cavitating bubbles, unleashing strong local mechanical effects. This enables engineered gas vesicles to serve as remotely actuated cell-killing and tissue-disrupting agents, and allows genetically engineered cells to lyse, release molecular payloads and produce local mechanical damage on command. We demonstrate the capabilities of biomolecular inertial cavitation in vitro, in cellulo and in vivo, including in a mouse model of tumour-homing probiotic therapy.

Development of the micro-patterned 3D neuronal-hydrogel model using soft-lithography for study a 3D neural network on a microelectrode array. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:1234-1237. [PMID: 34891510 DOI: 10.1109/embc46164.2021.9629822]

In vitro patterned neuronal models have been studied as one of the strategies to investigate the relationship between structural connectivity and functional activity of neural network. Despite the importance of three-dimensional (3D) cell models, most of these studies have been performed on two-dimensional models. In this study, we present a technique to construct the micro-pattern to 3D neuronal-hydrogel model using a micromolding in capillaries (MIMIC) technique on microelectrode array (MEA). Our technique was suitable to prevent the deformation of micro-patterned collagen model against the neuronal contracted tension during the network formation. The relationship between the growth directions of glial cells and micro-pattern direction was investigated. Lastly, we confirmed that our 3D model had synchronized activity among neurons in 3D. This model is expected to be used as a tool to study the relationship between structural connectivity and functional activity in the 3D environment.

CalTrack: High-Throughput Automated Calcium Transient Analysis in Cardiomyocytes

[Figure: see text].

Adult mouse dorsal root ganglia neurons form aberrant glutamatergic connections in dissociated cultures

Cultured sensory neurons can exhibit complex activity patterns following stimulation in terms of increased excitability and interconnected responses of multiple neurons. Although these complex activity patterns suggest a network-like configuration, research so far had little interest in synaptic network formation ability of the sensory neurons. To identify interaction profiles of Dorsal Root Ganglia (DRG) neurons and explore their putative connectivity, we developed an in vitro experimental approach. A double transgenic mouse model, expressing genetically encoded calcium indicator (GECI) in their glutamatergic neurons, was produced. Dissociated DRG cultures from adult mice were prepared with a serum-free protocol and no additional growth factors or cytokines were utilized for neuronal sensitization. DRG neurons were grown on microelectrode arrays (MEA) to induce stimulus-evoked activity with a modality-free stimulation strategy. With an almost single-cell level electrical stimulation, spontaneous and evoked activity of GCaMP6s expressing neurons were detected under confocal microscope. Typical responses were analyzed, and correlated calcium events were detected across individual DRG neurons. Next, correlated responses were successfully blocked by glutamatergic receptor antagonists, which indicated functional synaptic coupling. Immunostaining confirmed the presence of synapses mainly in the axonal terminals, axon-soma junctions and axon-axon intersection sites. Concisely, the results presented here illustrate a new type of neuron-to-neuron interaction in cultured DRG neurons conducted through synapses. The developed assay can be a valuable tool to analyze individual and collective responses of the cultured sensory neurons.

A simple Ca2+-imaging approach to neural network analyses in cultured neurons

BACKGROUND

Ca²⁺-imaging is a powerful tool to measure neuronal dynamics and network activity. To monitor network-level changes in cultured neurons, neuronal activity is often evoked by electrical or optogenetic stimulation and assessed using multi-electrode arrays or sophisticated imaging. Although such approaches allow detailed network analyses, multi-electrode arrays lack single-cell precision, whereas optical physiology generally requires advanced instrumentation that may not be universally available.

NEW METHOD

Here we developed a simple, stimulation-free protocol with associated Matlab algorithms that enables scalable analyses of spontaneous network activity in cultured human and mouse neurons. The approach allows analysis of the overall network activity and of single-neuron dynamics, and is amenable to screening purposes.

RESULTS

We validated the new protocol by assessing human neurons with a heterozygous conditional deletion of Munc18-1, and mouse neurons with a homozygous conditional deletion of neurexins. The approach described enabled identification of differential changes in these mutant neurons, allowing quantifications of the synchronous firing rate at the network level and of the amplitude and frequency of Ca²⁺-spikes at the single-neuron level. These results demonstrate the utility of the approach.

COMPARISION WITH EXISTING METHODS

Compared with current imaging platforms, our method is simple, scalable, accessible, and easy to implement. It enables quantification of more detailed parameters than multi-electrode arrays, but does not have the resolution and depth of more sophisticated yet labour-intensive methods, such as patch-clamp electrophysiology.

CONCLUSION

The method reported here is scalable for a rapid direct assessment of neuronal function in culture, and can be applied to both human and mouse neurons. Thus, the method can serve as a basis for phenotypical analysis of mutations and for drug discovery efforts.

Microglial Gi-dependent dynamics regulate brain network hyperexcitability

Microglial surveillance is a key feature of brain physiology and disease. We found that Gi-dependent microglial dynamics prevent neuronal network hyperexcitability. By generating Mg^PTX mice to genetically inhibit Gi in microglia, we showed that sustained reduction of microglia brain surveillance and directed process motility induced spontaneous seizures and increased hypersynchrony upon physiologically evoked neuronal activity in awake adult mice. Thus, Gi-dependent microglia dynamics may prevent hyperexcitability in neurological diseases.

NeuroCore formation during differentiation of neurospheres of mouse embryonic neural stem cells

Neural stem cells (NSCs) in the embryonic neocortex have the potential to generate a well-organized laminar architecture of the cerebral cortex through precise regulation of the proliferation, differentiation, and migration of neural cells. NSCs can be isolated in vitro and expanded as cell clusters, called neurospheres, which are primarily related to the proliferation ability of NSCs. Conversely, the tissue-organizing properties of NSCs via regulated differentiation and migration of the cells are not well understood. In this study, we established a three-dimensional (3D) differentiation model of neurospheres, which produce unique neuronal clusters, termed NeuroCore (NC). NC formation was initiated by the aggregation of young neurons. Upon maturation of the neurons and the establishment of radial glia-like structures, the initial organization of the NCs transformed into a glomeruli-like arrangement of cortical neurons. These neurons expressed multiple markers of upper and deep cortical neurons. Taken together, we propose that NSCs in vitro maintain some aspects of their original in vivo tissue-organizing properties, providing an alternative opportunity to explore the fundamental components of brain histogenesis in vitro.

Response Latency Tuning by Retinal Circuits Modulates Signal Efficiency

In the visual system, retinal ganglion cells (RGCs) of various subtypes encode preprocessed photoreceptor signals into a spike output which is then transmitted towards the brain through parallel feature pathways. Spike timing determines how each feature signal contributes to the output of downstream neurons in visual brain centers, thereby influencing efficiency in visual perception. In this study, we demonstrate a marked population-wide variability in RGC response latency that is independent of trial-to-trial variability and recording approach. RGC response latencies to simple visual stimuli vary considerably in a heterogenous cell population but remain reliable when RGCs of a single subtype are compared. This subtype specificity, however, vanishes when the retinal circuitry is bypassed via direct RGC electrical stimulation. This suggests that latency is primarily determined by the signaling speed through retinal pathways that provide subtype specific inputs to RGCs. In addition, response latency is significantly altered when GABA inhibition or gap junction signaling is disturbed, which further supports the key role of retinal microcircuits in latency tuning. Finally, modulation of stimulus parameters affects individual RGC response delays considerably. Based on these findings, we hypothesize that retinal microcircuits fine-tune RGC response latency, which in turn determines the context-dependent weighing of each signal and its contribution to visual perception.

CALIMA: The semi-automated open-source calcium imaging analyzer

BACKGROUND AND OBJECTIVE

Ever since its discovery, calcium imaging has proven its worth in discovering new insights into the mechanisms of cellular communication. Yet, the analysis of the data generated by calcium imaging experiments demands a large amount of time from researchers. Tools enabling automated and semi-automated analysis are available, but often they allow automating only a part of the data analysis process. Therefore, we developed CALIMA (https://aethelraed.nl/calima), a free and open-source standalone software tool that provides an opportunity to quickly detect cells, to obtain the calcium spikes, and to determine the underlying network structure of neuronal cell cultures.

METHODS

Owing to the difference of Gaussians algorithm applied for the cell detection, CALIMA is able to detect regions of interest (ROIs) quickly. The z-scoring algorithm provides a means to set the requirements for spike detection, and the neuronal connections can be reconstructed by analyzing the cross-correlation between the cellular activity. We evaluated CALIMA's reliability, speed, and functionality with a special focus on neuronal cell detection and network reconstruction. The evaluation was performed by using real-life data such as a known example dataset (cultured primary rat cortical neurons, University of Pennsylvania) and by analyzing video graphic footage of in vitro brain cell samples (SH-SY5Y neuroblastoma cultures, one sample with synchronous neuron firing). The obtained results were compared to the corresponding outcomes observed on same datasets for other similar software solutions. Moreover, we compared the results of segmentation and peak detection analysis, the ones obtained using CALIMA and those acquired manually.

RESULTS

CALIMA was able to detect the cells in the cultures within seconds. The average sensitivity was 82% across the datasets checked, comparing favorably with the alternative software solutions. Using the correct parameters, CALIMA's Ca-spikes detection sensitivity reached 96%. Lastly, neuronal networks were reconstructed by combining the data on the ROI's activity and the cell's positions, finding the most likely inter-cell connections.

CONCLUSIONS

We found that CALIMA proved to be a robust and fast tool to analyze the data of experiments for the digital reconstruction of the neuronal cellular network while being able to process the analysis steps with minimal user input required and in a time efficient manner.

SICT: automated detection and supervised inspection of fast Ca2+ transients

Recent advances in live Ca²⁺ imaging with increasing spatial and temporal resolution offer unprecedented opportunities, but also generate an unmet need for data processing. Here we developed SICT, a MATLAB program that automatically identifies rapid Ca²⁺ rises in time-lapse movies with low signal-to-noise ratios, using fluorescent indicators. A graphical user interface allows visual inspection of automatically detected events, reducing manual labour to less than 10% while maintaining quality control. The detection performance was tested using synthetic data with various signal-to-noise ratios. The event inspection phase was evaluated by four human observers. Reliability of the method was demonstrated in a direct comparison between manual and SICT-aided analysis. As a test case in cultured neurons, SICT detected an increase in frequency and duration of spontaneous Ca²⁺ transients in the presence of caffeine. This new method speeds up the analysis of elementary Ca²⁺ transients.

All-Optical Assay to Study Biological Neural Networks

We introduce a novel all-optical assay for functional studies of biological neural networks in vitro. We created a novel optogenetic construct named OptoCaMP which is a combination of a channelrhodopsin variant (CheRiff) and a red genetically encoded calcium indicator (GECI) (jRCaMP1b). It enables simultaneous optical stimulation and recording from large population of neurons with single-cell readout. Additionally, we have developed a spatio-temporal all-optical assay to simultaneously stimulate a sub-section of a neural network and record evoked calcium activity, in both stimulated and non-stimulated neurons, thus allowing the investigation of the spread of excitation through an interconnected network. Finally, we demonstrate the sensitivity of this assay to the change of neural network connectivity.

PeakCaller: an automated graphical interface for the quantification of intracellular calcium obtained by high-content screening

Background

Intracellular calcium is an important ion involved in the regulation and modulation of many neuronal functions. From regulating cell cycle and proliferation to initiating signaling cascades and regulating presynaptic neurotransmitter release, the concentration and timing of calcium activity governs the function and fate of neurons. Changes in calcium transients can be used in high-throughput screening applications as a basic measure of neuronal maturity, especially in developing or immature neuronal cultures derived from stem cells.

Results

Using human induced pluripotent stem cell derived neurons and dissociated mouse cortical neurons combined with the calcium indicator Fluo-4, we demonstrate that PeakCaller reduces type I and type II error in automated peak calling when compared to the oft-used PeakFinder algorithm under both basal and pharmacologically induced conditions.

Conclusion

Here we describe PeakCaller, a novel MATLAB script and graphical user interface for the quantification of intracellular calcium transients in neuronal cultures. PeakCaller allows the user to set peak parameters and smoothing algorithms to best fit their data set. This new analysis script will allow for automation of calcium measurements and is a powerful software tool for researchers interested in high-throughput measurements of intracellular calcium.

Influence of tetramethylenedisulfotetramine on synchronous calcium oscillations at distinct developmental stages of hippocampal neuronal cultures

The spatial and temporal patterns of spontaneous synchronous Ca²⁺ oscillations (SCOs) regulate physiological pathways that influence neuronal development, excitability, and health. Hippocampal neuronal cultures (HN) and neuron/glia co-cultures (HNG) produced from neonatal mice were loaded with Fluo-4/AM and SCOs recorded in real-time using a Fluorescence Imaging Plate Reader at different developmental stages in vitro. HNG showed an earlier onset of SCOs, with low amplitude and low frequency SCOs at 4days in vitro (DIV), whereas HN were quiescent at this point. SCO amplitude peaked at 9 DIV for both cultures. SCO network frequency peaked at 12 DIV in HN, whereas in HNG the frequency peaked at 6 DIV. SCO patterns were associated with the temporal development of neuronal networks and their ratio of glutamatergic to GABAergic markers of excitatory/inhibitory balance. HN and HNG exhibited differential responses to the convulsant tetramethylenedisulfotetramine (TETS) and were highly dependent on DIV. In HN, TETS triggered an acute rise of intracellular Ca²⁺ (Phase I response) only in 14 DIV and a sustained decrease of SCO frequency with increased amplitude (Phase II response) at all developmental stages. In HNG, TETS decreased the SCO frequency and increased the amplitude at 6 and 14 but not 9 DIV. There was no acute Ca²⁺ rise (Phase I response) in any age of HNG tested with TETS. These data demonstrated the importance of glia and developmental stage in modulating neuronal responses to TETS. Our results illustrate the applicability of the model for investigating how caged convulsants elicit abnormal network activity during the development of HN and HNG cultures in vitro.

Reconstruction of burst activity from calcium imaging of neuronal population via Lq minimization and interval screening

Calcium imaging is becoming an increasingly popular technology to indirectly measure activity patterns in local neuronal networks. Based on the dependence of calcium fluorescence on neuronal spiking, two-photon calcium imaging affords single-cell resolution of neuronal population activity. However, it is still difficult to reconstruct neuronal activity from complex calcium fluorescence traces, particularly for traces contaminated by noise. Here, we describe a robust and efficient neuronal-activity reconstruction method that utilizes Lq minimization and interval screening (IS), which we refer to as LqIS. The simulation results show that LqIS performs satisfactorily in terms of both accuracy and speed of reconstruction. Reconstruction of simulation and experimental data also shows that LqIS has advantages in terms of the recall rate, precision rate, and timing error. Finally, LqIS is demonstrated to effectively reconstruct neuronal burst activity from calcium fluorescence traces recorded from large-size neuronal population.