The genetic basis of sleep disorders.

The contribution of genes, environment and gene-environment interactions to sleep disorders is increasingly recognized. Well-documented familial and twin sleep disorder studies suggest an important influence of genetic factors. However, only few sleep disorders have an established genetic basis including four rare diseases that may result from a single gene mutation: fatal familial insomnia, familial advanced sleep-phase syndrome, chronic primary insomnia, and narcolepsy with cataplexy. However, most sleep disorders are complex in terms of their genetic susceptibility together with the variable expressivity of the phenotype even within a same family. Recent linkage, genome-wide and candidate gene association studies resulted in the identification of gene mutations, gene localizations, or evidence for susceptibility genes and/or loci in several sleep disorders. Molecular techniques including mainly genome-wide linkage and association studies are further required to identify the contribution of new genes. These identified susceptibility genetic determinants will provide clues to better understand pathogenesis of sleep disorders, to assess the risk for diseases and also to find new drug targets to treat and to prevent the underlying conditions. We reviewed here the role of genetic basis in most of key sleep disorders.

Microsatellite variation reveals low genetic subdivision in a chromosome race of Sorex araneus (Mammalia, Insectivora)

Two hundred and forty-five individuals of the common shrew (Sorex araneus, Insectivora, Mammalia) from 24 sampling localities situated in four different valleys of the western European Alps were genotyped for six microsatellite loci. Allelic variability ranged from 3 to 32 different alleles at a single locus and the average gene diversity over all loci was 0.69. An analysis for F and R statistics revealed weak genetic population subdivision (Fst = 0.032; Rst = 0.016). This suggests considerable gene flow and little phylogeographic structure within and between valleys. We tested whether a stepwise mutation model (SMM) better explained variation at the microsatellite loci than an infinite allele model (IAM). No trend in favor of either model was detected.

Genetic Variants of Serum Uric Acid and Gout: An Analysis of > 170,000 Individuals

Background/Purpose: Gout is a common and excruciatingly painful inflammatory arthritis caused by hyperuricemia. In addition to various lifestyle risk factors, a substantial genetic predisposition to gout has long been recognized. The Global Urate Genetics Consortium (GUGC) has aimed to comprehensively investigate the genetics of serum uric acid and gout using data from _ 140,000 individuals of European-ancestry, 8,340 individuals of Indian ancestry, 5,820 African-Americans, and 15,286 Japanese. Methods: We performed discovery GWAS meta-analyses of serum urate levels (n_110,347 individuals) followed by replication analyses (n_32,813 different individuals). Our gout analysis involved 3,151 cases and 68,350 controls, including 1,036 incident gout cases that met the American College of Rheumatology Criteria. We also examined the association of gout with fractional excretion of uric acid (n_6,799). A weighted genetic urate score was constructed based on the number of risk alleles across urate-associated loci, and their association with the risk of gout was evaluated. Furthermore, we examined implicated transcript expression in cis (expression quantitative trait loci databases) for potential insights into the gene underlying the association signal. Finally, in order to further identify urate-associated genomic regions, we performed functional network analyses that incorporated prior knowledge on molecular interactions in which the gene products of implicated genes operate. Results: We identified and replicated 28 genome-wide significant loci in association with serum urate (P 5_10_8), including all previously-reported loci as well as 18 novel genetic loci. Unlike the majority of previouslyidentified loci, none of the novel loci appeared to be obvious candidates for urate transport. Rather, they were mapped to genes that encode for purine production, transcription, or growth factors with broad downstream responses. Besides SLC2A9 and ABCG2, no additional regions contained SNPs that differed significantly (P _ 5_10_8) between sexes. Urateincreasing alleles were associated with an increased risk of gout for all loci. The urate genetic risk score (ranging from 10 to 45) was significantly associated with an increased odds of prevalent gout (OR per unit increase, 1.11; 95% CI, 1.09-1.14) and incident gout (OR, 1.10; 95% CI, 1.08-1.13). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. Detailed characterization of the loci revealed associations with transcript expression and the fractional excretion of urate. Network analyses implicated the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. Conclusion: The novel genetic candidates identified in this urate/gout consortium study, the largest to date, highlight the importance of metabolic control of urate production and urate excretion. The modulation by signaling processes that influence metabolic pathways such as glycolysis and the pentose phosphate pathway appear to be central mechanisms underpinned by the novel GWAS candidates. These findings may have implications for further research into urate-lowering drugs to treat and prevent gout.

