Cultured Mesenchymal Cells from Nasal Turbinate as a Cellular Model of the Neurodevelopmental Component of Schizophrenia Etiology.

The study of neurodevelopmental molecular mechanisms in schizophrenia requires the development of adequate biological models such as patient-derived cells and their derivatives. We previously utilized cell lines with neural progenitor properties (CNON) derived from the superior or middle turbinates of patients with schizophrenia and control groups to study schizophrenia-specific gene expression. In this study, we analyzed single-cell RNA seq data from two CNON cell lines (one derived from an individual with schizophrenia (SCZ) and the other from a control group) and two biopsy samples from the middle turbinate (MT) (also from an individual with SCZ and a control). We compared our data with previously published data regarding the olfactory neuroepithelium and demonstrated that CNON originated from a single cell type present both in middle turbinate and the olfactory neuroepithelium and expressed in multiple markers of mesenchymal cells. To define the relatedness of CNON to the developing human brain, we also compared CNON datasets with scRNA-seq data derived from an embryonic brain and found that the expression profile of the CNON closely matched the expression profile one of the cell types in the embryonic brain. Finally, we evaluated the differences between SCZ and control samples to assess the utility and potential benefits of using CNON single-cell RNA seq to study the etiology of schizophrenia.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Reconstructing genetic history of Siberian and Northeastern European populations.

Siberia and Northwestern Russia are home to over 40 culturally and linguistically diverse indigenous ethnic groups, yet genetic variation and histories of peoples from this region are largely uncharacterized. We present deep whole-genome sequencing data (∼38×) from 28 individuals belonging to 14 distinct indigenous populations from that region. We combined these data sets with additional 32 modern-day and 46 ancient human genomes to reconstruct genetic histories of several indigenous Northern Eurasian populations. We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ust'-Ishim individual who was previously believed to have no modern-day descendants. Western Siberians trace 57% of their ancestry to ancient North Eurasians, represented by the 24,000-yr-old Siberian Mal'ta boy MA-1. Eastern Siberian populations formed a distinct sublineage that separated from other East Asian populations ∼10,000 yr ago. In addition, we uncovered admixtures between Siberians and Eastern European hunter-gatherers from Samara, Karelia, Hungary, and Sweden (from 8000-6600 yr ago); Yamnaya people (5300-4700 yr ago); and modern-day Northeastern Europeans. Our results provide new insights into genetic histories of Siberian and Northeastern European populations and evidence of ancient gene flow from Siberia into Europe.

Regulation of synaptic nlg-1/neuroligin abundance by the skn-1/Nrf stress response pathway protects against oxidative stress.

The Nrf family of transcription factors mediates adaptive responses to stress and longevity, but the identities of the crucial Nrf targets, and the tissues in which they function in multicellular organisms to promote survival, are not known. Here, we use whole transcriptome RNA sequencing to identify 810 genes whose expression is controlled by the SKN-1/Nrf2 negative regulator WDR-23 in the nervous system of Caenorhabditis elegans. Among the genes identified is the synaptic cell adhesion molecule nlg-1/neuroligin. We find that the synaptic abundance of NLG-1 protein increases following pharmacological treatments that generate oxidative stress or by the genetic activation of skn-1. Increasing nlg-1 dosage correlates with increased survival in response to oxidative stress, whereas genetic inactivation of nlg-1 reduces survival and impairs skn-1-mediated stress resistance. We identify a canonical SKN-1 binding site in the nlg-1 promoter that binds to SKN-1 in vitro and is necessary for SKN-1 and toxin-mediated increases in nlg-1 expression in vivo. Together, our results suggest that SKN-1 activation in the nervous system can confer protection to organisms in response to stress by directly regulating nlg-1/neuroligin expression.

Assessing characteristics of RNA amplification methods for single cell RNA sequencing.

