Better together: Non-additive interactions between schizophrenia risk genes.

Zhang, Zhang, Forrest et al combine allele-specific open chromatin (ASoC) mapping and CRISPR-editing to evaluate the functional impact of schizophrenia risk variants on human neuronal gene expression, synaptic development, and function. In doing so, they uncover surprising non-additive effects between target genes regulated by the same risk variant.

Modeling gene × environment interactions in PTSD using human neurons reveals diagnosis-specific glucocorticoid-induced gene expression.

Post-traumatic stress disorder (PTSD) can develop following severe trauma, but the extent to which genetic and environmental risk factors contribute to individual clinical outcomes is unknown. Here, we compared transcriptional responses to hydrocortisone exposure in human induced pluripotent stem cell (hiPSC)-derived glutamatergic neurons and peripheral blood mononuclear cells (PBMCs) from combat veterans with PTSD (n = 19 hiPSC and n = 20 PBMC donors) and controls (n = 20 hiPSC and n = 20 PBMC donors). In neurons only, we observed diagnosis-specific glucocorticoid-induced changes in gene expression corresponding with PTSD-specific transcriptomic patterns found in human postmortem brains. We observed glucocorticoid hypersensitivity in PTSD neurons, and identified genes that contribute to this PTSD-dependent glucocorticoid response. We find evidence of a coregulated network of transcription factors that mediates glucocorticoid hyper-responsivity in PTSD. These findings suggest that induced neurons represent a platform for examining the molecular mechanisms underlying PTSD, identifying biomarkers of stress response, and conducting drug screening to identify new therapeutics.

Attenuated transcriptional response to pro-inflammatory cytokines in schizophrenia hiPSC-derived neural progenitor cells.

Maternal immune activation (MIA) during prenatal development is an environmental risk factor for psychiatric disorders including schizophrenia (SZ). Converging lines of evidence from human and animal model studies suggest that elevated cytokine levels in the maternal and fetal compartments are an important indication of the mechanisms driving this association. However, there is variability in susceptibility to the psychiatric risk conferred by MIA, likely influenced by genetic factors. How MIA interacts with a genetic profile susceptible to SZ is challenging to test in animal models. To address this gap, we examined whether differential gene expression responses occur in forebrain-lineage neural progenitor cells (NPCs) derived from human induced pluripotent stem cells (hiPSC) generated from three individuals with a diagnosis of schizophrenia and three healthy controls. Following acute (24 h) treatment with either interferon-gamma (IFNγ; 25 ng/µl) or interleukin (IL)-1ß (10 ng/µl), we identified, by RNA sequencing, 3380 differentially expressed genes (DEGs) in the IFNγ-treated control lines (compared to untreated controls), and 1980 DEGs in IFNγ-treated SZ lines (compared to untreated SZ lines). Out of 4137 genes that responded significantly to IFNγ across all lines, 1223 were common to both SZ and control lines. The 2914 genes that appeared to respond differentially to IFNγ treatment in SZ lines were subjected to a further test of significance (multiple testing correction applied to the interaction effect between IFNγ treatment and SZ diagnosis), yielding 359 genes that passed the significance threshold. There were no differentially expressed genes in the IL-1ß-treatment conditions after Benjamini-Hochberg correction. Gene set enrichment analysis however showed that IL-1ß impacts immune function and neuronal differentiation. Overall, our data suggest that a) SZ NPCs show an attenuated transcriptional response to IFNγ treatment compared to controls; b) Due to low IL-1ß receptor expression in NPCs, NPC cultures appear to be less responsive to IL-1ß than IFNγ; and c) the genes differentially regulated in SZ lines - in the face of a cytokine challenge - are primarily associated with mitochondrial, "loss-of-function", pre- and post-synaptic gene sets. Our findings particularly highlight the role of early synaptic development in the association between maternal immune activation and schizophrenia risk.

Neuroligin-3 and neuroligin-4X form nanoscopic clusters and regulate growth cone organization and size.

The cell-adhesion proteins neuroligin-3 and neuroligin-4X (NLGN3/4X) have well described roles in synapse formation. NLGN3/4X are also expressed highly during neurodevelopment. However, the role these proteins play during this period is unknown. Here we show that NLGN3/4X localized to the leading edge of growth cones where it promoted neuritogenesis in immature human neurons. Super-resolution microscopy revealed that NLGN3/4X clustering induced growth cone enlargement and influenced actin filament organization. Critically, these morphological effects were not induced by autism spectrum disorder (ASD)-associated NLGN3/4X variants. Finally, actin regulators p21-activated kinase 1 and cofilin were found to be activated by NLGN3/4X and involved in mediating the effects of these adhesion proteins on actin filaments, growth cones and neuritogenesis. These data reveal a novel role for NLGN3 and NLGN4X in the development of neuronal architecture, which may be altered in the presence of ASD-associated variants.

