Functional Effects of Schizophrenia-Linked Genetic Variants on Intrinsic Single-Neuron Excitability: A Modeling Study

Single-cell multi-cohort dissection of the schizophrenia transcriptome

The complexity and heterogeneity of schizophrenia have hindered mechanistic elucidation and the development of more effective therapies. Here, we performed single-cell dissection of schizophrenia-associated transcriptomic changes in the human prefrontal cortex across 140 individuals in two independent cohorts. Excitatory neurons were the most affected cell group, with transcriptional changes converging on neurodevelopment and synapse-related molecular pathways. Transcriptional alterations included known genetic risk factors, suggesting convergence of rare and common genomic variants on neuronal population-specific alterations in schizophrenia. Based on the magnitude of schizophrenia-associated transcriptional change, we identified two populations of individuals with schizophrenia marked by expression of specific excitatory and inhibitory neuronal cell states. This single-cell atlas links transcriptomic changes to etiological genetic risk factors, contextualizing established knowledge within the human cortical cytoarchitecture and facilitating mechanistic understanding of schizophrenia pathophysiology and heterogeneity.

Agranular frontal cortical microcircuit underlying cognitive control in macaques

The error-related negativity and an N2-component recorded over medial frontal cortex index core functions of cognitive control. While they are known to originate from agranular frontal areas, the underlying microcircuit mechanisms remain elusive. Most insights about microcircuit function have been derived from variations of the so-called canonical microcircuit model. These microcircuit architectures are based extensively on studies from granular sensory cortical areas in monkeys, cats, and rodents. However, evidence has shown striking cytoarchitectonic differences across species and differences in the functional relationships across cortical layers in agranular compared to granular sensory areas. In this minireview, we outline a tentative microcircuit model underlying cognitive control in the agranular frontal cortex of primates. The model incorporates the main GABAergic interneuron subclasses with specific laminar arrangements and target regions on pyramidal cells. We emphasize the role of layer 5 pyramidal cells in error and conflict detection. We offer several specific questions necessary for creating a specific intrinsic microcircuit model of the agranular frontal cortex.

Genetic mechanisms for impaired synaptic plasticity in schizophrenia revealed by computational modelling

Schizophrenia phenotypes are suggestive of impaired cortical plasticity in the disease, but the mechanisms of these deficits are unknown. Genomic association studies have implicated a large number of genes that regulate neuromodulation and plasticity, indicating that the plasticity deficits have a genetic origin. Here, we used biochemically detailed computational modelling of post-synaptic plasticity to investigate how schizophrenia-associated genes regulate long-term potentiation (LTP) and depression (LTD). We combined our model with data from post-mortem mRNA expression studies (CommonMind gene-expression datasets) to assess the consequences of altered expression of plasticity-regulating genes for the amplitude of LTP and LTD. Our results show that the expression alterations observed post mortem, especially those in anterior cingulate cortex, lead to impaired PKA-pathway-mediated LTP in synapses containing GluR1 receptors. We validated these findings using a genotyped EEG dataset where polygenic risk scores for synaptic and ion channel-encoding genes as well as modulation of visual evoked potentials (VEP) were determined for 286 healthy controls. Our results provide a possible genetic mechanism for plasticity impairments in schizophrenia, which can lead to improved understanding and, ultimately, treatment of the disorder.

Akkouh I, Hughes T, Frei E, A. Andreassen O, Djurovic S. Using iPSC Models to Understand the Role of Estrogen in Neuron-Glia Interactions in Schizophrenia and Bipolar Disorder

