Using sequences of life-events to predict human lives

Here we represent human lives in a way that shares structural similarity to language, and we exploit this similarity to adapt natural language processing techniques to examine the evolution and predictability of human lives based on detailed event sequences. We do this by drawing on a comprehensive registry dataset, which is available for Denmark across several years, and that includes information about life-events related to health, education, occupation, income, address and working hours, recorded with day-to-day resolution. We create embeddings of life-events in a single vector space, showing that this embedding space is robust and highly structured. Our models allow us to predict diverse outcomes ranging from early mortality to personality nuances, outperforming state-of-the-art models by a wide margin. Using methods for interpreting deep learning models, we probe the algorithm to understand the factors that enable our predictions. Our framework allows researchers to discover potential mechanisms that impact life outcomes as well as the associated possibilities for personalized interventions.

Uncovering Cortical Units of Processing From Multi-Layered Connectomes

Modern diffusion and functional magnetic resonance imaging (dMRI/fMRI) provide non-invasive high-resolution images from which multi-layered networks of whole-brain structural and functional connectivity can be derived. Unfortunately, the lack of observed correspondence between the connectivity profiles of the two modalities challenges the understanding of the relationship between the functional and structural connectome. Rather than focusing on correspondence at the level of connections we presently investigate correspondence in terms of modular organization according to shared canonical processing units. We use a stochastic block-model (SBM) as a data-driven approach for clustering high-resolution multi-layer whole-brain connectivity networks and use prediction to quantify the extent to which a given clustering accounts for the connectome within a modality. The employed SBM assumes a single underlying parcellation exists across modalities whilst permitting each modality to possess an independent connectivity structure between parcels thereby imposing concurrent functional and structural units but different structural and functional connectivity profiles. We contrast the joint processing units to their modality specific counterparts and find that even though data-driven structural and functional parcellations exhibit substantial differences, attributed to modality specific biases, the joint model is able to achieve a consensus representation that well accounts for both the functional and structural connectome providing improved representations of functional connectivity compared to using functional data alone. This implies that a representation persists in the consensus model that is shared by the individual modalities. We find additional support for this viewpoint when the anatomical correspondence between modalities is removed from the joint modeling. The resultant drop in predictive performance is in general substantial, confirming that the anatomical correspondence of processing units is indeed present between the two modalities. Our findings illustrate how multi-modal integration admits consensus representations well-characterizing each individual modality despite their biases and points to the importance of multi-layered connectomes as providing supplementary information regarding the brain's canonical processing units.

Noise-assisted variational quantum thermalization

Preparing thermal states on a quantum computer can have a variety of applications, from simulating many-body quantum systems to training machine learning models. Variational circuits have been proposed for this task on near-term quantum computers, but several challenges remain, such as finding a scalable cost-function, avoiding the need of purification, and mitigating noise effects. We propose a new algorithm for thermal state preparation that tackles those three challenges by exploiting the noise of quantum circuits. We consider a variational architecture containing a depolarizing channel after each unitary layer, with the ability to directly control the level of noise. We derive a closed-form approximation for the free-energy of such circuit and use it as a cost function for our variational algorithm. By evaluating our method on a variety of Hamiltonians and system sizes, we find several systems for which the thermal state can be approximated with a high fidelity. However, we also show that the ability for our algorithm to learn the thermal state strongly depends on the temperature: while a high fidelity can be obtained for high and low temperatures, we identify a specific range for which the problem becomes more challenging. We hope that this first study on noise-assisted thermal state preparation will inspire future research on exploiting noise in variational algorithms.

Real-Time Decoding of Attentional States Using Closed-Loop EEG Neurofeedback

