Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Fingermarks and Other Impressions Left by the Human Body - A Review (August 2007-July 2010)

The purpose of this paper is to review the scientific literature from August 2007 to July 2010. The review is focused on more than 420 published papers. The review will not cover information coming from international meetings available only in abstract form. Fingermarks constitute an important chapter with coverage of the identification process as well as detection techniques on various surfaces. We note that the research has been very dense both at exploring and understanding current detection methods as well as bringing groundbreaking techniques to increase the number of marks detected from various objects. The recent report from the US National Research Council (NRC) is a milestone that has promoted a critical discussion on the state of forensic science and its associated research. We can expect a surge of interest in research in relation to cognitive aspect of mark and print comparison, establishment of relevant forensic error rates and statistical modelling of the selectivity of marks' attributes. Other biometric means of forensic identification such as footmarks or earmarks are also covered in the report. Compared to previous years, we noted a decrease in the number of submission in these areas. No doubt that the NRC report has set the seed for further investigation of these fields as well.

Evaluating sampling strategies and logistic regression methods for modelling complex land cover changes

The role of land cover change as a significant component of global change has become increasingly recognized in recent decades. Large databases measuring land cover change, and the data which can potentially be used to explain the observed changes, are also becoming more commonly available. When developing statistical models to investigate observed changes, it is important to be aware that the chosen sampling strategy and modelling techniques can influence results. We present a comparison of three sampling strategies and two forms of grouped logistic regression models (multinomial and ordinal) in the investigation of patterns of successional change after agricultural land abandonment in Switzerland. Results indicated that both ordinal and nominal transitional change occurs in the landscape and that the use of different sampling regimes and modelling techniques as investigative tools yield different results. Synthesis and applications. Our multimodel inference identified successfully a set of consistently selected indicators of land cover change, which can be used to predict further change, including annual average temperature, the number of already overgrown neighbouring areas of land and distance to historically destructive avalanche sites. This allows for more reliable decision making and planning with respect to landscape management. Although both model approaches gave similar results, ordinal regression yielded more parsimonious models that identified the important predictors of land cover change more efficiently. Thus, this approach is favourable where land cover change pattern can be interpreted as an ordinal process. Otherwise, multinomial logistic regression is a viable alternative.

Techniques for CT and MR, post-processing, radiation

Brain perfusion can be assessed by CT and MR. For CT, two major techniques are used. First, Xenon CT is an equilibrium technique based on a freely diffusible tracer. First pass of iodinated contrast injected intravenously is a second method, more widely available. Both methods are proven to be robust and quantitative, thanks to the linear relationship between contrast concentration and x-ray attenuation. For the CT methods, concern regarding x-ray doses delivered to the patients need to be addressed. MR is also able to assess brain perfusion using the first pass of gadolinium based contrast agent injected intravenously. This method has to be considered as a semi-quantitative because of the non linear relationship between contrast concentration and MR signal changes. Arterial spin labeling is another MR method assessing brain perfusion without injection of contrast. In such case, the blood flow in the carotids is magnetically labelled by an external radiofrequency pulse and observed during its first pass through the brain. Each of this various CT and MR techniques have advantages and limits that will be illustrated and summarized.Learning Objectives:1. To understand and compare the different techniques for brain perfusion imaging.2. To learn about the methods of acquisition and post-processing of brain perfusion by first pass of contrast agent for CT and MR.3. To learn about non contrast MR methods (arterial spin labelling).

Mirror me : the effect of the voting advice application smartvote on voting preferences and behavior of Swiss voters

