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.

A random linear functional approach to efficiency bounds

This paper suggests a method for obtaining efficiency bounds in models containing either only infinite-dimensional parameters or both finite- and infinite-dimensional parameters (semiparametric models). The method is based on a theory of random linear functionals applied to the gradient of the log-likelihood functional and is illustrated by computing the lower bound for Cox's regression model

Random regression models to estimate genetic parameters for milk production of Guzerat cows using orthogonal Legendre polynomials

The objective of this work was to compare random regression models for the estimation of genetic parameters for Guzerat milk production, using orthogonal Legendre polynomials. Records (20,524) of test-day milk yield (TDMY) from 2,816 first-lactation Guzerat cows were used. TDMY grouped into 10-monthly classes were analyzed for additive genetic effect and for environmental and residual permanent effects (random effects), whereas the contemporary group, calving age (linear and quadratic effects) and mean lactation curve were analized as fixed effects. Trajectories for the additive genetic and permanent environmental effects were modeled by means of a covariance function employing orthogonal Legendre polynomials ranging from the second to the fifth order. Residual variances were considered in one, four, six, or ten variance classes. The best model had six residual variance classes. The heritability estimates for the TDMY records varied from 0.19 to 0.32. The random regression model that used a second-order Legendre polynomial for the additive genetic effect, and a fifth-order polynomial for the permanent environmental effect is adequate for comparison by the main employed criteria. The model with a second-order Legendre polynomial for the additive genetic effect, and that with a fourth-order for the permanent environmental effect could also be employed in these analyses.

Self-avoiding random walks and Olbers' paradox

In this paper, we prove that a self-avoiding walk of infinite length provides a structure that would resolve Olbers' paradox. That is, if the stars of a universe were distributed like the vertices of an infinite random walk with each segment length of about a parsec, then the night sky could be as dark as actually observed on the Earth. Self-avoiding random walk structure can therefore resolve the Olbers' paradox even in a static universe.

Moderate hypercapnia exerts beneficial effects on splanchnic energy metabolism during endotoxemia.

PURPOSE: Low tidal volume ventilation and permissive hypercapnia are required in patients with sepsis complicated by ARDS. The effects of hypercapnia on tissue oxidative metabolism in this setting are unknown. We therefore determined the effects of moderate hypercapnia on markers of systemic and splanchnic oxidative metabolism in an animal model of endotoxemia. METHODS: Anesthetized rats maintained at a PaCO(2) of 30, 40 or 60 mmHg were challenged with endotoxin. A control group (PaCO(2) 40 mmHg) received isotonic saline. Hemodynamic variables, arterial lactate, pyruvate, and ketone bodies were measured at baseline and after 4 h. Tissue adenosine triphosphate (ATP) and lactate were measured in the small intestine and the liver after 4 h. RESULTS: Endotoxin resulted in low cardiac output, increased lactate/pyruvate ratio and decreased ketone body ratio. These changes were not influenced by hypercapnia, but were more severe with hypocapnia. In the liver, ATP decreased and lactate increased independently from PaCO(2) after endotoxin. In contrast, the drop of ATP and the rise in lactate triggered by endotoxin in the intestine were prevented by hypercapnia. CONCLUSIONS: During endotoxemia in rats, moderate hypercapnia prevents the deterioration of tissue energetics in the intestine.

Predicting Non Return to Work after Orthopaedic Trauma: The Wallis Occupational Rehabilitation RisK (WORRK) Model.