Sustained attention is a cognitive ability to maintain task focus over extended periods of time (Mackworth, 1948; Chun, Golomb, & Turk-Browne, 2011). In this study, scalp electroencephalography (EEG) signals were processed in real time using a 32 dry-electrode system during a sustained visual attention task. An attention training paradigm was implemented, as designed in DeBettencourt, Cohen, Lee, Norman, and Turk-Browne (2015) in which the composition of a sequence of blended images is updated based on the participant's decoded attentional level to a primed image category. It was hypothesized that a single neurofeedback training session would improve sustained attention abilities. Twenty-two participants were trained on a single neurofeedback session with behavioral pretraining and posttraining sessions within three consecutive days. Half of the participants functioned as controls in a double-blinded design and received sham neurofeedback. During the neurofeedback session, attentional states to primed categories were decoded in real time and used to provide a continuous feedback signal customized to each participant in a closed-loop approach. We report a mean classifier decoding error rate of 34.3% (chance = 50%). Within the neurofeedback group, there was a greater level of task-relevant attentional information decoded in the participant's brain before making a correct behavioral response than before an incorrect response. This effect was not visible in the control group (interaction p=7.23e-4), which strongly indicates that we were able to achieve a meaningful measure of subjective attentional state in real time and control participants' behavior during the neurofeedback session. We do not provide conclusive evidence whether the single neurofeedback session per se provided lasting effects in sustained attention abilities. We developed a portable EEG neurofeedback system capable of decoding attentional states and predicting behavioral choices in the attention task at hand. The neurofeedback code framework is Python based and open source, and it allows users to actively engage in the development of neurofeedback tools for scientific and translational use.

Single-Trial Decoding of Scalp EEG under Natural Conditions

There is significant current interest in decoding mental states from electroencephalography (EEG) recordings. EEG signals are subject-specific, are sensitive to disturbances, and have a low signal-to-noise ratio, which has been mitigated by the use of laboratory-grade EEG acquisition equipment under highly controlled conditions. In the present study, we investigate single-trial decoding of natural, complex stimuli based on scalp EEG acquired with a portable, 32 dry-electrode sensor system in a typical office setting. We probe generalizability by a leave-one-subject-out cross-validation approach. We demonstrate that support vector machine (SVM) classifiers trained on a relatively small set of denoised (averaged) pseudotrials perform on par with classifiers trained on a large set of noisy single-trial samples. We propose a novel method for computing sensitivity maps of EEG-based SVM classifiers for visualization of EEG signatures exploited by the SVM classifiers. Moreover, we apply an NPAIRS resampling framework for estimation of map uncertainty, and thus show that effect sizes of sensitivity maps for classifiers trained on small samples of denoised data and large samples of noisy data are similar. Finally, we demonstrate that the average pseudotrial classifier can successfully predict the class of single trials from withheld subjects, which allows for fast classifier training, parameter optimization, and unbiased performance evaluation in machine learning approaches for brain decoding.

Sequence Classification Using Third-Order Moments

Model-based classification of sequence data using a set of hidden Markov models is a well-known technique. The involved score function, which is often based on the class-conditional likelihood, can, however, be computationally demanding, especially for long data sequences. Inspired by recent theoretical advances in spectral learning of hidden Markov models, we propose a score function based on third-order moments. In particular, we propose to use the Kullback-Leibler divergence between theoretical and empirical third-order moments for classification of sequence data with discrete observations. The proposed method provides lower computational complexity at classification time than the usual likelihood-based methods. In order to demonstrate the properties of the proposed method, we perform classification of both simulated data and empirical data from a human activity recognition study.

Spatio-temporal reconstruction of brain dynamics from EEG with a Markov prior

Electroencephalography (EEG) can capture brain dynamics in high temporal resolution. By projecting the scalp EEG signal back to its origin in the brain also high spatial resolution can be achieved. Source localized EEG therefore has potential to be a very powerful tool for understanding the functional dynamics of the brain. Solving the inverse problem of EEG is however highly ill-posed as there are many more potential locations of the EEG generators than EEG measurement points. Several well-known properties of brain dynamics can be exploited to alleviate this problem. More short ranging connections exist in the brain than long ranging, arguing for spatially focal sources. Additionally, recent work (Delorme et al., 2012) argues that EEG can be decomposed into components having sparse source distributions. On the temporal side both short and long term stationarity of brain activation are seen. We summarize these insights in an inverse solver, the so-called "Variational Garrote" (Kappen and Gómez, 2013). Using a Markov prior we can incorporate flexible degrees of temporal stationarity. Through spatial basis functions spatially smooth distributions are obtained. Sparsity of these are inherent to the Variational Garrote solver. We name our method the MarkoVG and demonstrate its ability to adapt to the temporal smoothness and spatial sparsity in simulated EEG data. Finally a benchmark EEG dataset is used to demonstrate MarkoVG's ability to recover non-stationary brain dynamics.