Schizophrenia (SCZ) and bipolar disorder (BIP) are severe mental disorders with a considerable disease burden worldwide due to early age of onset, chronicity, and lack of efficient treatments or prevention strategies. Whilst our current knowledge is that SCZ and BIP are highly heritable and share common pathophysiological mechanisms associated with cellular signaling, neurotransmission, energy metabolism, and neuroinflammation, the development of novel therapies has been hampered by the unavailability of appropriate models to identify novel targetable pathomechanisms. Recent data suggest that neuron-glia interactions are disturbed in SCZ and BIP, and are modulated by estrogen (E2). However, most of the knowledge we have so far on the neuromodulatory effects of E2 came from studies on animal models and human cell lines, and may not accurately reflect many processes occurring exclusively in the human brain. Thus, here we highlight the advantages of using induced pluripotent stem cell (iPSC) models to revisit studies of mechanisms underlying beneficial effects of E2 in human brain cells. A better understanding of these mechanisms opens the opportunity to identify putative targets of novel therapeutic agents for SCZ and BIP. In this review, we first summarize the literature on the molecular mechanisms involved in SCZ and BIP pathology and the beneficial effects of E2 on neuron-glia interactions. Then, we briefly present the most recent developments in the iPSC field, emphasizing the potential of using patient-derived iPSCs as more relevant models to study the effects of E2 on neuron-glia interactions.

Sex Differences in Biophysical Signatures across Molecularly Defined Medial Amygdala Neuronal Subpopulations

The medial amygdala (MeA) is essential for processing innate social and non-social behaviors, such as territorial aggression and mating, which display in a sex-specific manner. While sex differences in cell numbers and neuronal morphology in the MeA are well established, if and how these differences extend to the biophysical level remain unknown. Our previous studies revealed that expression of the transcription factors, Dbx1 and Foxp2, during embryogenesis defines separate progenitor pools destined to generate different subclasses of MEA inhibitory output neurons. We have also previously shown that Dbx1-lineage and Foxp2-lineage neurons display different responses to innate olfactory cues and in a sex-specific manner. To examine whether these neurons also possess sex-specific biophysical signatures, we conducted a multidimensional analysis of the intrinsic electrophysiological profiles of these transcription factor defined neurons in the male and female MeA. We observed striking differences in the action potential (AP) spiking patterns across lineages, and across sex within each lineage, properties known to be modified by different voltage-gated ion channels. To identify the potential mechanism underlying the observed lineage-specific and sex-specific differences in spiking adaptation, we conducted a phase plot analysis to narrow down putative ion channel candidates. Of these candidates, we found a subset expressed in a lineage-biased and/or sex-biased manner. Thus, our results uncover neuronal subpopulation and sex differences in the biophysical signatures of developmentally defined MeA output neurons, providing a potential physiological substrate for how the male and female MeA may process social and non-social cues that trigger innate behavioral responses.

The Role of Parvalbumin-positive Interneurons in Auditory Steady-State Response Deficits in Schizophrenia

Despite an increasing body of evidence demonstrating subcellular alterations in parvalbumin-positive (PV⁺) interneurons in schizophrenia, their functional consequences remain elusive. Since PV⁺ interneurons are involved in the generation of fast cortical rhythms, these changes have been hypothesized to contribute to well-established alterations of beta and gamma range oscillations in patients suffering from schizophrenia. However, the precise role of these alterations and the role of different subtypes of PV⁺ interneurons is still unclear. Here we used a computational model of auditory steady-state response (ASSR) deficits in schizophrenia. We investigated the differential effects of decelerated synaptic dynamics, caused by subcellular alterations at two subtypes of PV⁺ interneurons: basket cells and chandelier cells. Our simulations suggest that subcellular alterations at basket cell synapses rather than chandelier cell synapses are the main contributor to these deficits. Particularly, basket cells might serve as target for innovative therapeutic interventions aiming at reversing the oscillatory deficits.

Computational Modeling of Genetic Contributions to Excitability and Neural Coding in Layer V Pyramidal Cells: Applications to Schizophrenia Pathology