Maintenance of ancestral sex chromosomes in Palearctic tree frogs: direct evidence from Hyla orientalis.

Contrasting with the situation found in birds and mammals, sex chromosomes are generally homomorphic in poikilothermic vertebrates. This homomorphy was recently shown to result from occasional X-Y recombinations (not from turnovers) in several European species of tree frogs (Hyla arborea, H. intermedia and H. molleri). Because of recombination, however, alleles at sex-linked loci were rarely diagnostic at the population level; support for sex linkage had to rely on multilocus associations, combined with occasional sex differences in allelic frequencies. Here, we use direct evidence, obtained from anatomical and histological analyses of offspring with known pedigrees, to show that the Eastern tree frog (H. orientalis) shares the same pair of sex chromosomes, with identical patterns of male heterogamety and complete absence of X-Y recombination in males. Conservation of an ancestral pair of sex chromosomes, regularly rejuvenated via occasional X-Y recombination, seems thus a widespread pattern among Hyla species. Sibship analyses also identified discrepancies between genotypic and phenotypic sex among offspring, associated with abnormal gonadal development, suggesting a role for sexually antagonistic genes on the sex chromosomes.

ON THE ROLE OF PERIOD-2 IN THE CIRCADIAN AND HOMEOSTATIC REGULATION OF SLEEP

Humans spend one third of their life sleeping, then we could raise the basic question: Why do we sleep? Despite the fact that we still don't fully understand its function, we made much progress in understanding at different levels how sleep is regulated. One model suggests that sleep is regulated by two processes: a homeostatic process that tracks the need for sleep and by a circadian rhythm that determines the preferred time-of-day sleep occurs. At the molecular level circadian rhythms are a property of interlocking transcriptional regula-tors referred to as clock genes. The heterodimeric transcription factors BMAL1::CLOCK/NPAS2 drive the transcription of many target genes including the clock genes Cryptochome1 (Cry1), Cry2, Period1 (Per1), and Per2. The encoded CRY/PER proteins are transcriptional inhibitors of BMAL1::CLOCK/NPAS2 thereby providing negative feedback to their own transcription. These genes seem, however, also involved in sleep homeostasis because the brain expression of clock genes, es-pecially that of Per2, increase as a function of time-spent-awake and because mice lacking clock genes display altered sleep homeostasis. The aim of first part of my doctoral work has been to advance our understanding the link that exists between sleep homeostasis and circadian rhythms investigating a possible mechanism by which sleep deprivation could alter clock gene expression by quantifying DNA-binding of the core-clock genes BMAL1, CLOCK and NPAS2 to their target chromatin loci including the E-box enhancers of the Per2 promoter. We made use of chromatin immunoprecipitation (ChIP) and quantitative poly-merase chain reaction (qPCR) to show that DNA-binding of CLOCK and BMAL1 to their target genes changes as a function of time-of-day in both liver and cerebral cortex. We then performed a 6h sleep deprivation (SD) and observed a significant decrease in DNA-binding of CLOCK and BMAL1 to Dbp. This is consistent with a decrease in Dbp mRNA levels after SD. The DNA-binding of NPAS2 and BMAL1 to Per2 was similarly decreased following SD. However, SD has been previously shown to in-crease Per2 expression in the cortex which seems paradoxical. Our results demonstrate that sleep-wake history can affect the molecular clock machinery directly at the level of the chromatin thereby altering the cortical expression of Dbp and Per2, and likely other targets. However, the precise dy-namic relationship between DNA-binding and mRNA expression, especially for Per2, remains elusive. The second aim of my doctoral work has been to perform an in depth characterization of cir-cadian rhythmicity, sleep architecture, analyze the response to SD in full null-Per2 knock-out (Per2-/-) mice, and Per1-/- mice, as well as their double knock-out offspring (Per1,2-/-) and littermate wildtype (Wt) mice. The techniques used include locomotor activity recording by passive infrared (PIR) sen-sors, EEG/EMG surgery, recording, and analysis, and cerebral cortex extraction and quantification of mRNA levels by qPCR. Under standard LD12:12 conditions, we found that wakefulness onset, as well as the time courses of clock gene expression in the brain and corticosterone plasma levels were ad-vanced by about 2h in Per2-/- mice compared to Wt mice. When released under constant dark condi-tions almost all Per2-/- mice (97%) became arrhythmic immediately. From these observations, we conclude that while Per2-/- mice seem to be able to anticipate dark onset, this does not result from a self-sustained circadian clock. Our results suggest instead that the earlier onset of activity results from a labile, not-self sustained 22h rhythm linked to light onset suggesting the existence of a light-driven rhythm. Analyses of sleep under LD12:12 conditions revealed that in both Per2-/- and Per1,2-/- mice the same sleep phenotypes are observed compared to Wt mice: increased NREM sleep frag-mentation and inability to adequately compensate the loss of NREM sleep. That suggests a possible role of PER2 in sleep consolidation and recovery.