Voting Advice Applications (VAAs) have become a central component of election campaigns worldwide. Through matching political preferences of voters to parties and candidates, the web application grants voters a look into their political mirror and reveals the most suitable political choices to them in terms of policy congruence. Both the dense and concise information on the electoral offer and the comparative nature of the application make VAAs an unprecedented information source for electoral decision making. In times where electoral choices are found to be highly individualized and driven by political issue positions, an ever increasing number of voters turn to VAAs before casting their ballots. With VAAs in high demand, the question of their effects on voters has become a pressing research topic. In various countries, survey research has been used to proclaim an impact of VAAs on electoral behavior, yet practically all studies fail to provide the scientific evidence that would allow for making such claims. In this thesis, I set out to systematically establish the causal link between VAA use and electoral behavior, using various data sources and appropriate statistical techniques in doing so. The focus lies on the Swiss VAA smartvote, introduced in the forefront of the 2003 Swiss federal elections and meanwhile an integral part of the national election campaign, smartvote has produced over a million voting recommendations in the last Swiss federal elections to an active electorate of two million, potentially guiding a vast amount of voters in their choices on the ballot. In order to determine the effect of the VAA on electoral behavior, I analyze both voting preferences and choice among Swiss voters during two consecutive election periods. First, I introduce statistical techniques to adequately examine VAA effects in observational studies and use them to demonstrate that voters who used smartvote prior to the 2007 Swiss federal elections were significantly more likely to swing vote in the elections than non- users. Second, I analyze preference voting during the same election and show that the smartvote voting recommendation inclines politically knowledgeable voters to modify their ballots and cast candidate specific preference votes. Third, to further tackle the indication that smartvote use affects the preference structure of voters, I employ an experimental research design to demonstrate that voters who use the application tend to strengthen their vote propensities for their most preferred party and adapt their overall party preferences in a way that they consider more than one party as eligible vote options after engaging with the application. Finally, vote choice is examined for the 2011 Swiss federal election, showing once more that the VAA initiated a change of party choice among voters. In sum, this thesis presents empirical evidence for the transformative effect of the Swiss VAA smartvote on the electoral behavior.

Care complexity, mood, and quality of life in liver pre-transplant patients.

RATIONALE: This study was intended to document the frequency of care complexity in liver transplant candidates, and its association with mood disturbance and poor health-related quality of life (HRQoL). METHODS: Consecutive patients fulfilling inclusion criteria, recruited in three European hospitals, were assessed with INTERMED, a reliable and valid method for the early assessment of bio-psychosocial health risks and needs. Blind to the results, they were also assessed with the Hospital Anxiety and Depression Scale (HADS). HRQoL was documented with the EuroQol and the SF36. Statistical analysis included multivariate and multilevel techniques. RESULTS: Among patients fulfilling inclusion criteria, 60 patients (75.9%) completed the protocol and 38.3% of them were identified as "complex" by INTERMED, but significant between-center differences were found. In support of the working hypothesis, INTERMED scores were significantly associated with all measures of both the SF36 and the EuroQol, and also with the HADS. A one point increase in the INTERMED score results in a reduction in 0.93 points in EuroQol and a 20% increase in HADS score. CONCLUSIONS: INTERMED-measured case complexity is frequent in liver transplant candidates but varies widely between centers. The use of this method captures in one instrument multiple domains of patient status, including mood disturbances and reduced HRQoL.

Comparison of the site distribution of melanoma in New Zealand and Canada

A comparison of the site distribution of cutaneous malignant melanoma in New Zealand and Canada was performed. This series deals with 41,331 incident cases registered between 1968 and 1990 and is the largest to date to evaluate the influence of age and gender on the site distribution of melanoma. Site-specific, age-standardized rates per unit surface area and relative tumour density were assessed by gender and country and differences compared with statistical techniques adapted to this context. The age-standardized rates for all sites were higher in New Zealand than in Canada, the ratio being 3.2 for men and 3.8 for women. Occurrence of melanoma was denser for chronically than intermittently exposed sites in both New Zealand and Canada. The highest incidence rate per unit area was for the ears in men which was more than 5 times the rate for the entire body in each country. For each gender, melanomas were relatively commoner on the trunk and the face in Canada, and on the lower limbs in New Zealand. The variations in the site distribution were similar in each country and consistent with the effect of differential patterns of sun exposure between genders. Our results show that the levels of risk of melanoma between phenotypically comparable populations exposed to different amount of UV radiation vary in a site-specific manner, especially for intermittently exposed sites. This suggests that both environmental conditions and lifestyle factors influence the site distribution of melanoma in these two populations.

