Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Inflammatory responses in aggregating rat brain cell cultures subjected to different demyelinating conditions.

To study inflammatory reactions occurring in relation to demyelination, aggregating rat brain cell cultures were subjected to three different demyelinating insults, i.e., (i) lysophosphatidylcholine (LPC), (ii) interferon-gamma combined with lipopolysaccharide (IFN-gamma+LPS), and (iii) anti-MOG antibodies plus complement (alpha-MOG+C). Demyelination was assessed by measuring the expression of myelin basic protein (MBP) and myelin oligodendrocyte glycoprotein (MOG), and the activity of 2',3'-cyclic nucleotide 3'-phosphohydrolase (CNP). The accompanying inflammatory reactions were examined by the quantification of microglia-specific staining, by immunostaining for glial fibrillary acidic protein (GFAP), and by measuring the mRNA expression of a panel of inflammation-related genes. It was found that all three demyelinating insults decreased the expression of MBP and MOG, and induced microglial reactivity. LPC and alpha-MOG+C, but not IFN-gamma+LPS, decreased CNP activity; they also caused the appearance of macrophagic microglia, and increased GFAP staining indicating astrogliosis. LPC affected also the integrity of neurons and astrocytes. LPC and IFN-gamma+LPS upregulated the expression of the inflammation-related genes IL-6, TNF-alpha, Ccl5, Cxcl1, and iNOS, although to different degrees. Other inflammatory markers were upregulated by only one of the three insults, e.g., Cxcl2 by LPC; IL-1beta and IL-15 by IFN-gamma+LPS; and IFN-gamma by alpha-MOG+C. These findings indicate that each of the three demyelinating insults caused distinct patterns of demyelination and inflammatory reactivity, and that of the demyelinating agents tested only LPC exhibited general toxicity.

Thyrotropin-releasing hormone-stimulated thyrotropin expression involves islet-brain-1/c-Jun N-terminal kinase interacting protein-1.

Islet-brain-1 (IB1)/c-Jun N-terminal kinase interacting protein 1 (JIP-1) is a scaffold protein that is expressed at high levels in neurons and the endocrine pancreas. IB1/JIP-1 interacts with the c-Jun N-terminal kinase and mediates the specific physiological stimuli (such as cytokines). However, the potential role of the protein in the pituitary has not been evaluated. Herein, we examined expression of the gene encoding IB1/JIP-1 and its translated product in the anterior pituitary gland and a pituitary cell line, GH3. We then examined the potential role of IB1/JIP-1 in controlling TSH-beta gene expression. Exposure of GH3 cells to TRH stimulated the expression of IB1/JIP-1 protein levels, mRNA, and transcription of the promoter. The increase of IB1/JIP-1 content by transient transfection study of a vector encoding IB1/JIP-1 or by the stimulation of TRH stimulates TSH-beta promoter activity. This effect is not found in the presence of a mutated nonfunctional (IB1S59N) IB1/JIP-1 protein. Together, these facts point to a central role of the IB1/JIP-1 protein in the control of TRH-mediated TSH-beta stimulation.

Constructing brain functional networks from EEG: partial and unpartial correlations.

We consider electroencephalograms (EEGs) of healthy individuals and compare the properties of the brain functional networks found through two methods: unpartialized and partialized cross-correlations. The networks obtained by partial correlations are fundamentally different from those constructed through unpartial correlations in terms of graph metrics. In particular, they have completely different connection efficiency, clustering coefficient, assortativity, degree variability, and synchronization properties. Unpartial correlations are simple to compute and they can be easily applied to large-scale systems, yet they cannot prevent the prediction of non-direct edges. In contrast, partial correlations, which are often expensive to compute, reduce predicting such edges. We suggest combining these alternative methods in order to have complementary information on brain functional networks.

The role of energetic value in dynamic brain response adaptation during repeated food image viewing.

The repeated presentation of simple objects as well as biologically salient objects can cause the adaptation of behavioral and neural responses during the visual categorization of these objects. Mechanisms of response adaptation during repeated food viewing are of particular interest for better understanding food intake beyond energetic needs. Here, we measured visual evoked potentials (VEPs) and conducted neural source estimations to initial and repeated presentations of high-energy and low-energy foods as well as non-food images. The results of our study show that the behavioral and neural responses to food and food-related objects are not uniformly affected by repetition. While the repetition of images displaying low-energy foods and non-food modulated VEPs as well as their underlying neural sources and increased behavioral categorization accuracy, the responses to high-energy images remained largely invariant between initial and repeated encounters. Brain mechanisms when viewing images of high-energy foods thus appear less susceptible to repetition effects than responses to low-energy and non-food images. This finding is likely related to the superior reward value of high-energy foods and might be one reason why in particular high-energetic foods are indulged although potentially leading to detrimental health consequences.

