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

Inflammation-related pathology in the olfactory epithelium: its impact on the olfactory system in psychotic disorders

Smell deficits and neurobiological changes in the olfactory bulb (OB) and olfactory epithelium (OE) have been observed in schizophrenia and related disorders. The OE is the most peripheral olfactory system located outside the cranium, and is connected with the brain via direct neuronal projections to the OB. Nevertheless, it is unknown whether and how a disturbance of the OE affects the OB in schizophrenia and related disorders. Addressing this gap would be the first step in studying the impact of OE pathology in the disease pathophysiology in the brain. In this cross-species study, we observed that chronic, local OE inflammation with a set of upregulated genes in an inducible olfactory inflammation (IOI) mouse model led to a volume reduction, layer structure changes, and alterations of neuron functionality in the OB. Furthermore, IOI model also displayed behavioral deficits relevant to negative symptoms (avolition) in parallel to smell deficits. In first episode psychosis (FEP) patients, we observed a significant alteration in immune/inflammation-related molecular signatures in olfactory neuronal cells (ONCs) enriched from biopsied OE and a significant reduction in the OB volume, compared with those of healthy controls (HC). The increased expression of immune/inflammation-related molecules in ONCs was significantly correlated to the OB volume reduction in FEP patients, but no correlation was found in HCs. Moreover, the increased expression of human orthologues of the IOI genes in ONCs was significantly correlated with the OB volume reduction in FEP, but not in HCs. Together, our study implies a potential mechanism of the OE-OB pathology in patients with psychotic disorders (schizophrenia and related disorders). We hope that this mechanism may have a cross-disease implication, including COVID-19-elicited mental conditions that include smell deficits.

3D Genome Reconstruction from Partially Phased Hi-C Data

The 3-dimensional (3D) structure of the genome is of significant importance for many cellular processes. In this paper, we study the problem of reconstructing the 3D structure of chromosomes from Hi-C data of diploid organisms, which poses additional challenges compared to the better-studied haploid setting. With the help of techniques from algebraic geometry, we prove that a small amount of phased data is sufficient to ensure finite identifiability, both for noiseless and noisy data. In the light of these results, we propose a new 3D reconstruction method based on semidefinite programming, paired with numerical algebraic geometry and local optimization. The performance of this method is tested on several simulated datasets under different noise levels and with different amounts of phased data. We also apply it to a real dataset from mouse X chromosomes, and we are then able to recover previously known structural features.

Inflammation-related pathology in the olfactory epithelium: its impact on the olfactory system in psychotic disorders

Smell deficits and neurobiological changes in the olfactory bulb (OB) and olfactory epithelium (OE) have been observed in schizophrenia and related disorders. The OE is the most peripheral olfactory system located outside the cranium, and is connected with the brain via direct neuronal projections to the OB. Nevertheless, it is unknown whether and how a disturbance of the OE affects the OB in schizophrenia and related disorders. Addressing this gap would be the first step in studying the impact of OE pathology in the disease pathophysiology in the brain. In this cross-species study, we observed that chronic, local OE inflammation with a set of upregulated genes in an inducible olfactory inflammation (IOI) mouse model led to a volume reduction, layer structure changes, and alterations of neuron functionality in the OB. Furthermore, IOI model also displayed behavioral deficits relevant to negative symptoms (avolition) in parallel to smell deficits. In first episode psychosis (FEP) patients, we observed a significant alteration in immune/inflammation-related molecular signatures in olfactory neuronal cells (ONCs) enriched from biopsied OE and a significant reduction in the OB volume, compared with those of healthy controls (HC). The increased expression of immune/inflammation-related molecules in ONCs was significantly correlated to the OB volume reduction in FEP patients, but no correlation was found in HCs. Moreover, the increased expression of human orthologues of the IOI genes in ONCs was significantly correlated with the OB volume reduction in FEP, but not in HCs. Together, our study implies a potential mechanism of the OE-OB pathology in patients with psychotic disorders (schizophrenia and related disorders). We hope that this mechanism may have a cross-disease implication, including COVID-19-elicited mental conditions that include smell deficits.

Unraveling the Link between Olfactory Deficits and Neuropsychiatric Disorders

Smell loss has caught public attention during the recent COVID-19 pandemic. Research on olfactory function in health and disease gains new momentum. Smell deficits have long been recognized as an early clinical sign associated with neuropsychiatric disorders. Here we review research on the associations between olfactory deficits and neuropathological conditions, focusing on recent progress in four areas: (1) human clinical studies of the correlations between smell deficits and neuropsychiatric disorders; (2) development of olfactory mucosa-derived tissue and cell models for studying the molecular pathologic mechanisms; (3) recent findings in brain imaging studies of structural and functional connectivity changes in olfactory pathways in neuropsychiatric disorders; and (4) application of preclinical animal models to validate and extend the findings from human subjects. Together, these studies have provided strong evidence of the link between the olfactory system and neuropsychiatric disorders, highlighting the relevance of deepening our understanding of the role of the olfactory system in pathophysiological processes. Following the lead of studies reviewed here, future research in this field may open the door to the early detection of neuropsychiatric disorders, personalized treatment approaches, and potential therapeutic interventions through nasal administration techniques, such as nasal brush or nasal spray.

Olfactory neuronal cells as a promising tool to realize the "druggable genome" approach for drug discovery in neuropsychiatric disorders

