Genetic variants affect diurnal glucose levels throughout the day.

Circadian rhythms not only coordinate the timing of wake and sleep but also regulate homeostasis within the body, including glucose metabolism. However, the genetic variants that contribute to temporal control of glucose levels have not been previously examined. Using data from 420,000 individuals from the UK Biobank and replicating our findings in 100,000 individuals from the Estonian Biobank, we show that diurnal serum glucose is under genetic control. We discover a robust temporal association of glucose levels at the Melatonin receptor 1B (MTNR1B) (rs10830963, P = 1e-22) and a canonical circadian pacemaker gene Cryptochrome 2 (CRY2) loci (rs12419690, P = 1e-16). Furthermore, we show that sleep modulates serum glucose levels and the genetic variants have a separate mechanism of diurnal control. Finally, we show that these variants independently modulate risk of type 2 diabetes. Our findings, together with earlier genetic and epidemiological evidence, show a clear connection between sleep and metabolism and highlight variation at MTNR1B and CRY2 as temporal regulators for glucose levels.

Genetic contribution to disease-course severity and progression in the SUPER-Finland study, a cohort of 10,403 individuals with psychotic disorders.

Genetic factors contribute to the susceptibility of psychotic disorders, but less is known how they affect psychotic disease-course development. Utilizing polygenic scores (PGSs) in combination with longitudinal healthcare data with decades of follow-up we investigated the contributing genetics to psychotic disease-course severity and diagnostic shifts in the SUPER-Finland study, encompassing 10 403 genotyped individuals with a psychotic disorder. To longitudinally track the study participants' past disease-course severity, we created a psychiatric hospitalization burden metric using the full-coverage and nation-wide Finnish in-hospital registry (data from 1969 and onwards). Using a hierarchical model, ranking the psychotic diagnoses according to clinical severity, we show that high schizophrenia PGS (SZ-PGS) was associated with progression from lower ranked psychotic disorders to schizophrenia (OR = 1.32 [1.23-1.43], p = 1.26e-12). This development manifested already at psychotic illness onset as a higher psychiatric hospitalization burden, the proxy for disease-course severity. In schizophrenia (n = 5 479), both a high SZ-PGS and a low educational attainment PGS (EA-PGS) were associated with increased psychiatric hospitalization burden (p = 1.00e-04 and p = 4.53e-10). The SZ-PGS and the EA-PGS associated with distinct patterns of hospital usage. In individuals with high SZ-PGS, the increased hospitalization burden was composed of longer individual hospital stays, while low EA-PGS associated with shorter but more frequent hospital visits. The negative effect of a low EA-PGS was found to be partly mediated via substance use disorder, a major risk factor for hospitalizations. In conclusion, we show that high SZ-PGS and low EA-PGS both impacted psychotic disease-course development negatively but resulted in different disease-course trajectories.

Evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

Intellectual disability (ID) is a common disorder, yet there is a wide spectrum of impairment from mild to profoundly affected individuals. Mild ID is seen as the low extreme of the general distribution of intelligence, while severe ID is often seen as a monogenic disorder caused by rare, pathogenic, highly penetrant variants. To investigate the genetic factors influencing mild and severe ID, we evaluated rare and common variation in the Northern Finland Intellectual Disability cohort (n = 1096 ID patients), a cohort with a high percentage of mild ID (n = 550) and from a population bottleneck enriched in rare, damaging variation. Despite this enrichment, we found only a small percentage of ID was due to recessive Finnish-enriched variants (0.5%). A larger proportion was linked to dominant variation, with a significant burden of rare, damaging variation in both mild and severe ID. This rare variant burden was enriched in more severe ID (p = 2.4e-4), patients without a relative with ID (p = 4.76e-4), and in those with features associated with monogenic disorders. We also found a significant burden of common variants associated with decreased cognitive function, with no difference between mild and more severe ID. When we included common and rare variants in a joint model, the rare and common variants had additive effects in both mild and severe ID. A multimodel inference approach also found that common and rare variants together best explained ID status (ΔAIC = 16.8, ΔBIC = 10.2). Overall, we report evidence for the additivity of rare and common variant burden throughout the spectrum of intellectual disability.

