The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development.

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome.

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

The neuronal chromatin landscape in adult schizophrenia brains is linked to early fetal development.

Non-coding variants increase risk of neuropsychiatric disease. However, our understanding of the cell-type specific role of the non-coding genome in disease is incomplete. We performed population scale (N=1,393) chromatin accessibility profiling of neurons and non-neurons from two neocortical brain regions: the anterior cingulate cortex and dorsolateral prefrontal cortex. Across both regions, we observed notable differences in neuronal chromatin accessibility between schizophrenia cases and controls. A per-sample disease pseudotime was positively associated with genetic liability for schizophrenia. Organizing chromatin into cis- and trans-regulatory domains, identified a prominent neuronal trans-regulatory domain (TRD1) active in immature glutamatergic neurons during fetal development. Polygenic risk score analysis using genetic variants within chromatin accessibility of TRD1 successfully predicted susceptibility to schizophrenia in the Million Veteran Program cohort. Overall, we present the most extensive resource to date of chromatin accessibility in the human cortex, yielding insights into the cell-type specific etiology of schizophrenia.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level.

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease.

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis.

Population-level variation in enhancer expression identifies disease mechanisms in the human brain.

Identification of risk variants for neuropsychiatric diseases within enhancers underscores the importance of understanding population-level variation in enhancer function in the human brain. Besides regulating tissue-specific and cell-type-specific transcription of target genes, enhancers themselves can be transcribed. By jointly analyzing large-scale cell-type-specific transcriptome and regulome data, we cataloged 30,795 neuronal and 23,265 non-neuronal candidate transcribed enhancers. Examination of the transcriptome in 1,382 brain samples identified robust expression of transcribed enhancers. We explored gene-enhancer coordination and found that enhancer-linked genes are strongly implicated in neuropsychiatric disease. We identified expression quantitative trait loci (eQTLs) for both genes and enhancers and found that enhancer eQTLs mediate a substantial fraction of neuropsychiatric trait heritability. Inclusion of enhancer eQTLs in transcriptome-wide association studies enhanced functional interpretation of disease loci. Overall, our study characterizes the gene-enhancer regulome and genetic mechanisms in the human cortex in both healthy and diseased states.

Genetics of the human microglia regulome refines Alzheimer's disease risk loci.

Microglia are brain myeloid cells that play a critical role in neuroimmunity and the etiology of Alzheimer's disease (AD), yet our understanding of how the genetic regulatory landscape controls microglial function and contributes to AD is limited. Here, we performed transcriptome and chromatin accessibility profiling in primary human microglia from 150 donors to identify genetically driven variation and cell-specific enhancer-promoter (E-P) interactions. Integrative fine-mapping analysis identified putative regulatory mechanisms for 21 AD risk loci, of which 18 were refined to a single gene, including 3 new candidate risk genes (KCNN4, FIBP and LRRC25). Transcription factor regulatory networks captured AD risk variation and identified SPI1 as a key putative regulator of microglia expression and AD risk. This comprehensive resource capturing variation in the human microglia regulome provides insights into the etiology of neurodegenerative disease.

The Neuroepigenome: Implications of Chemical and Physical Modifications of Genomic DNA in Schizophrenia.

Schizophrenia is a chronic mental illness with a substantial genetic component. To unfold the complex etiology of schizophrenia, it is important to understand the interplay between genetic and nongenetic factors. Genetic factors involve variation in the DNA sequences of protein-coding genes, which directly contribute to phenotypic traits, and variation in noncoding sequences, which comprise 98% of the genome and contain DNA elements known to play a role in regulating gene expression. The epigenome refers to the chemical modifications on both DNA and the structural proteins that package DNA into the nucleus, which together regulate gene expression in specific cell types, conditions, and developmental stages. The dynamic nature of the epigenome makes it an ideal tool to investigate the relationship between inherited genetic mutations associated with schizophrenia and altered gene regulation throughout the course of brain development. In this review, we focus on the current understanding of the role of epigenetic marks and their three-dimensional nuclear organization in the developmental trajectory of distinct brain cell types to decipher the complex gene regulatory mechanisms that are disrupted in schizophrenia.

Chromatin domain alterations linked to 3D genome organization in a large cohort of schizophrenia and bipolar disorder brains.