"Druggable genome" is a novel concept that emphasizes the importance of using the information of genome-wide genetic studies for drug discovery and development. Successful precedents of "druggable genome" have recently emerged for some disorders by combining genomic and gene expression profiles with medical and pharmacological knowledge. One of the key premises for the success is the good access to disease-relevant tissues from "living" patients in which we may observe molecular expression changes in association with symptomatic alteration. Thus, given brain biopsies are ethically and practically difficult, the application of the "druggable genome" approach is challenging for neuropsychiatric disorders. Here, to fill this gap, we propose the use of olfactory neuronal cells (ONCs) biopsied and established via nasal biopsy from living subjects. By using candidate genes that were proposed in a study in which genetic information, postmortem brain expression profiles, and pharmacological knowledge were considered for cognition in the general population, we addressed the utility of ONCs in the "druggable genome" approach by using the clinical and cell resources of an established psychosis cohort in our group. Through this pilot effort, we underscored the chloride voltage-gated channel 2 (CLCN2) gene as a possible druggable candidate for early-stage psychosis. The CLCN2 gene expression was associated with verbal memory, but not with other dimensions in cognition, nor psychiatric manifestations (positive and negative symptoms). The association between this candidate molecule and verbal memory was also confirmed at the protein level. By using ONCs from living subjects, we now provide more specific information regarding molecular expression and clinical phenotypes. The use of ONCs also provides the opportunity of validating the relationship not only at the RNA level but also protein level, leading to the potential of functional assays in the future. Taken together, we now provide evidence that supports the utility of ONCs as a tool for the "druggable genome" approach in translational psychiatry.

The miR-124-AMPAR pathway connects polygenic risks with behavioral changes shared between schizophrenia and bipolar disorder

Schizophrenia (SZ) and bipolar disorder (BP) are highly heritable major psychiatric disorders that share a substantial portion of genetic risk as well as their clinical manifestations. This raises a fundamental question of whether, and how, common neurobiological pathways translate their shared polygenic risks into shared clinical manifestations. This study shows the miR-124-3p-AMPAR pathway as a key common neurobiological mediator that connects polygenic risks with behavioral changes shared between these two psychotic disorders. We discovered the upregulation of miR-124-3p in neuronal cells and the postmortem prefrontal cortex from both SZ and BP patients. Intriguingly, the upregulation is associated with the polygenic risks shared between these two disorders. Seeking mechanistic dissection, we generated a mouse model that upregulates miR-124-3p in the medial prefrontal cortex. We demonstrated that the upregulation of miR-124-3p increases GRIA2-lacking calcium-permeable AMPARs and perturbs AMPAR-mediated excitatory synaptic transmission, leading to deficits in the behavioral dimensions shared between SZ and BP.

A Comparison of Topologically Associating Domain Callers Based on Hi-C Data

Topologically associating domains (TADs) are local chromatin interaction domains, which have been shown to play an important role in gene expression regulation. TADs were originally discovered in the investigation of 3D genome organization based on High-throughput Chromosome Conformation Capture (Hi-C) data. Continuous considerable efforts have been dedicated to developing methods for detecting TADs from Hi-C data. Different computational methods for TADs identification vary in their assumptions and criteria in calling TADs. As a consequence, the TADs called by these methods differ in their similarities and biological features they are enriched in. In this work, we performed a systematic comparison of twenty-six TAD callers. We first compared the TADs and gaps between adjacent TADs across different methods, resolutions, and sequencing depths. We then assessed the quality of TADs and TAD boundaries according to three criteria: the decay of contact frequencies over the genomic distance, enrichment and depletion of regulatory elements around TAD boundaries, and reproducibility of TADs and TAD boundaries in replicate samples. Last, due to the lack of a gold standard of TADs, we also evaluated the performance of the methods on synthetic datasets. We discussed the key principles of TAD callers, and pinpointed current situation in the detection of TADs. We provide a concise, comprehensive, and systematic framework for evaluating the performance of TAD callers, and expect our work will provide useful guidance in choosing suitable approaches for the detection and evaluation of TADs.

Characterizing chromatin interactions of regulatory elements and nucleosome positions, using Hi-C, Micro-C, and promoter capture Micro-C

BACKGROUND

Regulatory elements such as promoters, enhancers, and insulators interact each other to mediate molecular processes. To capture chromatin interactions of regulatory elements, 3C-derived methods such as Hi-C and Micro-C are developed. Here, we generated and analyzed Hi-C, Micro-C, and promoter capture Micro-C datasets with different sequencing depths to study chromatin interactions of regulatory elements and nucleosome positions in human prostate cancer cells.

RESULTS

Compared to Hi-C, Micro-C identifies more high-resolution loops, including ones around structural variants. By evaluating the effect of sequencing depth, we revealed that more than 2 billion reads of Micro-C are needed to detect chromatin interactions at 1 kb resolution. Moreover, we found that deep-sequencing identifies additional long-range loops that are longer than 1 Mb in distance. Furthermore, we found that more than 50% of the loops are involved in insulators while less than 10% of the loops are promoter-enhancer loops. To comprehensively capture chromatin interactions that promoters are involved in, we performed promoter capture Micro-C. Promoter capture Micro-C identifies loops near promoters with a lower amount of sequencing reads. Sequencing of 160 million reads of promoter capture Micro-C resulted in reaching a plateau of identifying loops. However, there was still a subset of promoters that are not involved in loops even after deep-sequencing. By integrating Micro-C with NOMe-seq and ChIP-seq, we found that active promoters involved in loops have a more accessible region with lower levels of DNA methylation and more highly phased nucleosomes, compared to active promoters that are not involved in loops.

