A basic bifurcation structure from bursting to spiking of injured nerve fibers in a two-dimensional parameter space

Two different bifurcation scenarios of firing patterns with decreasing extracellular calcium concentrations were observed in identical sciatic nerve fibers of a chronic constriction injury (CCI) model when the extracellular 4-aminopyridine concentrations were fixed at two different levels. Both processes proceeded from period-1 bursting to period-1 spiking via complex or simple processes. Multiple typical experimental examples manifested dynamics closely matching those simulated in a recently proposed 4-dimensional model to describe the nonlinear dynamics of the CCI model, which included most cases of the bifurcation scenarios. As the extracellular 4-aminopyridine concentrations is increased, the structure of the bifurcation scenario becomes more complex. The results provide a basic framework for identifying the relationships between different neural firing patterns and different bifurcation scenarios and for revealing the complex nonlinear dynamics of neural firing patterns. The potential roles of the basic bifurcation structures in identifying the information process mechanism are discussed.     

Unstable dynamics of a periodically driven oscillator in the presence of noise

Two qualitatively different unstable dynamical behaviours are shown to arise from the application of a periodic input to a simple mathematical model of an oscillator in the presence of noise. Rhythms similar to quasiperiodic dynamics may arise when there is a low amplitude periodic input, while with high amplitude inputs, patterns with irregular skipped or intercalated beats are found. These two qualitatively different types of unstable dynamics are similar to those observed in the respiratory activity of mechanically ventilated cats. A number of numerical simulations are performed to illustrate the quantitative properties of the two unstable patterns and to show how the quantitative properties can be compared with experimental data.     

Reproducibility and discriminability of brain patterns of semantic categories enhanced by congruent audiovisual stimuli

One of the central questions in cognitive neuroscience is the precise neural representation, or brain pattern, associated with a semantic category. In this study, we explored the influence of audiovisual stimuli on the brain patterns of concepts or semantic categories through a functional magnetic resonance imaging (fMRI) experiment. We used a pattern search method to extract brain patterns corresponding to two semantic categories: "old people" and "young people." These brain patterns were elicited by semantically congruent audiovisual, semantically incongruent audiovisual, unimodal visual, and unimodal auditory stimuli belonging to the two semantic categories. We calculated the reproducibility index, which measures the similarity of the patterns within the same category. We also decoded the semantic categories from these brain patterns. The decoding accuracy reflects the discriminability of the brain patterns between two categories. The results showed that both the reproducibility index of brain patterns and the decoding accuracy were significantly higher for semantically congruent audiovisual stimuli than for unimodal visual and unimodal auditory stimuli, while the semantically incongruent stimuli did not elicit brain patterns with significantly higher reproducibility index or decoding accuracy. Thus, the semantically congruent audiovisual stimuli enhanced the within-class reproducibility of brain patterns and the between-class discriminability of brain patterns, and facilitate neural representations of semantic categories or concepts. Furthermore, we analyzed the brain activity in superior temporal sulcus and middle temporal gyrus (STS/MTG). The strength of the fMRI signal and the reproducibility index were enhanced by the semantically congruent audiovisual stimuli. Our results support the use of the reproducibility index as a potential tool to supplement the fMRI signal amplitude for evaluating multimodal integration.     

Variation in moisture duration as a driver of coexistence by the storage effect in desert annual plants

Temporal environmental variation is a leading hypothesis for the coexistence of desert annual plants. Environmental variation is hypothesized to cause species-specific patterns of variation in germination, which then generates the storage effect coexistence mechanism. However, it has never been shown how sufficient species differences in germination patterns for multispecies coexistence can arise from a shared fluctuating environment. Here we show that nonlinear germination responses to a single fluctuating physical environmental factor can lead to sufficient differences between species in germination pattern for the storage effect to yield coexistence of multiple species. We derive these nonlinear germination responses from experimental data on the effects of varying soil moisture duration. Although these nonlinearities lead to strong species asymmetries in germination patterns, the relative nonlinearity coexistence mechanism is minor compared with the storage effect. However, these asymmetries mean that the storage effect can be negative for some species, which then only persist in the face of interspecific competition through average fitness advantages. This work shows how a low dimensional physical environment can nevertheless stabilize multispecies coexistence when the species have different nonlinear responses to common conditions, as supported by our experimental data.     

Extracting representations of cognition across neuroimaging studies improves brain decoding

