Parkinson-Associated SNCA Enhancer Variants Revealed by Open Chromatin in Mouse Dopamine Neurons

An Emerging Role for Enhancer RNAs in Brain Disorders

Noncoding DNA undergoes widespread context-dependent transcription to produce noncoding RNAs. In recent decades, tremendous advances in genomics and transcriptomics have revealed important regulatory roles for noncoding DNA elements and the RNAs that they produce. Enhancers are one such element that are well-established drivers of gene expression changes in response to a variety of factors such as external stimuli, cellular responses, developmental cues, and disease states. They are known to act at long distances, interact with multiple target gene loci simultaneously, synergize with other enhancers, and associate with dynamic chromatin architectures to form a complex regulatory network. Recent advances in enhancer biology have revealed that upon activation, enhancers transcribe long noncoding RNAs, known as enhancer RNAs (eRNAs), that have been shown to play important roles in enhancer-mediated gene regulation and chromatin-modifying activities. In the brain, enhancer dysregulation and eRNA transcription has been reported in numerous disorders from acute injuries to chronic neurodegeneration. Because this is an emerging area, a comprehensive understanding of eRNA function has not yet been achieved in brain disorders; however, the findings to date have illuminated a role for eRNAs in activity-driven gene expression and phenotypic outcomes. In this review, we highlight the breadth of the current literature on eRNA biology in brain health and disease and discuss the challenges as well as focus areas and strategies for future in-depth research on eRNAs in brain health and disease.

Unravelling cell type-specific responses to Parkinson's Disease at single cell resolution

Parkinson's Disease (PD) is the second most common neurodegenerative disorder. The pathological hallmark of PD is loss of dopaminergic neurons and the presence of aggregated α-synuclein, primarily in the substantia nigra pars compacta (SNpc) of the midbrain. However, the molecular mechanisms that underlie the pathology in different cell types is not currently understood. Here, we present a single nucleus transcriptome analysis of human post-mortem SNpc obtained from 15 sporadic Parkinson's Disease (PD) cases and 14 Controls. Our dataset comprises ∼84K nuclei, representing all major cell types of the brain, allowing us to obtain a transcriptome-level characterization of these cell types. Importantly, we identify multiple subpopulations for each cell type and describe specific gene sets that provide insights into the differing roles of these subpopulations. Our findings reveal a significant decrease in neuronal cells in PD samples, accompanied by an increase in glial cells and T cells. Subpopulation analyses demonstrate a significant depletion of tyrosine hydroxylase (TH) enriched astrocyte, microglia and oligodendrocyte populations in PD samples, as well as TH enriched neurons, which are also depleted. Moreover, marker gene analysis of the depleted subpopulations identified 28 overlapping genes, including those associated with dopamine metabolism (e.g., ALDH1A1, SLC6A3 & SLC18A2). Overall, our study provides a valuable resource for understanding the molecular mechanisms involved in dopaminergic neuron degeneration and glial responses in PD, highlighting the existence of novel subpopulations and cell type-specific gene sets.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Activation of HCA2 regulates microglial responses to alleviate neurodegeneration in LPS-induced in vivo and in vitro models

BACKGROUND

Previous studies have shown a close association between an altered immune system and Parkinson's disease (PD). Neuroinflammation inhibition may be an effective measure to prevent PD. Recently, numerous reports have highlighted the potential of hydroxy-carboxylic acid receptor 2 (HCA2) in inflammation-related diseases. Notably, the role of HCA2 in neurodegenerative diseases is also becoming more widely known. However, its role and exact mechanism in PD remain to be investigated. Nicotinic acid (NA) is one of the crucial ligands of HCA2, activating it. Based on such findings, this study aimed to examine the effect of HCA2 on neuroinflammation and the role of NA-activated HCA2 in PD and its underlying mechanisms.

METHODS