Acromelic frontonasal dysostosis and ZSWIM6 mutation: phenotypic spectrum and mosaicism

Acromelic frontonasal dysostosis (AFND) is a distinctive and rare frontonasal malformation that presents in combination with brain and limb abnormalities. A single recurrent heterozygous missense substitution in ZSWIM6, encoding a protein of unknown function, was previously shown to underlie this disorder in four unrelated cases. Here we describe four additional individuals from three families, comprising two sporadic subjects (one of whom had no limb malformation) and a mildly affected female with a severely affected son. In the latter family we demonstrate parental mosaicism through deep sequencing of DNA isolated from a variety of tissues, which each contain different levels of mutation. This has important implications for genetic counselling.

Multiview Bayesian Correlated Component Analysis

Correlated component analysis as proposed by Dmochowski, Sajda, Dias, and Parra (2012) is a tool for investigating brain process similarity in the responses to multiple views of a given stimulus. Correlated components are identified under the assumption that the involved spatial networks are identical. Here we propose a hierarchical probabilistic model that can infer the level of universality in such multiview data, from completely unrelated representations, corresponding to canonical correlation analysis, to identical representations as in correlated component analysis. This new model, which we denote Bayesian correlated component analysis, evaluates favorably against three relevant algorithms in simulated data. A well-established benchmark EEG data set is used to further validate the new model and infer the variability of spatial representations across multiple subjects.

The Center for Integrated Molecular Brain Imaging (Cimbi) database

We here describe a multimodality neuroimaging containing data from healthy volunteers and patients, acquired within the Lundbeck Foundation Center for Integrated Molecular Brain Imaging (Cimbi) in Copenhagen, Denmark. The data is of particular relevance for neurobiological research questions related to the serotonergic transmitter system with its normative data on the serotonergic subtype receptors 5-HT1A, 5-HT1B, 5-HT2A, and 5-HT4 and the 5-HT transporter (5-HTT), but can easily serve other purposes. The Cimbi database and Cimbi biobank were formally established in 2008 with the purpose to store the wealth of Cimbi-acquired data in a highly structured and standardized manner in accordance with the regulations issued by the Danish Data Protection Agency as well as to provide a quality-controlled resource for future hypothesis-generating and hypothesis-driven studies. The Cimbi database currently comprises a total of 1100 PET and 1000 structural and functional MRI scans and it holds a multitude of additional data, such as genetic and biochemical data, and scores from 17 self-reported questionnaires and from 11 neuropsychological paper/computer tests. The database associated Cimbi biobank currently contains blood and in some instances saliva samples from about 500 healthy volunteers and 300 patients with e.g., major depression, dementia, substance abuse, obesity, and impulsive aggression. Data continue to be added to the Cimbi database and biobank.

How many separable sources? Model selection in independent components analysis

Unlike mixtures consisting solely of non-Gaussian sources, mixtures including two or more Gaussian components cannot be separated using standard independent components analysis methods that are based on higher order statistics and independent observations. The mixed Independent Components Analysis/Principal Components Analysis (mixed ICA/PCA) model described here accommodates one or more Gaussian components in the independent components analysis model and uses principal components analysis to characterize contributions from this inseparable Gaussian subspace. Information theory can then be used to select from among potential model categories with differing numbers of Gaussian components. Based on simulation studies, the assumptions and approximations underlying the Akaike Information Criterion do not hold in this setting, even with a very large number of observations. Cross-validation is a suitable, though computationally intensive alternative for model selection. Application of the algorithm is illustrated using Fisher's iris data set and Howells' craniometric data set. Mixed ICA/PCA is of potential interest in any field of scientific investigation where the authenticity of blindly separated non-Gaussian sources might otherwise be questionable. Failure of the Akaike Information Criterion in model selection also has relevance in traditional independent components analysis where all sources are assumed non-Gaussian.