Cognitive brain imaging is accumulating datasets about the neural substrate of many different mental processes. Yet, most studies are based on few subjects and have low statistical power. Analyzing data across studies could bring more statistical power; yet the current brain-imaging analytic framework cannot be used at scale as it requires casting all cognitive tasks in a unified theoretical framework. We introduce a new methodology to analyze brain responses across tasks without a joint model of the psychological processes. The method boosts statistical power in small studies with specific cognitive focus by analyzing them jointly with large studies that probe less focal mental processes. Our approach improves decoding performance for 80% of 35 widely-different functional-imaging studies. It finds commonalities across tasks in a data-driven way, via common brain representations that predict mental processes. These are brain networks tuned to psychological manipulations. They outline interpretable and plausible brain structures. The extracted networks have been made available; they can be readily reused in new neuro-imaging studies. We provide a multi-study decoding tool to adapt to new data.     

Decoding Multiple Sound Categories in the Human Temporal Cortex Using High Resolution fMRI

Perception of sound categories is an important aspect of auditory perception. The extent to which the brain’s representation of sound categories is encoded in specialized subregions or distributed across the auditory cortex remains unclear. Recent studies using multivariate pattern analysis (MVPA) of brain activations have provided important insights into how the brain decodes perceptual information. In the large existing literature on brain decoding using MVPA methods, relatively few studies have been conducted on multi-class categorization in the auditory domain. Here, we investigated the representation and processing of auditory categories within the human temporal cortex using high resolution fMRI and MVPA methods. More importantly, we considered decoding multiple sound categories simultaneously through multi-class support vector machine-recursive feature elimination (MSVM-RFE) as our MVPA tool. Results show that for all classifications the model MSVM-RFE was able to learn the functional relation between the multiple sound categories and the corresponding evoked spatial patterns and classify the unlabeled sound-evoked patterns significantly above chance. This indicates the feasibility of decoding multiple sound categories not only within but across subjects. However, the across-subject variation affects classification performance more than the within-subject variation, as the across-subject analysis has significantly lower classification accuracies. Sound category-selective brain maps were identified based on multi-class classification and revealed distributed patterns of brain activity in the superior temporal gyrus and the middle temporal gyrus. This is in accordance with previous studies, indicating that information in the spatially distributed patterns may reflect a more abstract perceptual level of representation of sound categories. Further, we show that the across-subject classification performance can be significantly improved by averaging the fMRI images over items, because the irrelevant variations between different items of the same sound category are reduced and in turn the proportion of signals relevant to sound categorization increases.     

Decoding Accuracy in Supplementary Motor Cortex Correlates with Perceptual Sensitivity to Tactile Roughness

Perceptual sensitivity to tactile roughness varies across individuals for the same degree of roughness. A number of neurophysiological studies have investigated the neural substrates of tactile roughness perception, but the neural processing underlying the strong individual differences in perceptual roughness sensitivity remains unknown. In this study, we explored the human brain activation patterns associated with the behavioral discriminability of surface texture roughness using functional magnetic resonance imaging (fMRI). First, a whole-brain searchlight multi-voxel pattern analysis (MVPA) was used to find brain regions from which we could decode roughness information. The searchlight MVPA revealed four brain regions showing significant decoding results: the supplementary motor area (SMA), contralateral postcentral gyrus (S1), and superior portion of the bilateral temporal pole (STP). Next, we evaluated the behavioral roughness discrimination sensitivity of each individual using the just-noticeable difference (JND) and correlated this with the decoding accuracy in each of the four regions. We found that only the SMA showed a significant correlation between neuronal decoding accuracy and JND across individuals; Participants with a smaller JND (i.e., better discrimination ability) exhibited higher decoding accuracy from their voxel response patterns in the SMA. Our findings suggest that multivariate voxel response patterns presented in the SMA represent individual perceptual sensitivity to tactile roughness and people with greater perceptual sensitivity to tactile roughness are likely to have more distinct neural representations of different roughness levels in their SMA.     

Three-dimensional static modeling of the lumbar spine