For in vivo studies, 10-week-old male C57BL/6 and HCA2^-/- mice were injected with LPS in the substantia nigra (SN) to construct a PD model. The motor behavior of mice was detected using open field, pole-climbing and rotor experiment. The damage to the mice's dopaminergic neurons was detected using immunohistochemical staining and western blotting methods. In vitro, inflammatory mediators (IL-6, TNF-α, iNOS and COX-2) and anti-inflammatory factors (Arg-1, Ym-1, CD206 and IL-10) were detected using RT-PCR, ELISA and immunofluorescence. Inflammatory pathways (AKT, PPARγ and NF-κB) were delineated by RT-PCR and western blotting. Neuronal damage was detected using CCK8, LDH, and flow cytometry assays.

RESULTS

HCA2^-/- increases mice susceptibility to dopaminergic neuronal injury, motor deficits, and inflammatory responses. Mechanistically, HCA2 activation in microglia promotes anti-inflammatory microglia and inhibits pro-inflammatory microglia by activating AKT/PPARγ and inhibiting NF-κB signaling pathways. Further, HCA2 activation in microglia attenuates microglial activation-mediated neuronal injury. Moreover, nicotinic acid (NA), a specific agonist of HCA2, alleviated dopaminergic neuronal injury and motor deficits in PD mice by activating HCA2 in microglia in vivo.

CONCLUSIONS

Niacin receptor HCA2 modulates microglial phenotype to inhibit neurodegeneration in LPS-induced in vivo and in vitro models.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate of disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo , mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable a proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Boosting tissue-specific prediction of active cis-regulatory regions through deep learning and Bayesian optimization techniques

BACKGROUND

Cis-regulatory regions (CRRs) are non-coding regions of the DNA that fine control the spatio-temporal pattern of transcription; they are involved in a wide range of pivotal processes such as the development of specific cell-lines/tissues and the dynamic cell response to physiological stimuli. Recent studies showed that genetic variants occurring in CRRs are strongly correlated with pathogenicity or deleteriousness. Considering the central role of CRRs in the regulation of physiological and pathological conditions, the correct identification of CRRs and of their tissue-specific activity status through Machine Learning methods plays a major role in dissecting the impact of genetic variants on human diseases. Unfortunately, the problem is still open, though some promising results have been already reported by (deep) machine-learning based methods that predict active promoters and enhancers in specific tissues or cell lines by encoding epigenetic or spectral features directly extracted from DNA sequences.

RESULTS

We present the experiments we performed to compare two Deep Neural Networks, a Feed-Forward Neural Network model working on epigenomic features, and a Convolutional Neural Network model working only on genomic sequence, targeted to the identification of enhancer- and promoter-activity in specific cell lines. While performing experiments to understand how the experimental setup influences the prediction performance of the methods, we particularly focused on (1) automatic model selection performed by Bayesian optimization and (2) exploring different data rebalancing setups for reducing negative unbalancing effects.

CONCLUSIONS

Results show that (1) automatic model selection by Bayesian optimization improves the quality of the learner; (2) data rebalancing considerably impacts the prediction performance of the models; test set rebalancing may provide over-optimistic results, and should therefore be cautiously applied; (3) despite working on sequence data, convolutional models obtain performance close to those of feed forward models working on epigenomic information, which suggests that also sequence data carries informative content for CRR-activity prediction. We therefore suggest combining both models/data types in future works.

Intronic enhancers of the human SNCA gene predominantly regulate its expression in brain in vivo

Evidence from patients with Parkinson's disease (PD) and our previously reported α-synuclein (SNCA) transgenic rat model support the idea that increased SNCA protein is a substantial risk factor of PD pathogenesis. However, little is known about the transcription control of the human SNCA gene in the brain in vivo. Here, we identified that the DYT6 gene product THAP1 (THAP domain-containing apoptosis-associated protein 1) and its interaction partner CTCF (CCCTC-binding factor) act as transcription regulators of SNCA. THAP1 controls SNCA intronic enhancers' activities, while CTCF regulates its enhancer-promoter loop formation. The SNCA intronic enhancers present neurodevelopment-dependent activities and form enhancer clusters similar to "super-enhancers" in the brain, in which the PD-associated single-nucleotide polymorphisms are enriched. Deletion of the SNCA intronic enhancer clusters prevents the release of paused RNA polymerase II from its promoter and subsequently reduces its expression drastically in the brain, which may provide new therapeutic approaches to prevent its accumulation and thus related neurodegenerative diseases defined as synucleinopathies.