Non-parametric Bayesian graph models reveal community structure in resting state fMRI

Modeling of resting state functional magnetic resonance imaging (rs-fMRI) data using network models is of increasing interest. It is often desirable to group nodes into clusters to interpret the communication patterns between nodes. In this study we consider three different nonparametric Bayesian models for node clustering in complex networks. In particular, we test their ability to predict unseen data and their ability to reproduce clustering across datasets. The three generative models considered are the Infinite Relational Model (IRM), Bayesian Community Detection (BCD), and the Infinite Diagonal Model (IDM). The models define probabilities of generating links within and between clusters and the difference between the models lies in the restrictions they impose upon the between-cluster link probabilities. IRM is the most flexible model with no restrictions on the probabilities of links between clusters. BCD restricts the between-cluster link probabilities to be strictly lower than within-cluster link probabilities to conform to the community structure typically seen in social networks. IDM only models a single between-cluster link probability, which can be interpreted as a background noise probability. These probabilistic models are compared against three other approaches for node clustering, namely Infomap, Louvain modularity, and hierarchical clustering. Using 3 different datasets comprising healthy volunteers' rs-fMRI we found that the BCD model was in general the most predictive and reproducible model. This suggests that rs-fMRI data exhibits community structure and furthermore points to the significance of modeling heterogeneous between-cluster link probabilities.

Chitin-calcium alginate composite fibers for wound care dressings spun from ionic liquid solution

Chitin-calcium alginate composite fibers were prepared from a solution of high molecular weight chitin extracted from shrimp shells and alginic acid in the ionic liquid 1-ethyl-3-methylimidazolium acetate by dry-jet wet spinning into an aqueous bath saturated with CaCO₃. The fibers exhibited a significant proportion of the individual properties of both calcium alginate and chitin. Ultimate stress values were close to values obtained for calcium alginate fibers, and the absorption capacities measured were consistent with those reported for current wound care dressings. Wound healing studies (rat model, histological evaluation) indicated that chitin-calcium alginate covered wound sites underwent normal wound healing with re-epithelialization and that coverage of the dermal fibrosis with hyperplastic epidermis was consistently complete after only 7 days of treatment. Using a single patch per wound per animal during the entire study, all rat wounds achieved 95-99% closure by day 10 with complete wound closure by day 14.

Frontal alpha oscillations distinguish leaders from followers: multivariate decoding of mutually interacting brains

Successful social interactions rely upon the abilities of two or more people to mutually exchange information in real-time, while simultaneously adapting to one another. The neural basis of social cognition has mostly been investigated in isolated individuals, and more recently using two-person paradigms to quantify the neuronal dynamics underlying social interaction. While several studies have shown the relevance of understanding complementary and mutually adaptive processes, the neural mechanisms underlying such coordinative behavioral patterns during joint action remain largely unknown. Here, we employed a synchronized finger-tapping task while measuring dual-EEG from pairs of human participants who either mutually adjusted to each other in an interactive task or followed a computer metronome. Neurophysiologically, the interactive condition was characterized by a stronger suppression of alpha and low-beta oscillations over motor and frontal areas in contrast to the non-interactive computer condition. A multivariate analysis of two-brain activity to classify interactive versus non-interactive trials revealed asymmetric patterns of the frontal alpha-suppression in each pair, during both task anticipation and execution, such that only one member showed the frontal component. Analysis of the behavioral data showed that this distinction coincided with the leader-follower relationship in 8/9 pairs, with the leaders characterized by the stronger frontal alpha-suppression. This suggests that leaders invest more resources in prospective planning and control. Hence our results show that the spontaneous emergence of leader-follower relationships in dyadic interactions can be predicted from EEG recordings of brain activity prior to and during interaction. Furthermore, this emphasizes the importance of investigating complementarity in joint action.

The smartphone brain scanner: a portable real-time neuroimaging system