Pyramidal cells in layer V of the neocortex are one of the most widely studied neuron types in the mammalian brain. Due to their role as integrators of feedforward and cortical feedback inputs, they are well-positioned to contribute to the symptoms and pathology in mental disorders-such as schizophrenia-that are characterized by a mismatch between the internal perception and external inputs. In this modeling study, we analyze the input/output properties of layer V pyramidal cells and their sensitivity to modeled genetic variants in schizophrenia-associated genes. We show that the excitability of layer V pyramidal cells and the way they integrate inputs in space and time are altered by many types of variants in ion-channel and Ca²⁺ transporter-encoding genes that have been identified as risk genes by recent genome-wide association studies. We also show that the variability in the output patterns of spiking and Ca²⁺ transients in layer V pyramidal cells is altered by these model variants. Importantly, we show that many of the predicted effects are robust to noise and qualitatively similar across different computational models of layer V pyramidal cells. Our modeling framework reveals several aspects of single-neuron excitability that can be linked to known schizophrenia-related phenotypes and existing hypotheses on disease mechanisms. In particular, our models predict that single-cell steady-state firing rate is positively correlated with the coding capacity of the neuron and negatively correlated with the amplitude of a prepulse-mediated adaptation and sensitivity to coincidence of stimuli in the apical dendrite and the perisomatic region of a layer V pyramidal cell. These results help to uncover the voltage-gated ion-channel and Ca²⁺ transporter-associated genetic underpinnings of schizophrenia phenotypes and biomarkers.

Machine learning in schizophrenia genomics, a case-control study using 5,090 exomes

Our hypothesis is that machine learning (ML) analysis of whole exome sequencing (WES) data can be used to identify individuals at high risk for schizophrenia (SCZ). This study applies ML to WES data from 2,545 individuals with SCZ and 2,545 unaffected individuals, accessed via the database of genotypes and phenotypes (dbGaP). Single nucleotide variants and small insertions and deletions were annotated by ANNOVAR using the reference genome hg19/GRCh37. Rare (predicted functional) variants with a minor allele frequency ≤1% and genotype quality ≥90 including missense, frameshift, stop gain, stop loss, intronic, and exonic splicing variants were selected. A file containing all cases and controls, the names of genes with variants meeting our criteria, and the number of variants per gene for each individual, was used for ML analysis. The supervised machine-learning algorithm used the patterns of variants observed in the different genes to determine which subset of genes can best predict that an individual is affected. Seventy percent of the data was used to train the algorithm and the remaining 30% of data (n = 1,526) was used to evaluate its efficiency. The supervised ML algorithm, gradient boosted trees with regularization (eXtreme Gradient Boosting implementation) was the best performing algorithm yielding promising results (accuracy: 85.7%, specificity: 86.6%, sensitivity: 84.9%, area under the receiver-operator characteristic curve: 0.95). The top 50 features (genes) of the algorithm were analyzed using bioinformatics resources for new insights about the pathophysiology of SCZ. This manuscript presents a novel predictor which could potentially enable studies exploring disease-modifying intervention in the early stages of the disease.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling

Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na⁺ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca²⁺ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.

A stepwise neuron model fitting procedure designed for recordings with high spatial resolution: Application to layer 5 pyramidal cells

BACKGROUND

Recent progress in electrophysiological and optical methods for neuronal recordings provides vast amounts of high-resolution data. In parallel, the development of computer technology has allowed simulation of ever-larger neuronal circuits. A challenge in taking advantage of these developments is the construction of single-cell and network models in a way that faithfully reproduces neuronal biophysics with subcellular level of details while keeping the simulation costs at an acceptable level.

NEW METHOD

In this work, we develop and apply an automated, stepwise method for fitting a neuron model to data with fine spatial resolution, such as that achievable with voltage sensitive dyes (VSDs) and Ca²⁺ imaging.

RESULT

We apply our method to simulated data from layer 5 pyramidal cells (L5PCs) and construct a model with reduced neuronal morphology. We connect the reduced-morphology neurons into a network and validate against simulated data from a high-resolution L5PC network model.

COMPARISON WITH EXISTING METHODS

Our approach combines features from several previously applied model-fitting strategies. The reduced-morphology neuron model obtained using our approach reliably reproduces the membrane-potential dynamics across the dendrites as predicted by the full-morphology model.

CONCLUSIONS

The network models produced using our method are cost-efficient and predict that interconnected L5PCs are able to amplify delta-range oscillatory inputs across a large range of network sizes and topologies, largely due to the medium after hyperpolarization mediated by the Ca²⁺-activated SK current.

