Spatiotemporal Activity Patterns Research Articles

Spatiotemporal resonance in mouse primary visual cortex

Human primary visual cortex (V1) responds more strongly, or resonates, when exposed to ∼10, ∼15-20, and ∼40-50Hz rhythmic flickering light. Full-field flicker also evokes the perception of hallucinatory geometric patterns, which mathematical models explain as standing-wave formations emerging from periodic forcing at resonant frequencies of the simulated neural network. However, empirical evidence for such flicker-induced standing waves in the visual cortex was missing. We recorded cortical responses to flicker in awake mice using high-spatial-resolution widefield imaging in combination with high-temporal-resolution glutamate-sensing fluorescent reporter (iGluSnFR). The temporal frequency tuningcurves in the mouse V1 were similar to those observed in humans, showing a banded structure with multiple resonance peaks (8, 15, and 33Hz). Spatially, all flicker frequencies evoked responses inV1 corresponding to retinotopic stimulus location, but some evoked additional peaks. These flicker-induced cortical patterns displayed standing-wave characteristics and matched linear wave equationsolutions in an area restricted to the visual cortex. Taken together, the interaction of periodic traveling waves with cortical area boundaries leads to spatiotemporal activity patterns that may affect perception.

Synanthropic rodents and urbanization processes: understanding the spatiotemporal pattern of rodent activity during urbanization works in a low-income neighborhood

Synanthropic rodents and urbanization processes: understanding the spatiotemporal pattern of rodent activity during urbanization works in a low-income neighborhood

Utilising activity patterns of a complex biophysical network model to optimise intra-striatal deep brain stimulation

A large-scale biophysical network model for the isolated striatal body is developed to optimise potential intrastriatal deep brain stimulation applied to, e.g. obsessive-compulsive disorder. The model is based on modified Hodgkin–Huxley equations with small-world connectivity, while the spatial information about the positions of the neurons is taken from a detailed human atlas. The model produces neuronal spatiotemporal activity patterns segregating healthy from pathological conditions. Three biomarkers were used for the optimisation of stimulation protocols regarding stimulation frequency, amplitude and localisation: the mean activity of the entire network, the frequency spectrum of the entire network (rhythmicity) and a combination of the above two. By minimising the deviation of the aforementioned biomarkers from the normal state, we compute the optimal deep brain stimulation parameters, regarding position, amplitude and frequency. Our results suggest that in the DBS optimisation process, there is a clear trade-off between frequency synchronisation and overall network activity, which has also been observed during in vivo studies.

Data-driven and equation-free methods for neurological disorders: analysis and control of the striatum network.

The striatum as part of the basal ganglia is central to both motor, and cognitive functions. Here, we propose a large-scale biophysical network for this part of the brain, using modified Hodgkin-Huxley dynamics to model neurons, and a connectivity informed by a detailed human atlas. The model shows different spatio-temporal activity patterns corresponding to lower (presumably normal) and increased cortico-striatal activation (as found in, e.g., obsessive-compulsive disorder), depending on the intensity of the cortical inputs. By applying equation-free methods, we are able to perform a macroscopic network analysis directly from microscale simulations. We identify the mean synaptic activity as the macroscopic variable of the system, which shows similarity with local field potentials. The equation-free approach results in a numerical bifurcation and stability analysis of the macroscopic dynamics of the striatal network. The different macroscopic states can be assigned to normal/healthy and pathological conditions, as known from neurological disorders. Finally, guided by the equation-free bifurcation analysis, we propose a therapeutic close loop control scheme for the striatal network.

Predator in proximity: how does a large carnivore respond to anthropogenic pressures at fine-scales? Implications for interface area management.