CONCLUSION

We determined the required sequencing depth for Micro-C and promoter capture Micro-C to generate high-resolution chromatin interaction maps and loops. We also investigated the effect of sequencing coverage of Hi-C, Micro-C, and promoter capture Micro-C on detecting chromatin loops. Our analyses suggest the presence of distinct regulatory element groups, which are differently involved in nucleosome positions and chromatin interactions. This study does not only provide valuable insights on understanding chromatin interactions of regulatory elements, but also present guidelines for designing research projects on chromatin interactions among regulatory elements.

Characterizing and Targeting Genes Regulated by Transcription Factor MYBL2 in Lung Adenocarcinoma Cells

Overexpression of MYBL2 is associated with poor survival of lung adenocarcinoma patients, but the molecular mechanism by which it regulates transcription and carcinogenesis has not yet been elucidated. In this study, we performed ChIP-seq using an MYBL2-targeted antibody and discovered that MYBL2 primarily binds to the promoters of highly expressed genes in lung adenocarcinoma cells. Using a knockdown experiment of MYBL2 and global transcriptome profiling, we identified that over a thousand genes are dysregulated by MYBL2, and MYBL2 acts as a transcriptional activator in lung adenocarcinoma cells. Moreover, we revealed that the binding sites of FOXM1 are largely shared with MYBL2 binding sites, and genes involved in cell cycle phase transitions are regulated by these transcription factors. We furthermore investigated the effect of a previously reported FOXM1 inhibitor, FDI-6, in lung adenocarcinoma cells. We demonstrated that FDI-6 decreases the proliferation of lung adenocarcinoma cells and inhibits the activities of FOXM1 as well as MYBL2. Moreover, we found that genes involved in cell death and cell cycle are inhibited by FDI-6. Overall, our findings suggest that MYBL2 and FOXM1 activate cell cycle genes together, acting as oncogenic transcription factors in lung adenocarcinoma cells, and they are potential treatment targets for the disease.

Functional variant rs2270363 on 16p13.3 confers schizophrenia risk by regulating NMRAL1

Recent genome-wide association studies have reported multiple schizophrenia risk loci, yet the functional variants and their roles in schizophrenia remain to be characterized. Here we identify a functional single nucleotide polymorphism (rs2270363: G>A) at the schizophrenia risk locus 16p13.3. rs2270363 lies in the E-box element of the promoter of NMRAL1 and disrupts binding of the basic helix-loop-helix leucine zipper family proteins, including USF1, MAX and MXI1. We validated the regulatory effects of rs2270363 using reporter gene assays and electrophoretic mobility shift assay. Besides, expression quantitative trait loci analysis showed that the risk allele (A) of rs2270363 was significantly associated with elevated NMRAL1 expression in the human brain. Transcription factors knockdown and CRISPR-Cas9-mediated editing further confirmed the regulatory effects of the genomic region containing rs2270363 on NMRAL1. Intriguingly, NMRAL1 was significantly downregulated in the brain of schizophrenia patients compared with healthy subjects, and knockdown of Nmral1 expression affected proliferation and differentiation of mouse neural stem cells, as well as genes and pathways associated with brain development and synaptic transmission. Of note, Nmral1 knockdown resulted in significant decrease of dendritic spine density, revealing the potential pathophysiological mechanisms of NMRAL1 in schizophrenia. Finally, we independently confirmed the association between rs2270363 and schizophrenia in the Chinese population and found that the risk allele of rs2270363 was the same in European and Chinese populations. These lines of evidence suggest that rs2270363 may confer schizophrenia risk by regulating NMRAL1, a gene whose expression dysregulation might be involved in the pathogenesis of schizophrenia by affecting neurodevelopment and synaptic plasticity.

Modeling common and rare genetic risk factors of neuropsychiatric disorders in human induced pluripotent stem cells

Recent genome-wide association studies (GWAS) and whole-exome sequencing of neuropsychiatric disorders, especially schizophrenia, have identified a plethora of common and rare disease risk variants/genes. Translating the mounting human genetic discoveries into novel disease biology and more tailored clinical treatments is tied to our ability to causally connect genetic risk variants to molecular and cellular phenotypes. When combined with the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated (Cas) nuclease-mediated genome editing system, human induced pluripotent stem cell (hiPSC)-derived neural cultures (both 2D and 3D organoids) provide a promising tractable cellular model for bridging the gap between genetic findings and disease biology. In this review, we first conceptualize the advances in understanding the disease polygenicity and convergence from the past decade of iPSC modeling of different types of genetic risk factors of neuropsychiatric disorders. We then discuss the major cell types and cellular phenotypes that are most relevant to neuropsychiatric disorders in iPSC modeling. Finally, we critically review the limitations of iPSC modeling of neuropsychiatric disorders and outline the need for implementing and developing novel methods to scale up the number of iPSC lines and disease risk variants in a systematic manner. Sufficiently scaled-up iPSC modeling and a better functional interpretation of genetic risk variants, in combination with cutting-edge CRISPR/Cas9 gene editing and single-cell multi-omics methods, will enable the field to identify the specific and convergent molecular and cellular phenotypes in precision for neuropsychiatric disorders.