Chromosomal organization, scaling from the 147-base pair (bp) nucleosome to megabase-ranging domains encompassing multiple transcriptional units, including heritability loci for psychiatric traits, remains largely unexplored in the human brain. In this study, we constructed promoter- and enhancer-enriched nucleosomal histone modification landscapes for adult prefrontal cortex from H3-lysine 27 acetylation and H3-lysine 4 trimethylation profiles, generated from 388 controls and 351 individuals diagnosed with schizophrenia (SCZ) or bipolar disorder (BD) (n = 739). We mapped thousands of cis-regulatory domains (CRDs), revealing fine-grained, 104-106-bp chromosomal organization, firmly integrated into Hi-C topologically associating domain stratification by open/repressive chromosomal environments and nuclear topography. Large clusters of hyper-acetylated CRDs were enriched for SCZ heritability, with prominent representation of regulatory sequences governing fetal development and glutamatergic neuron signaling. Therefore, SCZ and BD brains show coordinated dysregulation of risk-associated regulatory sequences assembled into kilobase- to megabase-scaling chromosomal domains.

Breast Tumor Classification in Ultrasound Images Using Combined Deep and Handcrafted Features.

This study aims to enable effective breast ultrasound image classification by combining deep features with conventional handcrafted features to classify the tumors. In particular, the deep features are extracted from a pre-trained convolutional neural network model, namely the VGG19 model, at six different extraction levels. The deep features extracted at each level are analyzed using a features selection algorithm to identify the deep feature combination that achieves the highest classification performance. Furthermore, the extracted deep features are combined with handcrafted texture and morphological features and processed using features selection to investigate the possibility of improving the classification performance. The cross-validation analysis, which is performed using 380 breast ultrasound images, shows that the best combination of deep features is obtained using a feature set, denoted by CONV features that include convolution features extracted from all convolution blocks of the VGG19 model. In particular, the CONV features achieved mean accuracy, sensitivity, and specificity values of 94.2%, 93.3%, and 94.9%, respectively. The analysis also shows that the performance of the CONV features degrades substantially when the features selection algorithm is not applied. The classification performance of the CONV features is improved by combining these features with handcrafted morphological features to achieve mean accuracy, sensitivity, and specificity values of 96.1%, 95.7%, and 96.3%, respectively. Furthermore, the cross-validation analysis demonstrates that the CONV features and the combined CONV and morphological features outperform the handcrafted texture and morphological features as well as the fine-tuned VGG19 model. The generalization performance of the CONV features and the combined CONV and morphological features is demonstrated by performing the training using the 380 breast ultrasound images and the testing using another dataset that includes 163 images. The results suggest that the combined CONV and morphological features can achieve effective breast ultrasound image classifications that increase the capability of detecting malignant tumors and reduce the potential of misclassifying benign tumors.

Genetic Variation in Long-Range Enhancers.

Cis-regulatory elements (CREs), including insulators, promoters, and enhancers, play critical roles in the establishment and maintenance of normal cellular function. Within each cell, the 3D structure of chromatin is arranged in specific patterns to expose the CREs required for optimal spatiotemporal regulation of gene expression. CREs can act over large distances along the linear genome, facilitated by looping of the intervening chromatin to allow direct interaction between distal regulatory elements and their target genes. A number of pathologies are associated with dysregulation of CRE function, including developmental disorders, cancers, and neuropsychiatric disease. A majority of known neuropsychiatric disease risk loci are noncoding, and increasing evidence suggests that they contribute to disease through disruption of CREs. As such, rather than directly altering the amino acid content of proteins, these variants are instead thought to affect where, when, and to what extent a given gene is expressed. The distances over which CREs can operate often render their target genes difficult to identify. Furthermore, as many risk loci contain multiple variants in high linkage disequilibrium, identification of the causative single nucleotide polymorphism(s) therein is not straightforward. Thus, deciphering the genetic etiology of complex neuropsychiatric disorders presents a significant challenge.

Spatial Organization of Single mRNPs at Different Stages of the Gene Expression Pathway.