Driven by habitat loss and fragmentation, large carnivores are increasingly navigating human-dominated landscapes, where their activity is restricted and their behaviour altered. This movement, however, raises significant concerns and costs for people living nearby. While intricately linked, studies often isolate human and carnivore impacts, hindering effective management efforts. Hence, in this study, we brought these two into a common framework, focusing on an interface area between the critical tiger habitat and the human-dominated multiple-use buffer area of a central Indian protected area. We employed a fine-scale camera trap survey complemented by GPS-collar movement data to understand spatio-temporal activity patterns and adjustments of tigers in response to anthropogenic pressures. We used an occupancy framework to evaluate space use, Bayesian circular GLMs to model temporal activity, and home range and step length analyses to assess the movement patterns of tigers. Further, we used predation-risk models to understand conflict patterns as a function of tiger presence and other habitat variables. Despite disturbance, a high proportion of the sampled area was occupied by 17 unique tigers (ψ = 0.76; CI [0.73-0.92]). The distance to villages (β ± SE = 0.63 ± 0.21) and the relative abundance of large-bodied wild prey (β ± SE = 0.72 ± 0.37) emerged as key predictors of tiger space use probability, indicating a preference for wild prey by tigers, while human influences constrained their habitat utilisation. Distance to villages was also identified as the most significant predictor of the tigers' temporal activity (μ ± σ = 3.03 ± 0.06 rad) that exhibited higher nocturnality near villages. A total of 11% of tiger home ranges were within village boundaries, accompanied by faster movement in these areas (displacement 40-82% higher). Livestock depredation probability by tigers increased with proximity to villages (P = 0.002) and highway (P = 0.003). Although tiger space use probability (P = 0.056) and wild prey abundance (P = 0.134) were non-significant at the 0.05 threshold, their presence in the best-fit predation-risk model suggests their contextual relevance for understanding conflict risk. The results highlight the importance of appropriately managing livestock near human infrastructures to effectively mitigate conflict. Shared space of carnivores and humans requires dynamic site-specific actions grounded in evidence-based decision-making. This study emphasises the importance of concurrently addressing the intricate interactions between humans and large carnivores, particularly the latter's behavioural adaptations and role in conflict dynamics. Such an integrated approach is essential to unravel cause-effect relationships and promote effective interface management in human-dominated landscapes.

STEERING FROM THE REAR: COORDINATION OF CENTRAL PATTERN GENERATORS UNDERLYING NAVIGATION BY ASCENDING INTERNEURONS.

Understanding how animals coordinate movements to achieve goals is a fundamental pursuit in neuroscience. Here we explore how neurons that reside in posterior lower-order regions of a locomotor system project to anterior higher-order regions to influence steering and navigation. We characterized the anatomy and functional role of a population of ascending interneurons in the ventral nerve cord of Drosophila larvae. Through electron microscopy reconstructions and light microscopy, we determined that the cholinergic 19f cells receive input primarily from premotor interneurons and synapse upon a diverse array of postsynaptic targets within the anterior segments including other 19f cells. Calcium imaging of 19f activity in isolated central nervous system (CNS) preparations in relation to motor neurons revealed that 19f neurons are recruited into most larval motor programmes. 19f activity lags behind motor neuron activity and as a population, the cells encode spatio-temporal patterns of locomotor activity in the larval CNS. Optogenetic manipulations of 19f cell activity in isolated CNS preparations revealed that they coordinate the activity of central pattern generators underlying exploratory headsweeps and forward locomotion in a context and location specific manner. In behaving animals, activating 19f cells suppressed exploratory headsweeps and slowed forward locomotion, while inhibition of 19f activity potentiated headsweeps, slowing forward movement. Inhibiting activity in 19f cells ultimately affected the ability of larvae to remain in the vicinity of an odor source during an olfactory navigation task. Overall, our findings provide insights into how ascending interneurons monitor motor activity and shape interactions amongst rhythm generators underlying complex navigational tasks.

TAD border deletion at the Kit locus causes tissue-specific ectopic activation of a neighboring gene

Topologically associated domains (TADs) restrict promoter-enhancer interactions, thereby maintaining the spatiotemporal pattern of gene activity. However, rearrangements of the TADs boundaries do not always lead to significant changes in the activity pattern. Here, we investigated the consequences of the TAD boundaries deletion on the expression of developmentally important genes encoding tyrosine kinase receptors: Kit, Kdr, Pdgfra. We used genome editing in mice to delete the TADs boundaries at the Kit locus and characterized chromatin folding and gene expression in pure cultures of fibroblasts, mast cells, and melanocytes. We found that although Kit is highly active in both mast cells and melanocytes, deletion of the TAD boundary between the Kit and Kdr genes results in ectopic activation only in melanocytes. Thus, the epigenetic landscape, namely the mutual arrangement of enhancers and actively transcribing genes, is important for predicting the consequences of the TAD boundaries removal. We also found that mice without a TAD border between the Kit and Kdr genes have a phenotypic manifestation of the mutation — a lighter coloration. Thus, the data obtained shed light on the principles of interaction between the 3D chromatin organization and epigenetic marks in the regulation of gene activity.