Délire religieux treize ans après traumatisme crânien [Religious delusion 13 years after brain injury].

We report here with a case of religious delusion in a 39 years old woman. She had suffered a head injury with right temporal concussion 13 years before but had no earlier history of psychiatric disorder. In view of the fact that this acute psychiatric state lasted for a short duration of time and that personality and affects were preserved, this incident is compared to the schizophreniform disorder of the type DSM-III-R. The hypothesis of an acquired predisposition due to head injury has been put forward as an explanation.

Comparison of T1 relaxation times of the neurochemical profile in rat brain at 9.4 tesla and 14.1 tesla.

Knowledge of T(1) relaxation times can be important for accurate relative and absolute quantification of brain metabolites, for sensitivity optimizations, for characterizing molecular dynamics, and for studying changes induced by various pathological conditions. (1)H T(1) relaxation times of a series of brain metabolites, including J-coupled ones, were determined using a progressive saturation (PS) technique that was validated with an adiabatic inversion-recovery (IR) method. The (1)H T(1) relaxation times of 16 functional groups of the neurochemical profile were measured at 14.1T and 9.4T. Overall, the T(1) relaxation times found at 14.1T were, within the experimental error, identical to those at 9.4T. The T(1)s of some coupled spin resonances of the neurochemical profile were measured for the first time (e.g., those of gamma-aminobutyrate [GABA], aspartate [Asp], alanine [Ala], phosphoethanolamine [PE], glutathione [GSH], N-acetylaspartylglutamate [NAAG], and glutamine [Gln]). Our results suggest that T(1) does not increase substantially beyond 9.4T. Furthermore, the similarity of T(1) among the metabolites (approximately 1.5 s) suggests that T(1) relaxation time corrections for metabolite quantification are likely to be similar when using rapid pulsing conditions. We therefore conclude that the putative T(1) increase of metabolites has a minimal impact on sensitivity when increasing B(0) beyond 9.4T.

Epidermal growth factor (EGF) stimulation of cultured brain cells. I. Enhancement of the developmental increase in glial enzymatic activity.

Serum-free aggregating cell cultures of fetal rat telencephalon grown in the presence of 3 ng/ml (5 X 10(-10) M) epidermal growth factor (EGF) until day 12 showed 2- to 3-fold increased activities in the two glial enzymes, glutamine synthetase (GLU-S) and 2',3'-cyclic nucleotide 3'-phosphohydrolase (CNPase). This effect was concentration-dependent, with maximal stimulation in cultures treated daily with 3 ng/ml EGF. Addition of EGF during the first 10 culture days was sufficient to produce a maximal stimulation of both GLU-S and CNPase on day 19, whereas treatments starting on day 12 were ineffective. The stimulation of GLU-S preceded that of CNPase. The EGF-induced increase in GLU-S activity was not directly dependent on the presence of insulin, triiodothyronine, or hydrocortisone in the medium, whereas insulin was required for the stimulation of CNPase. A single dose of 5 ng/ml EGF on day 2 caused a slight but significant decrease in DNA synthesis after day 6. The present results indicate that in serum-free aggregating cell cultures of fetal rat telencephalon EGF partially inhibits DNA synthesis, and stimulates an early step in glial differentiation.

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone.

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long-standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well-studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best-supported scenario was a long period of allopatric isolation during the first three-quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.

Hierarchical analyses of genetic differentiation in a hybrid zone of Sorex araneus (Insectivora : Soricidae)

Microsatellites are used to unravel the fine-scale genetic structure of a hybrid zone between chromosome races Valais and Cordon of the common shrew (Sorex araneus) located in the French Alps. A total of 269 individuals collected between 1992 and 1995 was typed for seven microsatellite loci. A modified version of the classical multiple correspondence analysis is carried out. This analysis clearly shows the dichotomy between the two races. Several approaches are used to study genetic structuring. Gene flow is clearly reduced between these chromosome races and is estimated at one migrant every two generations using X-statistics and one migrant per generation using F-statistics. Hierarchical F- and R-statistics are compared and their efficiency to detect inter- and intraracial patterns of divergence is discussed. Within-race genetic structuring is significant, but remains weak. F-ST displays similar values on both sides of the hybrid zone, although no environmental barriers are found on the Cordon side, whereas the Valais side is divided by several mountain rivers. We introduce the exact G-test to microsatellite data which proved to be a powerful test to detect genetic differentiation within as well as among races. The genetic background of karyotypic hybrids was compared with the genetic background of pure parental forms using a CRT-MCA. Our results indicate that, without knowledge of the karyotypes, we would not have been able to distinguish these hybrids from karyotypically pure samples.