The impact of rare protein coding genetic variation on adult cognitive function.

Compelling evidence suggests that human cognitive function is strongly influenced by genetics. Here, we conduct a large-scale exome study to examine whether rare protein-coding variants impact cognitive function in the adult population (n = 485,930). We identify eight genes (ADGRB2, KDM5B, GIGYF1, ANKRD12, SLC8A1, RC3H2, CACNA1A and BCAS3) that are associated with adult cognitive function through rare coding variants with large effects. Rare genetic architecture for cognitive function partially overlaps with that of neurodevelopmental disorders. In the case of KDM5B we show how the genetic dosage of one of these genes may determine the variability of cognitive, behavioral and molecular traits in mice and humans. We further provide evidence that rare and common variants overlap in association signals and contribute additively to cognitive function. Our study introduces the relevance of rare coding variants for cognitive function and unveils high-impact monogenic contributions to how cognitive function is distributed in the normal adult population.

Functional Characterization of Six SLCO1B1 (OATP1B1) Variants Observed in Finnish Individuals with a Psychotic Disorder.

Variants in the SLCO1B1 (solute carrier organic anion transporter family member 1B1) gene encoding the OATP1B1 (organic anion transporting polypeptide 1B1) protein are associated with altered transporter function that can predispose patients to adverse drug effects with statin treatment. We explored the effect of six rare SLCO1B1 single nucleotide variants (SNVs) occurring in Finnish individuals with a psychotic disorder on expression and functionality of the OATP1B1 protein. The SUPER-Finland study has performed exome sequencing on 9381 individuals with at least one psychotic episode during their lifetime. SLCO1B1 SNVs were annotated with PHRED-scaled combined annotation-dependent (CADD) scores and the Ensembl variant effect predictor. In vitro functionality studies were conducted for the SNVs with a PHRED-scaled CADD score of >10 and predicted to be missense. To estimate possible changes in transport activity caused by the variants, transport of 2',7'-dichlorofluorescein (DCF) in OATP1B1-expressing HEK293 cells was measured. According to the findings, additional tests with rosuvastatin and estrone sulfate were conducted. The amount of OATP1B1 in crude membrane fractions was quantified using a liquid chromatography tandem mass spectrometry-based quantitative targeted absolute proteomics analysis. Six rare missense variants of SLCO1B1 were identified in the study population, located in transmembrane helix 3: c.317T>C (p.106I>T), intracellular loop 2: c.629G>T (p.210G>V), c.633A>G (p.211I>M), c.639T>A (p.213N>L), transmembrane helix 6: 820A>G (p.274I>V), and the C-terminal end: 2005A>C (p.669N>H). Of these variants, SLCO1B1 c.629G>T (p.210G>V) resulted in the loss of in vitro function, abolishing the uptake of DCF, estrone sulfate, and rosuvastatin and reducing the membrane protein expression to 31% of reference OATP1B1. Of the six rare missense variants, SLCO1B1 c.629G>T (p.210G>V) causes a loss of function of OATP1B1 transport in vitro and severely decreases membrane protein abundance. Carriers of SLCO1B1 c.629G>T might be susceptible to altered pharmacokinetics of OATP1B1 substrate drugs and might have increased likelihood of adverse drug effects such as statin-associated musculoskeletal symptoms.

A novel variant in SMG9 causes intellectual disability, confirming a role for nonsense-mediated decay components in neurocognitive development.