BACKGROUND: Recently, measurement of RNA at single cell resolution has yielded surprising insights. Methods for single-cell RNA sequencing (scRNA-seq) have received considerable attention, but the broad reliability of single cell methods and the factors governing their performance are still poorly known. RESULTS: Here, we conducted a large-scale control experiment to assess the transfer function of three scRNA-seq methods and factors modulating the function. All three methods detected greater than 70% of the expected number of genes and had a 50% probability of detecting genes with abundance greater than 2 to 4 molecules. Despite the small number of molecules, sequencing depth significantly affected gene detection. While biases in detection and quantification were qualitatively similar across methods, the degree of bias differed, consistent with differences in molecular protocol. Measurement reliability increased with expression level for all methods and we conservatively estimate measurements to be quantitative at an expression level greater than ~5-10 molecules. CONCLUSIONS: Based on these extensive control studies, we propose that RNA-seq of single cells has come of age, yielding quantitative biological information.

Transcriptional Gene Silencing of the Autism-Associated Long Noncoding RNA MSNP1AS in Human Neural Progenitor Cells.

The long noncoding RNA MSNP1AS (moesin pseudogene 1, antisense) is a functional element that was previously associated with autism spectrum disorder (ASD) with genome-wide significance. Expression of MSNP1AS was increased 12-fold in the cerebral cortex of individuals with ASD and 22-fold in individuals with a genome-wide significantly associated ASD genetic marker on chromosome 5p14.1. Overexpression of MSNP1AS in human neuronal cells caused decreased expression of moesin protein, which is involved in neuronal process stability. In this study, we hypothesize that MSNP1AS knockdown impacts global transcriptome levels. We transfected the human neural progenitor cell line SK- N-SH with constructs that caused a 50% suppression of MSNP1AS expression. After 24 h, cells were harvested for total RNA isolation. Strand-specific RNA sequencing analysis indicated altered expression of 1,352 genes, including altered expression of 318 genes following correction for multiple comparisons. Expression of the OAS2 gene was increased >150-fold, a result that was validated by quantitative PCR. Gene ontology analysis of the 318 genes with altered expression following correction for multiple comparisons indicated that upregulated genes were significantly enriched for genes involved in immune response, and downregulated genes were significantly enriched for genes involved in chromatin remodeling. These data indicate multiple transcriptional and translational functions of MSNP1AS that impact ASD-relevant biological processes. Chromatin remodeling and immune response are biological processes implicated by genes with rare mutations associated with ASD. Our data suggest that the functional elements implicated by association of common genetic variants impact the same biological processes, suggesting a possible shared common molecular pathway of ASD.

The conserved SKN-1/Nrf2 stress response pathway regulates synaptic function in Caenorhabditis elegans.

The Nrf family of transcription factors plays a critical role in mediating adaptive responses to cellular stress and defends against neurodegeneration, aging, and cancer. Here, we report a novel role for the Caenorhabditis elegans Nrf homolog SKN-1 in regulating synaptic transmission at neuromuscular junctions (NMJs). Activation of SKN-1, either by acute pharmacological treatment with the mitochondrial toxin sodium arsenite or by mutations that cause constitutive SKN-1 activation, results in defects in neuromuscular function. Additionally, elimination of the conserved WD40 repeat protein WDR-23, a principal negative regulator of SKN-1, results in impaired locomotion and synaptic vesicle and neuropeptide release from cholinergic motor axons. Mutations that abolish skn-1 activity restore normal neuromuscular function to wdr-23 mutants and animals treated with toxin. We show that negative regulation of SKN-1 by WDR-23 in the intestine, but not at neuromuscular junctions, is necessary and sufficient for proper neuromuscular function. WDR-23 isoforms differentially localize to the outer membranes of mitochondria and to nuclei, and the effects of WDR-23 on neuromuscular function are dependent on its interaction with cullin E3 ubiquitin ligase. Finally, whole-transcriptome RNA sequencing of wdr-23 mutants reveals an increase in the expression of known SKN-1/Nrf2-regulated stress-response genes, as well as neurotransmission genes not previously implicated in SKN-1/Nrf2 responses. Together, our results indicate that SKN-1/Nrf2 activation may be a mechanism through which cellular stress, detected in one tissue, affects cellular function of a distal tissue through endocrine signaling. These results provide insight into how SKN-1/Nrf2 might protect the nervous system from damage in response to oxidative stress.