Combining low-cost wireless EEG sensors with smartphones offers novel opportunities for mobile brain imaging in an everyday context. Here we present the technical details and validation of a framework for building multi-platform, portable EEG applications with real-time 3D source reconstruction. The system – Smartphone Brain Scanner – combines an off-the-shelf neuroheadset or EEG cap with a smartphone or tablet, and as such represents the first fully portable system for real-time 3D EEG imaging. We discuss the benefits and challenges, including technical limitations as well as details of real-time reconstruction of 3D images of brain activity. We present examples of brain activity captured in a simple experiment involving imagined finger tapping, which shows that the acquired signal in a relevant brain region is similar to that obtained with standard EEG lab equipment. Although the quality of the signal in a mobile solution using an off-the-shelf consumer neuroheadset is lower than the signal obtained using high-density standard EEG equipment, we propose mobile application development may offset the disadvantages and provide completely new opportunities for neuroimaging in natural settings.

Smartphones as pocketable labs: visions for mobile brain imaging and neurofeedback

Mobile brain imaging solutions, such as the Smartphone Brain Scanner, which combines low cost wireless EEG sensors with open source software for real-time neuroimaging, may transform neuroscience experimental paradigms. Normally subject to the physical constraints in labs, neuroscience experimental paradigms can be transformed into dynamic environments allowing for the capturing of brain signals in everyday contexts. Using smartphones or tablets to access text or images may enable experimental design capable of tracing emotional responses when shopping or consuming media, incorporating sensorimotor responses reflecting our actions into brain machine interfaces, and facilitating neurofeedback training over extended periods. Even though the quality of consumer neuroheadsets is still lower than laboratory equipment and susceptible to environmental noise, we show that mobile neuroimaging solutions, like the Smartphone Brain Scanner, complemented by 3D reconstruction or source separation techniques may support a range of neuroimaging applications and thus become a valuable addition to high-end neuroimaging solutions.

An evaluation of EEG scanner's dependence on the imaging technique, forward model computation method, and array dimensionality

EEG source reconstruction involves solving an inverse problem that is highly ill-posed and dependent on a generally fixed forward propagation model. In this contribution we compare a low and high density EEG setup's dependence on correct forward modeling. Specifically, we examine how different forward models affect the source estimates obtained using four inverse solvers Minimum-Norm, LORETA, Minimum-Variance Adaptive Beamformer, and Sparse Bayesian Learning.

Specialized tools are needed when searching the web for rare disease diagnoses

In our recent paper, we study web search as an aid in the process of diagnosing rare diseases. To answer the question of how well Google Search and PubMed perform, we created an evaluation framework with 56 diagnostic cases and made our own specialized search engine, FindZebra (findzebra.com). FindZebra uses a set of publicly available curated sources on rare diseases and an open-source information retrieval system, Indri. Our evaluation and the feedback received after the publication of our paper both show that FindZebra outperforms Google Search and PubMed. In this paper, we summarize the original findings and the response to FindZebra, discuss why Google Search is not designed for specialized tasks and outline some of the current trends in using web resources and social media for medical diagnosis.

Nonlinear denoising and analysis of neuroimages with kernel principal component analysis and pre-image estimation

We investigate the use of kernel principal component analysis (PCA) and the inverse problem known as pre-image estimation in neuroimaging: i) We explore kernel PCA and pre-image estimation as a means for image denoising as part of the image preprocessing pipeline. Evaluation of the denoising procedure is performed within a data-driven split-half evaluation framework. ii) We introduce manifold navigation for exploration of a nonlinear data manifold, and illustrate how pre-image estimation can be used to generate brain maps in the continuum between experimentally defined brain states/classes. We base these illustrations on two fMRI BOLD data sets - one from a simple finger tapping experiment and the other from an experiment on object recognition in the ventral temporal lobe.

Model selection for Gaussian kernel PCA denoising

We propose kernel parallel analysis (kPA) for automatic kernel scale and model order selection in Gaussian kernel principal component analysis (KPCA). Parallel analysis is based on a permutation test for covariance and has previously been applied for model order selection in linear PCA, we here augment the procedure to also tune the Gaussian kernel scale of radial basis function based KPCA. We evaluate kPA for denoising of simulated data and the U.S. postal data set of handwritten digits. We find that kPA outperforms other heuristics to choose the model order and kernel scale in terms of signal-to-noise ratio of the denoised data.

