Chow–Liu trees are sufficient predictive models for reproducing key features of functional networks of periictal EEG time-series

Seizure freedom in patients suffering from pharmacoresistant epilepsies is still not achieved in 20–30% of all cases. Hence, current therapies need to be improved, based on a more complete understanding of ictogenesis. In this respect, the analysis of functional networks derived from intracranial electroencephalographic (iEEG) data has recently become a standard tool. Functional networks however are purely descriptive models and thus are conceptually unable to predict fundamental features of iEEG time-series, e.g., in the context of therapeutical brain stimulation. In this paper we present some first steps towards overcoming the limitations of functional network analysis, by showing that its results are implied by a simple predictive model of time-sliced iEEG time-series. More specifically, we learn distinct graphical models (so called Chow–Liu (CL) trees) as models for the spatial dependencies between iEEG signals. Bayesian inference is then applied to the CL trees, allowing for an analytic derivation/prediction of functional networks, based on thresholding of the absolute value Pearson correlation coefficient (CC) matrix. Using various measures, the thus obtained networks are then compared to those which were derived in the classical way from the empirical CC-matrix. In the high threshold limit we find (a) an excellent agreement between the two networks and (b) key features of periictal networks as they have previously been reported in the literature. Apart from functional networks, both matrices are also compared element-wise, showing that the CL approach leads to a sparse representation, by setting small correlations to values close to zero while preserving the larger ones. Overall, this paper shows the validity of CL-trees as simple, spatially predictive models for periictal iEEG data. Moreover, we suggest straightforward generalizations of the CL-approach for modeling also the temporal features of iEEG signals.

Identification of a new genomic hot spot of evolutionary diversification of protein function.

Establishment of phylogenetic relationships remains a challenging task because it is based on computational analysis of genomic hot spots that display species-specific sequence variations. Here, we identify a species-specific thymine-to-guanine sequence variation in the Glrb gene which gives rise to species-specific splice donor sites in the Glrb genes of mouse and bushbaby. The resulting splice insert in the receptor for the inhibitory neurotransmitter glycine (GlyR) conveys synaptic receptor clustering and specific association with a particular synaptic plasticity-related splice variant of the postsynaptic scaffold protein gephyrin. This study identifies a new genomic hot spot which contributes to phylogenetic diversification of protein function and advances our understanding of phylogenetic relationships.

Inference of evolutionary forces acting on human biological pathways

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald–Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

Mechanisms of intrinsic beating variability in cardiac cell cultures and model pacemaker networks

Heart rate variability (HRV) exhibits fluctuations characterized by a power law behavior of its power spectrum. The interpretation of this nonlinear HRV behavior, resulting from interactions between extracardiac regulatory mechanisms, could be clinically useful. However, the involvement of intrinsic variations of pacemaker rate in HRV has scarcely been investigated. We examined beating variability in spontaneously active incubating cultures of neonatal rat ventricular myocytes using microelectrode arrays. In networks of mathematical model pacemaker cells, we evaluated the variability induced by the stochastic gating of transmembrane currents and of calcium release channels and by the dynamic turnover of ion channels. In the cultures, spontaneous activity originated from a mobile focus. Both the beat-to-beat movement of the focus and beat rate variability exhibited a power law behavior. In the model networks, stochastic fluctuations in transmembrane currents and stochastic gating of calcium release channels did not reproduce the spatiotemporal patterns observed in vitro. In contrast, long-term correlations produced by the turnover of ion channels induced variability patterns with a power law behavior similar to those observed experimentally. Therefore, phenomena leading to long-term correlated variations in pacemaker cellular function may, in conjunction with extracardiac regulatory mechanisms, contribute to the nonlinear characteristics of HRV.

Cortical reorganisation of cerebral networks after childhood stroke: impact on outcome.