Spatially Resolved Transcriptomic Signatures of Hippocampal Subregions and Arc-Expressing Ensembles in Active Place Avoidance Memory.

The rodent hippocampus is a spatially organized neuronal network that supports the formation of spatial and episodic memories. We conducted bulk RNA sequencing and spatial transcriptomics experiments to measure gene expression changes in the dorsal hippocampus following the recall of active place avoidance (APA) memory. Through bulk RNA sequencing, we examined the gene expression changes following memory recall across the functionally distinct subregions of the dorsal hippocampus. We found that recall induced differentially expressed genes (DEGs) in the CA1 and CA3 hippocampal subregions were enriched with genes involved in synaptic transmission and synaptic plasticity, while DEGs in the dentate gyrus (DG) were enriched with genes involved in energy balance and ribosomal function. Through spatial transcriptomics, we examined gene expression changes following memory recall across an array of spots encompassing putative memory-associated neuronal ensembles marked by the expression of the IEGs Arc, Egr1, and c-Jun. Within samples from both trained and untrained mice, the subpopulations of spatial transcriptomic spots marked by these IEGs were transcriptomically and spatially distinct from one another. DEGs detected between Arc+ and Arc- spots exclusively in the trained mouse were enriched in several memory-related gene ontology terms, including "regulation of synaptic plasticity" and "memory." Our results suggest that APA memory recall is supported by regionalized transcriptomic profiles separating the CA1 and CA3 from the DG, transcriptionally and spatially distinct IEG expressing spatial transcriptomic spots, and biological processes related to synaptic plasticity as a defining the difference between Arc+ and Arc- spatial transcriptomic spots.

Mutations in LGI1 cause autosomal-dominant partial epilepsy with auditory features.

The epilepsies are a common, clinically heterogeneous group of disorders defined by recurrent unprovoked seizures. Here we describe identification of the causative gene in autosomal-dominant partial epilepsy with auditory features (ADPEAF, MIM 600512), a rare form of idiopathic lateral temporal lobe epilepsy characterized by partial seizures with auditory disturbances. We constructed a complete, 4.2-Mb physical map across the genetically implicated disease-gene region, identified 28 putative genes (Fig. 1) and resequenced all or part of 21 genes before identifying presumptive mutations in one copy of the leucine-rich, glioma-inactivated 1 gene (LGI1) in each of five families with ADPEAF. Previous studies have indicated that loss of both copies of LGI1 promotes glial tumor progression. We show that the expression pattern of mouse Lgi1 is predominantly neuronal and is consistent with the anatomic regions involved in temporal lobe epilepsy. Discovery of LGI1 as a cause of ADPEAF suggests new avenues for research on pathogenic mechanisms of idiopathic epilepsies.

Mutant small heat-shock protein 27 causes axonal Charcot-Marie-Tooth disease and distal hereditary motor neuropathy.

Charcot-Marie-Tooth disease (CMT) is the most common inherited neuromuscular disease and is characterized by considerable clinical and genetic heterogeneity. We previously reported a Russian family with autosomal dominant axonal CMT and assigned the locus underlying the disease (CMT2F; OMIM 606595) to chromosome 7q11-q21 (ref. 2). Here we report a missense mutation in the gene encoding 27-kDa small heat-shock protein B1 (HSPB1, also called HSP27) that segregates in the family with CMT2F. Screening for mutations in HSPB1 in 301 individuals with CMT and 115 individuals with distal hereditary motor neuropathies (distal HMNs) confirmed the previously observed mutation and identified four additional missense mutations. We observed the additional HSPB1 mutations in four families with distal HMN and in one individual with CMT neuropathy. Four mutations are located in the Hsp20-alpha-crystallin domain, and one mutation is in the C-terminal part of the HSP27 protein. Neuronal cells transfected with mutated HSPB1 were less viable than cells expressing the wild-type protein. Cotransfection of neurofilament light chain (NEFL) and mutant HSPB1 resulted in altered neurofilament assembly in cells devoid of cytoplasmic intermediate filaments.

Cell Type Catalog of Middle Turbinate Epithelium (Cell Catalog of Turbinate Epithelium).