The neuronal adhesion protein D2 in differentiating aggregates of brain cells.

The D2-protein is a high molecular weight protein involved in interneuronal adhesion. The concentration of D2-protein was measured both in aggregates of fetal rat telencephalic cells cultured in a chemically defined medium and in developing forebrain. Both the concentration of the D2-protein and the degree of sialylation were changed in the cultures in parallel with the corresponding values obtained from postnatal forebrain. In the cultures the highest specific concentration of D2-protein was observed after 12 days in culture. This value was 2.7 times higher than the average value of adult rat forebrain. Antibodies to D2-protein have previously been shown to inhibit fasciculation of neuritic fibers extending from cultured explants of sympathetic ganglia. We investigated the effect of such antibodies on the differentiation of aggregating telencephalic cells. By adding surplus antibodies to the cultures from day 11 to day 16 we were able to decrease the specific concentration of D2-protein on the neurons by 53% measured at day 19. The decrease was not compensated fully even after further 10 days in the culture. Although the concentration of D2-protein was decreased during the period of synaptogenesis no change was found in the specific concentration of a marker of mature synapses, the D3-protein. Thus, in this culture system synaptogenesis could proceed to an unimpaired extent in the presence of a decreased concentration of a putatively involved adhesion molecule. However, the specific concentration of two markers of myelination, 2',3'-cyclic nucleotide 3'-phosphodiesterase and myelin basic protein, were both increased, suggesting an antibody-induced stimulation of myelination in the cultured aggregates.

Severe traumatic brain injury in a high-income country: an epidemiological study.

This adult cohort determined the incidence and patients' short-term outcomes of severe traumatic brain injury (sTBI) in Switzerland and age-related differences. A prospective cohort study with a follow-up at 14 days was performed. Patients ≥16 years of age sustaining sTBI and admitted to 1 of 11 trauma centers were included. sTBI was defined by an Abbreviated Injury Scale of the head (HAIS) score >3. The centers participated from 6 months to 3 years. The results are presented as percentages, medians, and interquartile ranges (IQRs). Subgroup analyses were performed for patients ≤65 years (younger) and >65 (elderly). sTBI was observed in 921 patients (median age, 55 years; IQR, 33-71); 683 (74.2%) were male. Females were older (median age, 67 years; IQR, 42-80) than males (52; IQR, 31-67; p<0.00001). The estimated incidence was 10.58 per 100,000 inhabitants per year. Blunt trauma was observed in 879 patients (95.4%) and multiple trauma in 283 (30.7%). Median Glasgow Coma Score (GCS) on the scene was 9 (IQR 4-14; 8 in younger, 12 in elderly) and in emergency departments 5 (IQR, 3-14; 3 in younger, 8 in elderly). Trauma mechanisms included the following: 484 patients with falls (52.6%; younger, 242 patients [50.0%]; elderly, 242 [50.0%]), 291 with road traffic accidents (31.6%; younger, 237 patients [81.4%]; elderly, 54 [18.6%]), and 146 with others (15.8%). Mortality was 30.2% (24.5% in younger, 40.9% in elderly). Median GCS at 14 days was 15 (IQR, 14-15) without differences among subgroups. Estimated incidence of sTBI in Switzerland was low, age was high, and mortality considerable. The elderly had higher initial GCS and a higher death rate, but high GCS at 14 days.

Computational anatomy for studying use-dependant brain plasticity.

In this article we provide a comprehensive literature review on the in vivo assessment of use-dependant brain structure changes in humans using magnetic resonance imaging (MRI) and computational anatomy. We highlight the recent findings in this field that allow the uncovering of the basic principles behind brain plasticity in light of the existing theoretical models at various scales of observation. Given the current lack of in-depth understanding of the neurobiological basis of brain structure changes we emphasize the necessity of a paradigm shift in the investigation and interpretation of use-dependent brain plasticity. Novel quantitative MRI acquisition techniques provide access to brain tissue microstructural properties (e.g., myelin, iron, and water content) in-vivo, thereby allowing unprecedented specific insights into the mechanisms underlying brain plasticity. These quantitative MRI techniques require novel methods for image processing and analysis of longitudinal data allowing for straightforward interpretation and causality inferences.