BACKGROUND Recovery after arterial ischaemic stroke is known to largely depend on the plastic properties of the brain. The present study examines changes in the network topography of the developing brain after stroke. Effects of brain damage are best assessed by examining entire networks rather than single sites of structural lesions. Relating these changes to post-stroke neuropsychological variables and motor abilities will improve understanding of functional plasticity after stroke. Inclusion of healthy controls will provide additional insight into children's normal brain development. Resting state functional magnetic resonance imaging is a valid approach to topographically investigate the reorganisation of functional networks after a brain lesion. Transcranial magnetic stimulation provides complementary output information. This study will investigate functional reorganisation after paediatric arterial ischaemic stroke by means of resting state functional magnetic resonance imaging and transcranial magnetic stimulation in a cross-sectional plus longitudinal study design. The general aim of this study is to better understand neuroplasticity of the developing brain after stroke in order to develop more efficacious therapy and to improve the post-stroke functional outcome. METHODS The cross-sectional part of the study will investigate the functional cerebral networks of 35 children with chronic arterial ischaemic stroke (time of the lesion >2 years). In the longitudinal part, 15 children with acute arterial ischaemic stroke (shortly after the acute phase of the stroke) will be included and investigations will be performed 3 times within the subsequent 9 months. We will also recruit 50 healthy controls, matched for age and sex. The neuroimaging and neurophysiological data will be correlated with neuropsychological and neurological variables. DISCUSSION This study is the first to combine resting state functional magnetic resonance imaging and transcranial magnetic stimulation in a paediatric population diagnosed with arterial ischaemic stroke. Thus, this study has the potential to uniquely contribute to the understanding of neuronal plasticity in the brains of healthy children and those with acute or chronic brain injury. It is expected that the results will lead to the development of optimal interventions after arterial ischaemic stroke.

Differential activation of the lateral premotor cortex during action observation

Background Action observation leads to neural activation of the human premotor cortex. This study examined how the level of motor expertise (expert vs. novice) in ballroom dancing and the visual viewpoint (internal vs. external viewpoint) influence this activation within different parts of this area of the brain. Results Sixteen dance experts and 16 novices observed ballroom dance videos from internal or external viewpoints while lying in a functional magnetic resonance imaging scanner. A conjunction analysis of all observation conditions showed that action observation activated distinct networks of premotor, parietal, and cerebellar structures. Experts revealed increased activation in the ventral premotor cortex compared to novices. An internal viewpoint led to higher activation of the dorsal premotor cortex. Conclusions The present results suggest that the ventral and dorsal premotor cortex adopt differential roles during action observation depending on the level of motor expertise and the viewpoint.

Genetics of male nuptial colour divergence between sympatric sister species of a Lake Victoria cichlid fish

The hypothesis of sympatric speciation by sexual selection has been contentious. Several recent theoretical models of sympatric speciation by disruptive sexual selection were tailored to apply to African cichlids. Most of this work concludes that the genetic architecture of female preference and male trait is a key determinant of the likelihood of disruptive sexual selection to result in speciation. We investigated the genetic architecture controlling male nuptial colouration in a sympatric sibling species pair of cichlid fish from Lake Victoria, which differ conspicuously in male colouration and female mating preferences for these. We estimated that the difference between the species in male nuptial red colouration is controlled by a minimum number of two to four genes with significant epistasis and dominance effects. Yellow colouration appears to be controlled by one gene with complete dominance. The two colours appear to be epistatically linked. Knowledge on how male colouration segregates in hybrid generations and on the number of genes controlling differences between species can help us assess whether assumptions made in simulation models of sympatric speciation by sexual selection are realistic. In the particular case of the two sister species that we studied a small number of genes causing major differences in male colouration may have facilitated the divergence in male colouration associated with speciation.

Parallel amino acid replacements in the rhodopsins of the rockfishes (Sebastes spp.) associated with shifts in habitat depth

Among various groups of fishes, a shift in peak wavelength sensitivity has been correlated with changes in their photic environments. The genus Sebastes is a radiation of marine fish species that inhabit a wide range of depths from intertidal to over 600 m. We examined 32 species of Sebastes for evidence of adaptive amino acid substitution at the rhodopsin gene. Fourteen amino acid positions were variable among these species. Maximum likelihood analyses identify several of these to be targets of positive selection. None of these correspond to previously identified critical amino acid sites, yet they may in fact be functionally important. The occurrence of independent parallel changes at certain amino acid positions reinforces this idea. Reconstruction of habitat depths of ancestral nodes in the phylogeny suggests that shallow habitats have been colonized independently in different lineages. The evolution of rhodopsin appears to be associated with changes in depth, with accelerated evolution in lineages that have had large changes in depth.

Genome-wide investigation reveals high evolutionary rates in annual model plants