BACKGROUND: Workers with persistent disabilities after orthopaedic trauma may need occupational rehabilitation. Despite various risk profiles for non-return-to-work (non-RTW), there is no available predictive model. Moreover, injured workers may have various origins (immigrant workers), which may either affect their return to work or their eligibility for research purposes. The aim of this study was to develop and validate a predictive model that estimates the likelihood of non-RTW after occupational rehabilitation using predictors which do not rely on the worker's background. METHODS: Prospective cohort study (3177 participants, native (51%) and immigrant workers (49%)) with two samples: a) Development sample with patients from 2004 to 2007 with Full and Reduced Models, b) External validation of the Reduced Model with patients from 2008 to March 2010. We collected patients' data and biopsychosocial complexity with an observer rated interview (INTERMED). Non-RTW was assessed two years after discharge from the rehabilitation. Discrimination was assessed by the area under the receiver operating curve (AUC) and calibration was evaluated with a calibration plot. The model was reduced with random forests. RESULTS: At 2 years, the non-RTW status was known for 2462 patients (77.5% of the total sample). The prevalence of non-RTW was 50%. The full model (36 items) and the reduced model (19 items) had acceptable discrimination performance (AUC 0.75, 95% CI 0.72 to 0.78 and 0.74, 95% CI 0.71 to 0.76, respectively) and good calibration. For the validation model, the discrimination performance was acceptable (AUC 0.73; 95% CI 0.70 to 0.77) and calibration was also adequate. CONCLUSIONS: Non-RTW may be predicted with a simple model constructed with variables independent of the patient's education and language fluency. This model is useful for all kinds of trauma in order to adjust for case mix and it is applicable to vulnerable populations like immigrant workers.

Determinants of the price of allowance in the European Union emissions trading market

Teorian mukaan täydellisen kilpailun päästöoikeuskauppamarkkinoilla päästöoikeuden hinta muodostuu markkinoilla vallitsevan päästöjen vähentämisen rajakustannuksen perusteella. Euroopan päästökauppamarkkinoilla päästöjen vähentämisen kustannuksia nostavat suhteellisen lyhyet päästökauppajaksot ja epävarmuus järjestelmän jatkuvuudesta. Toisaalta päästökaupan osallistujien yhteenlaskettu päästöjen vähentämisen tarve lienee suhteellisen vähäinen ellei olematon ensimmäisellä päästökauppajaksolla. Euroopan päästökauppamarkkinoilla päästöjen vähentämisen tarve ja päästöjenvähentämisen kustannukset ovat osittain riippuvaisia muuttuvista tekijöistä. Päästöoikeuden hintaan voivat vaikuttaa päästökauppajakson aikana tapahtuva teollisuuden suhdannevaihtelu, polttoaineiden hintojen heilahtelut sekä säätilojen vaihtelu. Päästökaupan ensimmäisinä kuukausina päästöoikeuden hintakehityksellä on ollut yhteyksiä tekijöihin, joiden muutosten tulisikin vaikuttaa päästökauppamarkkinoiden tasapainoon. Näitä tekijöitä ovat esimerkiksi polttoainemarkkinoiden ja sähkömarkkinoiden hintakehitys sekä vaihtelut säätiloissa.

Clinical chemistry variables in normal elderly and healthy ambulatory populations: comparison with reference values.

Laboratory values of the most commonly assayed clinical chemistry variables were determined in selected elderly and healthy ambulatory populations. The upper and lower limits (2.5 and 97.5 fractiles) were compared with the adult reference values in use in university hospitals of Switzerland. The results suggest that conventional adult reference values can be used for most variables in the elderly and that these values are also useful in an ambulatory population.

The Persistence of Finnish Equity Fund Performance

Tutkimuksen tavoitteena on selvittää, esiintyykö suomeen sijoittavilla osakerahastoilla menestyksen pysyvyyttä. Tutkimusaineisto koostuu kaikista suomalaisista osakerahastoista, jotka toimivat ajanjaksolla 15.1.1998-13.1.2005. Aineisto on vapaa selviytymisvinoumasta. Suorituskyvyn mittareina käytetään CAPM-alfaa sekä kolmi- ja nelifaktori-alfaa. Empiirisessä osassa osakerahastojen menestyksen pysyvyyttä testataan Spearmanin järjestyskorrelaatiotestillä. Evidenssi menestyksen pysyvyydestä jäi vähäiseksi, vaikkakin sitä esiintyi satunnaisesti kaikilla menestysmittareilla joillakin ranking- ja sijoitusperiodin yhdistelmillä. CAPM-alfalla tarkasteltuna tilastollisesti merkitsevää menestyksen pysyvyyttä esiintyi selvästi useammin kuin muilla menestysmittareilla. Tulokset tukevat viimeaikaisia kansainvälisiä tutkimuksia, joiden mukaan menestyksen pysyvyys riippuu usein mittaustavasta. Menestysmittareina käytettyjen regressiomallien merkitsevyystestit osoittavat multifaktorimallien selittävän osakerahastojen tuottoja CAPM:a paremmin. Lisätyt muuttujat parantavat merkittävästi CAPM:n selitysvoimaa.

