Regional differences in gene expression and promoter usage in aged human brains

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Molecular tools for recording and intervention of neuronal activity

Observing the activity of neural networks is critical for the identification of learning and memory processes, as well as abnormal activities of neural circuits in disease, particularly for the purpose of tracking disease progression. Methodologies for describing the activity history of neural networks using molecular biology techniques first utilized genes expressed by active neurons, followed by the application of recently developed techniques including optogenetics and incorporation of insights garnered from other disciplines, including chemistry and physics. In this review, we will discuss ways in which molecular biological techniques used to describe the activity of neural networks have evolved along with the potential for future development.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders, yet the role of cell-type-specific splicing or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), and uncovered a substantial contribution of transcript-isoform diversity, regulated by RNA binding proteins, in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

One-Sentence Summary

A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease.

Structured Abstract

INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

Altered Hippocampal and Striatal Expression of Endothelial Markers and VIP/PACAP Neuropeptides in a Mouse Model of Systemic Lupus Erythematosus

Neuropsychiatric systemic lupus erythematosus (NPSLE) is one of the most common and severe manifestations of lupus; however, its pathogenesis is still poorly understood. While there is sparse evidence suggesting that the ongoing autoimmunity may trigger pathogenic changes to the central nervous system (CNS) microvasculature, culminating in inflammatory/ischemic damage, further evidence is still needed. In this study, we used the spontaneous mouse model of SLE (NZBWF1 mice) to investigate the expression of genes and proteins associated with endothelial (dys)function: tissue and urokinase plasminogen activators (tPA and uPA), intercellular and vascular adhesion molecules 1 (ICAM-1 and VCAM-1), brain derived neurotrophic factor (BDNF), endothelial nitric oxide synthase (eNOS) and Krüppel-like factor 4 (KLF4) and neuroprotection/immune modulation: pituitary adenylate cyclase-activating peptide (PACAP), vasoactive intestinal peptide (VIP), PACAP receptor (PAC1), VIP receptors 1 and 2 (VPAC1 and VPAC2). Analyses were carried out both in the hippocampus and striatum of SLE mice of two different age groups (2 and 7 months old), since age correlates with disease severity. In the hippocampus, we identified a gene/protein expression profile indicative of mild endothelial dysfunction, which increased in severity in aged SLE mice. These alterations were paralleled by moderate alterations in the expression of VIP, PACAP and related receptors. In contrast, we report a robust upregulation of endothelial activation markers in the striatum of both young and aged mice, concurrent with significant induction of the VIP/PACAP system. These data identify molecular signatures of endothelial alterations in the hippocampus and striatum of NZBWF1 mice, which are accompanied by a heightened expression of endogenous protective/immune-modulatory neuropeptides. Collectively, our results support the idea that NPSLE may cause alterations of the CNS micro-vascular compartment that cannot be effectively counteracted by the endogenous activity of the neuropeptides PACAP and VIP.

DNA Methylation Signature of Aging: Potential Impact on the Pathogenesis of Parkinson's Disease

Regulation of gene expression by epigenetic modifications means lasting and heritable changes in the function of genes without alterations in the DNA sequence. Of all epigenetic mechanisms identified thus far, DNA methylation has been of particular interest in both aging and age-related disease research over the last decade given the consistency of site-specific DNA methylation changes during aging that can predict future health and lifespan. An increasing line of evidence has implied the dynamic nature of DNA (de)methylation events that occur throughout the lifespan has a role in the pathophysiology of aging and age-associated neurodegenerative conditions, including Parkinson's disease (PD). In this regard, PD methylome shows, to some extent, similar genome-wide changes observed in the methylome of healthy individuals of matching age. In this review, we start by providing a brief overview of studies outlining global patterns of DNA methylation, then its mechanisms and regulation, within the context of aging and PD. Considering diverging lines of evidence from different experimental and animal models of neurodegeneration and how they combine to shape our current understanding of tissue-specific changes in DNA methylome in health and disease, we report a high-level comparison of the genomic methylation landscapes of brain, with an emphasis on dopaminergic neurons in PD and in natural aging. We believe this will be particularly useful for systematically dissecting overlapping genome-wide alterations in DNA methylation during PD and healthy aging, and for improving our knowledge of PD-specific changes in methylation patterns independent of aging process.

