Quantitative genetics of sleep in inbred mice.

The timing and the organization of sleep architecture are mainly controlled by the circadian system, while sleep need and intensity are regulated by a homeostatic process. How independent these two systems are in regulating sleep is not well understood. In contrast to the impressive progress in the molecular genetics of circadian rhythms, little is known about the molecular basis of sleep. Nevertheless, as summarized here, phenotypic dissection of sleep into its most basic aspects can be used to identify both the single major genes and small effect quantitative trait loci involved. Although experimental models such as the mouse are more readily amenable to genetic analysis of sleep, similar approaches can be applied to humans.

SYSGENET: a meeting report from a new European network for systems genetics.

The first scientific meeting of the newly established European SYSGENET network took place at the Helmholtz Centre for Infection Research (HZI) in Braunschweig, April 7-9, 2010. About 50 researchers working in the field of systems genetics using mouse genetic reference populations (GRP) participated in the meeting and exchanged their results, phenotyping approaches, and data analysis tools for studying systems genetics. In addition, the future of GRP resources and phenotyping in Europe was discussed.

GWAS of human bitter taste perception identifies new loci and reveals additional complexity of bitter taste genetics.

Human perception of bitterness displays pronounced interindividual variation. This phenotypic variation is mirrored by equally pronounced genetic variation in the family of bitter taste receptor genes. To better understand the effects of common genetic variations on human bitter taste perception, we conducted a genome-wide association study on a discovery panel of 504 subjects and a validation panel of 104 subjects from the general population of São Paulo in Brazil. Correction for general taste-sensitivity allowed us to identify a SNP in the cluster of bitter taste receptors on chr12 (10.88- 11.24 Mb, build 36.1) significantly associated (best SNP: rs2708377, P = 5.31 × 10(-13), r(2) = 8.9%, β = -0.12, s.e. = 0.016) with the perceived bitterness of caffeine. This association overlaps with-but is statistically distinct from-the previously identified SNP rs10772420 influencing the perception of quinine bitterness that falls in the same bitter taste cluster. We replicated this association to quinine perception (P = 4.97 × 10(-37), r(2) = 23.2%, β = 0.25, s.e. = 0.020) and additionally found the effect of this genetic locus to be concentration specific with a strong impact on the perception of low, but no impact on the perception of high concentrations of quinine. Our study, thus, furthers our understanding of the complex genetic architecture of bitter taste perception.

Landscape genetics of the Alpine newt (Mesotriton alpestris) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values < or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values < or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Conservation genetics and population monitoring of the Alpine Salamander (Salamandra Atra) :

