Detection of Optimal Models in Parameter Space with Support Vector Machines

The paper proposes an approach aimed at detecting optimal model parameter combinations to achieve the most representative description of uncertainty in the model performance. A classification problem is posed to find the regions of good fitting models according to the values of a cost function. Support Vector Machine (SVM) classification in the parameter space is applied to decide if a forward model simulation is to be computed for a particular generated model. SVM is particularly designed to tackle classification problems in high-dimensional space in a non-parametric and non-linear way. SVM decision boundaries determine the regions that are subject to the largest uncertainty in the cost function classification, and, therefore, provide guidelines for further iterative exploration of the model space. The proposed approach is illustrated by a synthetic example of fluid flow through porous media, which features highly variable response due to the parameter values' combination.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

The sandstone-hosted Beverley uranium Deposit, Lake Frome Basin, South Australia: Mineralogy, geochemistry, and a time-constrained model for its genesis

The sandstone-hosted Beverley uranium deposit is located in terrestrial sediments in the Lake Frome basin in the North Flinders Ranges, South Australia. The deposit is 13 km from the U-rich Mesoproterozoic basement of the Mount Painter inlier, which is being uplifted 100 to 200 m above the basin by neotectonic activity that probably initiated in the early Pliocene. The mineralization was deposited mainly in organic matter-poor Miocene lacustrine sands and partly in the underlying reductive strata comprising organic matter-rich clays and silts. The bulk of the mineralization consists of coffinite and/or uraninite nodules, growing around Co-rich pyrite with an S isotope composition (delta S-34 = 1.0 +/- 0.3 parts per thousand), suggestive of an early diagenetic lacustrine origin. In contrast, authigenic sulfides in the bulk of the sediments have a negative S isotope signature (delta S-34 ranges from -26.2 to -35.5 parts per thousand), indicative of an origin via bacterially mediated sulfate reduction. Minor amounts of Zn-bearing native copper and native lead also support the presence of specific, reducing microenvironments in the ore zone. Small amounts of carnotite are associated with the coffinite ore and also occur beneath a paleosoil horizon overlying the uranium deposit. Provenance studies suggest that the host Miocene sediments were derived from the reworking of Early Cretaceous glacial or glaciolacustrine sediments ultimately derived from Paleozoic terranes in eastern Australia. In contrast, the overlying Pliocene strata were in part derived from the Mesoproterozoic basement inlier. Mass-balance and geochemical data confirm that granites of the Mount Painter domain were the ultimate source of U and BEE at Beverley. U-Pb dating of coffinite and carnotite suggest that the U mineralization is Pliocene (6.7-3.4 Ma). The suitability of the Beverley deposit for efficient mining via in situ leaching, and hence its economic value, are determined by the nature of the hosting sand unit, which provides the permeability and low reactivity required for high fluid flow and low chemical consumption. These favorable sedimentologic and geometrical features result from a complex conjunction of factors, including deposition in lacustrine shore environment, reworking of angular sands of glacial origin, deep Pliocene weathering, and proximity to an active fault exposing extremely U rich rocks.

Emplacement of ultramafic rocks into the continental crust monitored by light and other trace elements: An example from the Geisspfad body (Swiss-Italian Alps)

In order to evaluate the influence of continental crustal rocks on trace element budgets of serpentinized peridotites incorporated into the continental crust, we have analyzed the chemical composition of whole rock samples and minerals of the Geisspfad ultramafic complex (Swiss-Italian Alps). This complex represents a relict oceanic succession composed of serpentinites, ophicarbonates and metabasic rocks, emplaced into crustal gneisses during Alpine collision. Following peak metamorphic amphibolite facies conditions, fluid flow modified some of the trace element contents of ophicarbonates and deformed serpentinites close to the contact with country rocks. The fluid originated from the surrounding continental crustal rocks as documented by the increase of Pb in the serpentinites, and by the strongly negative all) values (-112 parts per thousand) of some ultramafic rocks close to the contact with surrounding gneisses. Little or no modification of the fluid mobile elements Li, B or U was observed in the serpentinite. In-situ analysis of light elements of serpentinite minerals indicate redistribution of light elements coupled to changes of mineral modes towards the outer 100-150 m of the massif. In the centre of the massif, Li is preferentially concentrated in olivine, while Be and B are hosted by tremolite. In contrast, at the outer rim of the massif, Li and Be are preferentially incorporated into diopside, and B into antigorite. This redistribution of light elements among the different minerals is visible in the serpentinite, at a maximum distance of -100-150 m from the ophicarbonate-metabasite contact. Our results show that interaction of ultramafic rocks and crust-derived fluids can be easily detected by studies of Pb and partial derivative D in whole rocks. We argue that small ultramafic bodies potentially record an emplacement-related trace element signature, and that crustal light element values in ultramafic rocks are not necessarily derived from a subducting slab. (C) 2008 Elsevier B.V. All rights reserved.