Visualization of nonlinear kernel models in neuroimaging by sensitivity maps

There is significant current interest in decoding mental states from neuroimages. In this context kernel methods, e.g., support vector machines (SVM) are frequently adopted to learn statistical relations between patterns of brain activation and experimental conditions. In this paper we focus on visualization of such nonlinear kernel models. Specifically, we investigate the sensitivity map as a technique for generation of global summary maps of kernel classification models. We illustrate the performance of the sensitivity map on functional magnetic resonance (fMRI) data based on visual stimuli. We show that the performance of linear models is reduced for certain scan labelings/categorizations in this data set, while the nonlinear models provide more flexibility. We show that the sensitivity map can be used to visualize nonlinear versions of kernel logistic regression, the kernel Fisher discriminant, and the SVM, and conclude that the sensitivity map is a versatile and computationally efficient tool for visualization of nonlinear kernel models in neuroimaging.

Algorithms for sparse nonnegative Tucker decompositions

There is a increasing interest in analysis of large-scale multiway data. The concept of multiway data refers to arrays of data with more than two dimensions, that is, taking the form of tensors. To analyze such data, decomposition techniques are widely used. The two most common decompositions for tensors are the Tucker model and the more restricted PARAFAC model. Both models can be viewed as generalizations of the regular factor analysis to data of more than two modalities. Nonnegative matrix factorization (NMF), in conjunction with sparse coding, has recently been given much attention due to its part-based and easy interpretable representation. While NMF has been extended to the PARAFAC model, no such attempt has been done to extend NMF to the Tucker model. However, if the tensor data analyzed are nonnegative, it may well be relevant to consider purely additive (i.e., nonnegative) Tucker decompositions). To reduce ambiguities of this type of decomposition, we develop updates that can impose sparseness in any combination of modalities, hence, proposed algorithms for sparse nonnegative Tucker decompositions (SN-TUCKER). We demonstrate how the proposed algorithms are superior to existing algorithms for Tucker decompositions when the data and interactions can be considered nonnegative. We further illustrate how sparse coding can help identify what model (PARAFAC or Tucker) is more appropriate for the data as well as to select the number of components by turning off excess components. The algorithms for SN-TUCKER can be downloaded from Mørup (2007).

Biclique communities

We present a method for detecting communities in bipartite networks. Based on an extension of the k -clique community detection algorithm, we demonstrate how modular structure in bipartite networks presents itself as overlapping bicliques. If bipartite information is available, the biclique community detection algorithm retains all of the advantages of the k -clique algorithm, but avoids discarding important structural information when performing a one-mode projection of the network. Further, the biclique community detection algorithm provides a level of flexibility by incorporating independent clique thresholds for each of the nonoverlapping node sets in the bipartite network.

Akathisia--seen from the patients' perspective

Akathisia is a side-effect that continues to affect a large percentage of psychiatric patients though the focus over recent years has moved away from extrapyramidal side-effects to metabolic problems. As cognitive behavioural therapy has become an integrated part of the management plan for psychotic patients an increased understanding of patients' interpretation of their symptoms and side-effects are encouraged. This case series illustrates different views that patients take on medication induced side-effects.

Bayesian model comparison in nonlinear BOLD fMRI hemodynamics

Nonlinear hemodynamic models express the BOLD (blood oxygenation level dependent) signal as a nonlinear, parametric functional of the temporal sequence of local neural activity. Several models have been proposed for both the neural activity and the hemodynamics. We compare two such combined models: the original balloon model with a square-pulse neural model (Friston, Mechelli, Turner, & Price, 2000) and an extended balloon model with a more sophisticated neural model (Buxton, Uludag, Dubowitz, & Liu, 2004). We learn the parameters of both models using a Bayesian approach, where the distribution of the parameters conditioned on the data is estimated using Markov chain Monte Carlo techniques. Using a split-half resampling procedure (Strother, Anderson, & Hansen, 2002), we compare the generalization abilities of the models as well as their reproducibility, for both synthetic and real data, recorded from two different visual stimulation paradigms. The results show that the simple model is the better one for these data.