The Effects of Locus Coeruleus Optogenetic Stimulation on Global Spatiotemporal Patterns in Rats.

Whole-brain intrinsic activity as detected by resting-state fMRI can be summarized by three primary spatiotemporal patterns. These patterns have been shown to change with different brain states, especially arousal. The noradrenergic locus coeruleus (LC) is a key node in arousal circuits and has extensive projections throughout the brain, giving it neuromodulatory influence over the coordinated activity of structurally separated regions. In this study, we used optogenetic-fMRI in rats to investigate the impact of LC stimulation on the global signal and three primary spatiotemporal patterns. We report small, spatially specific changes in global signal distribution as a result of tonic LC stimulation, as well as regional changes in spatiotemporal patterns of activity at 5 Hz tonic and 15 Hz phasic stimulation. We also found that LC stimulation had little to no effect on the spatiotemporal patterns detected by complex principal component analysis. These results show that the effects of LC activity on the BOLD signal in rats may be small and regionally concentrated, as opposed to widespread and globally acting.

Anti-Hebbian plasticity drives sequence learning in striatum

Spatio-temporal activity patterns have been observed in a variety of brain areas in spontaneous activity, prior to or during action, or in response to stimuli. Biological mechanisms endowing neurons with the ability to distinguish between different sequences remain largely unknown. Learning sequences of spikes raises multiple challenges, such as maintaining in memory spike history and discriminating partially overlapping sequences. Here, we show that anti-Hebbian spike-timing dependent plasticity (STDP), as observed at cortico-striatal synapses, can naturally lead to learning spike sequences. We design a spiking model of the striatal output neuron receiving spike patterns defined as sequential input from a fixed set of cortical neurons. We use a simple synaptic plasticity rule that combines anti-Hebbian STDP and non-associative potentiation for a subset of the presented patterns called rewarded patterns. We study the ability of striatal output neurons to discriminate rewarded from non-rewarded patterns by firing only after the presentation of a rewarded pattern. In particular, we show that two biological properties of striatal networks, spiking latency and collateral inhibition, contribute to an increase in accuracy, by allowing a better discrimination of partially overlapping sequences. These results suggest that anti-Hebbian STDP may serve as a biological substrate for learning sequences of spikes.

Spatiotemporal Organization of Prefrontal Norepinephrine Influences Neuronal Activity.

Norepinephrine (NE), a neuromodulator released by locus ceruleus (LC) neurons throughout the cortex, influences arousal and learning through extrasynaptic vesicle exocytosis. While NE within cortical regions has been viewed as a homogenous field, recent studies have demonstrated heterogeneous axonal dynamics and advances in GPCR-based fluorescent sensors permit direct observation of the local dynamics of NE at cellular scale. To investigate how the spatiotemporal dynamics of NE release in the prefrontal cortex (PFC) affect neuronal firing, we employed in vivo two-photon imaging of layer 2/3 of the PFC in order to observe fine-scale neuronal calcium and NE dynamics concurrently. In this proof of principle study, we found that local and global NE fields can decouple from one another, providing a substrate for local NE spatiotemporal activity patterns. Optic flow analysis revealed putative release and reuptake events which can occur at the same location, albeit at different times, indicating the potential to create a heterogeneous NE field. Utilizing generalized linear models, we demonstrated that cellular Ca2+ fluctuations are influenced by both the local and global NE field. However, during periods of local/global NE field decoupling, the local field drives cell firing dynamics rather than the global field. These findings underscore the significance of localized, phasic NE fluctuations for structuring cell firing, which may provide local neuromodulatory control of cortical activity.

Migraine aura discrimination using machine learning: an fMRI study during ictal and interictal periods.

Functional magnetic resonance imaging (fMRI) studies on migraine with aura are challenging due to the rarity of patients with triggered cases. This study optimized methodologies to explore differences in ictal and interictal spatiotemporal activation patterns based on visual stimuli using fMRI in two patients with unique aura triggers. Both patients underwent separate fMRI sessions during the ictal and interictal periods. The Gaussian Process Classifier (GPC) was used to differentiate these periods by employing a machine learning temporal embedding approach and spatiotemporal activation patterns based on visual stimuli. When restricted to visual and occipital regions, GPC had an improved performance, with accuracy rates for patients A and B of roughly 86-90% and 77-81%, respectively (p < 0.01). The algorithm effectively differentiated visual stimulation and rest periods and identified times when aura symptoms manifested, as evident from the varying predicted probabilities in the GPC models. These findings contribute to our understanding of the role of visual processing and brain activity patterns in migraine with aura and the significance of temporal embedding techniques in examining aura phenomena. This finding has implications for diagnostic tools and therapeutic techniques, especially for patients suffering from aura symptoms.