mRNAs form ribonucleoprotein complexes (mRNPs) by association with proteins that are crucial for mRNA metabolism. While the mRNP proteome has been well characterized, little is known about mRNP organization. Using a single-molecule approach, we show that mRNA conformation changes depending on its cellular localization and translational state. Compared to nuclear mRNPs and lncRNPs, association with ribosomes decompacts individual mRNAs, while pharmacologically dissociating ribosomes or sequestering them into stress granules leads to increased compaction. Moreover, translating mRNAs rarely show co-localized 5' and 3' ends, indicating either that mRNAs are not translated in a closed-loop configuration, or that mRNA circularization is transient, suggesting that a stable closed-loop conformation is not a universal state for all translating mRNAs.

Measuring mRNA Decay in Budding Yeast Using Single Molecule FISH.

Cellular mRNA levels are determined by the rates of mRNA synthesis and mRNA decay. Typically, mRNA degradation kinetics are measured on a population of cells that are either chemically treated or genetically engineered to inhibit transcription. However, these manipulations can affect the mRNA decay process itself by inhibiting regulatory mechanisms that govern mRNA degradation, especially if they occur on short time-scales. Recently, single molecule fluorescent in situ hybridization (smFISH) approaches have been implemented to quantify mRNA decay rates in single, unperturbed cells. Here, we provide a step-by-step protocol that allows quantification of mRNA decay in single Saccharomyces cerevisiae using smFISH. Our approach relies on fluorescent labeling of single cytoplasmic mRNAs and nascent mRNAs found at active sites of transcription, coupled with mathematical modeling to derive mRNA half-lives. Commercially available, single-stranded smFISH DNA oligonucleotides (smFISH probes) are used to fluorescently label mRNAs followed by the quantification of cellular and nascent mRNAs using freely available spot detection algorithms. Our method enables quantification of mRNA decay of any mRNA in single, unperturbed yeast cells and can be implemented to quantify mRNA turnover in a variety of cell types as well as tissues.

Single-cell profiling reveals that eRNA accumulation at enhancer-promoter loops is not required to sustain transcription.

Enhancers are intergenic DNA elements that regulate the transcription of target genes in response to signaling pathways by interacting with promoters over large genomic distances. Recent studies have revealed that enhancers are bi-directionally transcribed into enhancer RNAs (eRNAs). Using single-molecule fluorescence in situ hybridization (smFISH), we investigated the eRNA-mediated regulation of transcription during estrogen induction in MCF-7 cells. We demonstrate that eRNAs are localized exclusively in the nucleus and are induced with similar kinetics as target mRNAs. However, eRNAs are mostly nascent at enhancers and their steady-state levels remain lower than those of their cognate mRNAs. Surprisingly, at the single-allele level, eRNAs are rarely co-expressed with their target loci, demonstrating that active gene transcription does not require the continuous transcription of eRNAs or their accumulation at enhancers. When co-expressed, sub-diffraction distance measurements between nascent mRNA and eRNA signals reveal that co-transcription of eRNAs and mRNAs rarely occurs within closed enhancer-promoter loops. Lastly, basal eRNA transcription at enhancers, but not E2-induced transcription, is maintained upon depletion of MLL1 and ERα, suggesting some degree of chromatin accessibility prior to signal-dependent activation of transcription. Together, our findings suggest that eRNA accumulation at enhancer-promoter loops is not required to sustain target gene transcription.

The nuclear basket mediates perinuclear mRNA scanning in budding yeast.

After synthesis and transit through the nucleus, messenger RNAs (mRNAs) are exported to the cytoplasm through the nuclear pore complex (NPC). At the NPC, messenger ribonucleoproteins (mRNPs) first encounter the nuclear basket where mRNP rearrangements are thought to allow access to the transport channel. Here, we use single mRNA resolution live cell microscopy and subdiffraction particle tracking to follow individual mRNAs on their path toward the cytoplasm. We show that when reaching the nuclear periphery, RNAs are not immediately exported but scan along the nuclear periphery, likely to find a nuclear pore allowing export. Deletion or mutation of the nuclear basket proteins MLP1/2 or the mRNA binding protein Nab2 changes the scanning behavior of mRNPs at the nuclear periphery, shortens residency time at nuclear pores, and results in frequent release of mRNAs back into the nucleoplasm. These observations suggest a role for the nuclear basket in providing an interaction platform that keeps RNAs at the periphery, possibly to allow mRNP rearrangements before export.

