Genome-wide CRISPR screen in a mouse model of tumor growth and metastasis.

Genetic screens are powerful tools for identifying genes responsible for diverse phenotypes. Here we describe a genome-wide CRISPR/Cas9-mediated loss-of-function screen in tumor growth and metastasis. We mutagenized a non-metastatic mouse cancer cell line using a genome-scale library with 67,405 single-guide RNAs (sgRNAs). The mutant cell pool rapidly generates metastases when transplanted into immunocompromised mice. Enriched sgRNAs in lung metastases and late-stage primary tumors were found to target a small set of genes, suggesting that specific loss-of-function mutations drive tumor growth and metastasis. Individual sgRNAs and a small pool of 624 sgRNAs targeting the top-scoring genes from the primary screen dramatically accelerate metastasis. In all of these experiments, the effect of mutations on primary tumor growth positively correlates with the development of metastases. Our study demonstrates Cas9-based screening as a robust method to systematically assay gene phenotypes in cancer evolution in vivo.

Heterozygous deletion of the autism-associated gene CHD8 impairs synaptic function through widespread changes in gene expression and chromatin compaction.

Whole-exome sequencing of autism spectrum disorder (ASD) probands and unaffected family members has identified many genes harboring de novo variants suspected to play a causal role in the disorder. Of these, chromodomain helicase DNA-binding protein 8 (CHD8) is the most recurrently mutated. Despite the prevalence of CHD8 mutations, we have little insight into how CHD8 loss affects genome organization or the functional consequences of these molecular alterations in neurons. Here, we engineered two isogenic human embryonic stem cell lines with CHD8 loss-of-function mutations and characterized differences in differentiated human cortical neurons. We identified hundreds of genes with altered expression, including many involved in neural development and excitatory synaptic transmission. Field recordings and single-cell electrophysiology revealed a 3-fold decrease in firing rates and synaptic activity in CHD8+/- neurons, as well as a similar firing-rate deficit in primary cortical neurons from Chd8+/- mice. These alterations in neuron and synapse function can be reversed by CHD8 overexpression. Moreover, CHD8+/- neurons displayed a large increase in open chromatin across the genome, where the greatest change in compaction was near autism susceptibility candidate 2 (AUTS2), which encodes a transcriptional regulator implicated in ASD. Genes with changes in chromatin accessibility and expression in CHD8+/- neurons have significant overlap with genes mutated in probands for ASD, intellectual disability, and schizophrenia but not with genes mutated in healthy controls or other disease cohorts. Overall, this study characterizes key molecular alterations in genome structure and expression in CHD8+/- neurons and links these changes to impaired neuronal and synaptic function.

A Potential Source of Bias in Group-Level EEG Microstate Analysis.

Microstate analysis is a promising technique for analyzing high-density electroencephalographic data, but there are multiple questions about methodological best practices. Between and within individuals, microstates can differ both in terms of characteristic topographies and temporal dynamics, which leads to analytic challenges as the measurement of microstate dynamics is dependent on assumptions about their topographies. Here we focus on the analysis of group differences, using simulations seeded on real data from healthy control subjects to compare approaches that derive separate sets of maps within subgroups versus a single set of maps applied uniformly to the entire dataset. In the absence of true group differences in either microstate maps or temporal metrics, we found that using separate subgroup maps resulted in substantially inflated type I error rates. On the other hand, when groups truly differed in their microstate maps, analyses based on a single set of maps confounded topographic effects with differences in other derived metrics. We propose an approach to alleviate both classes of bias, based on a paired analysis of all subgroup maps. We illustrate the qualitative and quantitative impact of these issues in real data by comparing waking versus non-rapid eye movement sleep microstates. Overall, our results suggest that even subtle chance differences in microstate topography can have profound effects on derived microstate metrics and that future studies using microstate analysis should take steps to mitigate this large source of error.

Study Protocol: Global Research Initiative on the Neurophysiology of Schizophrenia (GRINS) project.

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.

Analysing an allelic series of rare missense variants of CACNA1I in a Swedish schizophrenia cohort.