Flying by the river side: Survey of bat distributions and environmental contexts along a 1000‐mile river corridor, Green and Colorado Rivers, USA

AbstractAimEmerging research shows how bioindicators, specifically bats, can serve as a means for monitoring conservation and management of riparian corridors for multiple taxonomic groups. To track changes in the composition or abundance of bioindicator species, researchers must attain a baseline in species presence and relative activity. We examined the spatial and temporal patterns of bat community composition and activity along a 1000‐mile river corridor to determine species diversity trends by latitude and habitat.LocationColorado River Basin.MethodsHere we describe the results from an acoustic bat survey conducted opportunistically on the 2019 Sesquicentennial Colorado River Exploring Expedition. This broad, 1000‐mile survey provides a baseline for species distributions over a large geographic range.ResultsIn total, we collected 63 nights of acoustic data over 70 days and recorded over 59,000 files equating to 45,363 caLL files (≥2 pulses). 18,490 (41% of caLL files) were identified as species (n = 19 bat species). We applied non‐metric multidimensional scaling to characterize spatiotemporal patterns of activity between species, as well as compared bat activity among river features and local environmental conditions (i.e., temperature and time since sunset) using an information theoretic approach.ConclusionSpecies composition varied by physiographic region and adjacent river habitat, thus providing a quantifiable measure of determining habitat quality along this major river system and providing baseline information for using bats as bioindicators of habitat quality.

Human mobility reshaped? Deciphering the impacts of the Covid-19 pandemic on activity patterns, spatial habits, and schedule habits

Despite the historically documented regularity in human mobility patterns, the relaxation of spatial and temporal constraints, brought by the widespread adoption of telecommuting and e-commerce during the COVID-19 pandemic, as well as a growing desire for flexible work arrangements in a post-pandemic work, indicates a potential reshaping of these patterns. In this paper, we investigate the multifaceted impacts of relaxed spatio-temporal constraints on human mobility, using well-established metrics from the travel behavior literature. Further, we introduce a novel metric for schedule regularity, accounting for specific day-of-week characteristics that previous approaches overlooked. Building on the large body of literature on the impacts of COVID-19 on human mobility, we make use of passively tracked Point of Interest (POI) data for approximately 21,700 smartphone users in the US, and analyze data between January 2020 and September 2022 to answer two key questions: (1) has the COVID-19 pandemic and its associated relaxation of spatio-temporal activity patterns reshaped the different aspects of human mobility, and (2) have we achieved a state of stable post-pandemic “new normal”? We hypothesize that the relaxation of the spatiotemporal constraints around key activities will result in people exhibiting less regular schedules. Findings reveal a complex landscape: while some mobility indicators have reverted to pre-pandemic norms, such as trip frequency and travel distance, others, notably at-home dwell-time, persist at altered levels, suggesting a recalibration rather than a return to past behaviors. Most notably, our analysis reveals a paradox: despite the documented large-scale shift towards flexible work arrangements, schedule habits have strengthened rather than relaxed, defying our initial hypotheses and highlighting a desire for regularity. The study’s results contribute to a deeper understanding of the post-pandemic “new normal”, offering key insights on how multiple facets of travel behavior were reshaped, if at all, by the COVID-19 pandemic, and will help inform transportation planning in a post-pandemic world.

A Systematic Review of the Potential Influence of Urbanization on the Regional Thunderstorm Process and Lightning Activity