Neuroinflammation represents a common theme amongst genetic and environmental risk factors for Alzheimer and Parkinson diseases

Multifactorial diseases are characterized by inter-individual variation in etiology, age of onset, and penetrance. These diseases tend to be relatively common and arise from the combined action of genetic and environmental factors; however, parsing the convoluted mechanisms underlying these gene-by-environment interactions presents a significant challenge to their study and management. For neurodegenerative disorders, resolving this challenge is imperative, given the enormous health and societal burdens they impose. The mechanisms by which genetic and environmental effects may act in concert to destabilize homeostasis and elevate risk has become a major research focus in the study of common disease. Emphasis is further being placed on determining the extent to which a unifying biological principle may account for the progressively diminishing capacity of a system to buffer disease phenotypes, as risk for disease increases. Data emerging from studies of common, neurodegenerative diseases are providing insights to pragmatically connect mechanisms of genetic and environmental risk that previously seemed disparate. In this review, we discuss evidence positing inflammation as a unifying biological principle of homeostatic destabilization affecting the risk, onset, and progression of neurodegenerative diseases. Specifically, we discuss how genetic variation associated with Alzheimer disease and Parkinson disease may contribute to pro-inflammatory responses, how such underlying predisposition may be exacerbated by environmental insults, and how this common theme is being leveraged in the ongoing search for effective therapeutic interventions.

Comparative chromatin accessibility upon BDNF stimulation delineates neuronal regulatory elements

Neuronal stimulation induced by the brain-derived neurotrophic factor (BDNF) triggers gene expression, which is crucial for neuronal survival, differentiation, synaptic plasticity, memory formation, and neurocognitive health. However, its role in chromatin regulation is unclear. Here, using temporal profiling of chromatin accessibility and transcription in mouse primary cortical neurons upon either BDNF stimulation or depolarization (KCl), we identify features that define BDNF-specific chromatin-to-gene expression programs. Enhancer activation is an early event in the regulatory control of BDNF-treated neurons, where the bZIP motif-binding Fos protein pioneered chromatin opening and cooperated with co-regulatory transcription factors (Homeobox, EGRs, and CTCF) to induce transcription. Deleting cis-regulatory sequences affect BDNF-mediated Arc expression, a regulator of synaptic plasticity. BDNF-induced accessible regions are linked to preferential exon usage by neurodevelopmental disorder-related genes and the heritability of neuronal complex traits, which were validated in human iPSC-derived neurons. Thus, we provide a comprehensive view of BDNF-mediated genome regulatory features using comparative genomic approaches to dissect mammalian neuronal stimulation.

Ranking reprogramming factors for cell differentiation

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

The role of Tumour Necrosis Factor in neuroinflammation associated with Parkinson's disease and targeted therapies

Neurodegenerative disorders Parkinson's disease is a progressive neurodegenerative disorder associated with neuroinflammatory responses that lead to the neurodegeneration of the dopaminergic neurons. These neuroinflammatory mechanisms involve various cytokines produced by the activated glial cells. Tumour Necrosis factor α (TNF α) is one of the major mediators of the neuroinflammation associated with neurodegeneration. TNF α has a dual role of neuroprotection and neurotoxicity in the brain. The effective pathways of TNF involve various signalling pathways transduced by the receptors TNFR1 and TNFR2. Effective therapeutic strategies have been produced targeting the neurotoxic behaviour of the Tumour Necrosis Factor and the associated neurodegeneration which includes the use of Dominant Negative Tumour Necrosis Factor (DN-TNF) inhibitors like XENP 345 and XPro®1595 and peroxisome proliferator receptor gamma (PPAR-γ) agonists.