Using the dCas9-KRAB system to repress gene expression in hiPSC-derived NGN2 neurons.

We describe a CRISPR inhibition (CRISPRi) protocol to repress endogenous gene expression (e.g., ATP6V1A) in human induced pluripotent stem cell-derived NGN2-induced glutamatergic neurons. CRISPRi enables efficient and precise gene repression of one or multiple target genes via delivering gRNA(s) to direct a dCas9-KRAB fusion protein to the gene(s) of interest. This protocol can also be adapted for gene activation and high-throughput gene manipulation, allowing assessment of the transcriptomic and phenotypic impact of candidate gene(s) associated with neurodevelopment or brain disease. For complete details on the use and execution of this protocol, please refer to Ho et al. (2017) and Wang et al. (2021).

Common Genetic Variation in Humans Impacts In Vitro Susceptibility to SARS-CoV-2 Infection.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant interindividual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly afflict healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants influence vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single-nucleotide polymorphism (rs4702), common in the population and located in the 3' UTR of the protease FURIN, influences alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can have an impact on viral infection and thus contribute to clinical heterogeneity in COVID-19. Ongoing genetic studies will help to identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs.

Analysis framework and experimental design for evaluating synergy-driving gene expression.

The mechanisms by which genetic risk variants interact with each other, as well as environmental factors, to contribute to complex genetic disorders remain unclear. We describe in detail our recently published approach to resolve distinct additive and synergistic transcriptomic effects after combinatorial manipulation of genetic variants and/or chemical perturbagens. Although first developed for CRISPR-based perturbation studies of isogenic human induced pluripotent stem cell-derived neurons, our methodology can be broadly applied to any RNA sequencing dataset, provided that raw read counts are available. Whereas other differential expression analyses reveal the effect of individual perturbations, here we specifically query interactions between two or more perturbagens, resolving the extent of non-additive (synergistic) interactions between perturbations. We discuss the careful experimental design required to resolve synergistic effects and considerations of statistical power and how to quantify observed synergy between experiments. Additionally, we speculate on potential future applications and explore the obvious limitations of this approach. Overall, by interrogating the effect of independent factors, alone and in combination, our analytic framework and experimental design facilitate the discovery of convergence and synergy downstream of gene and/or treatment perturbations hypothesized to contribute to complex diseases. We think that this protocol can be successfully applied by any scientist with bioinformatic skills and basic proficiency in the R programming language. Our computational pipeline ( https://github.com/nadschro/synergy-analysis ) is straightforward, does not require supercomputing support and can be conducted in a single day upon completion of RNA sequencing experiments.

Common genetic variation in humans impacts in vitro susceptibility to SARS-CoV-2 infection.

The host response to SARS-CoV-2, the etiologic agent of the COVID-19 pandemic, demonstrates significant inter-individual variability. In addition to showing more disease in males, the elderly, and individuals with underlying comorbidities, SARS-CoV-2 can seemingly render healthy individuals with profound clinical complications. We hypothesize that, in addition to viral load and host antibody repertoire, host genetic variants also impact vulnerability to infection. Here we apply human induced pluripotent stem cell (hiPSC)-based models and CRISPR-engineering to explore the host genetics of SARS-CoV-2. We demonstrate that a single nucleotide polymorphism (rs4702), common in the population at large, and located in the 3'UTR of the protease FURIN, impacts alveolar and neuron infection by SARS-CoV-2 in vitro. Thus, we provide a proof-of-principle finding that common genetic variation can impact viral infection, and thus contribute to clinical heterogeneity in SARS-CoV-2. Ongoing genetic studies will help to better identify high-risk individuals, predict clinical complications, and facilitate the discovery of drugs that might treat disease.

Neuronal impact of patient-specific aberrant NRXN1α splicing.

NRXN1 undergoes extensive alternative splicing, and non-recurrent heterozygous deletions in NRXN1 are strongly associated with neuropsychiatric disorders. We establish that human induced pluripotent stem cell (hiPSC)-derived neurons well represent the diversity of NRXN1α alternative splicing observed in the human brain, cataloguing 123 high-confidence in-frame human NRXN1α isoforms. Patient-derived NRXN1+/- hiPSC-neurons show a greater than twofold reduction in half of the wild-type NRXN1α isoforms and express dozens of novel isoforms from the mutant allele. Reduced neuronal activity in patient-derived NRXN1+/- hiPSC-neurons is ameliorated by overexpression of individual control isoforms in a genotype-dependent manner, whereas individual mutant isoforms decrease neuronal activity levels in control hiPSC-neurons. In a genotype-dependent manner, the phenotypic impact of patient-specific NRXN1+/- mutations can occur through a reduction in wild-type NRXN1α isoform levels as well as the presence of mutant NRXN1α isoforms.

Synergistic effects of common schizophrenia risk variants.

The mechanisms by which common risk variants of small effect interact to contribute to complex genetic disorders are unclear. Here, we apply a genetic approach, using isogenic human induced pluripotent stem cells, to evaluate the effects of schizophrenia (SZ)-associated common variants predicted to function as SZ expression quantitative trait loci (eQTLs). By integrating CRISPR-mediated gene editing, activation and repression technologies to study one putative SZ eQTL (FURIN rs4702) and four top-ranked SZ eQTL genes (FURIN, SNAP91, TSNARE1 and CLCN3), our platform resolves pre- and postsynaptic neuronal deficits, recapitulates genotype-dependent gene expression differences and identifies convergence downstream of SZ eQTL gene perturbations. Our observations highlight the cell-type-specific effects of common variants and demonstrate a synergistic effect between SZ eQTL genes that converges on synaptic function. We propose that the links between rare and common variants implicated in psychiatric disease risk constitute a potentially generalizable phenomenon occurring more widely in complex genetic disorders.

Psychosis Risk Candidate ZNF804A Localizes to Synapses and Regulates Neurite Formation and Dendritic Spine Structure.

BACKGROUND: Variation in the gene encoding zinc finger binding protein 804A (ZNF804A) is associated with schizophrenia and bipolar disorder. Evidence suggests that ZNF804A is a regulator of gene transcription and is present in nuclear and extranuclear compartments. However, a detailed examination of ZNF804A distribution and its neuronal functions has yet to be performed. METHODS: The localization of ZNF804A protein was examined in neurons derived from human neural progenitor cells, human induced pluripotent stem cells, or in primary rat cortical neurons. In addition, small interfering RNA-mediated knockdown of ZNF804A was conducted to determine its role in neurite formation, maintenance of dendritic spine morphology, and responses to activity-dependent stimulations. RESULTS: Endogenous ZNF804A protein localized to somatodendritic compartments and colocalized with the putative synaptic markers in young neurons derived from human neural progenitor cells and human induced pluripotent stem cells. In mature rat neurons, Zfp804A, the homolog of ZNF804A, was present in a subset of dendritic spines and colocalized with synaptic proteins in specific nanodomains, as determined by super-resolution microscopy. Interestingly, knockdown of ZNF804A attenuated neurite outgrowth in young neurons, an effect potentially mediated by reduced neuroligin-4 expression. Furthermore, knockdown of ZNF804A in mature neurons resulted in the loss of dendritic spine density and impaired responses to activity-dependent stimulation. CONCLUSIONS: These data reveal a novel subcellular distribution for ZNF804A within somatodendritic compartments and a nanoscopic organization at excitatory synapses. Moreover, our results suggest that ZNF804A plays an active role in neurite formation, maintenance of dendritic spines, and activity-dependent structural plasticity.

Characterisation of neurons derived from a cortical human neural stem cell line CTX0E16.

INTRODUCTION: Conditionally immortalised human neural progenitor cells (hNPCs) represent a robust source of native neural cells to investigate physiological mechanisms in both health and disease. However, in order to recognise the utility of such cells, it is critical to determine whether they retain characteristics of their tissue of origin and generate appropriate neural cell types upon differentiation. To this end, we have characterised the conditionally immortalised, cortically-derived, human NPC line, CTX0E16, investigating the molecular and cellular phenotype of differentiated neurons to determine whether they possess characteristics of cortical glutamatergic neurons. METHODS: Differentiated CTX0E16 cells were characterised by assessing expression of several neural fates markers, and examination of developing neuronal morphology. Expression of neurotransmitter receptors, signalling proteins and related proteins were assessed by q- and RT-PCR and complemented by Ca(2+) imaging, electrophysiology and assessment of ERK signalling in response to neurotransmitter ligand application. Finally, differentiated neurons were assessed for their ability to form putative synapses and to respond to activity-dependent stimulation. RESULTS: Differentiation of CTX0E16 hNPCs predominately resulted in the generation of neurons expressing markers of cortical and glutamatergic (excitatory) fate, and with a typical polarized neuronal morphology. Gene expression analysis confirmed an upregulation in the expression of cortical, glutamatergic and signalling proteins following differentiation. CTX0E16 neurons demonstrated Ca(2+) and ERK1/2 responses following exogenous neurotransmitter application, and after 6 weeks displayed spontaneous Ca(2+) transients and electrophysiological properties consistent with that of immature neurons. Differentiated CTX0E16 neurons also expressed a range of pre- and post-synaptic proteins that co-localized along distal dendrites, and moreover, displayed structural plasticity in response to modulation of neuronal activity. CONCLUSIONS: Taken together, these findings demonstrate that the CTX0E16 hNPC line is a robust source of cortical neurons, which display functional properties consistent with a glutamatergic phenotype. Thus CTX0E16 neurons can be used to study cortical cell function, and furthermore, as these neurons express a range of disease-associated genes, they represent an ideal platform with which to investigate neurodevelopmental mechanisms in native human cells in health and disease.