Objective: Single-cell RNA-sequencing of middle turbinate mucosa was performed to create the first single-cell transcriptome catalog of this part of the human body. Study Design: Basic science research. Setting: Single center, tertiary care center. Methods: Samples were obtained from the head of the middle turbinate from a healthy volunteer. After the specimen was prepared per lab protocol, cells were dissociated, resuspended, and counted. Single-cell libraries were then prepared according to the 10x Genomics protocol and sequenced using NovaSeq 6000 (Illumina). Sequencing data were processed using Cell Ranger, and clustering and gene expression analysis was performed using Seurat. Cell types were annotated through expression profiling of single cells using known markers and data from other single-cell studies. Results: Fourteen unique cell types were identified, including serous, goblet, club, basal, ciliated, endothelial, and mesenchymal cells, as well as multiple types of blood cells. Conclusion: This catalog provides a comprehensive depiction of the cellular composition of middle turbinate mucosa. By uncovering the cellular stratification of gene expression profiles in the healthy middle turbinate epithelium, the groundwork has been laid for further investigation into the molecular pathogenesis and targeted therapy of sinonasal disease.

Inhibiting Phosphatidylcholine Remodeling in Adipose Tissue Increases Insulin Sensitivity.

Cell membrane phosphatidylcholine (PC) composition is regulated by lysophosphatidylcholine acyltransferase (LPCAT); changes in membrane PC saturation are implicated in metabolic disorders. Here, we identified LPCAT3 as the major isoform of LPCAT in adipose tissue and created adipocyte-specific Lpcat3-knockout mice to study adipose tissue lipid metabolism. Transcriptome sequencing and plasma adipokine profiling were used to investigate how LPCAT3 regulates adipose tissue insulin signaling. LPCAT3 deficiency reduced polyunsaturated PCs in adipocyte plasma membranes, increasing insulin sensitivity. LPCAT3 deficiency influenced membrane lipid rafts, which activated insulin receptors and AKT in adipose tissue, and attenuated diet-induced insulin resistance. Conversely, higher LPCAT3 activity in adipose tissue from ob/ob, db/db, and high-fat diet-fed mice reduced insulin signaling. Adding polyunsaturated PCs to mature human or mouse adipocytes in vitro worsened insulin signaling. We suggest that targeting LPCAT3 in adipose tissue to manipulate membrane phospholipid saturation is a new strategy to treat insulin resistance.

Cultured Mesenchymal Cells from Nasal Turbinate as a Cellular Model of the Neurodevelopmental Component of Schizophrenia Etiology.

Study of the neurodevelopmental molecular mechanisms of schizophrenia requires the development of adequate biological models such as patient-derived cells and their derivatives. We previously used cell lines with neural progenitor properties (CNON) derived from superior or middle turbinates of patients with schizophrenia and control groups to study gene expression specific to schizophrenia. In this study, we compared single cell-RNA seq data from two CNON cell lines, one derived from an individual with schizophrenia (SCZ) and the other from a control group, with two biopsy samples from the middle turbinate (MT), also from an individual with SCZ and a control. In addition, we compared our data with previously published data from olfactory neuroepithelium (1). Our data demonstrated that CNON originated from a single cell type which is present both in middle turbinate and olfactory neuroepithelium. CNON express multiple markers of mesenchymal cells. In order to define relatedness of CNON to the developing human brain, we also compared CNON datasets with scRNA-seq data of embryonic brain (2) and found that the expression profile of CNON very closely matched one of the cell types in the embryonic brain. Finally, we evaluated differences between SCZ and control samples to assess usability and potential benefits of using single cell RNA-seq of CNON to study etiology of schizophrenia.

RseqFlow: workflows for RNA-Seq data analysis.

SUMMARY: We have developed an RNA-Seq analysis workflow for single-ended Illumina reads, termed RseqFlow. This workflow includes a set of analytic functions, such as quality control for sequencing data, signal tracks of mapped reads, calculation of expression levels, identification of differentially expressed genes and coding SNPs calling. This workflow is formalized and managed by the Pegasus Workflow Management System, which maps the analysis modules onto available computational resources, automatically executes the steps in the appropriate order and supervises the whole running process. RseqFlow is available as a Virtual Machine with all the necessary software, which eliminates any complex configuration and installation steps. AVAILABILITY AND IMPLEMENTATION: http://genomics.isi.edu/rnaseq CONTACT: wangying@xmu.edu.cn; knowles@med.usc.edu; deelman@isi.edu; tingchen@usc.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