Segmentation of age-related white matter changes in a clinical multi-center study

Age-related white matter changes (WMC) are thought to be a marker of vascular pathology, and have been associated with motor and cognitive deficits. In the present study, an optimized artificial neural network was used as an automatic segmentation method to produce probabilistic maps of WMC in a clinical multi-center study. The neural network uses information from T1- and T2-weighted and fluid attenuation inversion recovery (FLAIR) magnetic resonance (MR) scans, neighboring voxels and spatial location. Generalizability of the neural network was optimized by including the Optimal Brain Damage (OBD) pruning method in the training stage. Six optimized neural networks were produced to investigate the impact of different input information on WMC segmentation. The automatic segmentation method was applied to MR scans of 362 non-demented elderly subjects from 11 centers in the European multi-center study Leukoaraiosis And Disability (LADIS). Semi-manually delineated WMC were used for validating the segmentation produced by the neural networks. The neural network segmentation demonstrated high consistency between subjects and centers, making it a promising technique for large studies. For WMC volumes less than 10 ml, an increasing discrepancy between semi-manual and neural network segmentation was observed using the similarity index (SI) measure. The use of all three image modalities significantly improved cross-center generalizability compared to neural networks using the FLAIR image only. Expert knowledge not available to the neural networks was a minor source of discrepancy, while variation in MR scan quality constituted the largest source of error.

Multivariate strategies in functional magnetic resonance imaging

We discuss aspects of multivariate fMRI modeling, including the statistical evaluation of multivariate models and means for dimensional reduction. In a case study we analyze linear and non-linear dimensional reduction tools in the context of a 'mind reading' predictive multivariate fMRI model.

Effect of spatial alignment transformations in PCA and ICA of functional neuroimages

It has been previously observed that independent component analysis (ICA), if applied to data pooled in a particular way, may lessen the need for spatial alignment of scans in a functional neuroimaging study. In this paper, we seek to determine analytically the conditions under which this observation is true, not only for spatial ICA, but also for temporal ICA and for principal component analysis (PCA). In each case, we find conditions that the spatial alignment operator must satisfy to ensure invariance of the results. We illustrate our findings using functional magnetic-resonance imaging (fMRI) data. Our analysis is applicable to both intersubject and intrasubject spatial normalization.

ERPWAVELAB

The open source toolbox 'ERPWAVELAB' is developed for multi-channel time-frequency analysis of event related activity of EEG and MEG data. The toolbox provides tools for data analysis and visualization of the most commonly used measures of time-frequency transformed event related data as well as data decomposition through non-negative matrix and multi-way (tensor) factorization. The decompositions provided can accommodate additional dimensions like subjects, conditions or repeats and as such they are perfected for group analysis. Furthermore, the toolbox enables tracking of phase locked activity from one channel-time-frequency instance to another as well as tools for artifact rejection in the time-frequency domain. Several other features are highlighted. ERPWAVELAB can freely be downloaded from www.erpwavelab.org, requires EEGLAB [Delorme A, Makeig S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Meth 2004;134:9-21] and runs under MATLAB (The Mathworks, Inc.).

Does a positive end-expiratory pressure-induced reduction in stroke volume indicate preload responsiveness? An experimental study

BACKGROUND

Increases in positive end-expiratory pressure (PEEP) are often associated with cardiovascular depression, responding to fluid loading. Therefore, we hypothesized that if stroke volume (SV) is reduced by an increase in PEEP this reduction is an indicator of hypovolemia or preload responsiveness, i.e. that SV would increase by fluid administration at zero end-expiratory pressure (ZEEP). The relationship between the cardiovascular response to different PEEP levels and fluid load as well as the relation between change in SV as a result of change in preload (Frank-Starling relationship) were evaluated in a porcine model. In addition, other measures of fluid status were assessed.

METHODS