Tau mRNA Metabolism in Neurodegenerative Diseases: A Tangle Journey

Tau proteins are known to be mainly involved in regulation of microtubule dynamics. Besides this function, which is critical for axonal transport and signal transduction, tau proteins also have other roles in neurons. Moreover, tau proteins are turned into aggregates and consequently trigger many neurodegenerative diseases termed tauopathies, of which Alzheimer’s disease (AD) is the figurehead. Such pathological aggregation processes are critical for the onset of these diseases. Among the various causes of tau protein pathogenicity, abnormal tau mRNA metabolism, expression and dysregulation of tau post-translational modifications are critical steps. Moreover, the relevance of tau function to general mRNA metabolism has been highlighted recently in tauopathies. In this review, we mainly focus on how mRNA metabolism impacts the onset and development of tauopathies. Thus, we intend to portray how mRNA metabolism of, or mediated by, tau is associated with neurodegenerative diseases.

Spatiotemporal Regulation of Transcript Isoform Expression in the Hippocampus

Proper development and plasticity of hippocampal neurons require specific RNA isoforms to be expressed in the right place at the right time. Precise spatiotemporal transcript regulation requires the incorporation of essential regulatory RNA sequences into expressed isoforms. In this review, we describe several RNA processing strategies utilized by hippocampal neurons to regulate the spatiotemporal expression of genes critical to development and plasticity. The works described here demonstrate how the hippocampus is an ideal investigative model for uncovering alternate isoform-specific mechanisms that restrict the expression of transcripts in space and time.

DNA methylation changes in clonally propagated oil palm

Several hypomethylated sites within the Karma region of EgDEF1 and hotspot regions in chromosomes 1, 2, 3, and 5 may be associated with mantling. One of the main challenges faced by the oil palm industry is fruit abnormalities, such as the "mantled" phenotype that can lead to reduced yields. This clonal abnormality is an epigenetic phenomenon and has been linked to the hypomethylation of a transposable element within the EgDEF1 gene. To understand the epigenome changes in clones, methylomes of clonal oil palms were compared to methylomes of seedling-derived oil palms. Whole-genome bisulfite sequencing data from seedlings, normal, and mantled clones were analyzed to determine and compare the context-specific DNA methylomes. In seedlings, coding and regulatory regions are generally hypomethylated while introns and repeats are extensively methylated. Genes with a low number of guanines and cytosines in the third position of codons (GC₃-poor genes) were increasingly methylated towards their 3' region, while GC₃-rich genes remain demethylated, similar to patterns in other eukaryotic species. Predicted promoter regions were generally hypomethylated in seedlings. In clones, CG, CHG, and CHH methylation levels generally decreased in functionally important regions, such as promoters, 5' UTRs, and coding regions. Although random regions were found to be hypomethylated in clonal genomes, hypomethylation of certain hotspot regions may be associated with the clonal mantling phenotype. Our findings, therefore, suggest other hypomethylated CHG sites within the Karma of EgDEF1 and hypomethylated hotspot regions in chromosomes 1, 2, 3 and 5, are associated with mantling.

Glucocerebrosidase regulators SCARB2 and TFEB are up-regulated in Lewy body disease brain