The effects of dominance, regular inbreeding and sampling design on Q(ST), an estimator of population differentiation for quantitative traits.

To test whether quantitative traits are under directional or homogenizing selection, it is common practice to compare population differentiation estimates at molecular markers (F(ST)) and quantitative traits (Q(ST)). If the trait is neutral and its determinism is additive, then theory predicts that Q(ST) = F(ST), while Q(ST) > F(ST) is predicted under directional selection for different local optima, and Q(ST) < F(ST) is predicted under homogenizing selection. However, nonadditive effects can alter these predictions. Here, we investigate the influence of dominance on the relation between Q(ST) and F(ST) for neutral traits. Using analytical results and computer simulations, we show that dominance generally deflates Q(ST) relative to F(ST). Under inbreeding, the effect of dominance vanishes, and we show that for selfing species, a better estimate of Q(ST) is obtained from selfed families than from half-sib families. We also compare several sampling designs and find that it is always best to sample many populations (>20) with few families (five) rather than few populations with many families. Provided that estimates of Q(ST) are derived from individuals originating from many populations, we conclude that the pattern Q(ST) > F(ST), and hence the inference of directional selection for different local optima, is robust to the effect of nonadditive gene actions.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Limited heterogeneity of autoantigens and T cells in autoimmune diseases?

For many induced and spontaneous autoimmune diseases, a predominant role for T cells in the organ-specific destruction process has been shown. In one of the induced models of autoimmunity, experimental allergic encephalomyelitis (EAE), a very small heterogeneity of T-cell receptor (TcR) molecules is expressed by the pathogenic T cells in both rats and mice. Contrary to induced autoimmune diseases, little is known about the autoantigens recognized by these autoimmune T cells and the heterogeneity of their TcR in spontaneous autoimmune diseases. The aim of this work was to establish a system which allows characterization of relevant autoantigens in spontaneous insulin-dependent diabetes mellitus (IDDM) in non-obese diabetic (NOD) mice. A completely different approach was taken to characterize the gene products of the minor lymphocyte stimulatory (Mls) loci. These gene products are responsible for the clonal elimination or the clonal stimulation of T cells expressing particular TcR V beta genes and therefore could be implicated in induction of autoimmune diseases by oligoclonal T-cell populations. The finding that Mls antigens are encoded by retroviral sequences leads to the hypothesis that viruses could be the inducing agents of autoimmune diseases.

Molecular Proxies for Climate Maladaptation in a Long-Lived Tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Common variants in the ATP2B1 gene are associated with susceptibility to hypertension: the Japanese Millennium Genome Project.

Hypertension is one of the most common complex genetic disorders. We have described previously 38 single nucleotide polymorphisms (SNPs) with suggestive association with hypertension in Japanese individuals. In this study we extend our previous findings by analyzing a large sample of Japanese individuals (n=14 105) for the most associated SNPs. We also conducted replication analyses in Japanese of susceptibility loci for hypertension identified recently from genome-wide association studies of European ancestries. Association analysis revealed significant association of the ATP2B1 rs2070759 polymorphism with hypertension (P=5.3×10(-5); allelic odds ratio: 1.17 [95% CI: 1.09 to 1.26]). Additional SNPs in ATP2B1 were subsequently genotyped, and the most significant association was with rs11105378 (odds ratio: 1.31 [95% CI: 1.21 to 1.42]; P=4.1×10(-11)). Association of rs11105378 with hypertension was cross-validated by replication analysis with the Global Blood Pressure Genetics consortium data set (odds ratio: 1.13 [95% CI: 1.05 to 1.21]; P=5.9×10(-4)). Mean adjusted systolic blood pressure was highly significantly associated with the same SNP in a meta-analysis with individuals of European descent (P=1.4×10(-18)). ATP2B1 mRNA expression levels in umbilical artery smooth muscle cells were found to be significantly different among rs11105378 genotypes. Seven SNPs discovered in published genome-wide association studies were also genotyped in the Japanese population. In the combined analysis with replicated 3 genes, FGF5 rs1458038, CYP17A1, rs1004467, and CSK rs1378942, odds ratio of the highest risk group was 2.27 (95% CI: 1.65 to 3.12; P=4.6×10(-7)) compared with the lower risk group. In summary, this study confirmed common genetic variation in ATP2B1, as well as FGF5, CYP17A1, and CSK, to be associated with blood pressure levels and risk of hypertension.