Single-cell analysis of gene expression in the substantia nigra pars compacta of a pesticide-induced mouse model of Parkinson's disease

Exposure to pesticides in humans increases the risk of Parkinson’s disease (PD), but the mechanisms remain poorly understood. To elucidate these pathways, we dosed C57BL/6J mice with a combination of the pesticides maneb and paraquat. Behavioral analysis revealed motor deficits consistent with PD. Single-cell RNA sequencing of substantia nigra pars compacta revealed both cell-type-specific genes and genes expressed differentially between pesticide and control, including Fam241b, Emx2os, Bivm, Gm1439, Prdm15, and Rai2. Neurons had the largest number of significant differentially expressed genes, but comparable numbers were found in astrocytes and less so in oligodendrocytes. In addition, network analysis revealed enrichment in functions related to the extracellular matrix. These findings emphasize the importance of support cells in pesticide-induced PD and refocus our attention away from neurons as the sole agent of this disorder.

Chromatin Landscapes of Human Lung Cells Predict Potentially Functional Chronic Obstructive Pulmonary Disease Genome-Wide Association Study Variants

Genome-wide association studies (GWASs) have identified dozens of loci associated with risk of chronic obstructive pulmonary disease (COPD). However, identifying the causal variants and their functional role in the appropriate cell type remains a major challenge. We aimed to identify putative causal variants in 82 GWAS loci associated with COPD susceptibility and predict the regulatory impact of these variants in lung-cell types. We used an integrated approach featuring statistical fine mapping, open chromatin profiling, and machine learning to identify functional variants. We generated chromatin accessibility data using the Assay for Transposase-Accessible Chromatin with High-Throughput Sequencing (ATAC-seq) for human primary lung-cell types implicated in COPD pathobiology. We then evaluated the enrichment of COPD risk variants in lung-specific open chromatin regions and generated cell type-specific regulatory predictions for >6,500 variants corresponding to 82 COPD GWAS loci. Integration of the fine-mapped variants with lung open chromatin regions helped prioritize 22 variants in putative regulatory elements with potential functional effects. Comparison with functional predictions from 222 Encyclopedia of DNA Elements (ENCODE) cell samples revealed cell type-specific regulatory effects of COPD variants in the lung epithelium, endothelium, and immune cells. We identified potential causal variants for COPD risk by integrating fine mapping in GWAS loci with cell-specific regulatory profiling, highlighting the importance of leveraging the chromatin status in relevant cell types to predict the molecular effects of risk variants in lung disease.

The landscape of chromatin accessibility in skeletal muscle during embryonic development in pigs

BACKGROUND

The development of skeletal muscle in pigs during the embryonic stage is precisely regulated by transcriptional mechanisms, which depend on chromatin accessibility. However, how chromatin accessibility plays a regulatory role during embryonic skeletal muscle development in pigs has not been reported. To gain insight into the landscape of chromatin accessibility and the associated genome-wide transcriptome during embryonic muscle development, we performed ATAC-seq and RNA-seq analyses of skeletal muscle from pig embryos at 45, 70 and 100 days post coitus (dpc).

RESULTS

In total, 21,638, 35,447 and 60,181 unique regions (or peaks) were found across the embryos at 45 dpc (LW45), 70 dpc (LW70) and 100 dpc (LW100), respectively. More than 91% of the peaks were annotated within - 1 kb to 100 bp of transcription start sites (TSSs). First, widespread increases in specific accessible chromatin regions (ACRs) from embryos at 45 to 100 dpc suggested that the regulatory mechanisms became increasingly complicated during embryonic development. Second, the findings from integrated ATAC-seq and RNA-seq analyses showed that not only the numbers but also the intensities of ACRs could control the expression of associated genes. Moreover, the motif screening of stage-specific ACRs revealed some transcription factors that regulate muscle development-related genes, such as MyoG, Mef2c, and Mef2d. Several potential transcriptional repressors, including E2F6, OTX2 and CTCF, were identified among the genes that exhibited different regulation trends between the ATAC-seq and RNA-seq data.

