Combined genome-wide linkage and targeted association analysis of head circumference in autism spectrum disorder families.

BACKGROUND: It has long been recognized that there is an association between enlarged head circumference (HC) and autism spectrum disorder (ASD), but the genetics of HC in ASD is not well understood. In order to investigate the genetic underpinning of HC in ASD, we undertook a genome-wide linkage study of HC followed by linkage signal targeted association among a sample of 67 extended pedigrees with ASD. METHODS: HC measurements on members of 67 multiplex ASD extended pedigrees were used as a quantitative trait in a genome-wide linkage analysis. The Illumina 6K SNP linkage panel was used, and analyses were carried out using the SOLAR implemented variance components model. Loci identified in this way formed the target for subsequent association analysis using the Illumina OmniExpress chip and imputed genotypes. A modification of the qTDT was used as implemented in SOLAR. RESULTS: We identified a linkage signal spanning 6p21.31 to 6p22.2 (maximum LOD = 3.4). Although targeted association did not find evidence of association with any SNP overall, in one family with the strongest evidence of linkage, there was evidence for association (rs17586672, p = 1.72E-07). CONCLUSIONS: Although this region does not overlap with ASD linkage signals in these same samples, it has been associated with other psychiatric risk, including ADHD, developmental dyslexia, schizophrenia, specific language impairment, and juvenile bipolar disorder. The genome-wide significant linkage signal represents the first reported observation of a potential quantitative trait locus for HC in ASD and may be relevant in the context of complex multivariate risk likely leading to ASD.

Voluntary wheel running reduces voluntary consumption of ethanol in mice: identification of candidate genes through striatal gene expression profiling.

Hedonic substitution, where wheel running reduces voluntary ethanol consumption, has been observed in prior studies. Here, we replicate and expand on previous work showing that mice decrease voluntary ethanol consumption and preference when given access to a running wheel. While earlier work has been limited mainly to behavioral studies, here we assess the underlying molecular mechanisms that may account for this interaction. From four groups of female C57BL/6J mice (control, access to two-bottle choice ethanol, access to a running wheel, and access to both two-bottle choice ethanol and a running wheel), mRNA-sequencing of the striatum identified differential gene expression. Many genes in ethanol preference quantitative trait loci were differentially expressed due to running. Furthermore, we conducted Weighted Gene Co-expression Network Analysis and identified gene networks corresponding to each effect behavioral group. Candidate genes for mediating the behavioral interaction between ethanol consumption and wheel running include multiple potassium channel genes, Oprm1, Prkcg, Stxbp1, Crhr1, Gabra3, Slc6a13, Stx1b, Pomc, Rassf5 and Camta2. After observing an overlap of many genes and functional groups previously identified in studies of initial sensitivity to ethanol, we hypothesized that wheel running may induce a change in sensitivity, thereby affecting ethanol consumption. A behavioral study examining Loss of Righting Reflex to ethanol following exercise trended toward supporting this hypothesis. These data provide a rich resource for future studies that may better characterize the observed transcriptional changes in gene networks in response to ethanol consumption and wheel running.

Identifying rare variants for genetic risk through a combined pedigree and phenotype approach: application to suicide and asthma.