Techniques for CT and MR, post-processing, radiation.

Brain perfusion can be assessed by CT and MR. For CT, two major techniquesare used. First, Xenon CT is an equilibrium technique based on a freely diffusibletracer. First pass of iodinated contrast injected intravenously is a second method,more widely available. Both methods are proven to be robust and quantitative,thanks to the linear relationship between contrast concentration and x-ray attenuation.For the CT methods, concern regarding x-ray doses delivered to the patientsneed to be addressed. MR is also able to assess brain perfusion using the firstpass of gadolinium based contrast agent injected intravenously. This method hasto be considered as a semi-quantitative because of the non linear relationshipbetween contrast concentration and MR signal changes. Arterial spin labelingis another MR method assessing brain perfusion without injection of contrast. Insuch case, the blood flow in the carotids is magnetically labelled by an externalradiofrequency pulse and observed during its first pass through the brain. Eachof this various CT and MR techniques have advantages and limits that will be illustratedand summarised.Learning Objectives:1. To understand and compare the different techniques for brain perfusionimaging.2. To learn about the methods of acquisition and post-processing of brainperfusion by first pass of contrast agent for CT and MR.3. To learn about non contrast MR methods (arterial spin labelling).

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

The development of decision limits for the GH-2000 detection methodology using additional insulin-like growth factor-I and amino-terminal pro-peptide of type III collagen assays.

The GH-2000 and GH-2004 projects have developed a method for detecting GH misuse based on measuring insulin-like growth factor-I (IGF-I) and the amino-terminal pro-peptide of type III collagen (P-III-NP). The objectives were to analyze more samples from elite athletes to improve the reliability of the decision limit estimates, to evaluate whether the existing decision limits needed revision, and to validate further non-radioisotopic assays for these markers. The study included 998 male and 931 female elite athletes. Blood samples were collected according to World Anti-Doping Agency (WADA) guidelines at various sporting events including the 2011 International Association of Athletics Federations (IAAF) World Athletics Championships in Daegu, South Korea. IGF-I was measured by the Immunotech A15729 IGF-I IRMA, the Immunodiagnostic Systems iSYS IGF-I assay and a recently developed mass spectrometry (LC-MS/MS) method. P-III-NP was measured by the Cisbio RIA-gnost P-III-P, Orion UniQ? PIIINP RIA and Siemens ADVIA Centaur P-III-NP assays. The GH-2000 score decision limits were developed using existing statistical techniques. Decision limits were determined using a specificity of 99.99% and an allowance for uncertainty because of the finite sample size. The revised Immunotech IGF-I - Orion P-III-NP assay combination decision limit did not change significantly following the addition of the new samples. The new decision limits are applied to currently available non-radioisotopic assays to measure IGF-I and P-III-NP in elite athletes, which should allow wider flexibility to implement the GH-2000 marker test for GH misuse while providing some resilience against manufacturer withdrawal or change of assays. Copyright © 2015 John Wiley & Sons, Ltd.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Augmentation cystoplasty using pedicled and de-epithelialized gastric patches in the mini-pig model.

PURPOSE: The most common methods of bladder augmentation are gastrocystoplasty and enterocystoplasty. Gastrocystoplasty is advantageous due to minimal mucous secretion and a well developed muscular wall as well as good urodynamic properties of the patch. However, the permanent contact of urine with the gastric mucosa is not free of complications. We report the urodynamic, macroscopic and histological outcomes of a pedicled de-epithelialized gastric patch incorporated in the bladder. We compared the results to those of our previous study, which sought to analyze these techniques of patch coverage using sigmoid patches. MATERIALS AND METHODS: We performed 20 augmentation cystoplasties in the mini-pig model using a pedicled de-epithelialized gastric patch and 5 techniques of patch coverage. RESULTS: Three months after surgery all bladders had an increase in volume except those in which the auto-augmentation technique was used. However, all gastric patches were smaller compared to preoperative size. Many had irregular fibrosed inner surfaces and histological evaluation revealed a fibrosed newly formed submucosal layer with a complete urothelial coverage in every patch. No gastric mucosal remnant was found. CONCLUSIONS: De-epithelialized gastrocystoplasty is an attractive procedure that can increase bladder capacity as well as provide a complete urothelial lining without mucosal remnants. However, the success of this procedure seems to be limited by increased morbidity and fibrotic changes, and decreased surface of the patch.