Biallelic loss-of-function variants in the SMG9 gene, encoding a regulatory subunit of the mRNA nonsense-mediated decay (NMD) machinery, are reported to cause heart and brain malformation syndrome. Here we report five patients from three unrelated families with intellectual disability (ID) and a novel pathogenic SMG9 c.551 T > C p.(Val184Ala) homozygous missense variant, identified using exome sequencing. Sanger sequencing confirmed recessive segregation in each family. SMG9 c.551T > C p.(Val184Ala) is most likely an autozygous variant identical by descent. Characteristic clinical findings in patients were mild to moderate ID, intention tremor, pyramidal signs, dyspraxia, and ocular manifestations. We used RNA sequencing of patients and age- and sex-matched healthy controls to assess the effect of the variant. RNA sequencing revealed that the SMG9 c.551T > C variant did not affect the splicing or expression level of SMG9 gene products, and allele-specific expression analysis did not provide evidence that the nonsense mRNA-induced NMD was affected. Differential gene expression analysis identified prevalent upregulation of genes in patients, including the genes SMOX, OSBP2, GPX3, and ZNF155. These findings suggest that normal SMG9 function may be involved in transcriptional regulation without affecting nonsense mRNA-induced NMD. In conclusion, we demonstrate that the SMG9 c.551T > C missense variant causes a neurodevelopmental disorder and impacts gene expression. NMD components have roles beyond aberrant mRNA degradation that are crucial for neurocognitive development.

Polygenic burden has broader impact on health, cognition, and socioeconomic outcomes than most rare and high-risk copy number variants.