In the context of global climate change, lightning disasters have emerged as a serious environmental factor that restricts the sustainable development of megacities. This paper provides a review of the research on the impact of urbanization on thunderstorm processes and lightning activity, exploring various aspects, such as aerosols, urban thermal effects, urban dynamic effects, and building morphology. Despite numerous significant achievements in the study of the impact of air pollutants on lightning activity, there is no consensus on whether aerosols serve to enhance or inhibit lightning activity. The temperature difference between the urban underlying surface and the natural underlying surface could sustain and promote the occurrence and development of convective systems, thus enhancing lightning activity. In terms of urban dynamics, the barrier effect has led to the maximum center of lightning appearing at the edge of a built-up area, which might be associated with factors, such as urban heat island (UHI) intensity, wind speed, synoptic background, and city size. Additionally, the size of a city and the height of the buildings was also an influencing factor on lightning activity. In summary, scholars have made progress in understanding the characteristics and drivers of urban lightning activity in recent years, but there are still some urgent problems that need to be solved: (1) How to analyze, comprehensively, the spatiotemporal patterns of urban lightning activity under different thunderstorm intensity backgrounds? (2) How to conduct analysis to investigate the influence of alterations in the boundary layer structure, water–heat energy balance, and water vapor circulation processes on urban lightning activity in the context of urbanization? (3) How to couple numerical models of different scales to enhance the understanding of the impact of complex underlying surfaces on urban lightning activity? Future studies could investigate the relationship between urbanization and thunderstorm/lightning activity using a combination of observational data, numerical modeling, and laboratory experiments, which holds promise for providing valuable theoretical insights and technical support to enhance the prediction, nowcasting, early warning, and risk assessment of thunderstorms and lightning in urban areas.

Alpha and Beta Oscillations Differentially Support Word Production in a Rule-Switching Task.

Research into the role of brain oscillations in basic perceptual and cognitive functions has suggested that the alpha rhythm reflects functional inhibition while the beta rhythm reflects neural ensemble (re)activation. However, little is known regarding the generalization of these proposed fundamental operations to linguistic processes, such as speech comprehension and production. Here, we recorded magnetoencephalography in participants performing a novel rule-switching paradigm. Specifically, Dutch native speakers had to produce an alternative exemplar from the same category or a feature of a given target word embedded in spoken sentences (e.g., for the word "tuna", an exemplar from the same category-"seafood"-would be "shrimp", and a feature would be "pink"). A cue indicated the task rule-exemplar or feature-either before (pre-cue) or after (retro-cue) listening to the sentence. Alpha power during the working memory delay was lower for retro-cue compared with that for pre-cue in the left hemispheric language-related regions. Critically, alpha power negatively correlated with reaction times, suggestive of alpha facilitating task performance by regulating inhibition in regions linked to lexical retrieval. Furthermore, we observed a different spatiotemporal pattern of beta activity for exemplars versus features in the right temporoparietal regions, in line with the proposed role of beta in recruiting neural networks for the encoding of distinct categories. Overall, our study provides evidence for the generalizability of the role of alpha and beta oscillations from perceptual to more "complex, linguistic processes" and offers a novel task to investigate links between rule-switching, working memory, and word production.

Putative Spatiotemporal Changes in Inhibitor Activity during Cold Stratification of Sapium sebiferum Seeds

Putative Spatiotemporal Changes in Inhibitor Activity during Cold Stratification of Sapium sebiferum Seeds

A review of bovine tuberculosis transmission risk in European wildlife communities

Abstract Understanding how disease moves through wildlife communities is essential to managing outbreaks of zoonotic diseases across the globe. Bovine tuberculosis is a disease caused by the bacterium Mycobacterium bovis that can threaten domestic and wildlife species. The mechanism by which Mycobacterium bovis is spread between species is still poorly understood. Previous reviews are limited in the breadth of species considered and are primarily concerned with transmission from wildlife to domestic species. We conducted a review and analysis of Mycobacterium bovis prevalence rates in European wildlife species to identify species of concern for the transmission of bovine tuberculosis in a wildlife community. We subsequently conducted a narrative review of these species assessing the risk of Mycobacterium bovis transmission in a wildlife community based on available literature. We calculated weighted mean disease prevalence rates to be highest in fallow deer (Dama dama, 20%), Eurasian badgers (Meles meles, 11%), wild boar (Sus scrofa, 9%) and red foxes (Vulpes vulpes, 4%). We considered these species to be of particular concern for the transmission of Mycobacterium bovis and selected them as the focus of our narrative review and risk assessment. Our risk assessment considered disease pathology, spatiotemporal activity patterns and animal behaviour as factors affecting the likelihood of Mycobacterium bovis transmission between wildlife species. We found that prior research has principally focused on a few individual species, but that Mycobacterium bovis transmission through a wildlife community is likely more complex. We determined that disease transmission between multiple species may compound the severity of an outbreak of bovine tuberculosis. Broad, multi‐species sampling campaigns and standardised Mycobacterium bovis testing protocols should be implemented in future studies. We also determined that an in‐depth analysis of spatiotemporal overlap between species was needed to better assess the risk of transmission between wildlife species.