CONCLUSIONS

This work indicates that chromatin accessibility plays an important regulatory role in the embryonic muscle development of pigs and regulates the temporal and spatial expression patterns of key genes in muscle development by influencing the binding of transcription factors. Our results contribute to a better understanding of the regulatory dynamics of genes involved in pig embryonic skeletal muscle development.

Alpha synuclein (SNCA) rs7684318 variant contributes to Parkinson's disease risk by altering transcription factor binding related with Notch and Wnt signaling

In view of inconsistencies in the association studies of alpha synuclein (SNCA) rs7684318 (chr4: 90655003 T > C) with Parkinson's disease (PD), we conducted a meta-analysis to establish the association of this variant with PD and examined changes in transcription factor binding. SNCA rs7684318 C-allele was identified as genetic risk factor for PD in fixed (OR: 1.53, 95 % CI: 1.40-1.68, p < 0.0001) and random effect (OR: 1.65, 95 % CI: 1.30-2.09, p = 0.0003) models. Heterogeneity was observed in association (Tau²: 0.0576, H: 2.32, I²: 0.815, Q: 21.64, p = 0.0002). Egger's test showed no evidence of publication bias (p = 0.37). Subgroup analysis showed that rs7684318 is contributing to PD risk in Japanese (OR: 1.46, 95 % CI: 1.30-1.64) and Indian (OR: 2.63, 95 % CI: 1.79-3.86) populations while showing no significant association in Chinese population (OR: 1.68, 95 % CI: 0.93-3.02). Sensitivity analysis showed that exclusion of any one of the studies has no significant impact on the association, which justifies the robustness of the analysis. Tissue-specific DNase foot print analysis revealed that this variant contributes to increased transcription factor binding in midbrain, putamen and caudate nucleus. The substitution of T > C increased binding of RBPJ and GATA-family transcription factors; and decreased binding of NKX2 family, SNAI2, SNAI3, DMRT1, HOXA13, HOXB13, HOXC13, HOXD13, WT1, POU4F1, POU4F2, POU4F3 transcriptional factors. TRANSFAC and DNA curvature analyses substantiate the association of this variant with increased binding of GATA1 that contribute to intensity of DNA curvature peaks and splitting pattern. These studies along with the meta-analysis strongly suggest that the rs7684318 variant contributes to the pathophysiology of PD by modulating binding of transcription factors related to Notch and Wnt signalling pathways that are likely to impair dopmanergic transmission.

Massively Parallel Reporter Assays: Defining Functional Psychiatric Genetic Variants Across Biological Contexts

Neuropsychiatric phenotypes have long been known to be influenced by heritable risk factors, directly confirmed by the past decade of genetic studies that have revealed specific genetic variants enriched in disease cohorts. However, the initial hope that a small set of genes would be responsible for a given disorder proved false. The more complex reality is that a given disorder may be influenced by myriad small-effect noncoding variants and/or by rare but severe coding variants, many de novo. Noncoding genomic sequences-for which molecular functions cannot usually be inferred-harbor a large portion of these variants, creating a substantial barrier to understanding higher-order molecular and biological systems of disease. Fortunately, novel genetic technologies-scalable oligonucleotide synthesis, RNA sequencing, and CRISPR (clustered regularly interspaced short palindromic repeats)-have opened novel avenues to experimentally identify biologically significant variants en masse. Massively parallel reporter assays (MPRAs) are an especially versatile technique resulting from such innovations. MPRAs are powerful molecular genetics tools that can be used to screen thousands of untranscribed or untranslated sequences and their variants for functional effects in a single experiment. This approach, though underutilized in psychiatric genetics, has several useful features for the field. We review methods for assaying putatively functional genetic variants and regions, emphasizing MPRAs and the opportunities they hold for dissection of psychiatric polygenicity. We discuss literature applying functional assays in neurogenetics, highlighting strengths, caveats, and design considerations-especially regarding disease-relevant variables (cell type, neurodevelopment, and sex), and we ultimately propose applications of MPRA to both computational and experimental neurogenetics of polygenic disease risk.