Copy number variants (CNVs) are associated with syndromic and severe neurological and psychiatric disorders (SNPDs), such as intellectual disability, epilepsy, schizophrenia, and bipolar disorder. Although considered high-impact, CNVs are also observed in the general population. This presents a diagnostic challenge in evaluating their clinical significance. To estimate the phenotypic differences between CNV carriers and non-carriers regarding general health and well-being, we compared the impact of SNPD-associated CNVs on health, cognition, and socioeconomic phenotypes to the impact of three genome-wide polygenic risk score (PRS) in two Finnish cohorts (FINRISK, n = 23,053 and NFBC1966, n = 4895). The focus was on CNV carriers and PRS extremes who do not have an SNPD diagnosis. We identified high-risk CNVs (DECIPHER CNVs, risk gene deletions, or large [>1 Mb] CNVs) in 744 study participants (2.66%), 36 (4.8%) of whom had a diagnosed SNPD. In the remaining 708 unaffected carriers, we observed lower educational attainment (EA; OR = 0.77 [95% CI 0.66-0.89]) and lower household income (OR = 0.77 [0.66-0.89]). Income-associated CNVs also lowered household income (OR = 0.50 [0.38-0.66]), and CNVs with medical consequences lowered subjective health (OR = 0.48 [0.32-0.72]). The impact of PRSs was broader. At the lowest extreme of PRS for EA, we observed lower EA (OR = 0.31 [0.26-0.37]), lower-income (OR = 0.66 [0.57-0.77]), lower subjective health (OR = 0.72 [0.61-0.83]), and increased mortality (Cox's HR = 1.55 [1.21-1.98]). PRS for intelligence had a similar impact, whereas PRS for schizophrenia did not affect these traits. We conclude that the majority of working-age individuals carrying high-risk CNVs without SNPD diagnosis have a modest impact on morbidity and mortality, as well as the limited impact on income and educational attainment, compared to individuals at the extreme end of common genetic variation. Our findings highlight that the contribution of traditional high-risk variants such as CNVs should be analyzed in a broader genetic context, rather than evaluated in isolation.

The bradykinin system in stress and anxiety in humans and mice.

Pharmacological research in mice and human genetic analyses suggest that the kallikrein-kinin system (KKS) may regulate anxiety. We examined the role of the KKS in anxiety and stress in both species. In human genetic association analysis, variants in genes for the bradykinin precursor (KNG1) and the bradykinin receptors (BDKRB1 and BDKRB2) were associated with anxiety disorders (p < 0.05). In mice, however, neither acute nor chronic stress affected B1 receptor gene or protein expression, and B1 receptor antagonists had no effect on anxiety tests measuring approach-avoidance conflict. We thus focused on the B2 receptor and found that mice injected with the B2 antagonist WIN 64338 had lowered levels of a physiological anxiety measure, the stress-induced hyperthermia (SIH), vs controls. In the brown adipose tissue, a major thermoregulator, WIN 64338 increased expression of the mitochondrial regulator Pgc1a and the bradykinin precursor gene Kng2 was upregulated after cold stress. Our data suggests that the bradykinin system modulates a variety of stress responses through B2 receptor-mediated effects, but systemic antagonists of the B2 receptor were not anxiolytic in mice. Genetic variants in the bradykinin receptor genes may predispose to anxiety disorders in humans by affecting their function.

Focused multidimensional scaling: interactive visualization for exploration of high-dimensional data.

BACKGROUND: Visualization is an important tool for generating meaning from scientific data, but the visualization of structures in high-dimensional data (such as from high-throughput assays) presents unique challenges. Dimension reduction methods are key in solving this challenge, but these methods can be misleading- especially when apparent clustering in the dimension-reducing representation is used as the basis for reasoning about relationships within the data. RESULTS: We present two interactive visualization tools, distnet and focusedMDS, that help in assessing the validity of a dimension-reducing plot and in interactively exploring relationships between objects in the data. The distnet tool is used to examine discrepancies between the placement of points in a two dimensional visualization and the points' actual similarities in feature space. The focusedMDS tool is an intuitive, interactive multidimensional scaling tool that is useful for exploring the relationships of one particular data point to the others, that might be useful in a personalized medicine framework. CONCLUSIONS: We introduce here two freely available tools for visually exploring and verifying the validity of dimension-reducing visualizations and biological information gained from these. The use of such tools can confirm that conclusions drawn from dimension-reducing visualizations are not simply artifacts of the visualization method, but are real biological insights.

Stress-free automatic sleep deprivation using air puffs.

BACKGROUND: Sleep deprivation via gentle handling is time-consuming and personnel-intensive. NEW METHOD: We present here an automated sleep deprivation system via air puffs. Implanted EMG and EEG electrodes were used to assess sleep/waking states in six male Sprague-Dawley rats. Blood samples were collected from an implanted intravenous catheter every 4h during the 12-h light cycle on baseline, 8h of sleep deprivation via air puffs, and 8h of sleep deprivation by gentle handling days. RESULTS: The automated system was capable of scoring sleep and waking states as accurately as our offline version (∼90% for sleep) and with sufficient speed to trigger a feedback response within an acceptable amount of time (1.76s). Manual state scoring confirmed normal sleep on the baseline day and sleep deprivation on the two manipulation days (68% decrease in non-REM, 63% decrease in REM, and 74% increase in waking). No significant differences in levels of ACTH and corticosterone (stress hormones indicative of HPA axis activity) were found at any time point between baseline sleep and sleep deprivation via air puffs. COMPARISON WITH EXISTING METHOD: There were no significant differences in ACTH or corticosterone concentrations between sleep deprivation by air puffs and gentle handling over the 8-h period. CONCLUSIONS: Our system accurately detects sleep and delivers air puffs to acutely deprive rats of sleep with sufficient temporal resolution during the critical 4-5h post learning sleep-dependent memory consolidation period. The system is stress-free and a viable alternative to existing sleep deprivation techniques.

Sleep alterations following exposure to stress predict fear-associated memory impairments in a rodent model of PTSD.

Sleep abnormalities, such as insomnia, nightmares, hyper-arousal, and difficulty initiating or maintaining sleep, are diagnostic criteria of posttraumatic stress disorder (PTSD). The vivid dream state, rapid eye movement (REM) sleep, has been implicated in processing emotional memories. We have hypothesized that REM sleep is maladaptive in those suffering from PTSD. However, the precise neurobiological mechanisms regulating sleep disturbances following trauma exposure are poorly understood. Using single prolonged stress (SPS), a well-validated rodent model of PTSD, we measured sleep alterations in response to stressor exposure and over a subsequent 7-day isolation period during which the PTSD-like phenotype develops. SPS resulted in acute increases in REM sleep and transition to REM sleep, and decreased waking in addition to alterations in sleep architecture. The severity of the PTSD-like phenotype was later assessed by measuring freezing levels on a fear-associated memory test. Interestingly, the change in REM sleep following SPS was significantly correlated with freezing behavior during extinction recall assessed more than a week later. Reductions in theta (4-10 Hz) and sigma (10-15 Hz) band power during transition to REM sleep also correlated with impaired fear-associated memory processing. These data reveal that changes in REM sleep, transition to REM sleep, waking, and theta and sigma power may serve as sleep biomarkers to identify individuals with increased susceptibility to PTSD following trauma exposure.