A comparison of topologically associating domain callers over mammals at high resolution

BACKGROUND

Topologically associating domains (TADs) are locally highly-interacting genome regions, which also play a critical role in regulating gene expression in the cell. TADs have been first identified while investigating the 3D genome structure over High-throughput Chromosome Conformation Capture (Hi-C) interaction dataset. Substantial degree of efforts have been devoted to develop techniques for inferring TADs from Hi-C interaction dataset. Many TAD-calling methods have been developed which differ in their criteria and assumptions in TAD inference. Correspondingly, TADs inferred via these callers vary in terms of both similarities and biological features they are enriched in.

RESULT

We have carried out a systematic comparison of 27 TAD-calling methods over mammals. We use Micro-C, a recent high-resolution variant of Hi-C, to compare TADs at a very high resolution, and classify the methods into 3 categories: feature-based methods, Clustering methods, Graph-partitioning methods. We have evaluated TAD boundaries, gaps between adjacent TADs, and quality of TADs across various criteria. We also found particularly CTCF and Cohesin proteins to be effective in formation of TADs with corner dots. We have also assessed the callers performance on simulated datasets since a gold standard for TADs is missing. TAD sizes and numbers change remarkably between TAD callers and dataset resolutions, indicating that TADs are hierarchically-organized domains, instead of disjoint regions. A core subset of feature-based TAD callers regularly perform the best while inferring reproducible domains, which are also enriched for TAD related biological properties.

CONCLUSION

We have analyzed the fundamental principles of TAD-calling methods, and identified the existing situation in TAD inference across high resolution Micro-C interaction datasets over mammals. We come up with a systematic, comprehensive, and concise framework to evaluate the TAD-calling methods performance across Micro-C datasets. Our research will be useful in selecting appropriate methods for TAD inference and evaluation based on available data, experimental design, and biological question of interest. We also introduce our analysis as a benchmarking tool with publicly available source code.

FOXC1 Binds Enhancers and Promotes Cisplatin Resistance in Bladder Cancer

Simple Summary

In bladder cancer, cisplatin remains the front-line therapy, but drug resistance is common. Previously, we showed that cancer cells can spontaneously convert to an aggressive drug-resistant phenotype without mutational events. In the current work, we explored the epigenetic mechanism behind the conversion to the drug-resistant phenotype. We discovered that drug-resistant cells have differentially accessible enhancers, which are bound by FOXC1, a transcription factor that is overexpressed in these cells. Accordingly, FOXC1 knockout significantly attenuates the emergence of the drug-resistant phenotype and reduces cell survival upon cisplatin treatment. These findings suggest that FOXC1 binding at accessible enhancers promotes cisplatin drug resistance in bladder cancer cells. Therefore, FOXC1 targeting may be a new therapeutic avenue to mitigate cisplatin resistance and improve treatment efficacy in bladder cancer.

Abstract

Chemotherapy resistance is traditionally attributed to DNA mutations that confer a survival advantage under drug selection pressure. However, in bladder cancer and other malignancies, we and others have previously reported that cancer cells can convert spontaneously to an aggressive drug-resistant phenotype without prior drug selection or mutational events. In the current work, we explored possible epigenetic mechanisms behind this phenotypic plasticity. Using Hoechst dye exclusion and flow cytometry, we isolated the aggressive drug-resistant cells and analyzed their chromatin accessibility at regulatory elements. Compared to the rest of the cancer cell population, the aggressive drug-resistant cells exhibited enhancer accessibility changes. In particular, we found that differentially accessible enhancers were enriched for the FOXC1 transcription factor motif, and that FOXC1 was the most significantly overexpressed gene in aggressive drug-resistant cells. ChIP-seq analysis revealed that differentially accessible enhancers in aggressive drug-resistant cells had a higher FOXC1 binding, which regulated the expression of adjacent cancer-relevant genes like ABCB1 and ID3. Accordingly, cisplatin treatment of bladder cancer cells led to an increased FOXC1 expression, which mediated cell survival and conversion to a drug-resistant phenotype. Collectively, these findings suggest that FOXC1 contributes to phenotypic plasticity by binding enhancers and promoting a mutation-independent shift towards cisplatin resistance in bladder cancer.

Comprehensive 3D epigenomic maps define limbal stem/progenitor cell function and identity

The insights into how genome topology couples with epigenetic states to govern the function and identity of the corneal epithelium are poorly understood. Here, we generate a high-resolution Hi-C interaction map of human limbal stem/progenitor cells (LSCs) and show that chromatin multi-hierarchical organisation is coupled to gene expression. By integrating Hi-C, epigenome and transcriptome data, we characterize the comprehensive 3D epigenomic landscapes of LSCs. We find that super-silencers mediate gene repression associated with corneal development, differentiation and disease via chromatin looping and/or proximity. Super-enhancer (SE) interaction analysis identified a set of SE interactive hubs that contribute to LSC-specific gene activation. These active and inactive element-anchored loop networks occur within the cohesin-occupied CTCF-CTCF loops. We further reveal a coordinated regulatory network of core transcription factors based on SE-promoter interactions. Our results provide detailed insights into the genome organization principle for epigenetic regulation of gene expression in stratified epithelia.

