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.

Linking genetic and species diversities in a community of freshwater gastropods

Summary Biodiversity is usually studied through species or genetic diversities. To date, these two levels of diversity have remained the independent .fields of investigations of community ecologists and population geneticists. However, recent joint analyses of species and genetic diversities have suggested that common processes may underlie the two levels. Positive correlations between species diversity and genetic diversity may arise when the effects of drift and migration overwhelm selective effects. The first goal of this thesis was to make a joint investigation of the patterns of species and genetic diversity in a community of freshwater gastropods living in a floodplain habitat. The second goal was to determine, as far as possible, the relative influences of the processes underlying the patterns observed at each level. In chapter 2 we investigate the relative influences of different evolutionary forces in shaping the genetic structure of Radix balthica populations. Results revealed that the structure inferred using quantitative traits was lower or equal to the one inferred using neutral molecular markers. Consequently, the pattern of structure observed could be only due to random drift, possibly to uniform selection, but definitely not to selection for local optima. In chapter 3, we analyze the temporal variation of species and genetic diversities in five localities. An extended period of drought occurred at the end of the study period leading to decay of both species and genetic diversities. This parallel loss of diversity following a natural perturbation highlighted the role sometimes predominant of random drift over selection on patterns of biodiversity in a floodplain habitat. In chapter 4, we compare the spatial genetic structures of two sympatric species: Radix balthica and Planorbis carinatus. We found that R. balthica populations are weakly structured and have moderate to high values of gene diversity. In contrast, P. carinatus populations are highly structured and poorly diverse. Then we measured correlations between various indices of species and genetic diversity using genetic data .from the two species. We found only one significant correlation: between species richness and gene diversity of P. carinatus. This result highlights the .need to use genetic date from more than one species to infer correlations between species and genetic diversities. Overall, this thesis provided new insights into the common processes underlying patterns of species and genetic diversity. Résumé La biodiversité est généralement étudiée au niveau de la diversité génétique ou spécifique. Ces deux niveaux sont restés jusqu'à maintenant les domaines d'investigation séparés des généticiens des populations et des écologistes des communautés. Cependant, des analyses conjointes des diversités génétique et spécifique ont récemment suggéré que des processus similaires pouvaient influencer ces deux niveaux. Des corrélations positives entre les diversités génétique et spécifique pourraient être dues aux effets de migration et de dérive qui dominent les effets sélectifs. Le premier but de cette thèse était de faire une étude conjointe des diversités génétique et spécifique dans une communauté de gastéropodes d'eau douce. Le second objectif était de déterminer les influences relatives des différents processus liés à chaque niveau de diversité. Dans le chapitre 2 nous cherchons à déterminer quelles forces évolutives influencent la structure génétique de quatre populations de Radix balthica. La structure mesurée sur des traits quantitatifs s'est révélée être plus faible ou égale à celle mesurée avec des marqueurs moléculaires neutres. La structure observée pourrait ainsi être due uniquement à la dérive génétique, potentiellement à la sélection uniforme, mais en aucun cas à la sélection locale pour différents optima. Dans le chapitre 3 nous analysons la variation temporelle des diversités génétique et spécifique dans cinq localités. Une récente période de sécheresse a causé une diminution parallèle des deux niveaux de diversité. Cette perturbation à mis en évidence le rôle parfois prépondérant de la dérive par rapport à celui de la sélection dans le déterminisme de la biodiversité dans un écosytème alluvial. Dans le chapitre 4, nous comparons la structure génétique spatiale de deux espèces vivant en sympatrie : Radix balthica et Planorbis carinatus. Les populations de R. balthica sont peu structurées et présentent un niveau de diversité relativement élevé alors que celles de P. carinatus sont fortement structurées et peu diversifiées. Nous avons ensuite mesuré différentes corrélations entre les diversités génétique et spécifique, mais la seule relation significative a été trouvée entre la richesse spécifique et la diversité génétique de P. carinatus. Ainsi, cette thèse a permis de découvrir de nouveaux aspects des processus qui influencent en parallèle la diversité aux niveaux génétique et spécifique.

Characterization of nontypable haemophilus influenzae isolates recovered from adult patients with underlying chronic lung disease reveals genotypic and phenotypic traits associated with persistent infection