This paper presents three-dimensional static modeling of the human lumbar spine to be used in the formation of anatomically-correct movement patterns for a fully cable-actuated robotic lumbar spine which can mimic in vivo human lumbar spine movements to provide better hands-on training for medical students. The mathematical model incorporates five lumbar vertebrae between the first lumbar vertebra and the sacrum, with dimensions of an average adult human spine. The vertebrae are connected to each other by elastic elements, torsional springs and a spherical joint located at the inferoposterior corner in the mid-sagittal plane of the vertebral body. Elastic elements represent the ligaments that surround the facet joints and the torsional springs represent the collective effect of intervertebral disc which plays a major role in balancing torsional load during upper body motion and the remaining ligaments that support the spinal column. The elastic elements and torsional springs are considered to be nonlinear. The nonlinear stiffness constants for six motion types were solved using a multiobjective optimization technique. The quantitative comparison between the angles of rotations predicted by the proposed model and in the experimental data confirmed that the model yields angles of rotation close to the experimental data. The main contribution is that the new model can be used for all motions while the experimental data was only obtained at discrete measurement points.     

A Simple Model for the Cold Requirement of the Tulip

A simple model is proposed that qualitatively describes someimportant features of the cold requirement of the tulip. Themodel may be used as a basis for quantitative experimental investigationsinto the response of tulip growth and development to differentperiods of bulb storage at low temperatures.     

A fully nonlinear viscohyperelastic model for the brain tissue applicable to dynamic rates

Understanding the mechanical response of the brain to external loadings is of critical importance in investigating the pathological conditions of this tissue during injurious conditions. Such injurious loadings may occur at high rates, for example among others, during road traffic or sport accidents, falls, or due to explosions. Hence, investigating the injury mechanism and design of protective devices for the brain requires constitutive modeling of this tissue at such rates. Accordingly, this paper is aimed at critically investigating the physical background for viscohyperelastic modeling of the brain tissue with scrutinizing the elastic fields pertinent to large, time dependent deformations, and developing a fully nonlinear multimode Maxwell model that can mathematically explain such deformations. The proposed model can be calibrated using the simple monotonic uniaxial deformation of the sample extracted from the tissue, and does not require additional information from relaxation or creep experiments. The performance of the proposed model is examined using the experimental results of two different studies, which reveals a desirable agreement. The usefulness, limitations, and future developments of the proposed model are discussed in this paper.     

Generalized optimal spatial filtering using a kernel approach with application to EEG classification

Common spatial patterns (CSP) has been widely used for finding the linear spatial filters which are able to extract the discriminative brain activities between two different mental tasks. However, the CSP is difficult to capture the nonlinearly clustered structure from the non-stationary EEG signals. To relax the presumption of strictly linear patterns in the CSP, in this paper, a generalized CSP (GCSP) based on generalized singular value decomposition (GSVD) and kernel method is proposed. Our method is able to find the nonlinear spatial filters which are formulated in the feature space defined by a nonlinear mapping through kernel functions. Furthermore, in order to overcome the overfitting problem, the regularized GCSP is developed by adding the regularized parameters. The experimental results demonstrate that our method is an effective nonlinear spatial filtering method.     

Live biospeckle laser imaging of root tissues

Live imaging is now a central component for the study of plant developmental processes. Currently, most techniques are extremely constraining: they rely on the marking of specific cellular structures which generally apply to model species because they require genetic transformations. The biospeckle laser (BSL) system was evaluated as an instrument to measure biological activity in plant tissues. The system allows collecting biospeckle patterns from roots which are grown in gels. Laser illumination has been optimized to obtain the images without undesirable specular reflections from the glass tube. Data on two different plant species were obtained and the ability of three different methods to analyze the biospeckle patterns are presented. The results showed that the biospeckle could provide quantitative indicators of the molecular activity from roots which are grown in gel substrate in tissue culture. We also presented a particular experimental configuration and the optimal approach to analyze the images. This may serve as a basis to further works on live BSL in order to study root development.     

Experimental study of the behavioural mechanisms underlying self-organization in human crowds

In animal societies as well as in human crowds, many observed collective behaviours result from self-organized processes based on local interactions among individuals. However, models of crowd dynamics are still lacking a systematic individual-level experimental verification, and the local mechanisms underlying the formation of collective patterns are not yet known in detail. We have conducted a set of well-controlled experiments with pedestrians performing simple avoidance tasks in order to determine the laws ruling their behaviour during interactions. The analysis of the large trajectory dataset was used to compute a behavioural map that describes the average change of the direction and speed of a pedestrian for various interaction distances and angles. The experimental results reveal features of the decision process when pedestrians choose the side on which they evade, and show a side preference that is amplified by mutual interactions. The predictions of a binary interaction model based on the above findings were then compared with bidirectional flows of people recorded in a crowded street. Simulations generate two asymmetric lanes with opposite directions of motion, in quantitative agreement with our empirical observations. The knowledge of pedestrian behavioural laws is an important step ahead in the understanding of the underlying dynamics of crowd behaviour and allows for reliable predictions of collective pedestrian movements under natural conditions.     