Benchmarking therapeutic drug monitoring software: A systematic evaluation of available computer tools

Introduction: Therapeutic drug monitoring (TDM) aims at optimizing treatment by individualizing dosage regimen based on measurement of blood concentrations. Maintaining concentrations within a target range requires pharmacokinetic and clinical capabilities. Bayesian calculation represents a gold standard in TDM approach but requires computing assistance. In the last decades computer programs have been developed to assist clinicians in this assignment. The aim of this benchmarking was to assess and compare computer tools designed to support TDM clinical activities.¦Method: Literature and Internet search was performed to identify software. All programs were tested on common personal computer. Each program was scored against a standardized grid covering pharmacokinetic relevance, user-friendliness, computing aspects, interfacing, and storage. A weighting factor was applied to each criterion of the grid to consider its relative importance. To assess the robustness of the software, six representative clinical vignettes were also processed through all of them.¦Results: 12 software tools were identified, tested and ranked. It represents a comprehensive review of the available software's characteristics. Numbers of drugs handled vary widely and 8 programs offer the ability to the user to add its own drug model. 10 computer programs are able to compute Bayesian dosage adaptation based on a blood concentration (a posteriori adjustment) while 9 are also able to suggest a priori dosage regimen (prior to any blood concentration measurement), based on individual patient covariates, such as age, gender, weight. Among those applying Bayesian analysis, one uses the non-parametric approach. The top 2 software emerging from this benchmark are MwPharm and TCIWorks. Other programs evaluated have also a good potential but are less sophisticated (e.g. in terms of storage or report generation) or less user-friendly.¦Conclusion: Whereas 2 integrated programs are at the top of the ranked listed, such complex tools would possibly not fit all institutions, and each software tool must be regarded with respect to individual needs of hospitals or clinicians. Interest in computing tool to support therapeutic monitoring is still growing. Although developers put efforts into it the last years, there is still room for improvement, especially in terms of institutional information system interfacing, user-friendliness, capacity of data storage and report generation.

Do seismic waves sense fracture connectivity?

A defining characteristic of fractured rocks is their very high level of seismic attenuation, which so far has been assumed to be mainly due to wave-induced fluid flow (WIFF) between the fractures and the pore space of the embedding matrix. Using oscillatory compressibility simulations based on the quasi-static poroelastic equations, we show that another important, and as of yet undocumented, manifestation of WIFF is at play in the presence of fracture connectivity. This additional energy loss is predominantly due to fluid flow within the connected fractures and is sensitive to their lengths, permeabilities, and intersection angles. Correspondingly, it contains key information on the governing hydraulic properties of fractured rock masses and hence should be accounted for whenever realistic seismic models of such media are needed.

Effects of colonization asymmetries on metapopulation persistence.

Ocean currents, prevailing winds, and the hierarchical structures of river networks are known to create asymmetries in re-colonization between habitat patches. The impacts of such asymmetries on metapopulation persistence are seldom considered, especially rarely in theoretical studies. Considering three classical models (the island, the stepping stone and the distance-dependent model), we explore how metapopulation persistence is affected by (i) asymmetry in dispersal strength, in which the colonization rate between two patches differs in direction, and (ii) asymmetry in connectivity, in which the overall colonization pattern displays asymmetry (circulating or dendritic networks). Viability can be drastically reduced when directional bias in dispersal strength is higher than 25%. Re-colonization patterns that allow for strong local connectivity provide the highest persistence compared to systems that allow circulation. Finally, asymmetry has relatively weak effects when metapopulations maintain strong general connectivity.