Mutations in the glucocerebrosidase (GCase) gene (GBA) and GCase deficiency are major risk factors for Lewy body diseases. Decreased GCase activity enhances alpha-synuclein aggregation and disease development. Lysosomal integral membrane protein type 2, encoded by SCARB2, binds GCase targeting it to lysosomes and transcription factor EB (Tfeb) regulates lysosomal proteostasis. Our aim was to find out if GCase deficiency in Lewy body diseases is accompanied by SCARB2 and TFEB deregulation at the transcriptional level involving alternative splicing as well. Relative mRNA expression of two SCARB2 and two TFEB transcripts was studied by real-time PCR in post-mortem brain samples of cases with pure Lewy body pathology (LBP), cases with concomitant LBP and Alzheimer disease-like pathology, and controls. TFEB expression was increased in the temporal cortex and caudate nucleus of LBP cases, and SCARB2 was differentially expressed. Female-gender associated overexpression of all transcripts was found in the caudate nucleus, and disease duration associated TFEB expression changes in the temporal cortex. SCARB2 and TFEB expression correlated negatively with GBA mRNA expression in the temporal cortex. Our findings show disease-specific deregulation of TFEB and SCARB2 expression affecting alternative promoter usage and alternative splicing in Lewy body diseases.

Factors other than hTau overexpression that contribute to tauopathy-like phenotype in rTg4510 mice

The tauopathy-like phenotype observed in the rTg4510 mouse line, in which human tau_P301L expression specifically within the forebrain can be temporally controlled, has largely been attributed to high overexpression of mutant human tau in the forebrain region. Unexpectedly, we found that in a different mouse line with a targeted-insertion of the same transgene driven by the same tetracycline-TransActivator (tTA) allele, but with even higher overexpression of tauP301L than rTg4510, atrophy and tau histopathology are delayed, and a different behavioral profile is observed. This suggests that it is not overexpression of mutant human tau alone that contributes to the phenotype in rTg4510 mice. Furthermore we show that the tauopathy-like phenotype seen in rTg4510 requires a ~70-copy tau-transgene insertion in a 244 kb deletion in Fgf14, a ~7-copy tTA-transgene insertion in a 508 kb deletion that disrupts another five genes, in addition to high transgene overexpression. We propose that these additional effects need to be accounted for in any studies using rTg4510.

CpG traffic lights are markers of regulatory regions in human genome

Background

DNA methylation is involved in the regulation of gene expression. Although bisulfite-sequencing based methods profile DNA methylation at a single CpG resolution, methylation levels are usually averaged over genomic regions in the downstream bioinformatic analysis.

Results

We demonstrate that on the genome level a single CpG methylation can serve as a more accurate predictor of gene expression than an average promoter / gene body methylation. We define CpG traffic lights (CpG TL) as CpG dinucleotides with a significant correlation between methylation and expression of a gene nearby. CpG TL are enriched in all regulatory regions. Among all promoters, CpG TL are especially enriched in poised ones, suggesting involvement of DNA methylation in their regulation. Yet, binding of only a handful of transcription factors, such as NRF1, ETS, STAT and IRF-family members, could be regulated by direct methylation of transcription factor binding sites (TFBS) or its close proximity. For the majority of TF, an alternative scenario is more likely: methylation and inactivation of the whole regulatory element indirectly represses functional TF binding with a CpG TL being a reliable marker of such inactivation.

Conclusions

CpG TL provide a promising insight into mechanisms of enhancer activity and gene regulation linking methylation of single CpG to gene expression. CpG TL methylation can be used as reliable markers of enhancer activity and gene expression in applications, e.g. in clinic where measuring DNA methylation is easier compared to directly measuring gene expression due to more stable nature of DNA.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5387-1) contains supplementary material, which is available to authorized users.

An alternative transcript of the Alzheimer's disease risk gene SORL1 encodes a truncated receptor

SORL1 encodes a 250-kDa protein named sorLA, a functional sorting receptor for the amyloid precursor protein (APP). Several single nucleotide polymorphisms of the gene SORL1, encoding sorLA, are genetically associated with Alzheimer's disease (AD). In the existing literature, SORL1 is insufficiently described at the transcriptional level, and there is very limited amount of functional data defining different transcripts. We have characterized a SORL1 transcript containing a novel exon 30B. The transcript is expressed in most brain regions with highest expression in the temporal lobe and hippocampus. Exon 30B is spliced to exon 31, leading to a mature transcript that encodes an 829 amino acid sorLA receptor. This receptor variant lacks the binding site for APP and is unlikely to function in APP sorting. This transcript is expressed in equal amounts in the cerebellum from AD and non-AD individuals. Our data describe a transcript that encodes a truncated sorLA receptor, suggesting novel neuronal functions for sorLA and that alternative transcription provides a mechanism for SORL1 activity regulation.