Evolving Concepts of Bacterial Species

The same evolutionary forces that cause diversification in sexual eukaryotes are expected to cause diversification in bacteria. However, in bacteria, the wider variety of mechanisms for gene exchange (or lack thereof) increases the range of expected diversity patterns compared to those of sexual organisms. Two parallel concepts for bacterial speciation have developed, based on ecological divergence or barriers to recombination in turn. Recent evidence from DNA sequence data shows that both processes can generate independently evolving groups that are equivalent to sexual species and that represent separate arenas within which recombination (when it occurs), selection and drift occur. It remains unclear, however, how often different processes act in concert to generate simple units of diversity, or whether a more complex model of diversity is required, specifying hierarchical levels at which different cohesive processes operate. We advocate an integrative approach that evaluates the effects of multiple evolutionary forces on diversity patterns. There is also great potential for laboratory studies of bacterial evolution that test evolutionary mechanisms inferred from population genetic analyses of multi-locus and genome sequence data.     

Encoding brain network response to free viewing of videos

A challenging goal for cognitive neuroscience researchers is to determine how mental representations are mapped onto the patterns of neural activity. To address this problem, functional magnetic resonance imaging (fMRI) researchers have developed a large number of encoding and decoding methods. However, previous studies typically used rather limited stimuli representation, like semantic labels and Wavelet Gabor filters, and largely focused on voxel-based brain patterns. Here, we present a new fMRI encoding model to predict the human brain’s responses to free viewing of video clips which aims to deal with this limitation. In this model, we represent the stimuli using a variety of representative visual features in the computer vision community, which can describe the global color distribution, local shape and spatial information and motion information contained in videos, and apply the functional connectivity to model the brain’s activity pattern evoked by these video clips. Our experimental results demonstrate that brain network responses during free viewing of videos can be robustly and accurately predicted across subjects by using visual features. Our study suggests the feasibility of exploring cognitive neuroscience studies by computational image/video analysis and provides a novel concept of using the brain encoding as a test-bed for evaluating visual feature extraction.     

Modular neuromuscular control of human locomotion by central pattern generator

The central pattern generators (CPG) in the spinal cord are thought to be responsible for producing the rhythmic motor patterns during rhythmic activities. For locomotor tasks, this involves much complexity, due to a redundant system of muscle actuators with a large number of highly nonlinear muscles. This study proposes a reduced neural control strategy for the CPG, based on modular organization of the co-active muscles, i.e., muscle synergies. Four synergies were extracted from the EMG data of the major leg muscles of two subjects, during two gait trials each, using non-negative matrix factorization algorithm. A Matsuoka׳s four-neuron CPG model with mutual inhibition, was utilized to generate the rhythmic activation patterns of the muscle synergies, using the hip flexion angle and foot contact force information from the sensory afferents as inputs. The model parameters were tuned using the experimental data of one gait trial, which resulted in a good fitting accuracy (RMSEs between 0.0491 and 0.1399) between the simulation and experimental synergy activations. The model׳s performance was then assessed by comparing its predictions for the activation patterns of the individual leg muscles during locomotion with the relevant EMG data. Results indicated that the characteristic features of the complex activation patterns of the muscles were well reproduced by the model for different gait trials and subjects. In general, the CPG- and muscle synergy-based model was promising in view of its simple architecture, yet extensive potentials for neuromuscular control, e.g., resolving redundancies, distributed and fast control, and modulation of locomotion by simple control signals.     

Making predictions of mangrove deforestation: a comparison of two methods in Kenya

