Search

Background: Aberrant functional connectivity is a hallmark of schizophrenia. The precise nature and mechanism of dysconnectivity in schizophrenia remains unclear, but evidence suggests that dysconnectivity is different in wake versus sleep. Microstate analysis uses electroencephalography (EEG) to investigate large-scale patterns of coordinated brain activity by clustering EEG data into a small set of recurring spatial patterns, or microstates. We hypothesized that this technique would allow us to probe connectivity between brain networks at a fine temporal resolution and uncover previously unknown sleep-specific dysconnectivity. Methods: We studied microstates during sleep in patients with schizophrenia by analyzing high-density EEG sleep data from 114 patients with schizophrenia and 79 control participants. We used a polarity-insensitive k-means analysis to extract a set of 6 microstate topographies. Results: These 6 states included 4 widely reported canonical microstates. In patients and control participants, falling asleep was characterized by a shift from microstates A, B, and C to microstates D, E, and F. Microstate F was decreased in patients during wake, and microstate E was decreased in patients during sleep. The complexity of microstate transitions was greater in patients than control participants during wake, but this reversed during sleep. Conclusions: Our findings reveal behavioral state–dependent patterns of cortical dysconnectivity in schizophrenia. Furthermore, these findings are largely unrelated to previous sleep-related EEG markers of schizophrenia such as decreased sleep spindles. Therefore, these findings are driven by previously undescribed sleep-related pathology in schizophrenia.

Abstract Background Objective and quantifiable markers are crucial for developing novel therapeutics for mental disorders by 1) stratifying clinically similar patients with different underlying neurobiological deficits and 2) objectively tracking disease trajectory and treatment response. Schizophrenia is often confounded with other psychiatric disorders, especially bipolar disorder, if based on cross-sectional symptoms. Awake and sleep EEG have shown promise in identifying neurophysiological differences as biomarkers for schizophrenia. However, most previous studies, while useful, were conducted in European and American populations, had small sample sizes, and utilized varying analytic methods, limiting comprehensive analyses or generalizability to diverse human populations. Furthermore, the extent to which wake and sleep neurophysiology metrics correlate with each other and with symptom severity or cognitive impairment remains unresolved. Moreover, how these neurophysiological markers compare across psychiatric conditions is not well characterized. The utility of biomarkers in clinical trials and practice would be significantly advanced by well-powered transdiagnostic studies. The Global Research Initiative on the Neurophysiology of Schizophrenia (GRINS) project aims to address these questions through a large, multi-center cohort study involving East Asian populations. To promote transparency and reproducibility, we describe the protocol for the GRINS project. Methods The research procedure consists of an initial screening interview followed by three subsequent sessions: an introductory interview, an evaluation visit, and an overnight neurophysiological recording session. Data from multiple domains, including demographic and clinical characteristics, behavioral performance (cognitive tasks, motor sequence tasks), and neurophysiological metrics (both awake and sleep electroencephalography), are collected by research groups specialized in each domain. Conclusion Pilot results from the GRINS project demonstrate the feasibility of this study protocol and highlight the importance of such research, as well as its potential to study a broader range of patients with psychiatric conditions. Through GRINS, we are generating a valuable dataset across multiple domains to identify neurophysiological markers of schizophrenia individually and in combination. By applying this protocol to related mental disorders often confounded with each other, we can gather information that offers insight into the neurophysiological characteristics and underlying mechanisms of these severe conditions, informing objective diagnosis, stratification for clinical research, and ultimately, the development of better-targeted treatment matching in the clinic.