Neural signatures of shared subjective affective engagement and disengagement during movie viewing.

When watching a negative emotional movie, we differ from person to person in the ease with which we engage and the difficulty with which we disengage throughout a temporally evolving narrative. We investigated neural responses of emotional processing, by considering inter-individual synchronization in subjective emotional engagement and disengagement. The neural underpinnings of these shared responses are ideally studied in naturalistic scenarios like movie viewing, wherein individuals emotionally engage and disengage at their own time and pace throughout the course of a narrative. Despite the rich data that naturalistic designs can bring to the study, there is a challenge in determining time-resolved behavioral markers of subjective engagement and disengagement and their underlying neural responses. We used a within-subject cross-over design instructing 22 subjects to watch clips of either neutral or sad content while undergoing functional magnetic resonance imaging (fMRI). Participants watched the same movies a second time while continuously annotating the perceived emotional intensity, thus enabling the mapping of brain activity and emotional experience. Our analyses revealed that between-participant similarity in waxing (engagement) and waning (disengagement) of emotional intensity was directly related to the between-participant similarity in spatiotemporal patterns of brain activation during the movie(s). Similar patterns of engagement reflected common activation in the bilateral ventromedial prefrontal cortex, regions often involved in self-referenced evaluation and generation of negative emotions. Similar patterns of disengagement reflected common activation in central executive and default mode network regions often involved in top-down emotion regulation. Together this work helps to better understand cognitive and neural mechanisms underpinning engagement and disengagement from emotionally evocative narratives.

LUR modeling of long-term average hourly concentrations of NO2 using hyperlocal mobile monitoring data

Mobile monitoring campaigns have effectively captured spatial hyperlocal variations in long-term average concentrations of regulated and unregulated air pollutants. However, their application in estimating spatiotemporally varying maps has rarely been investigated. Tackling this gap, we investigated whether mobile measurements can assess long-term average nitrogen dioxide (NO2) concentrations for each hour of the day. Using mobile NO2 data monitored for 10 months in Amsterdam, we examined the performance of two spatiotemporal land use regression (LUR) methods, Spatiotemporal-Kriging and GTWR (Geographical and Temporal Weighted Regression), alongside two classical spatial LUR models developed separately for each hour. We found that mobile measurements follow the general pattern of fixed-site measurements, but with considerable deviations (indicating collection uncertainty). Leveraging heterogeneous spatiotemporal autocorrelations, GTWR smoothed these deviations and achieved an overall performance of an R2 of 0.49 and a Mean Absolute Error of 6.33 μg/m3, validated by long-term fixed-site measurements (out-of-sample). The other models tested were more affected by the collection uncertainty. We highlighted that the spatiotemporal variations captured in mobile measurements can be used to reconstruct long-term average hourly air pollution maps. These maps facilitate dynamic exposure assessments considering spatiotemporal human activity patterns.

A method for detecting spatiotemporal patterns of cancer biomarkers-evoked activity using radial basis function network extracted time-domain features from calcium imaging data

BackgroundTwo-photon calcium imaging is widely used to study the odor-evoked glomerular activity in the dorsal olfactory bulb of macrosmatic animals. The nonstationary character of activated patterns sets a limit on the use of a traditional image processing approaches. New methodThe developed method makes it possible to automatically map cancer biomarkers-activated glomeruli in the rat dorsal olfactory bulb. We interpolated fluorescence intensity of calcium dynamics based on the Gaussian RBF network and synthesized the physiological fluorescence model of the receptive glomerular field. ResultsThe experiments on 5 rats confirmed the correctness of the developed approach. Patterns evoked by the 6-methyl-5-hepten-2-one (stomach cancer biomarker) and benzene (lung cancer biomarker) were correctly identified. Comparison with existing methodsThe proposed method was compared with the nonnegative matrix factorization method and with the method based on computer vision algorithms. The developed approach showed better accuracy in experiments and provided the mathematical models of the odor-evoked patterns synthesis. These models can be used to generate synthetic images of odor-evoked glomerular activity and thus to overcome the problem of small experimental data collected in calcium imaging. ConclusionsThe proposed method should be considered part of the toolkit for fully automatic analysis of calcium imaging-based studies. Currently available methodology is not able to use breath biomarkers to reliably discriminate between cancer patients and healthy controls. Nevertheless, the effective identification of the spatial patterns of cancer biomarkers-evoked glomerular activity can serve as the foundation for highly sensitive biohybrid systems for cancer screening.