Does NETosis Contribute to the Bacterial Pathoadaptation in Cystic Fibrosis?

Significant advances in our understanding of neutrophil biology were made in the past several years. The exciting discovery that neutrophils deploy neutrophil extracellular traps (NETs) to catch pathogens paved the way for a series of additional studies to define the molecular mechanisms of NET generation and the biological significance of NETosis in acute and chronic pathologic conditions. This review highlights the latest knowledge regarding NET structures, deployment, and function, with an emphasis on current understanding of NET proteomes, their conservation, and significance in the context of cystic fibrosis (CF), a condition characterized by excessive extracellular DNA/NET presence. We also discuss how our understanding of NETosis yields novel therapeutic approaches and their applicability to CF.

Distinct susceptibilities of corneal Pseudomonas aeruginosa clinical isolates to neutrophil extracellular trap-mediated immunity.

Ocular bacterial keratitis, often associated with Pseudomonas aeruginosa bacterial infection, commonly occurs in contact lens wearers and may lead to vision impairment. In this study, we analyzed the contribution of neutrophil extracellular traps (NETs) to the mediation of protection during ocular keratitis. Both invasive and cytotoxic P. aeruginosa clinical isolates induced NET release by neutrophils. NETs carried the characteristic histone proteins, elastase, lysozyme, myeloperoxidase, and metabolic enzymes. While the invasive P. aeruginosa strains PAO1 (serogroup O5) and 6294 (serogroup O6) were trapped by NETs, the cytotoxic P. aeruginosa strains 6077, 6206 (serogroup O11), and PA14 (serogroup 010) were less sensitive to NET capture. The mechanism of escape by the cytotoxic strains from adhesion to NETs involved the shedding of outer membrane vesicles (OMVs) that outcompeted the cytotoxic P. aeruginosa strains for NET binding. When ocular infection was caused by an invasive strain in vivo, NETs were released at the ocular surface to capture bacteria, limiting their spread. Treatment with MNase I had a dose-dependent effect, with low doses of MNase speeding up bacterial clearance and high doses of MNase having toxic consequences. Cumulatively, our data suggest that NET-mediated immunity is a two-step process. Initially, pathogens attach to NET fragments; subsequently, upon nuclease activity, active serine proteases, which proteolytically degrade NET-associated proteins and promote DNase activity, are released. Therefore, a balance between NET production and NET degradation is needed to achieve maximal NET immunity.

Single-molecule resolution fluorescent in situ hybridization (smFISH) in the yeast S. cerevisiae.

Regulating gene expression is a major task for all cellular systems. RNA production and degradation plays a critical role in this process and accurately measuring cellular mRNA levels is essential to understanding gene expression regulation. Classical biochemical assays that study gene expression rely on extracting RNAs from large populations of cells, taking them out of their native context and thereby losing spatial information as well as cell-to-cell variability. In this chapter, we describe a fluorescent in situ hybridization (FISH) technique that circumvents this problem by detecting single RNAs in single cells. The technique employs multiple single-stranded short DNA probes fluorescently labeled with organic dyes that hybridize to target RNAs in fixed cells, allowing quantification and localization of RNAs at the single-cell level and at single-molecule resolution. The protocol described here has been optimized for the yeast S. cerevisiae.

Bimodal expression of PHO84 is modulated by early termination of antisense transcription.

Many Saccharomyces cerevisiae genes encode antisense transcripts, some of which are unstable and degraded by the exosome component Rrp6. Loss of Rrp6 results in the accumulation of long PHO84 antisense (AS) RNAs and repression of sense transcription through PHO84 promoter deacetylation. We used single-molecule resolution fluorescent in situ hybridization (smFISH) to investigate antisense-mediated transcription regulation. We show that PHO84 AS RNA acts as a bimodal switch, in which continuous, low-frequency antisense transcription represses sense expression within individual cells. Surprisingly, antisense RNAs do not accumulate at the PHO84 gene but are exported to the cytoplasm. Furthermore, rather than stabilizing PHO84 AS RNA, the loss of Rrp6 favors its elongation by reducing early transcription termination by Nrd1-Nab3-Sen1. These observations suggest that PHO84 silencing results from antisense transcription through the promoter rather than the static accumulation of antisense RNAs at the repressed gene.