Link between genotype and antimicrobial resistance in bovine mastitis-related Staphylococcus aureus strains, determined by comparing Swiss and French isolates from the Rhône Valley.

Staphylococcus aureus is a major bovine mastitis pathogen. Although the reported antimicrobial resistance was generally low, the emergence of new genetic clusters in bovine mastitis requires examination of the link between antimicrobial resistance and genotypes. Here, amplified fragment length polymorphism (AFLP) profiles and standard antimicrobial resistance profiles were determined in order to characterize a total of 343 S. aureus cow mastitis isolates from two geographically close regions of Switzerland and France. AFLP profiles revealed similar population compositions in the two regions, with 4 major clusters (C8, C20, C97, and C151), but the proportions of isolates in each cluster significantly diverged between the two countries (P = 9.2 × 10⁻⁹). Antimicrobial resistance was overall low (< 5% resistance to all therapeutically relevant molecules), with the exception of penicillin resistance, which was detected in 26% of the isolates. Penicillin resistance proportions differed between clusters, with only 1 to 2% of resistance associated with C20 and C151 and up to 70% associated with bovine C97. The prevalence of C20 and C8 was unexpectedly high and requires further investigation into the mechanism of adaptation to the bovine host. The strong association of penicillin resistance with few clusters highlights the fact that the knowledge of local epidemiology is essential for rational choices of antimicrobial treatment in the absence of susceptibility testing. Taken together, these observations argue in favor of more routine scrutiny of antimicrobial resistance and antibiotic-resistant clones in cattle and the farm environment.

Diagnostic methods and treatment options for focal cortical dysplasia.

Our inability to adequately treat many patients with refractory epilepsy caused by focal cortical dysplasia (FCD), surgical inaccessibility and failures are significant clinical drawbacks. The targeting of physiologic features of epileptogenesis in FCD and colocalizing functionality has enhanced completeness of surgical resection, the main determinant of outcome. Electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) and magnetoencephalography are helpful in guiding electrode implantation and surgical treatment, and high-frequency oscillations help defining the extent of the epileptogenic dysplasia. Ultra high-field MRI has a role in understanding the laminar organization of the cortex, and fluorodeoxyglucose-positron emission tomography (FDG-PET) is highly sensitive for detecting FCD in MRI-negative cases. Multimodal imaging is clinically valuable, either by improving the rate of postoperative seizure freedom or by reducing postoperative deficits. However, there is no level 1 evidence that it improves outcomes. Proof for a specific effect of antiepileptic drugs (AEDs) in FCD is lacking. Pathogenic mutations recently described in mammalian target of rapamycin (mTOR) genes in FCD have yielded important insights into novel treatment options with mTOR inhibitors, which might represent an example of personalized treatment of epilepsy based on the known mechanisms of disease. The ketogenic diet (KD) has been demonstrated to be particularly effective in children with epilepsy caused by structural abnormalities, especially FCD. It attenuates epigenetic chromatin modifications, a master regulator for gene expression and functional adaptation of the cell, thereby modifying disease progression. This could imply lasting benefit of dietary manipulation. Neurostimulation techniques have produced variable clinical outcomes in FCD. In widespread dysplasias, vagus nerve stimulation (VNS) has achieved responder rates >50%; however, the efficacy of noninvasive cranial nerve stimulation modalities such as transcutaneous VNS (tVNS) and noninvasive (nVNS) requires further study. Although review of current strategies underscores the serious shortcomings of treatment-resistant cases, initial evidence from novel approaches suggests that future success is possible.

Vasopressin-dependent coupling between sodium transport and water flow in a mouse cortical collecting duct cell line.