Genome topology provides a structural basis for epigenome-mediated transcriptional regulation in eukaryotes. Here the authors characterized the 3D genome of stratified squamous epithelia. They generated a Hi-C map of human limbal stem/progenitor cells (LSCs) and integrated these data with epigenomics, transcription factor binding maps, and transcriptome data.

Sex-specific involvement of the Notch-JAG pathway in social recognition

Under the hypothesis that olfactory neural epithelium gene expression profiles may be useful to look for disease-relevant neuronal signatures, we examined microarray gene expression in olfactory neuronal cells and underscored Notch-JAG pathway molecules in association with schizophrenia (SZ). The microarray profiling study underscored JAG1 as the most promising candidate. Combined with further validation with real-time PCR, downregulation of NOTCH1 was statistically significant. Accordingly, we reverse-translated the significant finding from a surrogate tissue for neurons, and studied the behavioral profile of Notch1^+/- mice. We found a specific impairment in social novelty recognition, whereas other behaviors, such as sociability, novel object recognition and olfaction of social odors, were normal. This social novelty recognition deficit was male-specific and was rescued by rapamycin treatment. Based on the results from the animal model, we next tested whether patients with psychosis might have male-specific alterations in social cognition in association with the expression of NOTCH1 or JAG1. In our first episode psychosis cohort, we observed a specific correlation between the expression of JAG1 and a face processing measure only in male patients. The expression of JAG1 was not correlated with any other cognitive and symptomatic scales in all subjects. Together, although we acknowledge the pioneering and exploratory nature, the present work that combines both human and animal studies in a reciprocal manner suggests a novel role for the Notch-JAG pathway in a behavioral dimension(s) related to social cognition in psychotic disorders in a male-specific manner.

Integrating Multimorbidity into a Whole-Body Understanding of Disease Using Spatial Genomics

Multimorbidity is characterized by multidimensional complexity emerging from interactions between multiple diseases across levels of biological (including genetic) and environmental determinants and the complex array of interactions between and within cells, tissues and organ systems. Advances in spatial genomic research have led to an unprecedented expansion in our ability to link alterations in genome folding with changes that are associated with human disease. Studying disease-associated genetic variants in the context of the spatial genome has enabled the discovery of transcriptional regulatory programmes that potentially link dysregulated genes to disease development. However, the approaches that have been used have typically been applied to uncover pathological molecular mechanisms occurring in a specific disease-relevant tissue. These forms of reductionist, targeted investigations are not appropriate for the molecular dissection of multimorbidity that typically involves contributions from multiple tissues. In this perspective, we emphasize the importance of a whole-body understanding of multimorbidity and discuss how spatial genomics, when integrated with additional omic datasets, could provide novel insights into the molecular underpinnings of multimorbidity.

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

From population to neuron: exploring common mediators for metabolic problems and mental illnesses

Major mental illnesses such as schizophrenia (SZ) and bipolar disorder (BP) frequently accompany metabolic conditions, but their relationship is still unclear, in particular at the mechanistic level. We implemented an approach of "from population to neuron", combining population-based epidemiological analysis with neurobiological experiments using cell and animal models based on a hypothesis built from the epidemiological study. We characterized high-quality population data, olfactory neuronal cells biopsied from patients with SZ or BP, and healthy subjects, as well as mice genetically modified for insulin signaling. We accessed the Danish Registry and observed (1) a higher incidence of diabetes in people with SZ or BP and (2) higher incidence of major mental illnesses in people with diabetes in the same large cohort. These epidemiological data suggest the existence of common pathophysiological mediators in both diabetes and major mental illnesses. We hypothesized that molecules associated with insulin resistance might be such common mediators, and then validated the hypothesis by using two independent sets of olfactory neuronal cells biopsied from patients and healthy controls. In the first set, we confirmed an enrichment of insulin signaling-associated molecules among the genes that were significantly different between SZ patients and controls in unbiased expression profiling data. In the second set, olfactory neuronal cells from SZ and BP patients who were not pre-diabetic or diabetic showed reduced IRS2 tyrosine phosphorylation upon insulin stimulation, indicative of insulin resistance. These cells also displayed an upregulation of IRS1 protein phosphorylation at serine-312 at baseline (without insulin stimulation), further supporting the concept of insulin resistance in olfactory neuronal cells from SZ patients. Finally, Irs2 knockout mice showed an aberrant response to amphetamine, which is also observed in some patients with major mental illnesses. The bi-directional relationships between major mental illnesses and diabetes suggest that there may be common pathophysiological mediators associated with insulin resistance underlying these mental and physical conditions.

Three-dimensional chromatin organization in brain function and dysfunction

The three-dimensional (3D) organization of chromatin within the nucleus is now recognized as a bona fide epigenetic property influencing genome function, replication, and maintenance. In the recent years, several studies have revealed how 3D chromatin organization is associated with brain function and its emerging role in disorders of the brain. 3D chromatin organization plays a crucial role in the development of different cell types of the nervous system and some neuronal cell types have adapted unique modifications to this organization that deviates from all other cell types. In post-mitotic neurons, dynamic changes in chromatin interactions in response to neuronal activity underlie learning and memory formation. Finally, new evidence directly links 3D chromatin organization to several disorders of the brain. These recent findings position 3D chromatin organization as a fundamental regulatory mechanism poised to reveal the etiology of brain function and dysfunctions.