Deforestation of mangroves is of global concern given their importance for carbon storage, biogeochemical cycling and the provision of other ecosystem services, but the links between rates of loss and potential drivers or risk factors are rarely evaluated. Here, we identified key drivers of mangrove loss in Kenya and compared two different approaches to predicting risk. Risk factors tested included various possible predictors of anthropogenic deforestation, related to population, suitability for land use change and accessibility. Two approaches were taken to modelling risk; a quantitative statistical approach and a qualitative categorical ranking approach. A quantitative model linking rates of loss to risk factors was constructed based on generalized least squares regression and using mangrove loss data from 1992 to 2000. Population density, soil type and proximity to roads were the most important predictors. In order to validate this model it was used to generate a map of losses of Kenyan mangroves predicted to have occurred between 2000 and 2010. The qualitative categorical model was constructed using data from the same selection of variables, with the coincidence of different risk factors in particular mangrove areas used in an additive manner to create a relative risk index which was then mapped. Quantitative predictions of loss were significantly correlated with the actual loss of mangroves between 2000 and 2010 and the categorical risk index values were also highly correlated with the quantitative predictions. Hence, in this case the relatively simple categorical modelling approach was of similar predictive value to the more complex quantitative model of mangrove deforestation. The advantages and disadvantages of each approach are discussed, and the implications for mangroves are outlined.     

Estimation of neurophysiological parameters from the waking EEG using a biophysical model of brain dynamics

This paper presents the results from using electroencephalographic (EEG) data to estimate the values of key neurophysiological parameters using a detailed biophysical model of brain activity. The model incorporates spatial and temporal aspects of cortical function including axonal transmission delays, synapto-dendritic rates, range-dependent connectivities, excitatory and inhibitory neural populations, and intrathalamic, intracortical, corticocortical and corticothalamic pathways. Parameter estimates were obtained by fitting the model's theoretical spectrum to EEG spectra from each of 100 healthy human subjects. Statistical analysis was used to infer significant parameter variations occurring between eyes-closed and eyes-open states, and a correlation matrix was used to investigate links between the parameter variations and traditional measures of quantitative EEG (qEEG). Accurate fits to all experimental spectra were observed, and both inter-subject and between-state variability were accounted for by the variance in the fitted biophysical parameters, which were in turn consistent with known independent experimental and theoretical estimates. These values thus provide physiological information regarding the state. transitions (eyes-closed vs. eyes-open) and phenomena including cortical idling and alpha desynchronization. The parameters are also consistent with traditional qEEG, but are more informative, since they provide links to underlying physiological processes. To our knowledge, this is the first study where a detailed biophysical model of the brain is used to estimate neurophysiological parameters underlying the transitions in a broad range (0.25-50 Hz) of EEG spectra obtained from a large set of human data.     

Coalescent processes when the distribution of offspring number among individuals is highly skewed

We report a complex set of scaling relationships between mutation and reproduction in a simple model of a population. These follow from a consideration of patterns of genetic diversity in a sample of DNA sequences. Five different possible limit processes, each with a different scaled mutation parameter, can be used to describe genetic diversity in a large population. Only one of these corresponds to the usual population genetic model, and the others make drastically different predictions about genetic diversity. The complexity arises because individuals can potentially have very many offspring. To the extent that this occurs in a given species, our results imply that inferences from genetic data made under the usual assumptions are likely to be wrong. Our results also uncover a fundamental difference between populations in which generations are overlapping and those in which generations are discrete. We choose one of the five limit processes that appears to be appropriate for some marine organisms and use a sample of genetic data from a population of Pacific oysters to infer the parameters of the model. The data suggest the presence of rare reproduction events in which approximately 8% of the population is replaced by the offspring of a single individual.     

Finite Element Model to Reproduce the Effect of Pre-Stress and Needle Insertion Velocity During Infusions into Brain Phantom Gel

Convection-enhanced delivery (CED) is a technique to bypass the blood-brain barrier and deliver therapeutic agents into the brain. However, animal studies and preliminary clinical trials have reported reduced efficacy to transport drugs in specific regions, attributed mainly to backflow, in which an annular zone is formed outside the catheter and the fluid preferentially flows toward the surface of the brain rather than through the tissue toward the targeted area. In this study, a finite element model of backflow was updated by implementing the pre-stress generated during needle insertion, which allows considering the effect of needle insertion velocity during CED infusions in agarose gel. The nonlinear mechanical properties of the agarose solutions were obtained by fitting experimental data from stress-relaxation tests. Additional experimental measurements of backflow lengths were used to adjust the pre-stress model. The developed model was able to reproduce changes of backflow length under different insertions velocities and flow rates. These findings reveal the relevance of considering the pre-stress in the tissue located around the needle surface during CED infusions into the brain.