KCSD-python, reliable current source density estimation with quality control.
Ontology highlight
ABSTRACT:
SUBMITTER: Chintaluri C
PROVIDER: S-EPMC10965101 | biostudies-literature | 2024 Mar
REPOSITORIES: biostudies-literature
Publications
kCSD-python, reliable current source density estimation with quality control.
Chintaluri Chaitanya C Bejtka Marta M Średniawa Władysław W Czerwiński Michał M Dzik Jakub M JM Jędrzejewska-Szmek Joanna J Wójcik Daniel K DK
PLoS computational biology 20240314 3
Interpretation of extracellular recordings can be challenging due to the long range of electric field. This challenge can be mitigated by estimating the current source density (CSD). Here we introduce kCSD-python, an open Python package implementing Kernel Current Source Density (kCSD) method and related tools to facilitate CSD analysis of experimental data and the interpretation of results. We show how to counter the limitations imposed by noise and assumptions in the method itself. kCSD-python ...[more]
Similar Datasets
Project description:Custom-built microscopes often require control of multiple hardware devices and precise hardware coordination. It is also desirable to have a solution that is scalable to complex systems and that is translatable between components from different manufacturers. Here we report Python-Microscope, a free and open-source Python library for high-performance control of arbitrarily complex and scalable custom microscope systems. Python-Microscope offers simple to use Python-based tools, abstracting differences between physical devices by providing a defined interface for different device types. Concrete implementations are provided for a range of specific hardware, and a framework exists for further expansion. Python-Microscope supports the distribution of devices over multiple computers while maintaining synchronisation via highly precise hardware triggers. We discuss the architectural features of Python-Microscope that overcome the performance problems often raised against Python and demonstrate the different use cases that drove its design: integration with user-facing projects, namely the Microscope-Cockpit project; control of complex microscopes at high speed while using the Python programming language; and use as a microscope simulation tool for software development.
| S-EPMC8520730 | biostudies-literature
Project description:Electroencephalography (EEG) and source estimation can be used to identify brain areas activated during a task, which could offer greater insight on cortical dynamics. Source estimation requires knowledge of the locations of the EEG electrodes. This could be provided with a template or obtained by digitizing the EEG electrode locations. Operator skill and inherent uncertainties of a digitizing system likely produce a range of digitization reliabilities, which could affect source estimation and the interpretation of the estimated source locations. Here, we compared the reliabilities of five digitizing methods (ultrasound, structured-light 3D scan, infrared 3D scan, motion capture probe, and motion capture) and determined the relationship between digitization reliability and source estimation uncertainty, assuming other contributors to source estimation uncertainty were constant. We digitized a mannequin head using each method five times and quantified the reliability and validity of each method. We created five hundred sets of electrode locations based on our reliability results and applied a dipole fitting algorithm (DIPFIT) to perform source estimation. The motion capture method, which recorded the locations of markers placed directly on the electrodes had the best reliability with an average electrode variability of 0.001 cm. Then, in order of decreasing reliability were the method using a digitizing probe in the motion capture system, an infrared 3D scanner, a structured-light 3D scanner, and an ultrasound digitization system. Unsurprisingly, uncertainty of the estimated source locations increased with greater variability of EEG electrode locations and less reliable digitizing systems. If EEG electrode location variability was ∽1 cm, a single source could shift by as much as 2 cm. To help translate these distances into practical terms, we quantified Brodmann area accuracy for each digitizing method and found that the average Brodmann area accuracy for all digitizing methods was >80%. Using a template of electrode locations reduced the Brodmann area accuracy to ∽50%. Overall, more reliable digitizing methods can reduce source estimation uncertainty, but the significance of the source estimation uncertainty depends on the desired spatial resolution. For accurate Brodmann area identification, any of the digitizing methods tested can be used confidently.
| S-EPMC6856631 | biostudies-literature
Project description:Fitting complex models to epidemiological data is a challenging problem: methodologies can be inaccessible to all but specialists, there may be challenges in adequately describing uncertainty in model fitting, the complex models may take a long time to run, and it can be difficult to fully capture the heterogeneity in the data. We develop an adaptive approximate Bayesian computation scheme to fit a variety of epidemiologically relevant data with minimal hyper-parameter tuning by using an adaptive tolerance scheme. We implement a novel kernel density estimation scheme to capture both dispersed and multi-dimensional data, and directly compare this technique to standard Bayesian approaches. We then apply the procedure to a complex individual-based simulation of lymphatic filariasis, a human parasitic disease. The procedure and examples are released alongside this article as an open access library, with examples to aid researchers to rapidly fit models to data. This demonstrates that an adaptive ABC scheme with a general summary and distance metric is capable of performing model fitting for a variety of epidemiological data. It also does not require significant theoretical background to use and can be made accessible to the diverse epidemiological research community.
| S-EPMC6227249 | biostudies-literature
Project description:Electrocardiogram (ECG) and photoplethysmogram (PPG) are commonly used to determine the vital signs of heart rate, respiratory rate, and oxygen saturation in patient monitoring. In addition to simple observation of those summarized indexes, waveform signals can be analyzed to provide deeper insights into disease pathophysiology and support clinical decisions. Such data, generated from continuous patient monitoring from both conventional bedside and low-cost wearable monitors, are increasingly accessible. However, the recorded waveforms suffer from considerable noise and artifacts and, hence, are not necessarily used prior to certain quality control (QC) measures, especially by those with limited programming experience. Various signal quality indices (SQIs) have been proposed to indicate signal quality. To facilitate and harmonize a wider usage of SQIs in practice, we present a Python package, named vital_sqi, which provides a unified interface to the state-of-the-art SQIs for ECG and PPG signals. The vital_sqi package provides with seven different peak detectors and access to more than 70 SQIs by using different settings. The vital_sqi package is designed with pipelines and graphical user interfaces to enable users of various programming fluency to use the package. Multiple SQI extraction pipelines can take the PPG and ECG waveforms and generate a bespoke SQI table. As these SQI scores represent the signal features, they can be input in any quality classifier. The package provides functions to build simple rule-based decision systems for signal segment quality classification using user-defined SQI thresholds. An experiment with a carefully annotated PPG dataset suggests thresholds for relevant PPG SQIs.
| S-EPMC9692103 | biostudies-literature
Project description:Single cell whole-genome sequencing (scWGS) is providing novel insights into the nature of genetic heterogeneity in normal and diseased cells. However, the whole-genome amplification process required for scWGS introduces biases into the resulting sequencing that can confound downstream analysis. Here, we present a statistical method, with an accompanying package PaSD-qc (Power Spectral Density-qc), that evaluates the properties and quality of single cell libraries. It uses a modified power spectral density to assess amplification uniformity, amplicon size distribution, autocovariance and inter-sample consistency as well as to identify chromosomes with aberrant read-density profiles due either to copy alterations or poor amplification. These metrics provide a standard way to compare the quality of single cell samples as well as yield information necessary to improve variant calling strategies. We demonstrate the usefulness of this tool in comparing the properties of scWGS protocols, identifying potential chromosomal copy number variation, determining chromosomal and subchromosomal regions of poor amplification, and selecting high-quality libraries from low-coverage data for deep sequencing. The software is available free and open-source at https://github.com/parklab/PaSDqc.
| S-EPMC5829578 | biostudies-literature
Project description:We describe the first mapping of biological current in a live heart using ultrasound current source density imaging (UCSDI). Ablation procedures that treat severe heart arrhythmias require detailed maps of the cardiac activation wave. The conventional procedure is time-consuming and limited by its poor spatial resolution (5-10 mm). UCSDI can potentially improve on existing mapping procedures. It is based on a pressure-induced change in resistivity known as the acousto-electric (AE) effect, which is spatially confined to the ultrasound focus. Data from 2 experiments are presented. A 540 kHz ultrasonic transducer (f/# = 1, focal length = 90 mm, pulse repetition frequency = 1600 Hz) was scanned over an isolated rabbit heart perfused with an excitation-contraction decoupler to reduce motion significantly while retaining electric function. Tungsten electrodes inserted in the left ventricle recorded simultaneously the AE signal and the low-frequency electrocardiogram (ECG). UCSDI displayed spatial and temporal patterns consistent with the spreading activation wave. The propagation velocity estimated from UCSDI was 0.25 +/- 0.05 mm/ms, comparable to the values obtained with the ECG signals. The maximum AE signal-to-noise ratio after filtering was 18 dB, with an equivalent detection threshold of 0.1 mA/ cm(2). This study demonstrates that UCSDI is a potentially powerful technique for mapping current flow and biopotentials in the heart.
| S-EPMC2823813 | biostudies-literature
Project description:Despite the widespread use of current-source density (CSD) analysis of extracellular potential recordings in the brain, the physical mechanisms responsible for the generation of the signal are still debated. While the extracellular potential is thought to be exclusively generated by the transmembrane currents, recent studies suggest that extracellular diffusive, advective and displacement currents-traditionally neglected-may also contribute considerably toward extracellular potential recordings. Here, we first justify the application of the electro-quasistatic approximation of Maxwell's equations to describe the electromagnetic field of physiological origin. Subsequently, we perform spatial averaging of currents in neural tissue to arrive at the notion of the CSD and derive an equation relating it to the extracellular potential. We show that, in general, the extracellular potential is determined by the CSD of membrane currents as well as the gradients of the putative extracellular diffusion current. The diffusion current can contribute significantly to the extracellular potential at frequencies less than a few Hertz; in which case it must be subtracted to obtain correct CSD estimates. We also show that the advective and displacement currents in the extracellular space are negligible for physiological frequencies while, within cellular membrane, displacement current contributes toward the CSD as a capacitive current. Taken together, these findings elucidate the relationship between electric currents and the extracellular potential in brain tissue and form the necessary foundation for the analysis of extracellular recordings.
| S-EPMC5413824 | biostudies-literature
Project description:BackgroundDespite recent success in pharmacologic treatment of depression, the inability to predict individual treatment response remains a liability. This study replicates and extends findings relating pretreatment electroencephalographic (EEG) alpha to treatment outcomes for serotonergic medications.MethodsResting EEG (eyes-open and eyes-closed) was recorded from a 67-electrode montage in 41 unmedicated depressed patients and 41 healthy control subjects. Patients were tested before receiving antidepressants including a serotonergic mode of action (selective serotonin reuptake inhibitor [SSRI], serotonin and norepinephrine reuptake inhibitor, or SSRI plus norepinephrine and dopamine reuptake inhibitor). EEG was quantified by frequency principal components analysis of spectra derived from reference-free current source density (CSD) waveforms, which sharpens and simplifies EEG topographies, disentangles them from artifact, and yields measures that more closely represent underlying neuronal current generators.ResultsPatients who did not respond to treatment had significantly less alpha CSD compared with responders or healthy control subjects, localizable to well-defined posterior generators. The alpha difference between responders and nonresponders was greater for eyes-closed than eyes-open conditions and was present across alpha subbands. A classification criterion based on the median alpha for healthy control subjects showed good positive predictive value (93.3) and specificity (92.3). There was no evidence of differential value for predicting response to an SSRI alone or dual treatment targeting serotonergic plus other monoamine neurotransmitters.ConclusionsFindings confirm the value of EEG alpha amplitude as a viable predictor of antidepressant response and suggest that personalized treatments for depression may be identified using simple electrophysiologic CSD measures.
| S-EPMC3142299 | biostudies-literature
Project description:Computational models of cardiac electrophysiology provided insights into arrhythmogenesis and paved the way toward tailored therapies in the last years. To fully leverage in silico models in future research, these models need to be adapted to reflect pathologies, genetic alterations, or pharmacological effects, however. A common approach is to leave the structure of established models unaltered and estimate the values of a set of parameters. Today's high-throughput patch clamp data acquisition methods require robust, unsupervised algorithms that estimate parameters both accurately and reliably. In this work, two classes of optimization approaches are evaluated: gradient-based trust-region-reflective and derivative-free particle swarm algorithms. Using synthetic input data and different ion current formulations from the Courtemanche et al. electrophysiological model of human atrial myocytes, we show that neither of the two schemes alone succeeds to meet all requirements. Sequential combination of the two algorithms did improve the performance to some extent but not satisfactorily. Thus, we propose a novel hybrid approach coupling the two algorithms in each iteration. This hybrid approach yielded very accurate estimates with minimal dependency on the initial guess using synthetic input data for which a ground truth parameter set exists. When applied to measured data, the hybrid approach yielded the best fit, again with minimal variation. Using the proposed algorithm, a single run is sufficient to estimate the parameters. The degree of superiority over the other investigated algorithms in terms of accuracy and robustness depended on the type of current. In contrast to the non-hybrid approaches, the proposed method proved to be optimal for data of arbitrary signal to noise ratio. The hybrid algorithm proposed in this work provides an important tool to integrate experimental data into computational models both accurately and robustly allowing to assess the often non-intuitive consequences of ion channel-level changes on higher levels of integration.
| S-EPMC4710757 | biostudies-literature
Project description:Pharmacokinetic (PK) parameter estimation is a critical and complex step in the model-informed precision dosing (MIPD) approach. The mapbayr package was developed to perform maximum a posteriori Bayesian estimation (MAP-BE) in R from any population PK model coded in mrgsolve. The performances of mapbayr were assessed using two approaches. First, "test" models with different features were coded, for example, first-order and zero-order absorption, lag time, time-varying covariates, Michaelis-Menten elimination, combined and exponential residual error, parent drug and metabolite, and small or large inter-individual variability (IIV). A total of 4000 PK profiles (combining single/multiple dosing and rich/sparse sampling) were simulated from each test model, and MAP-BE of parameters was performed in both mapbayr and NONMEM. Second, a similar procedure was conducted with seven "real" previously published models to compare mapbayr and NONMEM on a PK outcome used in MIPD. For the test models, 98% of mapbayr estimations were identical to those given by NONMEM. Some discordances could be observed when dose-related parameters were estimated or when models with large IIV were used. The exploration of objective function values suggested that mapbayr might outdo NONMEM in specific cases. For the real models, a concordance close to 100% on PK outcomes was observed. The mapbayr package provides a reliable solution to perform MAP-BE of PK parameters in R. It also includes functions dedicated to data formatting and reporting and enables the creation of standalone Shiny web applications dedicated to MIPD, whatever the model or the clinical protocol and without additional software other than R.
| S-EPMC8520754 | biostudies-literature