Nontypable Haemophilus influenzae (NTHi) has emerged as an important opportunistic pathogen causing infection in adults suffering obstructive lung diseases. Existing evidence associates chronic infection by NTHi to the progression of the chronic respiratory disease, but specific features of NTHi associated with persistence have not been comprehensively addressed. To provide clues about adaptive strategies adopted by NTHi during persistent infection, we compared sequential persistent isolates with newly acquired isolates in sputa from six patients with chronic obstructive lung disease. Pulse field gel electrophoresis (PFGE) identified three patients with consecutive persistent strains and three with new strains. Phenotypic characterisation included infection of respiratory epithelial cells, bacterial self-aggregation, biofilm formation and resistance to antimicrobial peptides (AMP). Persistent isolates differed from new strains in showing low epithelial adhesion and inability to form biofilms when grown under continuous-flow culture conditions in microfermenters. Self-aggregation clustered the strains by patient, not by persistence. Increasing resistance to AMPs was observed for each series of persistent isolates; this was not associated with lipooligosaccharide decoration with phosphorylcholine or with lipid A acylation. Variation was further analyzed for the series of three persistent isolates recovered from patient 1. These isolates displayed comparable growth rate, natural transformation frequency and murine pulmonary infection. Genome sequencing of these three isolates revealed sequential acquisition of single-nucleotide variants in the AMP permease sapC, the heme acquisition systems hgpB, hgpC, hup and hxuC, the 3-deoxy-D-manno-octulosonic acid kinase kdkA, the long-chain fatty acid transporter ompP1, and the phosphoribosylamine glycine ligase purD. Collectively, we frame a range of pathogenic traits and a repertoire of genetic variants in the context of persistent infection by NTHi.

Genetic and environmental influences on blood pressure and physical activity: a study of nuclear families from Muzambinho, Brazil

Blood pressure (BP) and physical activity (PA) levels are inversely associated. Since genetic factors account for the observed variation in each of these traits, it is possible that part of their association may be related to common genetic and/or environmental influences. Thus, this study was designed to estimate the genetic and environmental correlations of BP and PA phenotypes in nuclear families from Muzambinho, Brazil. Families including 236 offspring (6 to 24 years) and their 82 fathers and 122 mothers (24 to 65 years) were evaluated. BP was measured, and total PA (TPA) was assessed by an interview (commuting, occupational, leisure time, and school time PA). Quantitative genetic modeling was used to estimate maximal heritability (h²), and genetic and environmental correlations. Heritability was significant for all phenotypes (systolic BP: h² = 0.37 ± 0.10, P < 0.05; diastolic BP: h² = 0.39 ± 0.09, P < 0.05; TPA: h² = 0.24 ± 0.09, P < 0.05). Significant genetic (r g) and environmental (r e) correlations were detected between systolic and diastolic BP (r g = 0.67 ± 0.12 and r e = 0.48 ± 0.08, P < 0.05). Genetic correlations between BP and TPA were not significant, while a tendency to an environmental cross-trait correlation was found between diastolic BP and TPA (r e = -0.18 ± 0.09, P = 0.057). In conclusion, BP and PA are under genetic influences. Systolic and diastolic BP share common genes and environmental influences. Diastolic BP and TPA are probably under similar environmental influences.

Genotypic and phenotypic diversity of a baculovirus population within an individual insect host

It is becoming increasingly apparent that many pathogen populations, including those of insects, show high levels of genotypic variation. Baculoviruses are known to be highly variable, with isolates collected from the same species in different geographical locations frequently showing genetic variation and differences in their biology. More recent Studies at smaller scales have also shown that virus DNA profiles from individual larvae can show polymorphisms within and between populations of the same species. Here, we investigate the genotypic and phenotypic variation of an insect baculovirus infection within a single insect host. Twenty four genotypically distinct nucleopolyhedrovirus (NPV) variants were isolated from an individual pine beauty moth, Panolis flammea, caterpillar by in vivo cloning techniques. No variant appeared to be dominant in the population. The Pafl NPV variants have been mapped using three restriction endonucleases and shown to contain three hypervariable regions containing insertions of 70-750 bp. Comparison of seven of these variants in an alternative host, Mamestra brassicae, demonstrated that the variants differed significantly in both pathogenicity and speed of kill. The generation and maintenance of pathogen heterogeneity are discussed. (c) 2005 Elsevier Inc. All rights reserved.

Molecular and phenotypic characterization of the alternative seasonal growth habit and flowering time in barley (Hordeum vulgare ssp. vulgare L.)

Barley can be classified into three major agronomic types, based on its seasonal growth habit (SGH): spring, winter and alternative. Winter varieties require exposure to vernalization to promote subsequent flowering and are autumn-sown. Spring varieties proceed to flowering in the absence of vernalization and are sown in the spring. The ‘alternative’ (also known as ‘facultative’) SGH is only loosely defined and can be sown in autumn or spring. Here, we investigate the molecular genetic basis of alternative barley. Analysis of the major barley vernalization (VRN-H1, VRN-H2) and photoperiod (PPD-H1, PPD-H2) response genes in a collection of 386 varieties found alternative SGH to be characterized by specific allelic combinations. Spring varieties possessed spring loci at one or both of the vernalization response loci, combined with long-day non-responsive ppd-H1 alleles and wild-type alleles at the short-day photoperiod response locus, PPD-H2. Winter varieties possessed winter alleles at both vernalization loci, in combination with the mutant ppd-H2 allele conferring delayed flowering under short-day photoperiods. In contrast, all alternative varieties investigated possessed a single spring allele (either at VRN-H1 or at VRN-H2) combined with mutant ppd-H2 alleles. This allelic combination is found only in alternative types and is diagnostic for alternative SGH in the collection studied. Analysis of flowering time under controlled environment found alternative varieties flowered later than spring control lines, with the difference most pronounced under short-day photoperiods. This work provides genetic characterization of the alternative SGH phenotype, allowing precise manipulation of SGH and flowering time within breeding programmes, and provides the molecular tools for classification of all three SGH categories within national variety registration processes.