Estimating sex-specific dispersal rates with autosomal markers in hierarchically structured populations.

A recent study suggests that sex-specific dispersal rates can be quantitatively estimated on the basis of sex- and state-specific (pre- vs. postdispersal) F-statistics. In the present paper, we extend this approach to account for the hierarchical structure of natural populations, and we validate it through individual-based simulations. The model is applied to an empirical data set consisting of 536 individuals (males, females, and predispersal juveniles) of greater white-toothed shrews (Crocidura russula), sampled according to a hierarchical design and typed for seven autosomal microsatellite loci. From this dataset, dispersal is significantly female biased at the local scale (breeding-group level), but not at the larger scale (among local populations). We argue that selective pressures on dispersal are likely to depend on the spatial scale considered, and that short-distance dispersal should mainly respond to kin interactions (inbreeding or kin competition avoidance), which exert differential pressure on males and females.

From the Apennines to the Alps: colonization genetics of the naturally expanding Italian wolf (Canis lupus) population

Wolves in Italy strongly declined in the past and were confined south of the Alps since the turn of the last century, reduced in the 1970s to approximately 100 individuals surviving in two fragmented subpopulations in the central-southern Apennines. The Italian wolves are presently expanding in the Apennines, and started to recolonize the western Alps in Italy, France and Switzerland about 16 years ago. In this study, we used a population genetic approach to elucidate some aspects of the wolf recolonization process. DNA extracted from 3068 tissue and scat samples collected in the Apennines (the source populations) and in the Alps (the colony), were genotyped at 12 microsatellite loci aiming to assess (i) the strength of the bottleneck and founder effects during the onset of colonization; (ii) the rates of gene flow between source and colony; and (iii) the minimum number of colonizers that are needed to explain the genetic variability observed in the colony. We identified a total of 435 distinct wolf genotypes, which showed that wolves in the Alps: (i) have significantly lower genetic diversity (heterozygosity, allelic richness, number of private alleles) than wolves in the Apennines; (ii) are genetically distinct using pairwise F(ST) values, population assignment test and Bayesian clustering; (iii) are not in genetic equilibrium (significant bottleneck test). Spatial autocorrelations are significant among samples separated up to c. 230 km, roughly correspondent to the apparent gap in permanent wolf presence between the Alps and north Apennines. The estimated number of first-generation migrants indicates that migration has been unidirectional and male-biased, from the Apennines to the Alps, and that wolves in southern Italy did not contribute to the Alpine population. These results suggest that: (i) the Alps were colonized by a few long-range migrating wolves originating in the north Apennine subpopulation; (ii) during the colonization process there has been a moderate bottleneck; and (iii) gene flow between sources and colonies was moderate (corresponding to 1.25-2.50 wolves per generation), despite high potential for dispersal. Bottleneck simulations showed that a total of c. 8-16 effective founders are needed to explain the genetic diversity observed in the Alps. Levels of genetic diversity in the expanding Alpine wolf population, and the permanence of genetic structuring, will depend on the future rates of gene flow among distinct wolf subpopulation fragments.

In vitro whole-genome analysis identifies a susceptibility locus for HIV-1.

Advances in large-scale analysis of human genomic variability provide unprecedented opportunities to study the genetic basis of susceptibility to infectious agents. We report here the use of an in vitro system for the identification of a locus on HSA8q24.3 associated with cellular susceptibility to HIV-1. This locus was mapped through quantitative linkage analysis using cell lines from multigeneration families, validated in vitro, and followed up by two independent association studies in HIV-positive individuals. Single nucleotide polymorphism rs2572886, which is associated with cellular susceptibility to HIV-1 in lymphoblastoid B cells and in primary T cells, was also associated with accelerated disease progression in one of two cohorts of HIV-1-infected patients. Biological analysis suggests a role of the rs2572886 region in the regulation of the LY6 family of glycosyl-phosphatidyl-inositol (GPI)-anchored proteins. Genetic analysis of in vitro cellular phenotypes provides an attractive approach for the discovery of susceptibility loci to infectious agents.

Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions.

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.