SUMMARYIn the context of the biodiversity crisis, amphibians are experiencing the most severe worldwide decline of all vertebrates and are in urgent need of better management. Efficient conservation strategies rely on sound knowledge of the species biology and of the genetic and demographic processes that might impair their welfare. Nonetheless, these processes are poorly understood in amphibians. Delineating population boundaries remains consequently problematic for these species, while it is of critical importance to define adequate management units for conservation. In this study, our attention focused on the alpine salamander (Salamandra atra), a species that deserves much interest in terms of both conservation biology and evolution. This endemic alpine species shows peculiar life-history traits (viviparity, reduced activity period, slow maturation) and has a slow population turnover, which might be problematic for its persistence in a changing environment. Due to its elusive behaviour (individuals spend most of their time underground and are unavailable for sampling), dynamic processes of gene and individuals were poorly understood for that species. Consequently, its conservation status could hardly be reliably assessed. Similarly the fire salamander (Salamandra salamandra) also poses special challenges for conservation, as no clear demarcation of geographical populations exists and dispersal patterns are poorly known. Through a phylogeographic analysis, we first studied the evolutionary history of the alpine salamander to better document the distribution of the genetic diversity along its geographical range. This study highlighted the presence of multiple divergent lineages in Italy together with a clear genetic divergence between populations from Northern and Dinaric Alps. These signs of cryptic genetic differentiation, which are not accounted for by the current taxonomy of the species, should not be neglected for further definition of conservation units. In addition, our data supported glacial survival of the species in northern peripheral glacial réfugia and nunataks, a pattern rarely documented for long-lived species. Then, we evaluated the level of gene flow between populations at the local scale and tested for asymmetries in male versus female dispersal using both field-based (mark-recapture) and genetic approaches. This study revealed high level of gene flow between populations, which stems mainly from male dispersal. This corroborated the idea that salamanders are much better dispersers than hitherto thought and provided a well- supported example of male-biased dispersal in amphibians. In a third step, based on a mark- recapture survey, we addressed the problem of sampling unavailability in alpine salamanders and evaluated its impact on two monitoring methods. We showed that about three quarters of individuals were unavailable for sampling during sampling sessions, a proportion that can vary with climatic conditions. If not taken into account, these complexities would result in false assumptions on population trends and misdirect conservation efforts. Finally, regarding the daunting task of delineating management units, our attention was drawn on the fire salamander. We conducted a local population genetic study that revealed high levels of gene flow among sampling sites. Management units for this species should consequently be large. Interestingly, despite the presence of several landscape features often reported to act as barriers, genetic breaks occurred at unexpected places. This suggests that landscape features may rather have idiosyncratic effects on population structure. In conclusion, this work brought new insights on both genetic and demographic processes occurring in salamanders. The results suggest that some biological paradigms should be taken with caution when particular species are in focus. Species- specific studies remain thus fundamental for a better understanding of species evolution and conservation, particularly in the context of current global changes.RESUMEDans le contexte de la crise de la biodiversité actuelle, les amphibiens subissent le déclin le plus important de tous les vertébrés et ont urgemment besoin d'une meilleure protection. L'établissement de stratégies de conservation efficaces repose sur des connaissances solides de la biologie des espèces et des processus génétiques et démographiques pouvant menacer leur survie. Ces processus sont néanmoins encore peu étudiés chez les amphibiens.Dans cette étude, notre attention s'est portée sur la salamandre noire (Salamandra atra), une espèce endémique des Alpes dont les traits d'histoire de vie atypiques (viviparité, phase d'activité réduite, lent turnover des populations) pourraient la rendre très vulnérable face aux changements environnementaux. Par ailleurs, en raison de son comportement cryptique (les individus passent la plupart de leur temps sous terre) la dynamique des gènes et des individus est mal comprise chez cette espèce. Il est donc difficile d'évaluer son statut de conservation de manière fiable. La salamandre tachetée {Salamandra salamandra), pour qui il n'existe aucune démarcation géographique apparente des populations, pose également des problèmes en termes de gestion. Dans un premier temps, nous avons étudié l'histoire évolutive de la salamandre noire afin de mieux décrire la distribution de sa diversité génétique au sein de son aire géographique. Cela a permis de mettre en évidence la présence de multiples lignées en Italie, ainsi qu'une nette divergence entre les populations du nord des Alpes et des Alpes dinariques. Ces résultats seront à prendre en compte lorsqu'il s'agira de définir des unités de conservation pour cette espèce. D'autre part, nos données soutiennent l'hypothèse d'une survie glaciaire dans des refuges nordiques périglaciaires ou dans des nunataks, fait rarement documenté pour une espèce longévive. Nous avons ensuite évalué la différentiation génétique des populations à l'échelle locale, ce qui a révélé d'important flux de gènes, ainsi qu'une asymétrie de dispersion en faveur des mâles. Ces résultats corroborent l'idée que les amphibiens dispersent mieux que ce que l'on pensait, et fournissent un exemple robuste de dispersion biaisée en faveur des mâles chez les amphibiens. Nous avons ensuite abordé le problème de Γ inaccessibilité des individus à la capture. Nous avons montré qu'environ trois quarts des individus sont inaccessibles lors des échantillonnages, une proportion qui peut varier en fonction des conditions climatiques. Ignoré, ce processus pourrait entraîner une mauvaise interprétation des fluctuations de populations ainsi qu'une mauvaise allocation des efforts de conservation. Concernant la définition d'unités de gestion pour la salamandre tachetée, nous avons pu mettre en évidence un flux de gènes important entre les sites échantillonnés. Les unités de gestion pour cette espèce devraient donc être étendues. Etonnamment, malgré la présence de nombreuses barrières potentielles au flux de gènes, les démarcations génétiques sont apparues à des endroits inattendus. En conclusion, ce travail a apporté une meilleure compréhension des processus génétiques et démographiques en action chez les salamandres. Les résultats suggèrent que certains paradigmes biologiques devraient être considérés avec précaution quand il s'agit de les appliquer à des espèces particulières. Les études spécifiques demeurent donc fondamentales pour une meilleure compréhension de l'évolution des espèces et leur conservation, tout particulièrement dans le contexte des changements globaux actuels.

THE USE OF SIMULATIONS IN EVOLUTIONARY POPULATION GENETICS: APPLICATIONS ON HUMANS, OWLS AND VIRTUAL ORGANISMS

Summary (in English) Computer simulations provide a practical way to address scientific questions that would be otherwise intractable. In evolutionary biology, and in population genetics in particular, the investigation of evolutionary processes frequently involves the implementation of complex models, making simulations a particularly valuable tool in the area. In this thesis work, I explored three questions involving the geographical range expansion of populations, taking advantage of spatially explicit simulations coupled with approximate Bayesian computation. First, the neutral evolutionary history of the human spread around the world was investigated, leading to a surprisingly simple model: A straightforward diffusion process of migrations from east Africa throughout a world map with homogeneous landmasses replicated to very large extent the complex patterns observed in real human populations, suggesting a more continuous (as opposed to structured) view of the distribution of modern human genetic diversity, which may play a better role as a base model for further studies. Second, the postglacial evolution of the European barn owl, with the formation of a remarkable coat-color cline, was inspected with two rounds of simulations: (i) determine the demographic background history and (ii) test the probability of a phenotypic cline, like the one observed in the natural populations, to appear without natural selection. We verified that the modern barn owl population originated from a single Iberian refugium and that they formed their color cline, not due to neutral evolution, but with the necessary participation of selection. The third and last part of this thesis refers to a simulation-only study inspired by the barn owl case above. In this chapter, we showed that selection is, indeed, effective during range expansions and that it leaves a distinguished signature, which can then be used to detect and measure natural selection in range-expanding populations. Résumé (en français) Les simulations fournissent un moyen pratique pour répondre à des questions scientifiques qui seraient inabordable autrement. En génétique des populations, l'étude des processus évolutifs implique souvent la mise en oeuvre de modèles complexes, et les simulations sont un outil particulièrement précieux dans ce domaine. Dans cette thèse, j'ai exploré trois questions en utilisant des simulations spatialement explicites dans un cadre de calculs Bayésiens approximés (approximate Bayesian computation : ABC). Tout d'abord, l'histoire de la colonisation humaine mondiale et de l'évolution de parties neutres du génome a été étudiée grâce à un modèle étonnement simple. Un processus de diffusion des migrants de l'Afrique orientale à travers un monde avec des masses terrestres homogènes a reproduit, dans une très large mesure, les signatures génétiques complexes observées dans les populations humaines réelles. Un tel modèle continu (opposé à un modèle structuré en populations) pourrait être très utile comme modèle de base dans l'étude de génétique humaine à l'avenir. Deuxièmement, l'évolution postglaciaire d'un gradient de couleur chez l'Effraie des clocher (Tyto alba) Européenne, a été examiné avec deux séries de simulations pour : (i) déterminer l'histoire démographique de base et (ii) tester la probabilité qu'un gradient phénotypique, tel qu'observé dans les populations naturelles puisse apparaître sans sélection naturelle. Nous avons montré que la population actuelle des chouettes est sortie d'un unique refuge ibérique et que le gradient de couleur ne peux pas s'être formé de manière neutre (sans l'action de la sélection naturelle). La troisième partie de cette thèse se réfère à une étude par simulations inspirée par l'étude de l'Effraie. Dans ce dernier chapitre, nous avons montré que la sélection est, en effet, aussi efficace dans les cas d'expansion d'aire de distribution et qu'elle laisse une signature unique, qui peut être utilisée pour la détecter et estimer sa force.

Molecular genetics of "PRPF31" dominant mutations linked to retinitis pigmentosa

ABSTRACT : The retina is one of the most important human sensory tissues since it detects and transmits all visual information from the outside world to the brain. Retinitis pigmentosa (RP) is the name given to a group of inherited diseases that affect specifically the photoreceptors present in the retina and in many instances lead to blindness. Dominant mutations in PRPF31, a gene that encodes for a pre-mRNA splicing factor, cause retinitis pigmentosa with reduced penetrance. We functionally investigated a novel mutation, identified in a large family with autosomal dominant RP, and 7 other mutations, substitutions and microdeletions, in 12 patients from 7 families with PRPF31-linked RP. Seven mutations lead to PRPF31 mRNA with premature stop codons and one to mRNA lacking the exon containing the initiation codon. Quantification of PRPF31 mRNA and protein levels revealed a significant reduction in cell lines derived from patients, compared to non carriers of mutations in PRPF31. Allelic quantification of PRPF31 mRNA indicated that the level of mutated mRNA is very low compared to wild-type mRNA. No mutant protein was detected and the subnuclear localization of wild-type PRPF31 remains the same in cell lines from patients and controls. Blocking nonsense-mediated mRNA decay in cell lines derived from patients partially restored PRPF31 mutated mRNA but derived proteins were still undetectable, even when protein degradation pathways were inhibited. Our results demonstrated that the vast majority of PRPF31 mutations result in null alleles, since they are subject to surveillance mechanisms that degrade mutated mRNA and possibly block its translation. Altogether, these data indicate that the likely cause of PRPF31-linked RP is haploinsufficiency, rather than a dominant negative effect. Penetrance of PRPF31 mutations has been previously demonstrated to be inversely correlated with the level of PRPF31 mRNA, since high expression of wild-type PRPF31 mRNA protects from the disease. Consequently, we have investigated the genetic modifiers that control the expression of PRPF31 by quantifying PRPF31 mRNA levels in cell lines derived from 200 individuals from 15 families representative of the general population. By linkage analyses we identified a 8.2Mb-region on chromosome 14q21-23 that contains a gene involved in the modulation of PRPF31 expression. We also assessed apreviously-mapped penetrance factor invariably located on the wild-type allele and linked to the PRPF31 locus in asymptomatic patients from different families with RP. We demonstrated that this modifier increases the expression of both PRPF31 alleles already at the pre-mRNA level. Finally, our data suggest that PRPF31 mRNA expression and consequently the penetrance of PRPF31 mutations is modulated by at least 2 diffusible compounds, which act on both PRPF31 alleles during their transcription.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Narcolepsy and familial advanced sleep-phase syndrome: molecular genetics of sleep disorders.

Sleep disorders are very prevalent and represent an emerging worldwide epidemic. However, research into the molecular genetics of sleep disorders remains surprisingly one of the least active fields. Nevertheless, rapid progress is being made in several prototypical disorders, leading recently to the identification of the molecular pathways underlying narcolepsy and familial advanced sleep-phase syndrome. Since the first reports of spontaneous and induced loss-of-function mutations leading to hypocretin deficiency in human and animal models of narcolepsy, the role of this novel neurotransmission pathway in sleep and several other behaviors has gained extensive interest. Also, very recent studies using an animal model of familial advanced sleep-phase syndrome shed new light on the regulation of circadian rhythms.

A modular approach for integrative analysis of large-scale gene-expression and drug-response data

High-throughput technologies are now used to generate more than one type of data from the same biological samples. To properly integrate such data, we propose using co-modules, which describe coherent patterns across paired data sets, and conceive several modular methods for their identification. We first test these methods using in silico data, demonstrating that the integrative scheme of our Ping-Pong Algorithm uncovers drug-gene associations more accurately when considering noisy or complex data. Second, we provide an extensive comparative study using the gene-expression and drug-response data from the NCI-60 cell lines. Using information from the DrugBank and the Connectivity Map databases we show that the Ping-Pong Algorithm predicts drug-gene associations significantly better than other methods. Co-modules provide insights into possible mechanisms of action for a wide range of drugs and suggest new targets for therapy

Extreme Precipitation Modelling Using Geostatistics and Machine Learning Algorithms

The paper presents an approach for mapping of precipitation data. The main goal is to perform spatial predictions and simulations of precipitation fields using geostatistical methods (ordinary kriging, kriging with external drift) as well as machine learning algorithms (neural networks). More practically, the objective is to reproduce simultaneously both the spatial patterns and the extreme values. This objective is best reached by models integrating geostatistics and machine learning algorithms. To demonstrate how such models work, two case studies have been considered: first, a 2-day accumulation of heavy precipitation and second, a 6-day accumulation of extreme orographic precipitation. The first example is used to compare the performance of two optimization algorithms (conjugate gradients and Levenberg-Marquardt) of a neural network for the reproduction of extreme values. Hybrid models, which combine geostatistical and machine learning algorithms, are also treated in this context. The second dataset is used to analyze the contribution of radar Doppler imagery when used as external drift or as input in the models (kriging with external drift and neural networks). Model assessment is carried out by comparing independent validation errors as well as analyzing data patterns.

Genetics of sleep.

Molecular and genetic approaches in several species have provided new insights into the mechanisms of rest-activity and sleep-wake regulation. Many of these discoveries are believed to support hypotheses about sleep functions, which nevertheless remain elusive. In this review we discuss the specific contribution of both mammalian and invertebrate models to our understanding of the molecular basis of sleep.