YPK9 and WHI2 Negatively Interact during Oxidative Stress.

Yeast PARK9 (YPK9) shares homology with human ATP13A2, which encodes a polyamine transporter implicated in juvenile forms of Parkinson's disease. We used YPK9 to gain insight into how ATP13A2 affects cell growth and sensitivity to oxidative stress. Surprisingly, the YPK9 deletion strain from the Saccharomyces cerevisiae deletion collection (YKO) in wildtype BY4741 (mating type a) grew faster and was more resistant to hydrogen peroxide than a commercial, putative parental BY4741 wildtype strain (BY4741COM). In contrast, deleting YPK9 from BY4741COM rendered it very sensitive to hydrogen peroxide, suggesting its background is different from that of the deletion collection. Whole-genome sequencing revealed that BY4741COM and BY4741COMypk9∆ contain a novel premature stop codon near the 3' end of WHI2 (WHI2G1324T), whereas the collection's YPK9 deletion strain contains WHI2, which encodes a 486 amino acid protein, Whi2p. Replacing full-length WHI2 with the sequence coding for the predicted truncation (Whi2pE442*) rendered strains more sensitive to hydrogen peroxide, whereas the converse replacement rendered them more resistant. The sequences of WHI2 in 20 randomly chosen strains from the collection encode the full-length protein, indicating that the putative parental BY4741 WHI2G1324T strain's genetic background differs from that of the deletion collection. Examination of WHI2 sequences in several commonly used wildtype S. cerevisiae strains and isolates revealed other Whi2p truncations that might yield altered phenotypes. Together, these results demonstrate a novel premature stop codon in WHI2 that renders yeast sensitive to hydrogen peroxide; they also reveal a negative genetic interaction between WHI2 and YPK9 in the presence of hydrogen peroxide in the BY4741 background.

Transcriptome data of temporal and cingulate cortex in the Rett syndrome brain.

Rett syndrome is an X-linked neurodevelopmental disorder caused by mutation in the methyl-CpG-binding protein 2 gene (MECP2) in the majority of cases. We describe an RNA sequencing dataset of postmortem brain tissue samples from four females clinically diagnosed with Rett syndrome and four age-matched female donors. The dataset contains 16 transcriptomes, including two brain regions, temporal and cingulate cortex, for each individual. We compared our dataset with published transcriptomic analyses of postmortem brain tissue from Rett syndrome and found consistent gene expression alterations among regions of the cerebral cortex. Our data provide a valuable resource to explore the biology of the human brain in Rett syndrome.

Robust RNA-Seq of aRNA-amplified single cell material collected by patch clamp.

Most single cell RNA sequencing protocols start with single cells dispersed from intact tissue. High-throughput processing of the separated cells is enabled using microfluidics platforms. However, dissociation of tissue results in loss of information about cell location and morphology and potentially alters the transcriptome. An alternative approach for collecting RNA from single cells is to re-purpose the electrophysiological technique of patch clamp recording. A hollow patch pipette is attached to individual cells, enabling the recording of electrical activity, after which the cytoplasm may be extracted for single cell RNA-Seq ("Patch-Seq"). Since the tissue is not disaggregated, the location of cells is readily determined, and the morphology of the cells is maintained, making possible the correlation of single cell transcriptomes with cell location, morphology and electrophysiology. Recent Patch-Seq studies utilizes PCR amplification to increase amount of nucleic acid material to the level required for current sequencing technologies. PCR is prone to create biased libraries - especially with the extremely high degrees of exponential amplification required for single cell amounts of RNA. We compared a PCR-based approach with linear amplifications and demonstrate that aRNA amplification (in vitro transcription, IVT) is more sensitive and robust for single cell RNA collected by a patch clamp pipette.