Long Non-Coding RNAs in Neuronal Aging

The expansion of long non-coding RNAs (lncRNAs) in organismal genomes has been associated with the emergence of sophisticated regulatory networks that may have contributed to more complex neuronal processes, such as higher-order cognition. In line with the important roles of lncRNAs in the normal functioning of the human brain, dysregulation of lncRNA expression has been implicated in aging and age-related neurodegenerative disorders. In this paper, we discuss the function and expression of known neuronal-associated lncRNAs, their impact on epigenetic changes, the contribution of transposable elements to lncRNA expression, and the implication of lncRNAs in maintaining the 3D nuclear architecture in neurons. Moreover, we discuss how the complex molecular processes that are orchestrated by lncRNAs in the aged brain may contribute to neuronal pathogenesis by promoting protein aggregation and neurodegeneration. Finally, this review explores the possibility that age-related disturbances of lncRNA expression change the genomic and epigenetic regulatory landscape of neurons, which may affect neuronal processes such as neurogenesis and synaptic plasticity.

Alternative promoter usage generates novel shorter MAPT mRNA transcripts in Alzheimer's disease and progressive supranuclear palsy brains

Alternative promoter usage is an important mechanism for transcriptome diversity and the regulation of gene expression. Indeed, this alternative usage may influence tissue/subcellular specificity, protein translation and function of the proteins. The existence of an alternative promoter for MAPT gene was considered for a long time to explain differential tissue specificity and differential response to transcription and growth factors between mRNA transcripts. The alternative promoter usage could explain partly the different tau proteins expression patterns observed in tauopathies. Here, we report on our discovery of a functional alternative promoter for MAPT, located upstream of the gene's second exon (exon 1). By analyzing genome databases and brain tissue from control individuals and patients with Alzheimer's disease or progressive supranuclear palsy, we identified novel shorter transcripts derived from this alternative promoter. These transcripts are increased in patients' brain tissue as assessed by 5'RACE-PCR and qPCR. We suggest that these new MAPT isoforms can be translated into normal or amino-terminal-truncated tau proteins. We further suggest that activation of MAPT's alternative promoter under pathological conditions leads to the production of truncated proteins, changes in protein localization and function, and thus neurodegeneration.

Comprehensive promoter level expression quantitative trait loci analysis of the human frontal lobe

BACKGROUND

Expression quantitative trait loci (eQTL) analysis is a powerful method to detect correlations between gene expression and genomic variants and is widely used to interpret the biological mechanism underlying identified genome wide association studies (GWAS) risk loci. Numerous eQTL studies have been performed on different cell types and tissues of which the majority has been based on microarray technology.

METHODS

We present here an eQTL analysis based on cap analysis gene expression sequencing (CAGEseq) data created from human postmortem frontal lobe tissue combined with genotypes obtained through genotyping arrays, exome sequencing, and CAGEseq. Using CAGEseq as an expression profiling technique combined with these different genotyping techniques allows measurement of the molecular effect of variants on individual transcription start sites and increases the resolution of eQTL analysis by also including the non-annotated parts of the genome.

RESULTS

We identified 2410 eQTLs and show that non-coding transcripts are more likely to contain an eQTL than coding transcripts, in particular antisense transcripts. We provide evidence for how previously identified GWAS loci for schizophrenia (NRGN), Parkinson's disease, and Alzheimer's disease (PARK16 and MAPT loci) could increase the risk for disease at a molecular level. Furthermore, we demonstrate that CAGEseq improves eQTL analysis because variants obtained from CAGEseq are highly enriched for having a functional effect and thus are an efficient method towards the identification of causal variants.