The cis-Regulatory Element of SNCA Intron 4 Modulates Susceptibility to Parkinson's Disease in Han Chinese

Objective: A novel functional cis-regulatory element (CRE) located at SNCA intron 4 has recently been identified in association with Parkinson's disease (PD) risk in European descendants. We aimed to investigate whether this CRE is associated with PD in Han Chinese ethnicity. Methods: A Chinese cohort comprising 513 sporadic PD patients and 517 controls was recruited. CRE variants were identified by sequencing and then analyzed. Results: A total of nine variants were detected, namely eight single nucleotide variants and one new insertion variant. Two variants, rs17016188 and rs7684892, had minor allele frequency greater than 5%. A difference of rs17016188 was observed in males with the C allele serving as a recessive risk factor (p = 0.001, OR = 2.349, 95% CI = 1.414-3.901) following Bonferroni correction. Haplotypes of rs17016188 and rs7684892 showed distribution differences in the total and the male populations (p = 0.002 and 4.08 × 10^-5, respectively). Among the haplotypes, rs17016188/T-rs7684892/G was associated with a reduced risk for PD (p = 4.8 × 10^-4, OR = 0.731, 95% CI = 0.614-0.872). Conclusions: Our results provide insight into how the SNCA intron 4 CRE harbors variants and its contribution to PD risk in Chinese ethnicity.

Enhancer Predictions and Genome-Wide Regulatory Circuits

Spatiotemporal control of gene expression during development requires orchestrated activities of numerous enhancers, which are cis-regulatory DNA sequences that, when bound by transcription factors, support selective activation or repression of associated genes. Proper activation of enhancers is critical during embryonic development, adult tissue homeostasis, and regeneration, and inappropriate enhancer activity is often associated with pathological conditions such as cancer. Multiple consortia [e.g., the Encyclopedia of DNA Elements (ENCODE) Consortium and National Institutes of Health Roadmap Epigenomics Mapping Consortium] and independent investigators have mapped putative regulatory regions in a large number of cell types and tissues, but the sequence determinants of cell-specific enhancers are not yet fully understood. Machine learning approaches trained on large sets of these regulatory regions can identify core transcription factor binding sites and generate quantitative predictions of enhancer activity and the impact of sequence variants on activity. Here, we review these computational methods in the context of enhancer prediction and gene regulatory network models specifying cell fate.

Leveraging mouse chromatin data for heritability enrichment informs common disease architecture and reveals cortical layer contributions to schizophrenia

Genome-wide association studies have implicated thousands of noncoding variants across common human phenotypes. However, they cannot directly inform the cellular context in which disease-associated variants act. Here, we use open chromatin profiles from discrete mouse cell populations to address this challenge. We applied stratified linkage disequilibrium score regression and evaluated heritability enrichment in 64 genome-wide association studies, emphasizing schizophrenia. We provide evidence that mouse-derived human open chromatin profiles can serve as powerful proxies for difficult to obtain human cell populations, facilitating the illumination of common disease heritability enrichment across an array of human phenotypes. We demonstrate that signatures from discrete subpopulations of cortical excitatory and inhibitory neurons are significantly enriched for schizophrenia heritability with maximal enrichment in cortical layer V excitatory neurons. We also show that differences between schizophrenia and bipolar disorder are concentrated in excitatory neurons in cortical layers II-III, IV, and V, as well as the dentate gyrus. Finally, we leverage these data to fine-map variants in 177 schizophrenia loci nominating variants in 104/177. We integrate these data with transcription factor binding site, chromatin interaction, and validated enhancer data, placing variants in the cellular context where they may modulate risk.