Metacode: One code to rule them all

The code of codes or metacode is a microcosm where biological layers, as well as their codes, interact together allowing the continuity of information flow in organisms by increasing biological entities' complexity. Through this novel organic code, biological systems scale towards niches with higher informatic freedom building structures that increase the entropy in the universe. Code biology has developed a novel informational framework where biological entities strive themselves through the information flow carried out through organic codes consisting of two molecular or functional landscapes intertwined through arbitrary linkages via an adaptor whose nature is autonomous from molecular determinism. Here we will integrate genomic and epigenomic codes according to the evidence released in ENCODE (phase 3), psychENCODE and GTEx project, outlining the principles of the metacode, to address the continuous nature of biological systems and their inter-layered information flow. This novel complex metacode maps from very constrained sets of elements (i.e., regulation sites modulating gene expression) to new ones with greater freedom of decoding (i.e., a continuous cell phenotypic space). This leads to a new domain in code biology where biological systems are informatic attractors that navigate an energy metaspace through a complexity-noise balance, stalling in emergent niches where organic codes take meaning.

Characterisation of Neurospheres-Derived Cells from Human Olfactory Epithelium

A major problem in psychiatric research is a deficit of relevant cell material of neuronal origin, especially in large quantities from living individuals. One of the promising options is cells from the olfactory neuroepithelium, which contains neuronal progenitors that ensure the regeneration of olfactory receptors. These cells are easy to obtain with nasal biopsies and it is possible to grow and cultivate them in vitro. In this work, we used RNAseq expression profiling and immunofluorescence microscopy to characterise neurospheres-derived cells (NDC), that simply and reliably grow from neurospheres (NS) obtained from nasal biopsies. We utilized differential expression analysis to explore the molecular changes that occur during transition from NS to NDC. We found that processes associated with neuronal and vascular cells are downregulated in NDC. A comparison with public transcriptomes revealed a depletion of neuronal and glial components in NDC. We also discovered that NDC have several metabolic features specific to neuronal progenitors treated with the fungicide maneb. Thus, while NDC retain some neuronal/glial identity, additional protocol alterations are needed to use NDC for mass sample collection in psychiatric research.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

PsychENCODE and beyond: transcriptomics and epigenomics of brain development and organoids

Crucial decisions involving cell fate and connectivity that shape the distinctive development of the human brain occur in the embryonic and fetal stages-stages that are difficult to access and investigate in humans. The last decade has seen an impressive increase in resources-from atlases and databases to biological models-that is progressively lifting the curtain on this critical period. In this review, we describe the current state of genomic, transcriptomic, and epigenomic datasets charting the development of normal human brain with a particular focus on recent single-cell technologies. We discuss the emergence of brain organoids generated from pluripotent stem cells as a model to compensate for the limited availability of fetal tissue. Indeed, comparisons of neural lineages, transcriptional dynamics, and noncoding element activity between fetal brain and organoids have helped identify gene regulatory networks functioning at early stages of brain development. Altogether, we argue that large multi-omics investigations have pushed brain development into the "big data" era, and that current and future transversal approaches needed to leverage both fetal brain and organoid resources promise to answer major questions of brain biology and psychiatry.

On the in vivo origin of human nasal mesenchymal stem cell cultures

OBJECTIVES

Mesenchymal stem cells (MSCs), classically expanded in culture from bone marrow, are of broad interest to the regenerative medicine community. Human nasal turbinate mesenchymal-like stem cell cultures have also been described, defined by an in vitro phenotype similar to bone marrow MSCs. Nonetheless, the identity in vivo of the cells that give rise to nasal MSC-like cultures remains unclear, and these cells are often suggested to be related to olfactory lineages. Here, we sought to define the in vivo phenotype of human nasal MSC-like cells.

METHODS

Human turbinate tissue samples were used for RNA and immunohistochemical analysis. We also analyzed a recently published single cell RNA-sequencing dataset from adult human olfactory and respiratory mucosa samples from our lab, to focus on cell populations expressing MSC markers. Immunochemistry was performed to stain turbinate sections and nasal MSC cultures for selected markers.

RESULTS

While there is no single MSC-specific gene, we identified a human nasal mucosal cell population in vivo that uniquely expressed transcripts characteristic of typical MSC cultures, including ENG (CD105), NES, and CD34, and lacked expression of other transcripts associated with surface epithelia. The expression of transcription factors such as SOX17, EBF1, and FOXP1 suggests cells in the MSC-like cluster maintain an ability to direct cell fate, consistent with the behavior of nasal MSC-like cells in vitro. SOX17 was found to be uniformly expressed by nasal MSC cultures, consistent with the in vivo data. Immunohistochemistry of human nasal tissue samples indicated that ENG, CD34, and SOX17 expression localized selectively to cells surrounding blood vessels in the lamina propria.

CONCLUSION

Our findings provide evidence that the in vivo origin of nasal MSC-like cultures is likely a vascular or pericyte population, rather than cells related to the olfactory neuronal lineage.

LEVEL OF EVIDENCE

NA.

Rare Variant Burden Analysis within Enhancers Identifies CAV1 as an ALS Risk Gene