CONCLUSION

Our data contain both coding and non-coding transcripts and has the added value that we have identified eQTLs for variants directly adjacent to TSS. Future eQTL studies would benefit from combining CAGEseq with RNA sequencing for a more complete interpretation of the transcriptome and increased understanding of eQTL signals.

Preservation of methylated CpG dinucleotides in human CpG islands

BACKGROUND

CpG dinucleotides are extensively underrepresented in mammalian genomes. It is widely accepted that genome-wide CpG depletion is predominantly caused by an elevated CpG > TpG mutation rate due to frequent cytosine methylation in the CpG context. Meanwhile the CpG content in genomic regions called CpG islands (CGIs) is noticeably higher. This observation is usually explained by lower CpG > TpG substitution rates within CGIs due to reduced cytosine methylation levels.

RESULTS

By combining genome-wide data on substitutions and methylation levels in several human cell types we have shown that cytosine methylation in human sperm cells was strongly and consistently associated with increased CpG > TpG substitution rates. In contrast, this correlation was not observed for embryonic stem cells or fibroblasts. Surprisingly, the decreased sperm CpG methylation level was insufficient to explain the reduced CpG > TpG substitution rates in CGIs.

CONCLUSIONS

While cytosine methylation in human sperm cells is strongly associated with increased CpG > TpG substitution rates, substitution rates are significantly reduced within CGIs even after sperm CpG methylation levels and local GC content are controlled for. Our findings are consistent with strong negative selection preserving methylated CpGs within CGIs including intergenic ones.

Gene expression in the aging human brain: an overview

PURPOSE OF REVIEW

The review aims to provide a summary of recent developments in the study of gene expression in the aging human brain.

RECENT FINDINGS

Profiling differentially expressed genes or 'transcripts' in the human brain over the course of normal aging has provided valuable insights into the biological pathways that appear activated or suppressed in late life. Genes mediating neuroinflammation and immune system activation in particular, show significant age-related upregulation creating a state of vulnerability to neurodegenerative and neuropsychiatric disease in the aging brain. Cellular ionic dyshomeostasis and age-related decline in a host of molecular influences on synaptic efficacy may underlie neurocognitive decline in later life. Critically, these investigations have also shed light on the mobilization of protective genetic responses within the aging human brain that help determine health and disease trajectories in older age. There is growing interest in the study of pre and posttranscriptional regulators of gene expression, and the role of noncoding RNAs in particular, as mediators of the phenotypic diversity that characterizes human brain aging.

SUMMARY

Gene expression studies in healthy brain aging offer an opportunity to unravel the intricately regulated cellular underpinnings of neurocognitive aging as well as disease risk and resiliency in late life. In doing so, new avenues for early intervention in age-related neurodegenerative disease could be investigated with potentially significant implications for the development of disease-modifying therapies.

HOCOMOCO: expansion and enhancement of the collection of transcription factor binding sites models