Missing heritability in Parkinson's disease: the emerging role of non-coding genetic variation

Parkinson’s disease (PD) is a neurodegenerative disorder caused by a complex interplay of genetic and environmental factors. For the stratification of PD patients and the development of advanced clinical trials, including causative treatments, a better understanding of the underlying genetic architecture of PD is required. Despite substantial efforts, genome-wide association studies have not been able to explain most of the observed heritability. The majority of PD-associated genetic variants are located in non-coding regions of the genome. A systematic assessment of their functional role is hampered by our incomplete understanding of genotype–phenotype correlations, for example through differential regulation of gene expression. Here, the recent progress and remaining challenges for the elucidation of the role of non-coding genetic variants is reviewed with a focus on PD as a complex disease with multifactorial origins. The function of gene regulatory elements and the impact of non-coding variants on them, and the means to map these elements on a genome-wide level, will be delineated. Moreover, examples of how the integration of functional genomic annotations can serve to identify disease-associated pathways and to prioritize disease- and cell type-specific regulatory variants will be given. Finally, strategies for functional validation and considerations for suitable model systems are outlined. Together this emphasizes the contribution of rare and common genetic variants to the complex pathogenesis of PD and points to remaining challenges for the dissection of genetic complexity that may allow for better stratification, improved diagnostics and more targeted treatments for PD in the future.

Characterization of the chromatin accessibility in an Alzheimer's disease (AD) mouse model

Background

The pathological hallmarks of Alzheimer’s disease (AD) involve alterations in the expression of numerous genes associated with transcriptional levels, which are determined by chromatin accessibility. Here, the landscape of chromatin accessibility was studied to understand the outline of the transcription and expression of AD-associated metabolism genes in an AD mouse model.

Methods

The assay for transposase-accessible chromatin by sequencing (ATAC-seq) was used to investigate the AD-associated chromatin reshaping in the APPswe/PS1dE9 (APP/PS1) mouse model. ATAC-seq data in the hippocampus of 8-month-old APP/PS1 mice were generated, and the relationship between chromatin accessibility and gene expression was analyzed in combination with RNA sequencing. Gene ontology (GO) analysis was applied to elucidate biological processes and signaling pathways altered in APP/PS1 mice. Critical transcription factors were identified; alterations in chromatin accessibility were further confirmed using chromatin immunoprecipitation assays.

Results

We identified 1690 increased AD-associated chromatin-accessible regions in the hippocampal tissues of APP/PS1 mice. These regions were enriched in genes related to diverse signaling pathways, including the PI3K-Akt, Hippo, TGF-β, and Jak-Stat signaling pathways, which play essential roles in regulating cell proliferation, apoptosis, and inflammatory responses. A total of 1003 decreased chromatin-accessible regions were considered to be related with declined AD-associated biological processes including cellular response to hyperoxia and insulin stimulus, synaptic transmission, and positive regulation of autophagy. In the APP/PS1 hippocampus, 1090 genes were found to be upregulated and 1081 downregulated. Interestingly, enhanced ATAC-seq signal was found in approximately 740 genes, with 43 exhibiting upregulated mRNA levels. Several genes involved in AD development were found to have a significantly increased expression in APP/PS1 mice compared to controls, including Sele, Clec7a, Cst7, and Ccr6. The signatures of numerous transcription factors, including Olig2, NeuroD1, TCF4, and NeuroG2, were found enriched in the AD-associated accessible chromatin regions. The transcription-activating marks of H3K4me3 and H3K27ac were also found increased in the promoters of these genes. These results indicate that the mechanism for the upregulation of genes could be attributed to the enrichment of open chromatin regions with transcription factors motifs and the histone marks H3K4me3 and H3K27ac.

Conclusion

Our study reveals that alterations in chromatin accessibility may be an initial mechanism in AD pathogenesis.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s13195-020-00598-2.