Water balance is achieved through the ability of the kidney to control water reabsorption in the connecting tubule and the collecting duct. In a mouse cortical collecting duct cell line (mCCD(c11)), physiological concentrations of arginine vasopressin increased both electrogenic, amiloride-sensitive, epithelial sodium channel (ENaC)-mediated sodium transport measured by the short-circuit current (Isc) method and water flow (Jv apical to basal) measured by gravimetry with similar activation coefficient K(1/2) (6 and 12 pM, respectively). Jv increased linearly according to the osmotic gradient across the monolayer. A small but highly significant Jv was also measured under isoosmotic conditions. To test the coupling between sodium reabsorption and water flow, mCCD(c11) cells were treated for 24 h under isoosmotic condition with either diluent, amiloride, vasopressin or vasopressin and amiloride. Isc, Jv, and net chemical sodium fluxes were measured across the same monolayers. Around 30% of baseline and 50% of vasopressin-induced water flow is coupled to an amiloride-sensitive, ENaC-mediated, electrogenic sodium transport, whereas the remaining flow is coupled to an amiloride-insensitive, nonelectrogenic sodium transport mediated by an unknown electroneutral transporter. The mCCD(c11) cell line is a first example of a mammalian tight epithelium allowing quantitative study of the coupling between sodium and water transport. Our data are consistent with the 'near isoosmotic' fluid transport model.

Computational fluid dynamics of the right ventricular outflow tract and of the pulmonary artery: a bench model of flow dynamics.

OBJECTIVES: The reconstruction of the right ventricular outflow tract (RVOT) with valved conduits remains a challenge. The reoperation rate at 5 years can be as high as 25% and depends on age, type of conduit, conduit diameter and principal heart malformation. The aim of this study is to provide a bench model with computer fluid dynamics to analyse the haemodynamics of the RVOT, pulmonary artery, its bifurcation, and left and right pulmonary arteries that in the future may serve as a tool for analysis and prediction of outcome following RVOT reconstruction. METHODS: Pressure, flow and diameter at the RVOT, pulmonary artery, bifurcation of the pulmonary artery, and left and right pulmonary arteries were measured in five normal pigs with a mean weight of 24.6 ± 0.89 kg. Data obtained were used for a 3D computer fluid-dynamics simulation of flow conditions, focusing on the pressure, flow and shear stress profile of the pulmonary trunk to the level of the left and right pulmonary arteries. RESULTS: Three inlet steady flow profiles were obtained at 0.2, 0.29 and 0.36 m/s that correspond to the flow rates of 1.5, 2.0 and 2.5 l/min flow at the RVOT. The flow velocity profile was constant at the RVOT down to the bifurcation and decreased at the left and right pulmonary arteries. In all three inlet velocity profiles, low sheer stress and low-velocity areas were detected along the left wall of the pulmonary artery, at the pulmonary artery bifurcation and at the ostia of both pulmonary arteries. CONCLUSIONS: This computed fluid real-time model provides us with a realistic picture of fluid dynamics in the pulmonary tract area. Deep shear stress areas correspond to a turbulent flow profile that is a predictive factor for the development of vessel wall arteriosclerosis. We believe that this bench model may be a useful tool for further evaluation of RVOT pathology following surgical reconstructions.

Importance of ENaC-mediated sodium transport in alveolar fluid clearance using genetically-engineered mice.

The lung possesses specific transport systems that intra- and extracellularly maintain salt and fluid balance necessary for its function. At birth, the lungs rapidly transform into a fluid (Na(+))-absorbing organ to enable efficient gas exchange. Alveolar fluid clearance, which mainly depends on sodium transport in alveolar epithelial cells, is an important mechanism by which excess water in the alveoli is reabsorbed during the resolution of pulmonary edema. In this review, we will focus and summarize on the role of ENaC in alveolar lung liquid clearance and discuss recent data from mouse models with altered activity of epithelial sodium channel function in the lung, and more specifically in alveolar fluid clearance. Recent data studying mice with hyperactivity of ENaC or mice with reduced ENaC activity clearly illustrate the impaired lung fluid clearance in these adult mice. Further understanding of the physiological role of ENaC and its regulatory proteins implicated in salt and water balance in the alveolar cells may therefore help to develop new therapeutic strategies to improve gas exchange in pulmonary edema.

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.