Genetic evidence for role of integration of fast and slow neurotransmission in schizophrenia

The most recent genome-wide association studies (GWAS) of schizophrenia (SCZ) identified hundreds of risk variants potentially implicated in the disease. Further, novel statistical methodology designed for polygenic architecture revealed more potential risk variants. This can provide a link between individual genetic factors and the mechanistic underpinnings of SCZ. Intriguingly, a large number of genes coding for ionotropic and metabotropic receptors for various neurotransmitters-glutamate, γ-aminobutyric acid (GABA), dopamine, serotonin, acetylcholine and opioids-and numerous ion channels were associated with SCZ. Here, we review these findings from the standpoint of classical neurobiological knowledge of neuronal synaptic transmission and regulation of electrical excitability. We show that a substantial proportion of the identified genes are involved in intracellular cascades known to integrate 'slow' (G-protein-coupled receptors) and 'fast' (ionotropic receptors) neurotransmission converging on the protein DARPP-32. Inspection of the Human Brain Transcriptome Project database confirms that that these genes are indeed expressed in the brain, with the expression profile following specific developmental trajectories, underscoring their relevance to brain organization and function. These findings extend the existing pathophysiology hypothesis by suggesting a unifying role of dysregulation in neuronal excitability and synaptic integration in SCZ. This emergent model supports the concept of SCZ as an 'associative' disorder-a breakdown in the communication across different slow and fast neurotransmitter systems through intracellular signaling pathways-and may unify a number of currently competing hypotheses of SCZ pathophysiology.

Genetic tests in major psychiatric disorders-integrating molecular medicine with clinical psychiatry-why is it so difficult?

With the advent of post-genomic era, new technologies create extraordinary possibilities for diagnostics and personalized therapy, transforming todays' medicine. Rooted in both medical genetics and clinical psychiatry, the paper is designed as an integrated source of information of the current and potential future application of emerging genomic technologies as diagnostic tools in psychiatry, moving beyond the classical concept of patient approach. Selected approaches are presented, starting from currently used technologies (next-generation sequencing (NGS) and microarrays), followed by newer options (reverse phenotyping). Next, we describe an old concept in a new light (endophenotypes), subsequently coming up with a sophisticated and complex approach (gene networks) ending by a nascent field (computational psychiatry). The challenges and barriers that exist to translate genomic research to real-world patient assessment are further discussed. We emphasize the view that only a paradigm shift can bring a fundamental change in psychiatric practice, allowing to disentangle the intricacies of mental diseases. All the diagnostic methods, as described, are directed at uncovering the integrity of the system including many types of relations within a complex structure. The integrative system approach offers new opportunity to connect genetic background with specific diseases entities, or concurrently, with symptoms regardless of a diagnosis. To advance the field, we propose concerted cross-disciplinary effort to provide a diagnostic platform operating at the general level of genetic pathogenesis of complex-trait psychiatric disorders rather than at the individual level of a specific disease.

Multifactorial Modeling of Impairment of Evoked Gamma Range Oscillations in Schizophrenia

Despite a significant increase in efforts to identify biomarkers and endophenotypic measures of psychiatric illnesses, only a very limited amount of computational models of these markers and measures has been implemented so far. Moreover, existing computational models dealing with biomarkers typically only examine one possible mechanism in isolation, disregarding the possibility that other combinations of model parameters might produce the same network behavior (what has been termed "multifactoriality"). In this study we describe a step toward a computational instantiation of an endophenotypic finding for schizophrenia, namely the impairment of evoked auditory gamma and beta oscillations in schizophrenia. We explore the multifactorial nature of this impairment using an established model of primary auditory cortex, by performing an extensive search of the parameter space. We find that single network parameters contain only little information about whether the network will show impaired gamma entrainment and that different regions in the parameter space yield similar network level oscillation abnormalities. These regions in the parameter space, however, show strong differences in the underlying network dynamics. To sum up, we present a first step toward an in silico instantiation of an important biomarker of schizophrenia, which has great potential for the identification and study of disease mechanisms and for understanding of existing treatments and development of novel ones.