CACNA1I is implicated in the susceptibility to schizophrenia by large-scale genetic association studies of single nucleotide polymorphisms. However, the channelopathy of CACNA1I in schizophrenia is unknown. CACNA1I encodes CaV3.3, a neuronal voltage-gated calcium channel that underlies a subtype of T-type current that is important for neuronal excitability in the thalamic reticular nucleus and other regions of the brain. Here, we present an extensive functional characterization of 57 naturally occurring rare and common missense variants of CACNA1I derived from a Swedish schizophrenia cohort of more than 10 000 individuals. Our analysis of this allelic series of coding CACNA1I variants revealed that reduced CaV3.3 channel current density was the dominant phenotype associated with rare CACNA1I coding alleles derived from control subjects, whereas rare CACNA1I alleles from schizophrenia patients encoded CaV3.3 channels with altered responses to voltages. CACNA1I variants associated with altered current density primarily impact the ionic channel pore and those associated with altered responses to voltage impact the voltage-sensing domain. CaV3.3 variants associated with altered voltage dependence of the CaV3.3 channel and those associated with peak current density deficits were significantly segregated across affected and unaffected groups (Fisher's exact test, P = 0.034). Our results, together with recent data from the SCHEMA (Schizophrenia Exome Sequencing Meta-Analysis) cohort, suggest that reduced CaV3.3 function may protect against schizophrenia risk in rare cases. We subsequently modelled the effect of the biophysical properties of CaV3.3 channel variants on thalamic reticular nucleus excitability and found that compared with common variants, ultrarare CaV3.3-coding variants derived from control subjects significantly decreased thalamic reticular nucleus excitability (P = 0.011). When all rare variants were analysed, there was a non-significant trend between variants that reduced thalamic reticular nucleus excitability and variants that either had no effect or increased thalamic reticular nucleus excitability across disease status. Taken together, the results of our functional analysis of an allelic series of >50 CACNA1I variants in a schizophrenia cohort reveal that loss of function of CaV3.3 is a molecular phenotype associated with reduced disease risk burden, and our approach may serve as a template strategy for channelopathies in polygenic disorders.

Author Correction: Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

The original version of this paper contained an incorrect primer sequence. In the Methods subsection "Rampage libraries," the text for modification 3 stated that the reverse primer used for library indexing was 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGT-3'. The correct sequence of the oligonucleotide used is 5'-CAAGCAGAAGACGGCATACGAGATXXXXXXXXGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT-3'. This error has been corrected in the PDF and HTML versions of the paper.

Comprehensive comparative analysis of 5'-end RNA-sequencing methods.

Specialized RNA-seq methods are required to identify the 5' ends of transcripts, which are critical for studies of gene regulation, but these methods have not been systematically benchmarked. We directly compared six such methods, including the performance of five methods on a single human cellular RNA sample and a new spike-in RNA assay that helps circumvent challenges resulting from uncertainties in annotation and RNA processing. We found that the 'cap analysis of gene expression' (CAGE) method performed best for mRNA and that most of its unannotated peaks were supported by evidence from other genomic methods. We applied CAGE to eight brain-related samples and determined sample-specific transcription start site (TSS) usage, as well as a transcriptome-wide shift in TSS usage between fetal and adult brain.

Overexpression of NEUROG2 and NEUROG1 in human embryonic stem cells produces a network of excitatory and inhibitory neurons.

Overexpression of mouse neurogenin ( Neurog) 2 alone or in combination with mouse Neurog2/1 in human embryonic stem cells (hESCs) and human induced pluripotent stem cells (hiPSCs) can rapidly produce high-yield excitatory neurons. Here, we report a detailed characterization of human neuronal networks induced by the expression of human NEUROG2 together with human NEUROG2/1 in hESCs using molecular, cellular, and electrophysiological measurements over 60 d after induction. Both excitatory synaptic transmission and network firing activity increased over time. Strikingly, inhibitory synaptic transmission and GABAergic cells were identified from NEUROG2/1 induced neurons (iNs). To illustrate the application of such iNs, we demonstrated that the heterozygous knock out of SCN2A, whose loss-of-function mutation is strongly implicated in autism risk, led to a dramatic reduction in network activity in the NEUROG2/1 iNs. Our findings not only extend our understanding of the NEUROG2/1-induced human neuronal network but also substantiate NEUROG2/1 iNs as an in vitro system for modeling neuronal and functional deficits on a human genetic background.-Lu, C., Shi, X., Allen, A., Baez-Nieto, D., Nikish, A., Sanjana, N. E., Pan, J. Q. Overexpression of NEUROG2 and NEUROG1 in human embryonic stem cells produces a network of excitatory and inhibitory neurons.

Biological concepts in human sodium channel epilepsies and their relevance in clinical practice.

OBJECTIVE: Voltage-gated sodium channels (SCNs) share similar amino acid sequence, structure, and function. Genetic variants in the four human brain-expressed SCN genes SCN1A/2A/3A/8A have been associated with heterogeneous epilepsy phenotypes and neurodevelopmental disorders. To better understand the biology of seizure susceptibility in SCN-related epilepsies, our aim was to determine similarities and differences between sodium channel disorders, allowing us to develop a broader perspective on precision treatment than on an individual gene level alone. METHODS: We analyzed genotype-phenotype correlations in large SCN-patient cohorts and applied variant constraint analysis to identify severe sodium channel disease. We examined temporal patterns of human SCN expression and correlated functional data from in vitro studies with clinical phenotypes across different sodium channel disorders. RESULTS: Comparing 865 epilepsy patients (504 SCN1A, 140 SCN2A, 171 SCN8A, four SCN3A, 46 copy number variation [CNV] cases) and analysis of 114 functional studies allowed us to identify common patterns of presentation. All four epilepsy-associated SCN genes demonstrated significant constraint in both protein truncating and missense variation when compared to other SCN genes. We observed that age at seizure onset is related to SCN gene expression over time. Individuals with gain-of-function SCN2A/3A/8A missense variants or CNV duplications share similar characteristics, most frequently present with early onset epilepsy (<3 months), and demonstrate good response to sodium channel blockers (SCBs). Direct comparison of corresponding SCN variants across different SCN subtypes illustrates that the functional effects of variants in corresponding channel locations are similar; however, their clinical manifestation differs, depending on their role in different types of neurons in which they are expressed. SIGNIFICANCE: Variant function and location within one channel can serve as a surrogate for variant effects across related sodium channels. Taking a broader view on precision treatment suggests that in those patients with a suspected underlying genetic epilepsy presenting with neonatal or early onset seizures (<3 months), SCBs should be considered.

Tracks through the genome to physiological events.

NEW FINDINGS: What is the topic of this review? We discuss tools available to access genome-wide data sets that harbour cell-specific, brain region-specific and tissue-specific information on exon usage for several species, including humans. In this Review, we demonstrate how to access this information in genome databases and its enormous value to physiology. What advances does it highlight? The sheer scale of protein diversity that is possible from complex genes, including those that encode voltage-gated ion channels, is vast. But this choice is critical for a complete understanding of protein function in the most physiologically relevant context. Many proteins of great interest to physiologists and neuroscientists are structurally complex and located in specialized subcellular domains, such as neuronal synapses and transverse tubules of muscle. Genes that encode these critical signalling molecules (receptors, ion channels, transporters, enzymes, cell adhesion molecules, cell-cell interaction proteins and cytoskeletal proteins) are similarly complex. Typically, these genes are large; human Dystrophin (DMD) encodes a cytoskeletal protein of muscle and it is the largest naturally occurring gene at a staggering 2.3 Mb. Large genes contain many non-coding introns and coding exons; human Titin (TTN), which encodes a protein essential for the assembly and functioning of vertebrate striated muscles, has over 350 exons and consequently has an enormous capacity to generate different forms of Titin mRNAs that have unique exon combinations. Functional and pharmacological differences among protein isoforms originating from the same gene may be subtle but nonetheless of critical physiological significance. Standard functional, immunological and pharmacological approaches, so useful for characterizing proteins encoded by different genes, typically fail to discriminate among splice isoforms of individual genes. Tools are now available to access genome-wide data sets that harbour cell-specific, brain region-specific and tissue-specific information on exon usage for several species, including humans. In this Review, we demonstrate how to access this information in genome databases and its enormous value to physiology.

Electroencephalographic Microstates During Sleep and Wake in Schizophrenia.

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.

EEG microstates are stereotyped patterns of scalp voltage topography that provide information about brain activity at millisecond-level temporal resolution. We used this method to study brain activity in schizophrenia during sleep and wake. We found state-dependent case-control differences in EEG microstates that were unrelated to the results of classic EEG analyses. These differences reflect aberrant neural functioning during sleep in patients with schizophrenia.

Differential functional consequences of GRIN2A mutations associated with schizophrenia and neurodevelopmental disorders.

Human genetic studies have revealed rare missense and protein-truncating variants in GRIN2A, encoding for the GluN2A subunit of the NMDA receptors, that confer significant risk for schizophrenia (SCZ). Mutations in GRIN2A are also associated with epilepsy and developmental delay/intellectual disability (DD/ID). However, it remains enigmatic how alterations to the same protein can result in diverse clinical phenotypes. Here, we performed functional characterization of human GluN1/GluN2A heteromeric NMDA receptors that contain SCZ-linked GluN2A variants, and compared them to NMDA receptors with GluN2A variants associated with epilepsy or DD/ID. Our findings demonstrate that SCZ-associated GRIN2A variants were predominantly loss-of-function (LoF), whereas epilepsy and DD/ID-associated variants resulted in both gain- and loss-of-function phenotypes. We additionally show that M653I and S809R, LoF GRIN2A variants associated with DD/ID, exert a dominant-negative effect when co-expressed with a wild-type GluN2A, whereas E58Ter and Y698C, SCZ-linked LoF variants, and A727T, an epilepsy-linked LoF variant, do not. These data offer a potential mechanism by which SCZ/epilepsy and DD/ID-linked variants can cause different effects on receptor function and therefore result in divergent pathological outcomes.

TRPA5 encodes a thermosensitive ankyrin ion channel receptor in a triatomine insect.

As ectotherms, insects need heat-sensitive receptors to monitor environmental temperatures and facilitate thermoregulation. We show that TRPA5, a class of ankyrin transient receptor potential (TRP) channels absent in dipteran genomes, may function as insect heat receptors. In the triatomine bug Rhodnius prolixus (order: Hemiptera), a vector of Chagas disease, the channel RpTRPA5B displays a uniquely high thermosensitivity, with biophysical determinants including a large channel activation enthalpy change (72 kcal/mol), a high temperature coefficient (Q10 = 25), and in vitro temperature-induced currents from 53°C to 68°C (T0.5 = 58.6°C), similar to noxious TRPV receptors in mammals. Monomeric and tetrameric ion channel structure predictions show reliable parallels with fruit fly dTRPA1, with structural uniqueness in ankyrin repeat domains, the channel selectivity filter, and potential TRP functional modulator regions. Overall, the finding of a member of TRPA5 as a temperature-activated receptor illustrates the diversity of insect molecular heat detectors.

A spectrum of altered non-rapid eye movement sleep in schizophrenia.

Multiple facets of sleep neurophysiology, including electroencephalography (EEG) metrics such as non-rapid eye movement (NREM) spindles and slow oscillations, are altered in individuals with schizophrenia (SCZ). However, beyond group-level analyses, the extent to which NREM deficits vary among patients is unclear, as are their relationships to other sources of heterogeneity including clinical factors, ageing, cognitive profiles and medication regimens. Using newly collected high-density sleep EEG data on 103 individuals with SCZ and 68 controls, we first sought to replicate our previously reported group-level differences between patients and controls (original N=130) during N2 stage. Then in the combined sample (N=301 including 175 patients), we characterized patient-to-patient variability. We replicated all group-level mean differences and confirmed the high accuracy of our predictive model (AUC=0.93 for diagnosis). Compared to controls, patients showed significantly increased between-individual variability across many (26%) sleep metrics. Although multiple clinical and cognitive factors were associated with NREM metrics, collectively they did not account for much of the general increase in patient-to-patient variability. Medication regimen was a greater contributor to variability. Some sleep metrics including fast spindle density showed exaggerated age-related effects in SCZ, and patients exhibited older predicted biological ages based on the sleep EEG; further, among patients, certain medications exacerbated these effects, in particular olanzapine. Collectively, our results point to a spectrum of N2 sleep deficits among SCZ patients that can be measured objectively and at scale, with relevance to both the etiological heterogeneity of SCZ as well as potential iatrogenic effects of antipsychotic medication.

IntelliSleepScorer, a software package with a graphic user interface for automated sleep stage scoring in mice based on a light gradient boosting machine algorithm.

Machine learning has been applied in recent years to categorize sleep stages (NREM, REM, and wake) using electroencephalogram (EEG) recordings; however, a well-validated sleep scoring automatic pipeline in rodent research is still not publicly available. Here, we present IntelliSleepScorer, a software package with a graphic user interface to score sleep stages automatically in mice. IntelliSleepScorer uses the light gradient boosting machine (LightGBM) to score sleep stages for each epoch of recordings. We developed LightGBM models using a large cohort of data, which consisted of 5776 h of sleep EEG and electromyogram (EMG) signals across 519 unique recordings from 124 mice. The LightGBM model achieved an overall accuracy of 95.2% and a Cohen's kappa of 0.91, which outperforms the baseline models such as the logistic regression model (accuracy = 93.3%, kappa = 0.88) and the random forest model (accuracy = 94.3%, kappa = 0.89). The overall performance of the LightGBM model as well as the performance across different sleep stages are on par with that of the human experts. Most importantly, we validated the generalizability of the LightGBM models: (1) The LightGBM model performed well on two publicly available, independent datasets (kappa > = 0.80), which have different sampling frequency and epoch lengths; (2) The LightGBM model performed well on data recorded at a lower sampling frequency (kappa = 0.90); (3) The performance of the LightGBM model is not affected by the light/dark cycle; and (4) A modified LightGBM model performed well on data containing only one EEG and one EMG electrode (kappa > = 0.89). Taken together, the LightGBM models offer state-of-the-art performance for automatic sleep stage scoring in mice. Last, we implemented the IntelliSleepScorer software package based on the validated model to provide an out-of-box solution to sleep researchers (available for download at https://sites.broadinstitute.org/pan-lab/resources ).

Mouse mutants in schizophrenia risk genes GRIN2A and AKAP11 show EEG abnormalities in common with schizophrenia patients.

Schizophrenia is a heterogeneous psychiatric disorder with a strong genetic basis, whose etiology and pathophysiology remain poorly understood. Exome sequencing studies have uncovered rare, loss-of-function variants that greatly increase risk of schizophrenia [1], including loss-of-function mutations in GRIN2A (aka GluN2A or NR2A, encoding the NMDA receptor subunit 2A) and AKAP11 (A-Kinase Anchoring Protein 11). AKAP11 and GRIN2A mutations are also associated with bipolar disorder [2], and epilepsy and developmental delay/intellectual disability [1, 3, 4], respectively. Accessible in both humans and rodents, electroencephalogram (EEG) recordings offer a window into brain activity and display abnormal features in schizophrenia patients. Does loss of Grin2a or Akap11 in mice also result in EEG abnormalities? We monitored EEG in heterozygous and homozygous knockout Grin2a and Akap11 mutant mice compared with their wild-type littermates, at 3- and 6-months of age, across the sleep/wake cycle and during auditory stimulation protocols. Grin2a and Akap11 mutants exhibited increased resting gamma power, attenuated auditory steady-state responses (ASSR) at gamma frequencies, and reduced responses to unexpected auditory stimuli during mismatch negativity (MMN) tests. Sleep spindle density was reduced in a gene dose-dependent manner in Akap11 mutants, whereas Grin2a mutants showed increased sleep spindle density. The EEG phenotypes of Grin2a and Akap11 mutant mice show a variety of abnormal features that overlap considerably with human schizophrenia patients, reflecting systems-level changes caused by Grin2a and Akap11 deficiency. These neurophysiologic findings further substantiate Grin2a and Akap11 mutants as genetic models of schizophrenia and identify potential biomarkers for stratification of schizophrenia patients.

Scanning mutagenesis of the voltage-gated sodium channel NaV1.2 using base editing.

It is challenging to apply traditional mutational scanning to voltage-gated sodium channels (NaVs) and functionally annotate the large number of coding variants in these genes. Using a cytosine base editor and a pooled viability assay, we screen a library of 368 guide RNAs (gRNAs) tiling NaV1.2 to identify more than 100 gRNAs that change NaV1.2 function. We sequence base edits made by a subset of these gRNAs to confirm specific variants that drive changes in channel function. Electrophysiological characterization of these channel variants validates the screen results and provides functional mechanisms of channel perturbation. Most of the changes caused by these gRNAs are classifiable as loss of function along with two missense mutations that lead to gain of function in NaV1.2 channels. This two-tiered strategy to functionally characterize ion channel protein variants at scale identifies a large set of loss-of-function mutations in NaV1.2.

Brain-region-specific changes in neurons and glia and dysregulation of dopamine signaling in Grin2a mutant mice.

A genetically valid animal model could transform our understanding of schizophrenia (SCZ) disease mechanisms. Rare heterozygous loss-of-function (LoF) mutations in GRIN2A, encoding a subunit of the NMDA receptor, greatly increase the risk of SCZ. By transcriptomic, proteomic, and behavioral analyses, we report that heterozygous Grin2a mutant mice show (1) large-scale gene expression changes across multiple brain regions and in neuronal (excitatory and inhibitory) and non-neuronal cells (astrocytes and oligodendrocytes), (2) evidence of hypoactivity in the prefrontal cortex (PFC) and hyperactivity in the hippocampus and striatum, (3) an elevated dopamine signaling in the striatum and hypersensitivity to amphetamine-induced hyperlocomotion (AIH), (4) altered cholesterol biosynthesis in astrocytes, (5) a reduction in glutamatergic receptor signaling proteins in the synapse, and (6) an aberrant locomotor pattern opposite of that induced by antipsychotic drugs. These findings reveal potential pathophysiologic mechanisms, provide support for both the "hypo-glutamate" and "hyper-dopamine" hypotheses of SCZ, and underscore the utility of Grin2a-deficient mice as a genetic model of SCZ.

A spectrum of altered non-rapid eye movement sleep in schizophrenia.

Background: Multiple facets of sleep neurophysiology, including electroencephalography (EEG) metrics such as non-rapid eye movement (NREM) spindles and slow oscillations (SO), are altered in individuals with schizophrenia (SCZ). However, beyond group-level analyses which treat all patients as a unitary set, the extent to which NREM deficits vary among patients is unclear, as are their relationships to other sources of heterogeneity including clinical factors, illness duration and ageing, cognitive profiles and medication regimens. Using newly collected high density sleep EEG data on 103 individuals with SCZ and 68 controls, we first sought to replicate our previously reported (Kozhemiako et. al, 2022) group-level mean differences between patients and controls (original N=130). Then in the combined sample (N=301 including 175 patients), we characterized patient-to-patient variability in NREM neurophysiology. Results: We replicated all group-level mean differences and confirmed the high accuracy of our predictive model (Area Under the ROC Curve, AUC = 0.93 for diagnosis). Compared to controls, patients showed significantly increased between-individual variability across many (26%) sleep metrics, with patterns only partially recapitulating those for group-level mean differences. Although multiple clinical and cognitive factors were associated with NREM metrics including spindle density, collectively they did not account for much of the general increase in patient-to-patient variability. Medication regimen was a greater (albeit still partial) contributor to variability, although original group mean differences persisted after controlling for medications. Some sleep metrics including fast spindle density showed exaggerated age-related effects in SCZ, and patients exhibited older predicted biological ages based on an independent model of ageing and the sleep EEG. Conclusion: We demonstrated robust and replicable alterations in sleep neurophysiology in individuals with SCZ and highlighted distinct patterns of effects contrasting between-group means versus within-group variances. We further documented and controlled for a major effect of medication use, and pointed to greater age-related change in NREM sleep in patients. That increased NREM heterogeneity was not explained by standard clinical or cognitive patient assessments suggests the sleep EEG provides novel, nonredundant information to support the goals of personalized medicine. Collectively, our results point to a spectrum of NREM sleep deficits among SCZ patients that can be measured objectively and at scale, and that may offer a unique window on the etiological and genetic diversity that underlies SCZ risk, treatment response and prognosis.

Sources of Variation in the Spectral Slope of the Sleep EEG.

The 1/f spectral slope of the electroencephalogram (EEG) estimated in the γ frequency range has been proposed as an arousal marker that differentiates wake, nonrapid eye movement (NREM) sleep, and rapid eye movement (REM) sleep. Here, we sought to replicate and extend these findings in a large sample, providing a comprehensive characterization of how slope changes with age, sex, and its test-retest reliability as well as potential confounds that could affect the slope estimation. We used 10,255 whole-night polysomnograms (PSGs) from the National Sleep Research Resource (NSRR). All preprocessing steps were performed using an open-source Luna package and the spectral slope was estimated by fitting log-log linear regression models on the absolute power from 30 to 45 Hz separately for wake, NREM, and REM stages. We confirmed that the mean spectral slope grows steeper going from wake to NREM to REM sleep. We found that the choice of mastoid referencing scheme modulated the extent to which electromyogenic, or electrocardiographic artifacts were likely to bias 30- to 45-Hz slope estimates, as well as other sources of technical, device-specific bias. Nonetheless, within individuals, slope estimates were relatively stable over time. Both cross-sectionally and longitudinal, slopes tended to become shallower with increasing age, particularly for REM sleep; males tended to show flatter slopes than females across all states. Our findings support that spectral slope can be a valuable arousal marker for both clinical and research endeavors but also underscore the importance of considering interindividual variation and multiple methodological aspects related to its estimation.