Genetic and environmental factors affecting ultrasound measures of longissimus muscle area and backfat thickness in Nelore cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genetic and environmental heterogeneity of residual variance of weight traits in Nellore beef cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic and environmental effects on sexual precocity traits in Nellore cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Genetic and environmental factors that affect the best time of Thoroughbred horses in Brazil

On the basis of information provided by the Brazilian Association of Race Horse Breeders, we analysed the racing performance of 947 Thoroughbred horses in races held from 1985 to 1992. The performance was evaluated using the best, time of the animals. The variance component was obtained by the derivative-free restricted maximum likelihood method, and the model used contained fixed effects of the racing month and year, sex, race track, track condition, animal age, number of competitors in race, and distance, and the random animal effect. The low heritability estimate obtained (0.12) indicates that selection based on animal phenotypic value must induce small genetic changes in this trail.

Virulence attenuation and phenotypic variation of Paracoccidioides brasiliensis isolates obtained from armadillos and patients

Paracoccidioides brasiliensis is the etiological agent of paracoccidioidomycosis, the most important systemic mycosis in Latin America. The virulence profiles of five isolates of P. brasiliensis were studied in two different moments and correlated with some colonial phenotypic aspects. We observed a significant decrease in the virulence and an intense phenotypic variation in the mycelial colony. The recognition of all ranges of phenotypic and virulence variation of P. brasiliensis, as well as its physiological and genetic basis, will be important for a better comprehension of its pathogenic and epidemiological features.

Historical biogeography and conservation of the golden-striped salamander (Chioglossa lusitanica) in northwestern Iberia: integrating ecological, phenotypic and phylogeographic data

The golden-striped salamander (Chioglossa lusitanica) is an endemic species inhabiting stream-side habitats in mountainous areas in the northwestern Iberian Peninsula. This salamandrid is listed in the IUCN Red Data Book as a threatened species. The combination of bioclimatic modeling of the species distribution and multivariate analysis of genetic and phenotypic data strengthens previous hypotheses concerning the historical biogeography of C. lusitanica: the Pleistocene subdivision of the species' range and a process of postglacial recolonization. Discrepancies between bioclimatic modeling predictions and the present-day distribution suggest that the species may still be expanding its range northwards. We propose the identification of two distinct units for the conservation of the species and suggest that this information should be taken into account in defining key areas for conservation in the Iberian Peninsula.

Univariate and bivariate distribution of growth traits in beef buffaloes from Brazil

The aim of this study was to analyze the weight at birth (BW) and adjusted at 205 (W205), 365 (W365) and 550 (W55O) days in beef buffaloes from Brazil, using two approaches: parametric, by normal distribution, and non-parametric, by kernel function, and thus estimating the genetic, environmental and phenotypic correlation among traits. Information of 5,169 animals at birth (BW), 3,792 at 205 days (W205), 3.883 at 365 days (W365) and 1,524 at 550 days of age (W550) were used. The birth weight distribution presented an evident discrepancy in relation to the normal distribution. However, W205, W365 and W550 presented normal distributions. The birth weight presented weak genetic, environmental, and phenotypic associations with the other weight measurements. on the other hand, the weight traits at 205, 365, 550 days of age showed a high genetic correlation.

Fungemia by yeasts in Brazil: occurrence and phenotypic study of strains isolated at the Public Hospital, Botucatu, São Paulo

Objective. - In this study strains of yeasts isolated from the blood of human patients were analyzed taxonomically, their virulence factors were determined and compared, and phenotypic markers were used to compare the samples with respect to phenotypic differences across the range of patients as well as between samples isotated from the same patient.Methods. - the study involved a total of 75 strains of yeast isolated from the blood of in-patients of the Public Hospital, Botucatu, S (a) over tildeo Paulo, Brazil, with a clinical profile of fungemia. The hospital wards with the largest number of fungemias were neonatal intensive care units (ICUs) (32%) followed by gastric surgery (13.4%) and pediatric wards (10.7%). After identification, the samples were analyzed for the production of phospholipase and proteinase enzymes, and biotyped according to their susceptibility to killer toxins.Results. - the most frequent species found was Candida albicans (38.7%) followed by C. parapsilosis (30.7%). In terms of enzyme production, 98.7% of the 75 samples of yeast presented a strongly positive activity for proteinase; however, 78.7% did not present any phospholipasic activity. Six different biotypes were identified, the most frequent being 511 and 888.Conclusion. In association with phenotypic methods, genetic analyses should also be made of the samples under study to help in the rational development of a wider range of preventive measures and better control of hospital-contracted infections. (c) 2005 Published by Elsevier SAS.

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18. mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate. © 2011 Elsevier Inc.