Gene Expression in Patient-Derived Neural Progenitors Implicates WNT5A Signaling in the Etiology of Schizophrenia.

BACKGROUND: Genome-wide association studies of schizophrenia have demonstrated that variations in noncoding regions are responsible for most of the common variation heritability of the disease. It is hypothesized that these risk variants alter gene expression. Therefore, studying alterations in gene expression in schizophrenia may provide a direct approach to understanding the etiology of the disease. In this study we use cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells) as a genetically unaltered cellular model to elucidate the neurodevelopmental aspects of schizophrenia. METHODS: We performed a gene expression study using RNA sequencing of CNON cells from 111 control subjects and 144 individuals with schizophrenia. Differentially expressed genes were identified with DESeq2 software, using covariates to correct for sex, age, library batches, and 1 surrogate variable component. RESULTS: A total of 80 genes were differentially expressed (false discovery rate < 10%), showing enrichment in cell migration, cell adhesion, developmental process, synapse assembly, cell proliferation, and related Gene Ontology categories. Cadherin and Wnt signaling pathways were positive in overrepresentation test, and, in addition, many genes were specifically involved in WNT5A signaling. The differentially expressed genes were modestly, but significantly, enriched in the genes overlapping single nucleotide polymorphisms with genome-wide significant association from the Psychiatric Genomics Consortium genome-wide association study of schizophrenia. We also found substantial overlap with genes associated with other psychiatric disorders or brain development, enrichment in the same Gene Ontology categories as genes with mutations de novo in schizophrenia, and studies of induced pluripotent stem cell-derived neural progenitor cells. CONCLUSIONS: CNON cells are a good model of the neurodevelopmental aspects of schizophrenia and can be used to elucidate the etiology of the disorder.

Association and linkage analysis of candidate genes GRP, GRPR, CRHR1, and TACR1 in panic disorder.

Panic disorder (PD) is a debilitating anxiety disorder, characterized by recurrent episodes of intense fear that are accompanied by autonomic and psychological symptoms leading to behavioral impairment. Basic research implicates neuropeptide-signaling genes in the modulation of anxiety and stress. The genes encoding corticotropin releasing hormone receptor 1 (CRHR1), tachykinin receptor 1 (TACR1), gastrin releasing peptide (GRP), and gastrin releasing peptide receptor (GRPR) were selected as candidates for PD based on their biology. Linkage and association analysis in 120 multiplex U.S. PD pedigrees was performed using 21 single nucleotide polymorphisms (SNPs). Parametric and non-parametric linkage tests in pedigrees, for single point and multipoint analysis, revealed limited support for genetic linkage to TACR1 (parametric and non-parametric lod scores approximately 1). The family-based association test (FBAT) generated nominal support for allelic association in TACR1 (P = 0.02), and GRP (P = 0.02), findings which must be considered in the light of multiple comparisons. Further exploration of the GRP and TACR1 findings in large case-control PD samples may provide more definitive evidence implicating these loci in the genetic etiology of PD.

Deconvolution of transcriptional networks identifies TCF4 as a master regulator in schizophrenia.

Applying tissue-specific deconvolution of transcriptional networks to identify their master regulators (MRs) in neuropsychiatric disorders has been largely unexplored. Here, using two schizophrenia (SCZ) case-control RNA-seq datasets, one on postmortem dorsolateral prefrontal cortex (DLPFC) and another on cultured olfactory neuroepithelium, we deconvolved the transcriptional networks and identified TCF4 as a top candidate MR that may be dysregulated in SCZ. We validated TCF4 as a MR through enrichment analysis of TCF4-binding sites in induced pluripotent stem cell (hiPSC)-derived neurons and in neuroblastoma cells. We further validated the predicted TCF4 targets by knocking down TCF4 in hiPSC-derived neural progenitor cells (NPCs) and glutamatergic neurons (Glut_Ns). The perturbed TCF4 gene network in NPCs was more enriched for pathways involved in neuronal activity and SCZ-associated risk genes, compared to Glut_Ns. Our results suggest that TCF4 may serve as a MR of a gene network dysregulated in SCZ at early stages of neurodevelopment.