Suicidal behavior is a complex disorder, with evidence for genetic risk independent of other genetic risk factors including psychiatric disorders. Since 1996, over 3000 DNA samples from Utah suicide decedents have been collected and banked for research use through the Utah Medical Examiner. In addition, over 12,000 Utah suicides were identified through examination of death certificates back to 1904. By linking this data with the Utah Population Database, we have identified multiple extended pedigrees with increased risk for suicide completion. A number of medical conditions co-occur with suicide, including asthma, and this study was undertaken to identify genetic risk common to asthma and suicide. This study tests the hypothesis that a particular comorbid condition may identify a more homogeneous genetic subgroup, facilitating the identification of specific genetic risk factors in that group. From pedigrees at increased risk for suicide, we identified three pedigrees also at significantly increased familial risk for asthma. Five suicide decedents from each of these pedigrees, plus an additional three decedents not from these pedigrees with diagnosed asthma, and 10 decedents with close relatives with asthma were genotyped. Results were compared with 183 publicly available unaffected control exomes from 1000 Genomes and CEPH (Centre d'etude du polymorphisme humain) samples genotyped on the same platform. A further 432 suicide decedents were also genotyped as non-asthma suicide controls. Genotyping was done using the Infinium HumanExome BeadChip. For analysis, we used the pedigree extension of Variant Annotation, Analysis and Search Tool (pVAAST) to calculate the disease burden of each gene. The Phenotype Driven Variant Ontological Re-ranking tool (Phevor) then re-ranked our pVAAST results in context of the phenotype. Using asthma as a seed phenotype, Phevor traversed biomedical ontologies and identified genes with similar biological properties to those known to result in asthma. Our top associated genes included those related to neurodevelopment or neural signaling (brain-derived neurotrophic factor (BDNF), neutral sphingomyelinase 2 (SMPD2), homeobox b2 (HOXB2), neural cell adhesion molecule (NCAM2), heterogeneous nuclear ribonucleoprotein A0 (HNRNPA0)), inflammation (free fatty acid receptor 2 (FFAR2)) and inflammation with additional evidence of neuronal involvement (oxidized low density lipoprotein receptor 1 (OLR1), toll-like receptor 3 (TLR3)). Of particular interest, BDNF has been previously implicated in both psychiatric disorders and asthma. Our results demonstrate the utility of combining pedigree and co-occurring phenotypes to identify rare variants associated with suicide risk in conjunction with specific co-occurring conditions.

Genetic risk factors in two Utah pedigrees at high risk for suicide.

We have used unique population-based data resources to identify 22 high-risk extended pedigrees that show clustering of suicide over twice that expected from demographically adjusted incidence rates. In this initial study of genetic risk factors, we focused on two high-risk pedigrees. In the first of these (pedigree 12), 10/19 (53%) of the related suicides were female, and the average age at death was 30.95. In the second (pedigree 5), 7/51 (14%) of the suicides were female and the average age at death was 36.90. Six decedents in pedigree 12 and nine in pedigree 5 were genotyped with the Illumina HumanExome BeadChip. Genotypes were analyzed using the Variant Annotation, Analysis, and Search program package that computes likelihoods of risk variants using the functional impact of the DNA variation, aggregative scoring of multiple variants across each gene and pedigree structure. We prioritized variants that were: (1) shared across pedigree members, (2) rare in other Utah suicides not related to these pedigrees, (3) < or = 5% in genotyping data from 398 other Utah population controls and (4) < or = 5% frequency in publicly available sequence data from 1358 controls and/or in dbSNP. Results included several membrane protein genes (ANO5, and TMEM141 for pedigree 12 and FAM38A and HRCT1 for pedigree 5). Other genes with known neuronal involvement and/or previous associations with psychiatric conditions were also identified, including NFKB1, CASP9, PLXNB1 and PDE11A in pedigree 12, and THOC1, and AUTS2 in pedigree 5. Although the study is limited to variants included on the HumanExome BeadChip, these findings warrant further exploration, and demonstrate the utility of this high-risk pedigree resource to identify potential genes or gene pathways for future development of targeted interventions.

Transcriptome analysis of Inbred Long Sleep and Inbred Short Sleep mice.

Many studies have utilized the Inbred Long Sleep and Inbred Short Sleep mouse strains to model the genetic influence on initial sensitivity to ethanol. The mechanisms underlying this divergent phenotype are still not completely understood. In this study, we attempt to identify genes that are differentially expressed between these two strains and to identify baseline networks of co-expressed genes, which may provide insight regarding their phenotypic differences. We examined the whole brain and striatal transcriptomes of both strains, using next generation RNA sequencing techniques. Many genes were differentially expressed between strains, including several in chromosomal regions previously shown to influence initial sensitivity to ethanol. These results are in concordance with a similar sample of striatal transcriptomes measured using microarrays. In addition to the higher dynamic range, RNA-Seq is not hindered by high background noise or polymorphisms in probesets as with microarray technology, and we are able to analyze exome sequence of abundant genes. Furthermore, utilizing Weighted Gene Co-expression Network Analysis, we identified several modules of co-expressed genes corresponding to strain differences. Several candidate genes were identified, including protein phosphatase 1 regulatory unit 1b (Ppp1r1b), prodynorphin (Pdyn), proenkephalin (Penk), ras association (RalGDS/AF-6) domain family member 2 (Rassf2), myosin 1d (Myo1d) and transthyretin (Ttr). In addition, we propose a role for potassium channel activity as well as map kinase signaling in the observed phenotypic differences between the two strains.

Regulation of the cycling of timeless (tim) RNA.

Circadian rhythms in Drosophila depend upon expression of the timeless (tim) and period (per) genes, which encode interacting components of the endogenous clock. These two clock genes show a robust circadian oscillation in transcription rate as well as RNA and protein levels. Transcriptional activation of both genes requires the basic helix-loop-helix (bHLH) PAS transcription factors dCLOCK (dCLK) and CYCLE (CYC), which bind E-box elements. We investigated the role of E-box elements in regulating behavioral rhythmicity and tim gene expression. We show that mutation of the upstream E-box in the tim gene prevents the rescue by tim cDNA sequences of the arrhythmic tim(01) phenotype. RNA encoded by this mutated tim transgene fails to cycle and is expressed at low levels. While a tim transgene carrying a wild-type E-box restores behavioral rhythms, tim RNA levels are intermediate to those of the mutant E-box transgenic lines and wild type, and do not display high amplitude cycling. On the other hand, high-amplitude RNA cycling was consistently obtained with a tim transgene that contains genomic, rather than cDNA, sequences. To identify additional sequences that may be required for tim cycling, we investigated the role of an E-box in the first intron of the tim gene through cell culture experiments. In these experiments, the presence of this intron did not have any effect on the activation of tim transcription by dCLK/CYC. As the upstream E-box was implicated in activation by dCLK/CYC in cell culture, we assayed sequences containing this E-box for association with proteins in fly head extracts. These studies provide the first biochemical evidence for an in vivo complex containing dCLK and CYC that binds the tim upstream sequence and is detected at all times of day. Together, these data highlight molecular mechanisms that are critical for behavioral rhythms.

The period E-box is sufficient to drive circadian oscillation of transcription in vivo.

The minimum element from the Drosophila period promoter capable of driving in vivo cycling mRNA is the 69 bp circadian regulatory sequence (CRS). In cell culture, an 18 bp E-box element from the period promoter is regulated by five genes that are involved in the regulation of circadian expression in flies. This E-box is a target for transcriptional activation by bHLH-PAS proteins dCLOCK (dCLK) and CYCLE (CYC), this activation is inhibited by PERIOD (PER) and TIMELESS (TIM) together, and inhibition of dCLK/CYC by PER and TIM is blocked by CRYPTOCHROME (CRY) in the presence of light. Here, the same 18 bp E-box region generated rhythmic expression of luciferase in flies under both light-dark cycling and constant conditions. Flies heterozygous for the Clke(jrk) mutation maintained rhythmic expression from the E-box although at a lower level than wild type. Homozygous mutant Clk(jrk) animals had drastically lowered and arrhythmic expression. In a per01 background, expression from the E-box was high and not rhythmic. Transcription mediated by the per E-box was restricted to the same spatial pattern as the CRS. The per E-box DNA element and cognate binding proteins can confer per-like temporal and spatial expression. This demonstrates in vivo that the known circadian genes that form the core of the circadian oscillator in Drosophila integrate their activities at a single DNA element.

Specific sequences outside the E-box are required for proper per expression and behavioral rescue.

A 69 bp circadian regulatory sequence (CRS) upstream of the per gene is sufficient to drive circadian transcription, mediate proper spatial expression, and rescue behavioral rhythmicity in per01 flies. Within the CRS, an E-box is required for transcriptional activation by two basic-helix-loop-helix (bHLH) PERARNT-SIM (PAS) transcription factors, dCLOCK (dCLK) and CYCLE (CYC). To define sequences within the CRS that are required for spatial expression, circadian expression, and behavioral rhythmicity, a series of mutants that alter blocks of 3 to 12 nucleotides across the entire CRS were used to drive lacZ or per expression in vivo. As expected, the E-box within the CRS is necessary for high-level expression and behavioral rhythmicity, but sequences outside the E-box are also required for transcriptional activation, proper spatial expression, and behavioral rhythmicity. These results indicate that the dCLK-CYC target site extends beyond the E-box and that factors other than dCLK and CYC modulate per transcription.

Light-dependent sequestration of TIMELESS by CRYPTOCHROME.

Most organisms have circadian clocks consisting of negative feedback loops of gene regulation that facilitate adaptation to cycles of light and darkness. In this study, CRYPTOCHROME (CRY), a protein involved in circadian photoperception in Drosophila, is shown to block the function of PERIOD/TIMELESS (PER/TIM) heterodimeric complexes in a light-dependent fashion. TIM degradation does not occur under these conditions; thus, TIM degradation is uncoupled from abrogation of its function by light. CRY and TIM are part of the same complex and directly interact in yeast in a light-dependent fashion. PER/TIM and CRY influence the subcellular distribution of these protein complexes, which reside primarily in the nucleus after the perception of a light signal. Thus, CRY acts as a circadian photoreceptor by directly interacting with core components of the circadian clock.

Closing the circadian loop: CLOCK-induced transcription of its own inhibitors per and tim.

The circadian oscillator generates a rhythmic output with a period of about 24 hours. Despite extensive studies in several model systems, the biochemical mode of action has not yet been demonstrated for any of its components. Here, the Drosophila CLOCK protein was shown to induce transcription of the circadian rhythm genes period and timeless. dCLOCK functioned as a heterodimer with a Drosophila homolog of BMAL1. These proteins acted through an E-box sequence in the period promoter. The timeless promoter contains an 18-base pair element encompassing an E-box, which was sufficient to confer dCLOCK responsiveness to a reporter gene. PERIOD and TIMELESS proteins blocked dCLOCK's ability to transactivate their promoters via the E-box. Thus, dCLOCK drives expression of period and timeless, which in turn inhibit dCLOCK's activity and close the circadian loop.

The automated modified Russell viper venom time test for the lupus anticoagulant.

Antiphospholipid antibodies, both anticardiolipin antibody and lupus anticoagulant, are associated with a hypercoagulable state, manifested as thrombotic events and pregnancy losses. Multiple coagulation tests are available for the lupus anticoagulant, but few are in wide use. The modified Russell viper venom test (RVVT) is one of several excellent tests for the lupus anticoagulant, with high sensitivity and specificity. We present an automated method for measuring the RVVT and demonstrate its reproducibility.