Models of transcription factor (TF) binding sites provide a basis for a wide spectrum of studies in regulatory genomics, from reconstruction of regulatory networks to functional annotation of transcripts and sequence variants. While TFs may recognize different sequence patterns in different conditions, it is pragmatic to have a single generic model for each particular TF as a baseline for practical applications. Here we present the expanded and enhanced version of HOCOMOCO (http://hocomoco.autosome.ru and http://www.cbrc.kaust.edu.sa/hocomoco10), the collection of models of DNA patterns, recognized by transcription factors. HOCOMOCO now provides position weight matrix (PWM) models for binding sites of 601 human TFs and, in addition, PWMs for 396 mouse TFs. Furthermore, we introduce the largest up to date collection of dinucleotide PWM models for 86 (52) human (mouse) TFs. The update is based on the analysis of massive ChIP-Seq and HT-SELEX datasets, with the validation of the resulting models on in vivo data. To facilitate a practical application, all HOCOMOCO models are linked to gene and protein databases (Entrez Gene, HGNC, UniProt) and accompanied by precomputed score thresholds. Finally, we provide command-line tools for PWM and diPWM threshold estimation and motif finding in nucleotide sequences.

EpiFactors: a comprehensive database of human epigenetic factors and complexes

Epigenetics refers to stable and long-term alterations of cellular traits that are not caused by changes in the DNA sequence per se. Rather, covalent modifications of DNA and histones affect gene expression and genome stability via proteins that recognize and act upon such modifications. Many enzymes that catalyse epigenetic modifications or are critical for enzymatic complexes have been discovered, and this is encouraging investigators to study the role of these proteins in diverse normal and pathological processes. Rapidly growing knowledge in the area has resulted in the need for a resource that compiles, organizes and presents curated information to the researchers in an easily accessible and user-friendly form. Here we present EpiFactors, a manually curated database providing information about epigenetic regulators, their complexes, targets and products. EpiFactors contains information on 815 proteins, including 95 histones and protamines. For 789 of these genes, we include expressions values across several samples, in particular a collection of 458 human primary cell samples (for approximately 200 cell types, in many cases from three individual donors), covering most mammalian cell steady states, 255 different cancer cell lines (representing approximately 150 cancer subtypes) and 134 human postmortem tissues. Expression values were obtained by the FANTOM5 consortium using Cap Analysis of Gene Expression technique. EpiFactors also contains information on 69 protein complexes that are involved in epigenetic regulation. The resource is practical for a wide range of users, including biologists, pharmacologists and clinicians.

Database URL: http://epifactors.autosome.ru

Insights into the Transcriptional Architecture of Behavioral Plasticity in the Honey Bee Apis mellifera

Honey bee colonies exhibit an age-related division of labor, with worker bees performing discrete sets of behaviors throughout their lifespan. These behavioral states are associated with distinct brain transcriptomic states, yet little is known about the regulatory mechanisms governing them. We used CAGEscan (a variant of the Cap Analysis of Gene Expression technique) for the first time to characterize the promoter regions of differentially expressed brain genes during two behavioral states (brood care (aka “nursing”) and foraging) and identified transcription factors (TFs) that may govern their expression. More than half of the differentially expressed TFs were associated with motifs enriched in the promoter regions of differentially expressed genes (DEGs), suggesting they are regulators of behavioral state. Strikingly, five TFs (nf-kb, egr, pax6, hairy, and clockwork orange) were predicted to co-regulate nearly half of the genes that were upregulated in foragers. Finally, differences in alternative TSS usage between nurses and foragers were detected upstream of 646 genes, whose functional analysis revealed enrichment for Gene Ontology terms associated with neural function and plasticity. This demonstrates for the first time that alternative TSSs are associated with stable differences in behavior, suggesting they may play a role in organizing behavioral state.

Expression analysis of the long non-coding RNA antisense to Uchl1 (AS Uchl1) during dopaminergic cells' differentiation in vitro and in neurochemical models of Parkinson's disease

Antisense (AS) transcripts are RNA molecules that are transcribed from the opposite strand to sense (S) genes forming S/AS pairs. The most prominent configuration is when a lncRNA is antisense to a protein coding gene. Increasing evidences prove that antisense transcription may control sense gene expression acting at distinct regulatory levels. However, its contribution to brain function and neurodegenerative diseases remains unclear. We have recently identified AS Uchl1 as an antisense to the mouse Ubiquitin carboxy-terminal hydrolase L1 (Uchl1) gene (AS Uchl1), the synthenic locus of UCHL1/PARK5. This is mutated in rare cases of early-onset familial Parkinson's Disease (PD) and loss of UCHL1 activity has been reported in many neurodegenerative diseases. Importantly, manipulation of UchL1 expression has been proposed as tool for therapeutic intervention. AS Uchl1 induces UchL1 expression by increasing its translation. It is the representative member of SINEUPs (SINEB2 sequence to UP-regulate translation), a new functional class of natural antisense lncRNAs that activate translation of their sense genes. Here we take advantage of FANTOM5 dataset to identify the transcription start sites associated to S/AS pair at Uchl1 locus. We show that AS Uchl1 expression is under the regulation of Nurr1, a major transcription factor involved in dopaminergic cells' differentiation and maintenance. Furthermore, AS Uch1 RNA levels are strongly down-regulated in neurochemical models of PD in vitro and in vivo. This work positions AS Uchl1 RNA as a component of Nurr1-dependent gene network and target of cellular stress extending our understanding on the role of antisense transcription in the brain.

Promoter analysis reveals globally differential regulation of human long non-coding RNA and protein-coding genes

Transcriptional regulation of protein-coding genes is increasingly well-understood on a global scale, yet no comparable information exists for long non-coding RNA (lncRNA) genes, which were recently recognized to be as numerous as protein-coding genes in mammalian genomes. We performed a genome-wide comparative analysis of the promoters of human lncRNA and protein-coding genes, finding global differences in specific genetic and epigenetic features relevant to transcriptional regulation. These two groups of genes are hence subject to separate transcriptional regulatory programs, including distinct transcription factor (TF) proteins that significantly favor lncRNA, rather than coding-gene, promoters. We report a specific signature of promoter-proximal transcriptional regulation of lncRNA genes, including several distinct transcription factor binding sites (TFBS). Experimental DNase I hypersensitive site profiles are consistent with active configurations of these lncRNA TFBS sets in diverse human cell types. TFBS ChIP-seq datasets confirm the binding events that we predicted using computational approaches for a subset of factors. For several TFs known to be directly regulated by lncRNAs, we find that their putative TFBSs are enriched at lncRNA promoters, suggesting that the TFs and the lncRNAs may participate in a bidirectional feedback loop regulatory network. Accordingly, cells may be able to modulate lncRNA expression levels independently of mRNA levels via distinct regulatory pathways. Our results also raise the possibility that, given the historical reliance on protein-coding gene catalogs to define the chromatin states of active promoters, a revision of these chromatin signature profiles to incorporate expressed lncRNA genes is warranted in the future.

Effects of cytosine methylation on transcription factor binding sites

BACKGROUND

DNA methylation in promoters is closely linked to downstream gene repression. However, whether DNA methylation is a cause or a consequence of gene repression remains an open question. If it is a cause, then DNA methylation may affect the affinity of transcription factors (TFs) for their binding sites (TFBSs). If it is a consequence, then gene repression caused by chromatin modification may be stabilized by DNA methylation. Until now, these two possibilities have been supported only by non-systematic evidence and they have not been tested on a wide range of TFs. An average promoter methylation is usually used in studies, whereas recent results suggested that methylation of individual cytosines can also be important.

RESULTS

We found that the methylation profiles of 16.6% of cytosines and the expression profiles of neighboring transcriptional start sites (TSSs) were significantly negatively correlated. We called the CpGs corresponding to such cytosines "traffic lights". We observed a strong selection against CpG "traffic lights" within TFBSs. The negative selection was stronger for transcriptional repressors as compared with transcriptional activators or multifunctional TFs as well as for core TFBS positions as compared with flanking TFBS positions.

CONCLUSIONS

Our results indicate that direct and selective methylation of certain TFBS that prevents TF binding is restricted to special cases and cannot be considered as a general regulatory mechanism of transcription.

Measuring cell-type specific differential methylation in human brain tissue

The behavior of epigenetic mechanisms in the brain is obscured by tissue heterogeneity and disease-related histological changes. Not accounting for these confounders leads to biased results. We develop a statistical methodology that estimates and adjusts for celltype composition by decomposing neuronal and non-neuronal differential signal. This method provides a conceptual framework for deconvolving heterogeneous epigenetic data from postmortem brain studies. We apply it to find cell-specific differentially methylated regions between prefrontal cortex and hippocampus. We demonstrate the utility of the method on both Infinium 450k and CHARM data.