Background Excessive daytime sleepiness (EDS), experienced in 10% to 20% of the population, has been associated with cardiovascular disease and death. However, the condition is heterogeneous and is prevalent in individuals having short and long sleep duration. We sought to clarify the relationship between sleep duration subtypes of EDS with cardiovascular outcomes, accounting for these subtypes. Methods and Results We defined 3 sleep duration subtypes of excessive daytime sleepiness: normal (6–9 hours), short (9 hours), and compared these with a nonsleepy, normal‐sleep‐duration reference group. We analyzed their associations with incident myocardial infarction (MI) and stroke using medical records of 355 901 UK Biobank participants and performed 2‐sample Mendelian randomization for each outcome. Compared with healthy sleep, long‐sleep EDS was associated with an 83% increased rate of MI (hazard ratio, 1.83 [95% CI, 1.21–2.77]) during 8.2‐year median follow‐up, adjusting for multiple health and sociodemographic factors. Mendelian randomization analysis provided supporting evidence of a causal role for a genetic long‐sleep EDS subtype in MI (inverse‐variance weighted β=1.995, P=0.001). In contrast, we did not find evidence that other subtypes of EDS were associated with incident MI or any associations with stroke (P>0.05). Conclusions Our study suggests the previous evidence linking EDS with increased cardiovascular disease risk may be primarily driven by the effect of its long‐sleep subtype on higher risk of MI. Underlying mechanisms remain to be investigated but may involve sleep irregularity and circadian disruption, suggesting a need for novel interventions in this population.

Summary: Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

Abstract Obstructive sleep apnea (OSA) is a common disorder associated with increased risk of cardiovascular disease and mortality. Iron and heme metabolism, implicated in ventilatory control and OSA comorbidities, was associated with OSA phenotypes in recent admixture mapping and gene enrichment analyses. However, its causal contribution was unclear. In this study, we performed pathway-level transcriptional Mendelian randomization (MR) analysis to investigate the causal relationships between iron and heme related pathways and OSA. In primary analysis, we examined the expression level of four iron/heme Reactome pathways as exposures and four OSA traits as outcomes using cross-tissue cis-eQTLs from the Genotype-Tissue Expression portal and published genome-wide summary statistics of OSA. We identify a significant putative causal association between up-regulated heme biosynthesis pathway with higher sleep time percentage of hypoxemia (p = 6.14 × 10–3). This association is supported by consistency of point estimates in one-sample MR in the Multi-Ethnic Study of Atherosclerosis using high coverage DNA and RNA sequencing data generated by the Trans-Omics for Precision Medicine project. Secondary analysis for 37 additional iron/heme Gene Ontology pathways did not reveal any significant causal associations. This study suggests a causal association between increased heme biosynthesis and OSA severity.

Abstract Background Sleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing. Methods The study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation

Joint analysis of multiple traits can increase power and provide insights into shared genetic architecture. Here, Nguyen et al. develop multi-trait TADA (mTADA), an extension of TADA (transmission and de novo association test) that jointly analyses de novo mutations of traits for improved risk-gene identification power.

A main symptom of chronic insufficient sleep is excessive daytime sleepiness. Here, Wang et al. report 42 genome-wide significant loci for self-reported daytime sleepiness in 452,071 individuals from the UK Biobank that cluster into two biological subtypes of either sleep propensity or sleep fragmentation.

Sleep is essential for homeostasis and insufficient or excessive sleep are associated with adverse outcomes. Here, the authors perform GWAS for self-reported habitual sleep duration in adults, supported by accelerometer-derived measures, and identify genetic correlation with psychiatric and metabolic traits