Definition of key variables for the induction of optimal NY-ESO-1-specific T cells in HLA transgene mice.

An attractive treatment of cancer consists in inducing tumor-eradicating CD8(+) CTL specific for tumor-associated Ags, such as NY-ESO-1 (ESO), a strongly immunogenic cancer germ line gene-encoded tumor-associated Ag, widely expressed on diverse tumors. To establish optimal priming of ESO-specific CTL and to define critical vaccine variables and mechanisms, we used HLA-A2/DR1 H-2(-/-) transgenic mice and sequential immunization with immunodominant DR1- and A2-restricted ESO peptides. Immunization of mice first with the DR1-restricted ESO(123-137) peptide and subsequently with mature dendritic cells (DCs) presenting this and the A2-restriced ESO(157-165) epitope generated abundant, circulating, high-avidity primary and memory CD8(+) T cells that efficiently killed A2/ESO(157-165)(+) tumor cells. This prime boost regimen was superior to other vaccine regimes and required strong Th1 cell responses, copresentation of MHC class I and MHC class II peptides by the same DC, and resulted in upregulation of sphingosine 1-phosphate receptor 1, and thus egress of freshly primed CD8(+) T cells from the draining lymph nodes into circulation. This well-defined system allowed detailed mechanistic analysis, which revealed that 1) the Th1 cytokines IFN-gamma and IL-2 played key roles in CTL priming, namely by upregulating on naive CD8(+) T cells the chemokine receptor CCR5; 2) the inflammatory chemokines CCL4 (MIP-1beta) and CCL3 (MIP-1alpha) chemoattracted primed CD4(+) T cells to mature DCs and activated, naive CD8(+) T cells to DC-CD4 conjugates, respectively; and 3) blockade of these chemokines or their common receptor CCR5 ablated priming of CD8(+) T cells and upregulation of sphingosine 1-phosphate receptor 1. These findings provide new opportunities for improving T cell cancer vaccines.

Three essays on the determinants of output, inflation and interest rates

Résumé: Output, inflation and interest rates are key macroeconomic variables, in particular for monetary policy. In modern macroeconomic models they are driven by random shocks which feed through the economy in various ways. Models differ in the nature of shocks and their transmission mechanisms. This is the common theme underlying the three essays of this thesis. Each essay takes a different perspective on the subject: First, the thesis shows empirically how different shocks lead to different behavior of interest rates over the business cycle. For commonly analyzed shocks (technology and monetary policy errors), the patterns square with standard models. The big unknown are sources of inflation persistence. Then the thesis presents a theory of monetary policy, when the central bank can better observe structural shocks than the public. The public will then seek to infer the bank's extra knowledge from its policy actions and expectation management becomes a key factor of optimal policy. In a simple New Keynesian model, monetary policy becomes more concerned with inflation persistence than otherwise. Finally, the thesis points to the huge uncertainties involved in estimating the responses to structural shocks with permanent effects.

Determinació de l’estat de seguretat alimentària de dues zones d’Etiòpia i de les variables socioeconòmiques que les afecten

L’objectiu d’aquest estudi ha sigut generar dades especifiques locals de dos zones agroecològiques diferents d’Etiòpia ( el woreda Wukro i el woreda Abaya), mesurar l’estat de seguretat alimentària de la seva població, i els determinants socioeconòmics que els afecten, així com identificar els principals problemes metodològics, interpretatius que afecten a aquestes mesures i, en conseqüència, les intervencions que se’n poden derivar.