Genome-wide association analysis identifies novel loci for chronotype in 100,420 individuals from the UK Biobank

Our sleep timing preference, or chronotype, is a manifestation of our internal biological clock. Variation in chronotype has been linked to sleep disorders, cognitive and physical performance, and chronic disease. Here we perform a genome-wide association study of self-reported chronotype within the UK Biobank cohort (n=100,420). We identify 12 new genetic loci that implicate known components of the circadian clock machinery and point to previously unstudied genetic variants and candidate genes that might modulate core circadian rhythms or light-sensing pathways. Pathway analyses highlight central nervous and ocular systems and fear-response-related processes. Genetic correlation analysis suggests chronotype shares underlying genetic pathways with schizophrenia, educational attainment and possibly BMI. Further, Mendelian randomization suggests that evening chronotype relates to higher educational attainment. These results not only expand our knowledge of the circadian system in humans but also expose the influence of circadian characteristics over human health and life-history variables such as educational attainment.

Genome-wide screening for DNA variants associated with reading and language traits

Reading and language abilities are heritable traits that are likely to share some genetic influences with each other. To identify pleiotropic genetic variants affecting these traits, we first performed a genome-wide association scan (GWAS) meta-analysis using three richly characterized datasets comprising individuals with histories of reading or language problems, and their siblings. GWAS was performed in a total of 1862 participants using the first principal component computed from several quantitative measures of reading- and language-related abilities, both before and after adjustment for performance IQ. We identified novel suggestive associations at the SNPs rs59197085 and rs5995177 (uncorrected P≈10 for each SNP), located respectively at the CCDC136/FLNC and RBFOX2 genes. Each of these SNPs then showed evidence for effects across multiple reading and language traits in univariate association testing against the individual traits. FLNC encodes a structural protein involved in cytoskeleton remodelling, while RBFOX2 is an important regulator of alternative splicing in neurons. The CCDC136/FLNC locus showed association with a comparable reading/language measure in an independent sample of 6434 participants from the general population, although involving distinct alleles of the associated SNP. Our datasets will form an important part of on-going international efforts to identify genes contributing to reading and language skills. Genome-wide association scan meta-analysis for reading and language ability. © 2014 John Wiley & Sons Ltd and International Behavioural and Neural Genetics Society.

PCSK6 is associated with handedness in individuals with dyslexia

Approximately 90% of humans are right-handed. Handedness is a heritable trait, yet the genetic basis is not well understood. Here we report a genome-wide association study for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)]. The most highly associated marker, rs11855415 (P = 4.7 × 10-7), is located within PCSK6. Two independent cohorts with RD show the same trend, with the minor allele conferring greater relative right-hand skill. Meta-analysis of all three RD samples is genome-wide significant (n = 744, P = 2.0 × 10-8). Conversely, in the general population (n = 2666), we observe a trend towards reduced laterality of hand skill for the minor allele (P = 0.0020). These results provide molecular evidence that cerebral asymmetry and dyslexia are linked. Furthermore, PCSK6 is a protease that cleaves the left–right axis determining protein NODAL. Functional studies of PCSK6 promise insights into mechanisms underlying cerebral lateralization and dyslexia.

Common variants in left/right asymmetry genes and pathways are associated with relative hand skill

Humans display structural and functional asymmetries in brain organization, strikingly with respect to language and handedness. The molecular basis of these asymmetries is unknown. We report a genome-wide association study meta-analysis for a quantitative measure of relative hand skill in individuals with dyslexia [reading disability (RD)] (n = 728). The most strongly associated variant, rs7182874 (P = 8.68×10-9), is located in PCSK6, further supporting an association we previously reported. We also confirmed the specificity of this association in individuals with RD; the same locus was not associated with relative hand skill in a general population cohort (n = 2,666). As PCSK6 is known to regulate NODAL in the development of left/right (LR) asymmetry in mice, we developed a novel approach to GWAS pathway analysis, using gene-set enrichment to test for an over-representation of highly associated variants within the orthologs of genes whose disruption in mice yields LR asymmetry phenotypes. Four out of 15 LR asymmetry phenotypes showed an over-representation (FDR≤5%). We replicated three of these phenotypes; situs inversus, heterotaxia, and double outlet right ventricle, in the general population cohort (FDR≤5%). Our findings lead us to propose that handedness is a polygenic trait controlled in part by the molecular mechanisms that establish LR body asymmetry early in development. © 2013 Brandler et al.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10(-8)) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left-right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Further evidence for a parent-of-origin effect at the NOP9 locus on language-related phenotypes

Background - Specific language impairment (SLI) is a common neurodevelopmental disorder, observed in 5–10 % of children. Family and twin studies suggest a strong genetic component, but relatively few candidate genes have been reported to date. A recent genome-wide association study (GWAS) described the first statistically significant association specifically for a SLI cohort between a missense variant (rs4280164) in the NOP9 gene and language-related phenotypes under a parent-of-origin model. Replications of these findings are particularly challenging because the availability of parental DNA is required. Methods - We used two independent family-based cohorts characterised with reading- and language-related traits: a longitudinal cohort (n = 106 informative families) including children with language and reading difficulties and a nuclear family cohort (n = 264 families) selected for dyslexia. Results - We observed association with language-related measures when modelling for parent-of-origin effects at the NOP9 locus in both cohorts: minimum P = 0.001 for phonological awareness with a paternal effect in the first cohort and minimum P = 0.0004 for irregular word reading with a maternal effect in the second cohort. Allelic and parental trends were not consistent when compared to the original study. Conclusions - A parent-of-origin effect at this locus was detected in both cohorts, albeit with different trends. These findings contribute in interpreting the original GWAS report and support further investigations of the NOP9 locus and its role in language-related traits. A systematic evaluation of parent-of-origin effects in genetic association studies has the potential to reveal novel mechanisms underlying complex traits.

