Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium

Integrative functional genomic analysis of human brain development and neuropsychiatric risks

INTRODUCTION The brain is responsible for cognition, behavior, and much of what makes us uniquely human. The development of the brain is a highly complex process, and this process is reliant on precise regulation of molecular and cellular events grounded in the spatiotemporal regulation of the transcriptome. Disruption of this regulation can lead to neuropsychiatric disorders. RATIONALE The regulatory, epigenomic, and transcriptomic features of the human brain have not been comprehensively compiled across time, regions, or cell types. Understanding the etiology of neuropsychiatric disorders requires knowledge not just of endpoint differences between healthy and diseased brains but also of the developmental and cellular contexts in which these differences arise. Moreover, an emerging body of research indicates that many aspects of the development and physiology of the human brain are not well recapitulated in model organisms, and therefore it is necessary that neuropsychiatric disorders be understood in the broader context of the developing and adult human brain. RESULTS Here we describe the generation and analysis of a variety of genomic data modalities at the tissue and single-cell levels, including transcriptome, DNA methylation, and histone modifications across multiple brain regions ranging in age from embryonic development through adulthood. We observed a widespread transcriptomic transition beginning during late fetal development and consisting of sharply decreased regional differences. This reduction coincided with increases in the transcriptional signatures of mature neurons and the expression of genes associated with dendrite development, synapse development, and neuronal activity, all of which were temporally synchronous across neocortical areas, as well as myelination and oligodendrocytes, which were asynchronous. Moreover, genes including MEF2C, SATB2, and TCF4, with genetic associations to multiple brain-related traits and disorders, converged in a small number of modules exhibiting spatial or spatiotemporal specificity. CONCLUSION We generated and applied our dataset to document transcriptomic and epigenetic changes across human development and then related those changes to major neuropsychiatric disorders. These data allowed us to identify genes, cell types, gene coexpression modules, and spatiotemporal loci where disease risk might converge, demonstrating the utility of the dataset and providing new insights into human development and disease

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

Impact of the Autism-Associated Long Noncoding RNA MSNP1AS on Neuronal Architecture and Gene Expression in Human Neural Progenitor Cells

We previously identified the long noncoding RNA (lncRNA) MSNP1AS (moesin pseudogene 1, antisense) as a functional element revealed by genome wide significant association with autism spectrum disorder (ASD). MSNP1AS expression was increased in the postmortem cerebral cortex of individuals with ASD and particularly in individuals with the ASD-associated genetic markers on chromosome 5p14.1. Here, we mimicked the overexpression of MSNP1AS observed in postmortem ASD cerebral cortex in human neural progenitor cell lines to determine the impact on neurite complexity and gene expression. ReNcell CX and SK-N-SH were transfected with an overexpression vector containing full-length MSNP1AS. Neuronal complexity was determined by the number and length of neuronal processes. Gene expression was determined by strand-specific RNA sequencing. MSNP1AS overexpression decreased neurite number and neurite length in both human neural progenitor cell lines. RNA sequencing revealed changes in gene expression in proteins involved in two biological processes: protein synthesis and chromatin remodeling. These data indicate that overexpression of the ASD-associated lncRNA MSNP1AS alters the number and length of neuronal processes. The mechanisms by which MSNP1AS overexpression impacts neuronal differentiation may involve protein synthesis and chromatin structure. These same biological processes are also implicated by rare mutations associated with ASD, suggesting convergent mechanisms

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

Abstract 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

A New Variant of Charcot-Marie-Tooth Disease Type 2 Is Probably the Result of a Mutation in the Neurofilament-Light Gene

Charcot-Marie-Tooth (CMT) disease is the most common inherited motor and sensory neuropathy. The axonal form of the disease is designated as “CMT type 2” (CMT2). Although four loci known to be implicated in autosomal dominant CMT2 have been mapped thus far (on 1p35-p36, 3q13.1, 3q13-q22, and 7p14), no one causative gene is yet known. A large Russian family with CMT2 was found in the Mordovian Republic (Russia). Affected members had the typical CMT2 phenotype. Additionally, several patients suffered from hyperkeratosis, although the association, if any, between the two disorders is not clear. Linkage with the CMT loci already known (CMT1A, CMT1B, CMT2A, CMT2B, CMT2D, and a number of other CMT-related loci) was excluded. Genomewide screening pinpointed the disease locus in this family to chromosome 8p21, within a 16-cM interval between markers D8S136 and D8S1769. A maximum two-point LOD score of 5.93 was yielded by a microsatellite from the 5′ region of the neurofilament-light gene (NF-L). Neurofilament proteins play an important role in axonal structure and are implicated in several neuronal disorders. Screening of affected family members for mutations in the NF-L gene and in the tightly linked neurofilament-medium gene (NF-M) revealed the only DNA alteration linked with the disease: a A998C transversion in the first exon of NF-L, which converts a conserved Gln333 amino acid to proline. This alteration was not found in 180 normal chromosomes. Twenty unrelated CMT2 patients, as well as 26 others with an undetermined form of CMT, also were screened for mutations in NF-L, but no additional mutations were found. It is suggested that Gln333Pro represents a rare disease-causing mutation, which results in the CMT2 phenotype

Analysis of Gene Expression Variance in Schizophrenia Using Structural Equation Modeling

Schizophrenia (SCZ) is a psychiatric disorder of unknown etiology. There is evidence suggesting that aberrations in neurodevelopment are a significant attribute of schizophrenia pathogenesis and progression. To identify biologically relevant molecular abnormalities affecting neurodevelopment in SCZ we used cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells). Here, we tested the hypothesis that variance in gene expression differs between individuals from SCZ and control groups. In CNON cells, variance in gene expression was significantly higher in SCZ samples in comparison with control samples. Variance in gene expression was enriched in five molecular pathways: serine biosynthesis, PI3K-Akt, MAPK, neurotrophin and focal adhesion. More than 14% of variance in disease status was explained within the logistic regression model (C-value = 0.70) by predictors accounting for gene expression in 69 genes from these five pathways. Structural equation modeling (SEM) was applied to explore how the structure of these five pathways was altered between SCZ patients and controls. Four out of five pathways showed differences in the estimated relationships among genes: between KRAS and NF1, and KRAS and SOS1 in the MAPK pathway; between PSPH and SHMT2 in serine biosynthesis; between AKT3 and TSC2 in the PI3K-Akt signaling pathway; and between CRK and RAPGEF1 in the focal adhesion pathway. Our analysis provides evidence that variance in gene expression is an important characteristic of SCZ, and SEM is a promising method for uncovering altered relationships between specific genes thus suggesting affected gene regulation associated with the disease. We identified altered gene-gene interactions in pathways enriched for genes with increased variance in expression in SCZ. These pathways and loci were previously implicated in SCZ, providing further support for the hypothesis that gene expression variance plays important role in the etiology of SCZ

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