Abstract Background Integrating rare variation from trio family and case–control studies has successfully implicated specific genes contributing to risk of neurodevelopmental disorders (NDDs) including autism spectrum disorders (ASD), intellectual disability (ID), developmental disorders (DDs), and epilepsy (EPI). For schizophrenia (SCZ), however, while sets of genes have been implicated through the study of rare variation, only two risk genes have been identified. Methods We used hierarchical Bayesian modeling of rare-variant genetic architecture to estimate mean effect sizes and risk-gene proportions, analyzing the largest available collection of whole exome sequence data for SCZ (1,077 trios, 6,699 cases, and 13,028 controls), and data for four NDDs (ASD, ID, DD, and EPI; total 10,792 trios, and 4,058 cases and controls). Results For SCZ, we estimate there are 1,551 risk genes. There are more risk genes and they have weaker effects than for NDDs. We provide power analyses to predict the number of risk-gene discoveries as more data become available. We confirm and augment prior risk gene and gene set enrichment results for SCZ and NDDs. In particular, we detected 98 new DD risk genes at FDR < 0.05. Correlations of risk-gene posterior probabilities are high across four NDDs (ρ>0.55), but low between SCZ and the NDDs (ρ

Arc is an activity-regulated neuronal protein, but little is known about its interactions, assembly into multiprotein complexes, and role in human disease and cognition. We applied an integrated proteomic and genetic strategy by targeting a tandem affinity purification (TAP) tag and Venus fluorescent protein into the endogenous Arc gene in mice. This allowed biochemical and proteomic characterization of native complexes in wild-type and knockout mice. We identified many Arc-interacting proteins, of which PSD95 was the most abundant. PSD95 was essential for Arc assembly into 1.5-MDa complexes and activity-dependent recruitment to excitatory synapses. Integrating human genetic data with proteomic data showed that Arc-PSD95 complexes are enriched in schizophrenia, intellectual disability, autism, and epilepsy mutations and normal variants in intelligence. We propose that Arc-PSD95 postsynaptic complexes potentially affect human cognitive function.

A large portion of common variant loci associated with genetic risk for schizophrenia reside within noncoding sequence of unknown function. Here, we demonstrate promoter and enhancer enrichment in schizophrenia variants associated with expression quantitative trait loci (eQTL). The enrichment is greater when functional annotations derived from the human brain are used relative to peripheral tissues. Regulatory trait concordance analysis ranked genes within schizophrenia genome-wide significant loci for a potential functional role, based on colocalization of a risk SNP, eQTL, and regulatory element sequence. We identified potential physical interactions of noncontiguous proximal and distal regulatory elements. This was verified in prefrontal cortex and -induced pluripotent stem cell–derived neurons for the L-type calcium channel (CACNA1C) risk locus. Our findings point to a functional link between schizophrenia-associated noncoding SNPs and 3D genome architecture associated with chromosomal loopings and transcriptional regulation in the brain.

Sleep is an important indicator of a person's health, and its accurate and cost-effective quantification is of great value in healthcare. The gold standard for sleep assessment and the clinical diagnosis of sleep disorders is polysomnography (PSG). However, PSG requires an overnight clinic visit and trained technicians to score the obtained multimodality data. Wrist-worn consumer devices, such as smartwatches, are a promising alternative to PSG because of their small form factor, continuous monitoring capability, and popularity. Unlike PSG, however, wearables-derived data are noisier and far less information-rich because of the fewer number of modalities and less accurate measurements due to their small form factor. Given these challenges, most consumer devices perform two-stage (i.e., sleep-wake) classification, which is inadequate for deep insights into a person's sleep health. The challenging multi-class (three, four, or five-class) staging of sleep using data from wrist-worn wearables remains unresolved. The difference in the data quality between consumer-grade wearables and lab-grade clinical equipment is the motivation behind this study. In this paper, we present an artificial intelligence (AI) technique termed sequence-to-sequence LSTM for automated mobile sleep staging (SLAMSS), which can perform three-class (wake, NREM, REM) and four-class (wake, light, deep, REM) sleep classification from activity (i.e., wrist-accelerometry-derived locomotion) and two coarse heart rate measures-both of which can be reliably obtained from a consumer-grade wrist-wearable device. Our method relies on raw time-series datasets and obviates the need for manual feature selection. We validated our model using actigraphy and coarse heart rate data from two independent study populations: the Multi-Ethnic Study of Atherosclerosis (MESA; N = 808) cohort and the Osteoporotic Fractures in Men (MrOS; N = 817) cohort. SLAMSS achieves an overall accuracy of 79%, weighted F1 score of 0.80, 77% sensitivity, and 89% specificity for three-class sleep staging and an overall accuracy of 70-72%, weighted F1 score of 0.72-0.73, 64-66% sensitivity, and 89-90% specificity for four-class sleep staging in the MESA cohort. It yielded an overall accuracy of 77%, weighted F1 score of 0.77, 74% sensitivity, and 88% specificity for three-class sleep staging and an overall accuracy of 68-69%, weighted F1 score of 0.68-0.69, 60-63% sensitivity, and 88-89% specificity for four-class sleep staging in the MrOS cohort. These results were achieved with feature-poor inputs with a low temporal resolution. In addition, we extended our three-class staging model to an unrelated Apple Watch dataset. Importantly, SLAMSS predicts the duration of each sleep stage with high accuracy. This is especially significant for four-class sleep staging, where deep sleep is severely underrepresented. We show that, by appropriately choosing the loss function to address the inherent class imbalance, our method can accurately estimate deep sleep time (SLAMSS/MESA: 0.61±0.69 hours, PSG/MESA ground truth: 0.60±0.60 hours; SLAMSS/MrOS: 0.53±0.66 hours, PSG/MrOS ground truth: 0.55±0.57 hours;). Deep sleep quality and quantity are vital metrics and early indicators for a number of diseases. Our method, which enables accurate deep sleep estimation from wearables-derived data, is therefore promising for a variety of clinical applications requiring long-term deep sleep monitoring.

Motivated by the potential of objective neurophysiological markers to index thalamocortical function in patients with severe psychiatric illnesses, we comprehensively characterized key non-rapid eye movement (NREM) sleep parameters across multiple domains, their interdependencies, and their relationship to waking event-related potentials and symptom severity. In 72 schizophrenia (SCZ) patients and 58 controls, we confirmed a marked reduction in sleep spindle density in SCZ and extended these findings to show that fast and slow spindle properties were largely uncorrelated. We also describe a novel measure of slow oscillation and spindle interaction that was attenuated in SCZ. The main sleep findings were replicated in a demographically distinct sample, and a joint model, based on multiple NREM components, statistically predicted disease status in the replication cohort. Although also altered in patients, auditory event-related potentials elicited during wake were unrelated to NREM metrics. Consistent with a growing literature implicating thalamocortical dysfunction in SCZ, our characterization identifies independent NREM and wake EEG biomarkers that may index distinct aspects of SCZ pathophysiology and point to multiple neural mechanisms underlying disease heterogeneity. This study lays the groundwork for evaluating these neurophysiological markers, individually or in combination, to guide efforts at treatment and prevention as well as identifying individuals most likely to benefit from specific interventions.

Sleep disordered breathing (SDB)-related overnight hypoxemia is associated with cardiometabolic disease and other comorbidities. Understanding the genetic bases for variations in nocturnal hypoxemia may help understand mechanisms influencing oxygenation and SDB-related mortality. We conducted genome-wide association tests across 10 cohorts and 4 populations to identify genetic variants associated with three correlated measures of overnight oxyhemoglobin saturation: average and minimum oxyhemoglobin saturation during sleep and the percent of sleep with oxyhemoglobin saturation under 90%. The discovery sample consisted of 8,326 individuals. Variants with p < 1 × 10(-6) were analyzed in a replication group of 14,410 individuals. We identified 3 significantly associated regions, including 2 regions in multi-ethnic analyses (2q12, 10q22). SNPs in the 2q12 region associated with minimum SpO2 (rs78136548 p = 2.70 × 10(-10)). SNPs at 10q22 were associated with all three traits including average SpO2 (rs72805692 p = 4.58 × 10(-8)). SNPs in both regions were associated in over 20,000 individuals and are supported by prior associations or functional evidence. Four additional significant regions were detected in secondary sex-stratified and combined discovery and replication analyses, including a region overlapping Reelin, a known marker of respiratory complex neurons.These are the first genome-wide significant findings reported for oxyhemoglobin saturation during sleep, a phenotype of high clinical interest. Our replicated associations with HK1 and IL18R1 suggest that variants in inflammatory pathways, such as the biologically-plausible NLRP3 inflammasome, may contribute to nocturnal hypoxemia.

Background: Alzheimer’s disease (AD) and bipolar disorder (BIP) are complex traits influenced by numerous common genetic variants, most of which remain to be detected. Clinical and epidemiological evidence suggest that AD and BIP are related. However, it is not established if this relation is of genetic origin. Here, we applied statistical methods based on the conditional false discovery rate (FDR) framework to detect genetic overlap between AD and BIP and utilized this overlap to increase the power to identify common genetic variants associated with either or both traits.Methods: We obtained genome wide association studies data from the International Genomics of Alzheimer’s Project part 1 (17,008 AD cases and 37,154 controls) and the Psychiatric Genetic Consortium Bipolar Disorder Working Group (20,352 BIP cases and 31,358 controls). We used conditional QQ-plots to assess overlap in common genetic variants between AD and BIP. We exploited the genetic overlap to re-rank test-statistics for AD and BIP and improve detection of genetic variants using the conditional FDR framework.Results: Conditional QQ-plots demonstrated a polygenic overlap between AD and BIP. Using conditional FDR, we identified one novel genomic locus associated with AD, and nine novel loci associated with BIP. Further, we identified two novel loci jointly associated with AD and BIP implicating the MARK2 gene (lead SNP rs10792421, conjunctional FDR = 0.030, same direction of effect) and the VAC14 gene (lead SNP rs11649476, conjunctional FDR = 0.022, opposite direction of effect).Conclusion: We found polygenic overlap between AD and BIP and identified novel loci for each trait and two jointly associated loci. Further studies should examine if the shared loci implicating the MARK2 and VAC14 genes could explain parts of the shared and distinct features of AD and BIP.

Technological advances make it possible to use high-throughput sequencing as a primary discovery tool of medical genetics, specifically for assaying rare variation. Still this approach faces the analytic challenge that the influence of very rare variants can only be evaluated effectively as a group. A further complication is that any given rare variant could have no effect, could increase risk, or could be protective. We propose here the C-alpha test statistic as a novel approach for testing for the presence of this mixture of effects across a set of rare variants. Unlike existing burden tests, C-alpha, by testing the variance rather than the mean, maintains consistent power when the target set contains both risk and protective variants. Through simulations and analysis of case/control data, we demonstrate good power relative to existing methods that assess the burden of rare variants in individuals.

Patterns of genetic diversity have previously been shown to mirror geography on a global scale and within continents and individual countries. Using genome-wide SNP data on 5174 Swedes with extensive geographical coverage, we analyzed the genetic structure of the Swedish population. We observed strong differences between the far northern counties and the remaining counties. The population of Dalarna county, in north middle Sweden, which borders southern Norway, also appears to differ markedly from other counties, possibly due to this county having more individuals with remote Finnish or Norwegian ancestry than other counties. An analysis of genetic differentiation (based on pairwise F(st)) indicated that the population of Sweden's southernmost counties are genetically closer to the HapMap CEU samples of Northern European ancestry than to the populations of Sweden's northernmost counties. In a comparison of extended homozygous segments, we detected a clear divide between southern and northern Sweden with small differences between the southern counties and considerably more segments in northern Sweden. Both the increased degree of homozygosity in the north and the large genetic differences between the south and the north may have arisen due to a small population in the north and the vast geographical distances between towns and villages in the north, in contrast to the more densely settled southern parts of Sweden. Our findings have implications for future genome-wide association studies (GWAS) with respect to the matching of cases and controls and the need for within-county matching. We have shown that genetic differences within a single country may be substantial, even when viewed on a European scale. Thus, population stratification needs to be accounted for, even within a country like Sweden, which is often perceived to be relatively homogenous and a favourable resource for genetic mapping, otherwise inferences based on genetic data may lead to false conclusions.

Translating a set of disease regions into insight about pathogenic mechanisms requires not only the ability to identify the key disease genes within them, but also the biological relationships among those key genes. Here we describe a statistical method, Gene Relationships Among Implicated Loci (GRAIL), that takes a list of disease regions and automatically assesses the degree of relatedness of implicated genes using 250,000 PubMed abstracts. We first evaluated GRAIL by assessing its ability to identify subsets of highly related genes in common pathways from validated lipid and height SNP associations from recent genome-wide studies. We then tested GRAIL, by assessing its ability to separate true disease regions from many false positive disease regions in two separate practical applications in human genetics. First, we took 74 nominally associated Crohn's disease SNPs and applied GRAIL to identify a subset of 13 SNPs with highly related genes. Of these, ten convincingly validated in follow-up genotyping; genotyping results for the remaining three were inconclusive. Next, we applied GRAIL to 165 rare deletion events seen in schizophrenia cases (less than one-third of which are contributing to disease risk). We demonstrate that GRAIL is able to identify a subset of 16 deletions containing highly related genes; many of these genes are expressed in the central nervous system and play a role in neuronal synapses. GRAIL offers a statistically robust approach to identifying functionally related genes from across multiple disease regions--that likely represent key disease pathways. An online version of this method is available for public use (http://www.broad.mit.edu/mpg/grail/).