Amyotrophic lateral sclerosis (ALS) is an incurable neurodegenerative disease. CAV1 and CAV2 organize membrane lipid rafts (MLRs) important for cell signaling and neuronal survival, and overexpression of CAV1 ameliorates ALS phenotypes in vivo. Genome-wide association studies localize a large proportion of ALS risk variants within the non-coding genome, but further characterization has been limited by lack of appropriate tools. By designing and applying a pipeline to identify pathogenic genetic variation within enhancer elements responsible for regulating gene expression, we identify disease-associated variation within CAV1/CAV2 enhancers, which replicate in an independent cohort. Discovered enhancer mutations reduce CAV1/CAV2 expression and disrupt MLRs in patient-derived cells, and CRISPR-Cas9 perturbation proximate to a patient mutation is sufficient to reduce CAV1/CAV2 expression in neurons. Additional enrichment of ALS-associated mutations within CAV1 exons positions CAV1 as an ALS risk gene. We propose CAV1/CAV2 overexpression as a personalized medicine target for ALS.

TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

Use of the epigenetic toolbox to contextualize common variants associated with schizophrenia risk

Schizophrenia is a debilitating psychiatric disorder with a complex genetic architecture and limited understanding of its neuropathology, reflected by the lack of diagnostic measures and effective pharmacological treatments. Geneticists have recently identified more than 145 risk loci comprising hundreds of common variants of small effect sizes, most of which lie in noncoding genomic regions. This review will discuss how the epigenetic toolbox can be applied to contextualize genetic findings in schizophrenia. Progress in next-generation sequencing, along with increasing methodological complexity, has led to the compilation of genome-wide maps of DNA methylation, histone modifications, RNA expression, and more. Integration of chromatin conformation datasets is one of the latest efforts in deciphering schizophrenia risk, allowing the identification of genes in contact with regulatory variants across 100s of kilobases. Large-scale multiomics studies will facilitate the prioritization of putative causal risk variants and gene networks that contribute to schizophrenia etiology, informing clinical diagnostics and treatment downstream.
.

A Colorimetric Artificial Olfactory System for Airborne Improvised Explosive Identification

The detection of ultralow or nonvolatile target analytes remains a significant challenge for artificial olfactory systems even after decades of development, which severely limits their widespread application. To overcome this challenge, an artificial olfactory system based on a colorimetric hydrogel array is constructed for the first time as a universal representative. As an effective extension of conventional artificial olfactory systems that integrates the merits of its predecessors, the proposed system accurately mimics olfactory mucosa and specific odorant binding proteins using hydrogels endowed with specific colorimetric reagents for the detection of hypochlorite, chlorate, perchlorate, urea, and nitrate. Therefore, the proposed system is capable of detecting and discriminating between these five airborne improvised explosive microparticulates with a detection limit as low as 39.4 pg. Additionally, the system demonstrates good reusability over ten cycles, rapid response time of ≈0.2 s, and excellent discrimination properties, despite significant variation. This proof-of-concept study on colorimetric artificial olfactory systems yields a novel strategy for the direct and discriminative detection of nonvolatile airborne microparticulates.

A chromosomal connectome for psychiatric and metabolic risk variants in adult dopaminergic neurons

BACKGROUND

Midbrain dopaminergic neurons (MDN) represent 0.0005% of the brain's neuronal population and mediate cognition, food intake, and metabolism. MDN are also posited to underlay the neurobiological dysfunction of schizophrenia (SCZ), a severe neuropsychiatric disorder that is characterized by psychosis as well as multifactorial medical co-morbidities, including metabolic disease, contributing to markedly increased morbidity and mortality. Paradoxically, however, the genetic risk sequences of psychosis and traits associated with metabolic disease, such as body mass, show very limited overlap.

METHODS

We investigated the genomic interaction of SCZ with medical conditions and traits, including body mass index (BMI), by exploring the MDN's "spatial genome," including chromosomal contact landscapes as a critical layer of cell type-specific epigenomic regulation. Low-input Hi-C protocols were applied to 5-10 × 103 dopaminergic and other cell-specific nuclei collected by fluorescence-activated nuclei sorting from the adult human midbrain.

RESULTS

The Hi-C-reconstructed MDN spatial genome revealed 11 "Euclidean hot spots" of clustered chromatin domains harboring risk sequences for SCZ and elevated BMI. Inter- and intra-chromosomal contacts interconnecting SCZ and BMI risk sequences showed massive enrichment for brain-specific expression quantitative trait loci (eQTL), with gene ontologies, regulatory motifs and proteomic interactions related to adipogenesis and lipid regulation, dopaminergic neurogenesis and neuronal connectivity, and reward- and addiction-related pathways.

CONCLUSIONS

We uncovered shared nuclear topographies of cognitive and metabolic risk variants. More broadly, our PsychENCODE sponsored Hi-C study offers a novel genomic approach for the study of psychiatric and medical co-morbidities constrained by limited overlap of their respective genetic risk architectures on the linear genome.

Spatial genome exploration in the context of cognitive and neurological disease