Background: ;Rates of molecular evolution vary widely among species. While significant deviations from molecular clock have been found in many taxa, effects of life histories on molecular evolution are not fully understood. In plants, annual/perennial life history traits have long been suspected to influence the evolutionary rates at the molecular level. To date, however, the number of genes investigated on this subject is limited and the conclusions are mixed. To evaluate the possible heterogeneity in evolutionary rates between annual and perennial plants at the genomic level, we investigated 85 nuclear housekeeping genes, 10 non-housekeeping families, and 34 chloroplast;genes using the genomic data from model plants including Arabidopsis thaliana and Medicago truncatula for annuals and grape (Vitis vinifera) and popular (Populus trichocarpa) for perennials.;Results: ;According to the cross-comparisons among the four species, 74-82% of the nuclear genes and 71-97% of the chloroplast genes suggested higher rates of molecular evolution in the two annuals than those in the two perennials. The significant heterogeneity in evolutionary rate between annuals and perennials was consistently found both in nonsynonymous sites and synonymous sites. While a linear correlation of evolutionary rates in orthologous genes between species was observed in nonsynonymous sites, the correlation was weak or invisible in synonymous sites. This tendency was clearer in nuclear genes than in chloroplast genes, in which the overall;evolutionary rate was small. The slope of the regression line was consistently lower than unity, further confirming the higher evolutionary rate in annuals at the genomic level.;Conclusions: ;The higher evolutionary rate in annuals than in perennials appears to be a universal phenomenon both in nuclear and chloroplast genomes in the four dicot model plants we investigated. Therefore, such heterogeneity in evolutionary rate should result from factors that have genome-wide influence, most likely those associated with annual/perennial life history. Although we acknowledge current limitations of this kind of study, mainly due to a small sample size available and a distant taxonomic relationship of the model organisms, our results indicate that the genome-wide survey is a promising approach toward further understanding of the;mechanism determining the molecular evolutionary rate at the genomic level.

Investigating human audio-visual object perception with a combination of hypothesis-generating and hypothesis-testing fMRI analysis tools

Primate multisensory object perception involves distributed brain regions. To investigate the network character of these regions of the human brain, we applied data-driven group spatial independent component analysis (ICA) to a functional magnetic resonance imaging (fMRI) data set acquired during a passive audio-visual (AV) experiment with common object stimuli. We labeled three group-level independent component (IC) maps as auditory (A), visual (V), and AV, based on their spatial layouts and activation time courses. The overlap between these IC maps served as definition of a distributed network of multisensory candidate regions including superior temporal, ventral occipito-temporal, posterior parietal and prefrontal regions. During an independent second fMRI experiment, we explicitly tested their involvement in AV integration. Activations in nine out of these twelve regions met the max-criterion (A < AV > V) for multisensory integration. Comparison of this approach with a general linear model-based region-of-interest definition revealed its complementary value for multisensory neuroimaging. In conclusion, we estimated functional networks of uni- and multisensory functional connectivity from one dataset and validated their functional roles in an independent dataset. These findings demonstrate the particular value of ICA for multisensory neuroimaging research and using independent datasets to test hypotheses generated from a data-driven analysis.

Repeatability and heritability of behavioural types in a social cichlid

Aim. The quantitative genetics underlying correlated behavioural traits (‘‘animal personality’’) have hitherto been studied mainly in domesticated animals. Here we report the repeatability ( ) and heritability ( ) of behavioural types in the highly social cichlid fish Neolamprologus pulcher. Methods. We tested 1779 individuals repeatedly and calculated the of behavioural types by variance components estimation (GLMM REML), using 1327 offspring from 162 broods from 74 pairs. Results. Repeatability of behavioural types was significant and considerable (0.546), but declined from 0.83 between tests conducted on the same day, to 0.19 on tests conducted up to 1201 days apart. All estimates were significant but low (e.g., pair identity SE). Additionally, we found significant variation between broods nested within the parent(s), but these were not related to several environmental factors tested. Conclusions. We conclude that despite a considerable , in this cichlid species is low, and variability in behavioural type appears to be strongly affected by other (non)genetic effects.

Community Genetics Reveal Elevated Levels of Sympatric Gene Flow among Morphologically Similar but Not among Morphologically Dissimilar Species of Lake Victoria Cichlid Fish

We examined genetic structure among five species of Lake Victoria haplochromine cichlids in four island communities, using a full factorial sampling design that compared genetic differentiation between pairs of species and populations of varying morphological similarity and geographical proximity. We found that allopatric conspecific populations were on average significantly more strongly differentiated than sympatric heterospecific populations of morphologically similar species. Allopatric heterospecific populations of morphologically dissimilar species were most differentiated. Our work demonstrates that phenotypic divergence can be maintained and perhaps even evolve in sympatry despite considerable gene flow between species. Conversely, phenotypic resemblance among conspecific populations can be maintained despite geographical isolation. Additionally we show that anthropogenically increased hybridization does not affect all sympatric species evenly but predominantly affects morphologically similar and closely related species. This has important implications for the evolution of reproductive isolation between species These findings are also consistent with the hypothesis of speciation reversal due to weakening of divergent selection and reproductive isolation as a consequence of habitat homogenization and offers an evolutionary mechanistic explanation for the observation that species poor assemblages in turbid areas of the lake are characterized by just one or two species in each of a few morphologically distinct genera.