Synergistic epistasis and alternative hypotheses.

Inbreeding generally results in deleterious shifts in mean fitness. If the fitness response to increasing inbreeding coefficient is non-linear, this suggests a contribution of epistasis to inbreeding depression. In a cross-breeding experiment, Salathe & Ebert (2003. J. Evol. Biol. 16: 976-985) tested and found the presence of this non-linearity in Daphnia magna. They argue that epistatic interactions cause this non-linearity. We argue here that their experimental protocol does not allow disentangling the effect of synergistic epistasis from two alternative hypotheses, namely hybrid vigour and statistical non-independence of data.

Predicting smoking cessation and its relapse in HIV-infected patients: the Swiss HIV Cohort Study.

OBJECTIVES: The aim of the study was to assess whether prospective follow-up data within the Swiss HIV Cohort Study can be used to predict patients who stop smoking; or among smokers who stop, those who start smoking again. METHODS: We built prediction models first using clinical reasoning ('clinical models') and then by selecting from numerous candidate predictors using advanced statistical methods ('statistical models'). Our clinical models were based on literature that suggests that motivation drives smoking cessation, while dependence drives relapse in those attempting to stop. Our statistical models were based on automatic variable selection using additive logistic regression with component-wise gradient boosting. RESULTS: Of 4833 smokers, 26% stopped smoking, at least temporarily; because among those who stopped, 48% started smoking again. The predictive performance of our clinical and statistical models was modest. A basic clinical model for cessation, with patients classified into three motivational groups, was nearly as discriminatory as a constrained statistical model with just the most important predictors (the ratio of nonsmoking visits to total visits, alcohol or drug dependence, psychiatric comorbidities, recent hospitalization and age). A basic clinical model for relapse, based on the maximum number of cigarettes per day prior to stopping, was not as discriminatory as a constrained statistical model with just the ratio of nonsmoking visits to total visits. CONCLUSIONS: Predicting smoking cessation and relapse is difficult, so that simple models are nearly as discriminatory as complex ones. Patients with a history of attempting to stop and those known to have stopped recently are the best candidates for an intervention.

Lipid composition of fingermark residue and donor classification using GC/MS

Lipids available in fingermark residue represent important targets for enhancement and dating techniques. While it is well known that lipid composition varies among fingermarks of the same donor (intra-variability) and between fingermarks of different donors (inter-variability), the extent of this variability remains uncharacterised. Thus, this worked aimed at studying qualitatively and quantitatively the initial lipid composition of fingermark residue of 25 different donors. Among the 104 detected lipids, 43 were reported for the first time in the literature. Furthermore, palmitic acid, squalene, cholesterol, myristyl myristate and myristyl myristoleate were quantified and their correlation within fingermark residue was highlighted. Ten compounds were then selected and further studied as potential targets for dating or enhancement techniques. It was shown that their relative standard deviation was significantly lower for the intra-variability than for the inter-variability. Moreover, the use of data pretreatments could significantly reduce this variability. Based on these observations, an objective donor classification model was proposed. Hierarchical cluster analysis was conducted on the pre-treated data and the fingermarks of the 25 donors were classified into two main groups, corresponding to "poor" and "rich" lipid donors. The robustness of this classification was tested using fingermark replicates of selected donors. 86% of these replicates were correctly classified, showing the potential of such a donor classification model for research purposes in order to select representative donors based on compounds of interest.