The 'non-linear' genome, or the spatial proximity of non-contiguous sequences, emerges as an important regulatory layer for genome organization and function, including transcriptional regulation. Here, we review recent genome-scale chromosome conformation mappings ('Hi-C') in developing and adult human and mouse brain. Neural differentiation is associated with widespread remodeling of the chromosomal contact map, reflecting dynamic changes in cell-type-specific gene expression programs, with a massive (estimated 20-50%) net loss of chromosomal contacts that is specific for the neuronal lineage. Hi-C datasets provided an unexpected link between locus-specific abnormal expansion of repeat sequences positioned at the boundaries of self-associating topological chromatin domains, and monogenic neurodevelopmental and neurodegenerative disease. Furthermore, integrative cell-type-specific Hi-C and transcriptomic analysis uncovered an expanded genomic risk space for sequences conferring liability for schizophrenia and other cognitive disease. We predict that spatial genome exploration will deliver radically new insights into the brain nucleome in health and disease.

Druggable Transcriptional Networks in the Human Neurogenic Epigenome

Chromosome conformation capture methods have revealed the dynamics of genome architecture which is spatially organized into topologically associated domains, with gene regulation mediated by enhancer-promoter pairs in chromatin space. New evidence shows that endogenous hormones and several xenobiotics act within circumscribed topological domains of the spatial genome, impacting subsets of the chromatin contacts of enhancer-gene promoter pairs in cis and trans Results from the National Institutes of Health-funded PsychENCODE project and the study of chromatin remodeling complexes have converged to provide a clearer understanding of the organization of the neurogenic epigenome in humans. Neuropsychiatric diseases, including schizophrenia, bipolar spectrum disorder, autism spectrum disorder, attention deficit hyperactivity disorder, and other neuropsychiatric disorders are significantly associated with mutations in neurogenic transcriptional networks. In this review, we have reanalyzed the results from publications of the PsychENCODE Consortium using pharmacoinformatics network analysis to better understand druggable targets that control neurogenic transcriptional networks. We found that valproic acid and other psychotropic drugs directly alter these networks, including chromatin remodeling complexes, transcription factors, and other epigenetic modifiers. We envision a new generation of CNS therapeutics targeted at neurogenic transcriptional control networks, including druggable parts of chromatin remodeling complexes and master transcription factor-controlled pharmacogenomic networks. This may provide a route to the modification of interconnected gene pathways impacted by disease in patients with neuropsychiatric and neurodegenerative disorders. Direct and indirect therapeutic strategies to modify the master regulators of neurogenic transcriptional control networks may ultimately help extend the life span of CNS neurons impacted by disease.

A high-resolution 3D epigenomic map reveals insights into the creation of the prostate cancer transcriptome

To better understand the impact of chromatin structure on regulation of the prostate cancer transcriptome, we develop high-resolution chromatin interaction maps in normal and prostate cancer cells using in situ Hi-C. By combining the in situ Hi-C data with active and repressive histone marks, CTCF binding sites, nucleosome-depleted regions, and transcriptome profiling, we identify topologically associating domains (TADs) that change in size and epigenetic states between normal and prostate cancer cells. Moreover, we identify normal and prostate cancer-specific enhancer-promoter loops and involved transcription factors. For example, we show that FOXA1 is enriched in prostate cancer-specific enhancer-promoter loop anchors. We also find that the chromatin structure surrounding the androgen receptor (AR) locus is altered in the prostate cancer cells with many cancer-specific enhancer-promoter loops. This creation of 3D epigenomic maps enables a better understanding of prostate cancer biology and mechanisms of gene regulation.

In prostate cancer, chromatin structure can impact the transcriptome. Here, the authors develop high resolution chromatin interaction maps in prostate cancer cells using in situ Hi-C, revealing prostate cancer-specific TADs and enhancer-promoter loops surrounding the androgen receptor (AR) locus.

Using Openly Accessible Resources to Strengthen Causal Inference in Epigenetic Epidemiology of Neurodevelopment and Mental Health

The recent focus on the role of epigenetic mechanisms in mental health has led to several studies examining the association of epigenetic processes with psychiatric conditions and neurodevelopmental traits. Some studies suggest that epigenetic changes might be causal in the development of the psychiatric condition under investigation. However, other scenarios are possible, e.g., statistical confounding or reverse causation, making it particularly challenging to derive conclusions on causality. In the present review, we examine the evidence from human population studies for a possible role of epigenetic mechanisms in neurodevelopment and mental health and discuss methodological approaches on how to strengthen causal inference, including the need for replication, (quasi-)experimental approaches and Mendelian randomization. We signpost openly accessible resources (e.g., "MR-Base" "EWAS catalog" as well as tissue-specific methylation and gene expression databases) to aid the application of these approaches.

Transcriptome and epigenome landscape of human cortical development modeled in organoids

Genes implicated in neuropsychiatric disorders are active in human fetal brain, yet difficult to study in a longitudinal fashion. We demonstrate that organoids from human pluripotent cells model cerebral cortical development on the molecular level before 16 weeks postconception. A multiomics analysis revealed differentially active genes and enhancers, with the greatest changes occurring at the transition from stem cells to progenitors. Networks of converging gene and enhancer modules were assembled into six and four global patterns of expression and activity across time. A pattern with progressive down-regulation was enriched with human-gained enhancers, suggesting their importance in early human brain development. A few convergent gene and enhancer modules were enriched in autism-associated genes and genomic variants in autistic children. The organoid model helps identify functional elements that may drive disease onset.