Genome-wide analysis identifies a general requirement for polarity proteins in endocytic traffic

In a genome-wide RNA-mediated interference screen for genes required in membrane traffic - including endocytic uptake, recycling from endosomes to the plasma membrane, and secretion - we identified 168 candidate endocytosis regulators and 100 candidate secretion regulators. Many of these candidates are highly conserved among metazoans but have not been previously implicated in these processes. Among the positives from the screen, we identified PAR-3, PAR-6, PKC-3 and CDC-42, proteins that are well known for their importance in the generation of embryonic and epithelial-cell polarity. Further analysis showed that endocytic transport in Caenorhabditis elegans coelomocytes and human HeLa cells was also compromised after perturbation of CDC-42/Cdc42 or PAR-6/Par6 function, indicating a general requirement for these proteins in regulating endocytic traffic. Consistent with these results, we found that tagged CDC-42/Cdc42 is enriched on recycling endosomes in C. elegans and mammalian cells, suggesting a direct function in the regulation of transport.

Expression profiling of type 2 diabetes susceptibility genes in the pancreatic islets, adipose tissue and liver of obese mice

Type 2 diabetes (T2D) is characterized by impaired beta cell function and insulin resistance. T2D susceptibility genes identified by Genome-wide association studies (GWAS) are likely to have roles in both impaired insulin secretion from the beta cell as well as insulin resistance. The aim of this study was to use gene expression profiling to assess the effect of the diabetic milieu on the expression of genes involved in both insulin secretion and insulin resistance. We measured the expression of 43 T2D susceptibility genes in the islets, adipose and liver of leptin-deficient Ob/Ob mice compared with Ob/+ littermates. The same panel of genes were also profiled in cultured rodent adipocytes, hepatocytes and beta cells in response to high glucose conditions, to distinguish expression effects due to elevated glycemia from those on the causal pathway to diabetes or induced by other factors in the diabetic microenviroment. We found widespread deregulation of these genes in tissues from Ob/Ob mice, with differential regulation of 23 genes in adipose, 18 genes in liver and one gene (Tcf7l2) in islets of diabetic animals (Ob/Ob) compared to control (Ob/+) animals. However, these expression changes were in most cases not noted in glucose-treated adipocyte, hepatocyte or beta cell lines, indicating that they may not be an effect of hyperglycemia alone. This study indicates that expression changes are apparent with diabetes in both the insulin producing beta cells, but also in peripheral tissues involved in insulin resistance. This suggests that incidence or progression of diabetic phenotypes in a mouse model of diabetes is driven by both secretory and peripheral defects. © J. A. Barth Verlag in Georg Thieme Verlag KG Stuttgart New York.

Independent modulation of engagement and connectivity of the facial network during affect processing by CACNA1C and ANK3 risk genes for bipolar disorder

IMPORTANCE Genome-wide association studies (GWASs) indicate that single-nucleotide polymorphisms in the CACNA1C and ANK3 genes increase the risk for bipolar disorder (BD). The genes influence neuronal firing by modulating calcium and sodium channel functions, respectively. Both genes modulate ?-aminobutyric acid-transmitting interneuron function and can thus affect brain regional activation and interregional connectivity. OBJECTIVE To determine whether the genetic risk for BD associated with 2 GWAS-supported risk single-nucleotide polymorphisms at CACNA1C rs1006737 and ANK3 rs10994336 is mediated through changes in regional activation and interregional connectivity of the facial affect-processing network. DESIGN, SETTING, AND PARTICIPANTS Cross-sectional functional magnetic resonance imaging study at a research institute of 41 euthymic patients with BD and 46 healthy participants, all of British white descent. MAIN OUTCOMES AND MEASURES Blood oxygen level-dependent signal and effective connectivity measures during the facial affect-processing task. RESULTS In healthy carriers, both genetic risk variants were independently associated with increased regional engagement throughout the facial affect-processing network and increased effective connectivity between the visual and ventral prefrontal cortical regions. In contrast, BD carriers of either genetic risk variant exhibited pronounced reduction in ventral prefrontal cortical activation and visual-prefrontal effective connectivity. CONCLUSIONS AND RELEVANCE Our data demonstrate that the effect of CACNA1C rs1006737 and ANK3 rs10994336 (or genetic variants in linkage disequilibrium) on the brain converges on the neural circuitry involved in affect processing and provides a mechanism linking BD to genome-wide genetic risk variants.

The impact of the CACNA1C gene polymorphism on frontolimbic function in bipolar disorder

Genome-wide association studies in bipolar disorder (BD)1 have implicated a single-nucleotide polymorphism (rs1006737, G right arrow A) in the CACNA1C gene, which encodes for the alpha 1c (CAV1.2) subunit of the voltage-gated, L-type calcium channel. Neuroimaging studies of healthy individuals report that this risk allele modulates brain function within limbic (amygdala, anterior cingulate gyrus) and hippocampal regions during tasks of reward processing2, 3 and episodic memory. Moreover, animal studies suggest that the CaV1.2 L-type calcium channels influence emotional behaviour through enhanced neurotransmission via the lateral amygdala pathway. On the basis of this evidence, we tested the hypotheses that the CACNA1C rs1006737 risk allele will modulate neural responses within predefined prefrontal and subcortical regions of interest during emotional face processing and that this effect would be amplified in BD patients.