Eight, 20-22 kg, anesthetized, mechanically ventilated pigs were subjected to 0, 10, and 20 cm H(2)O PEEP at 10% (of estimated blood volume) hypovolemia, normo- and 10% hypervolemia, and to ZEEP at 20% hypervolemia. SV, cardiac output, intrathoracic blood volume and airway, esophageal, vascular pressures, stroke volume variations, left ventricular end-diastolic and end-systolic areas and respiratory variations in the diameter of the inferior vena cava were obtained.

RESULTS

At hypovolemia and normovolemia, 10 cm H(2)O PEEP induced a significant decrease in SV, while no change occurred at 10% hypervolemia. SV measured at ZEEP increased from hypovolemia to normovolemia and 10% hypervolemia, while no change was found between 10% and 20% hypervolemia. The sensitivity and specificity decrease in SV by PEEP indicating an increase in SV by fluids was 60-88% and 67%, respectively, depending on the volemic (preload) levels.

CONCLUSION

Although the overall results suggest that a change in SV by PEEP might predict preload responsiveness, the individual response of SV by 10 cm H(2)O PEEP and of the successive fluid administration seemed to be dependent on where on the Frank-Starling curve the heart function was located.

Model Selection for Convolutive ICA with an Application to Spatiotemporal Analysis of EEG

We present a new algorithm for maximum likelihood convolutive independent component analysis (ICA) in which components are unmixed using stable autoregressive filters determined implicitly by estimating a convolutive model of the mixing process. By introducing a convolutive mixing model for the components, we show how the order of the filters in the model can be correctly detected using Bayesian model selection. We demonstrate a framework for deconvolving a subspace of independent components in electroencephalography (EEG). Initial results suggest that in some cases, convolutive mixing may be a more realistic model for EEG signals than the instantaneous ICA model.

Proprioceptive evoked gamma oscillations

A proprioceptive stimulus consisting of a weight change of a handheld load has recently been shown to elicit an evoked potential. Previously, somatosensory gamma oscillations have only been evoked by electrical stimuli. We conjectured that a natural proprioceptive stimulus also would be able to evoke gamma oscillations. EEG was recorded using 64 channels in 14 healthy subjects. In each of three runs a stimulus of 100 g load increment in each hand was presented in 120 trials. Data were wavelet transformed and runs collapsed. Inter-trial phase coherence (ITPC) was computed as the best measure of precision of evoked oscillations. The window of interest was 25-75 Hz and 0-200 ms. ITPC maxima were determined by visual inspection and these results were compared to results of multi-way factor analyses. The predicted 40 Hz activity was observed 60 ms after stimulus onset in the parietal region contralateral to stimulus side and additionally an unexpected 20 Hz activity was observed slightly lateralized in the frontal central region. The gamma phase locking may be a manifestation of early somatosensory feature integration. The analyses suggest that the high frequency activity consists of two distinct underlying oscillatory processes, in need of further investigation.

Growth-rate regulated genes have profound impact on interpretation of transcriptome profiling in Saccharomyces cerevisiae

BACKGROUND

Growth rate is central to the development of cells in all organisms. However, little is known about the impact of changing growth rates. We used continuous cultures to control growth rate and studied the transcriptional program of the model eukaryote Saccharomyces cerevisiae, with generation times varying between 2 and 35 hours.

RESULTS

A total of 5930 transcripts were identified at the different growth rates studied. Consensus clustering of these revealed that half of all yeast genes are affected by the specific growth rate, and that the changes are similar to those found when cells are exposed to different types of stress (>80% overlap). Genes with decreased transcript levels in response to faster growth are largely of unknown function (>50%) whereas genes with increased transcript levels are involved in macromolecular biosynthesis such as those that encode ribosomal proteins. This group also covers most targets of the transcriptional activator RAP1, which is also known to be involved in replication. A positive correlation between the location of replication origins and the location of growth-regulated genes suggests a role for replication in growth rate regulation.

CONCLUSION

Our data show that the cellular growth rate has great influence on transcriptional regulation. This, in turn, implies that one should be cautious when comparing mutants with different growth rates. Our findings also indicate that much of the regulation is coordinated via the chromosomal location of the affected genes, which may be valuable information for the control of heterologous gene expression in metabolic engineering.