TDP-43 Epigenetic Facets and Their Neurodegenerative Implications.

Since its initial involvement in numerous neurodegenerative pathologies in 2006, either as a principal actor or as a cofactor, new pathologies implicating transactive response (TAR) DNA-binding protein 43 (TDP-43) are regularly emerging also beyond the neuronal system. This reflects the fact that TDP-43 functions are particularly complex and broad in a great variety of human cells. In neurodegenerative diseases, this protein is often pathologically delocalized to the cytoplasm, where it irreversibly aggregates and is subjected to various post-translational modifications such as phosphorylation, polyubiquitination, and cleavage. Until a few years ago, the research emphasis has been focused particularly on the impacts of this aggregation and/or on its widely described role in complex RNA splicing, whether related to loss- or gain-of-function mechanisms. Interestingly, recent studies have strengthened the knowledge of TDP-43 activity at the chromatin level and its implication in the regulation of DNA transcription and stability. These discoveries have highlighted new features regarding its own transcriptional regulation and suggested additional mechanistic and disease models for the effects of TPD-43. In this review, we aim to give a comprehensive view of the potential epigenetic (de)regulations driven by (and driving) this multitask DNA/RNA-binding protein.

Synaptic and transcriptomic features of cortical and amygdala pyramidal neurons predict inefficient fear extinction.

Fear-related disorders arise from inefficient fear extinction and have immeasurable social and economic costs. Here, we characterize mouse phenotypes that spontaneously show fear-independent behavioral traits predicting adaptive or maladaptive fear extinction. We find that, already before fear conditioning, specific morphological, electrophysiological, and transcriptomic patterns of cortical and amygdala pyramidal neurons predispose to fear-related disorders. Finally, by using an optogenetic approach, we show the possibility to rescue inefficient fear extinction by activating infralimbic pyramidal neurons and to impair fear extinction by activating prelimbic pyramidal neurons.

The Ying and Yang of Hydrogen Sulfide as a Paracrine/Autocrine Agent in Neurodegeneration: Focus on Amyotrophic Lateral Sclerosis.

Ever since its presence was reported in the brain, the nature and role of hydrogen sulfide (H2S) in the Central Nervous System (CNS) have changed. Consequently, H2S has been elected as the third gas transmitter, along with carbon monoxide and nitric oxide, and a number of studies have focused on its neuromodulatory and protectant functions in physiological conditions. The research on H2S has highlighted its many facets in the periphery and in the CNS, and its role as a double-faced compound, switching from protective to toxic depending on its concentration. In this review, we will focus on the bell-shaped nature of H2S as an angiogenic factor and as a molecule released by glial cells (mainly astrocytes) and non-neuronal cells acting on the surrounding environment (paracrine) or on the releasing cells themselves (autocrine). Finally, we will discuss its role in Amyotrophic Lateral Sclerosis, a paradigm of a neurodegenerative disease.

D4Z4 Methylation Levels Combined with a Machine Learning Pipeline Highlight Single CpG Sites as Discriminating Biomarkers for FSHD Patients.

The study describes a protocol for methylation analysis integrated with Machine Learning (ML) algorithms developed to classify Facio-Scapulo-Humeral Dystrophy (FSHD) subjects. The DNA methylation levels of two D4Z4 regions (DR1 and DUX4-PAS) were assessed by an in-house protocol based on bisulfite sequencing and capillary electrophoresis, followed by statistical and ML analyses. The study involved two independent cohorts, namely a training group of 133 patients with clinical signs of FSHD and 150 healthy controls (CTRL) and a testing set of 27 FSHD patients and 25 CTRL. As expected, FSHD patients showed significantly reduced methylation levels compared to CTRL. We utilized single CpG sites to develop a ML pipeline able to discriminate FSHD subjects. The model identified four CpGs sites as the most relevant for the discrimination of FSHD subjects and showed high metrics values (accuracy: 0.94, sensitivity: 0.93, specificity: 0.96). Two additional models were developed to differentiate patients with lower D4Z4 size and patients who might carry pathogenic variants in FSHD genes, respectively. Overall, the present model enables an accurate classification of FSHD patients, providing additional evidence for DNA methylation as a powerful disease biomarker that could be employed for prioritizing subjects to be tested for FSHD.

Transcriptomic and Network Analyses Reveal Immune Modulation by Endocannabinoids in Approach/Avoidance Traits.

Approach and avoidance (A/A) tendencies are stable behavioral traits in responding to rewarding and fearful stimuli. They represent the superordinate division of emotion, and individual differences in such traits are associated with disease susceptibility. The neural circuitry underlying A/A traits is retained to be the cortico-limbic pathway including the amygdala, the central hub for the emotional processing. Furthermore, A/A-specific individual differences are associated with the activity of the endocannabinoid system (ECS) and especially of CB1 receptors whose density and functionality in amygdala differ according to A/A traits. ECS markedly interacts with the immune system (IS). However, how the interplay between ECS and IS is associated with A/A individual differences is still ill-defined. To fill this gap, here we analyzed the interaction between the gene expression of ECS and immune system (IS) in relation to individual differences. To unveil the deep architecture of ECS-IS interaction, we performed cell-specific transcriptomics analysis. Differential gene expression profiling, functional enrichment, and protein-protein interaction network analyses were performed in amygdala pyramidal neurons of mice showing different A/A behavioral tendencies. Several altered pro-inflammatory pathways were identified as associated with individual differences in A/A traits, indicating the chronic activation of the adaptive immune response sustained by the interplay between endocannabinoids and the IS. Furthermore, results showed that the interaction between the two systems modulates synaptic plasticity and neuronal metabolism in individual difference-specific manner. Deepening our knowledge about ECS/IS interaction may provide useful targets for treatment and prevention of psychopathology associated with A/A traits.

Malat-1-PRC2-EZH1 interaction supports adaptive oxidative stress dependent epigenome remodeling in skeletal myotubes.

PRC2-mediated epigenetic function involves the interaction with long non-coding RNAs (lncRNAs). Although the identity of some of these RNAs has been elucidated in the context of developmental programs, their counterparts in postmitotic adult tissue homeostasis remain uncharacterized. To this aim, we used terminally differentiated postmitotic skeletal muscle cells in which oxidative stress induces the dynamic activation of PRC2-Ezh1 through Embryonic Ectoderm Develpment (EED) shuttling to the nucleus. We identify lncRNA Malat-1 as a necessary partner for PRC2-Ezh1-dependent response to oxidative stress. We show that in this pathway, PRC2-EZH1 dynamic assembly, and in turn stress induced skeletal muscle targeted genes repression, depends specifically on Malat-1. Our study reports about PRC2-RNA interactions in the physiological context of adaptive oxidative stress response and identifies the first lncRNA involved in PRC2-Ezh1 function.

Optogenetic Stimulation of Prelimbic Pyramidal Neurons Maintains Fear Memories and Modulates Amygdala Pyramidal Neuron Transcriptome.

Fear extinction requires coordinated neural activity within the amygdala and medial prefrontal cortex (mPFC). Any behavior has a transcriptomic signature that is modified by environmental experiences, and specific genes are involved in functional plasticity and synaptic wiring during fear extinction. Here, we investigated the effects of optogenetic manipulations of prelimbic (PrL) pyramidal neurons and amygdala gene expression to analyze the specific transcriptional pathways associated to adaptive and maladaptive fear extinction. To this aim, transgenic mice were (or not) fear-conditioned and during the extinction phase they received optogenetic (or sham) stimulations over photo-activable PrL pyramidal neurons. At the end of behavioral testing, electrophysiological (neural cellular excitability and Excitatory Post-Synaptic Currents) and morphological (spinogenesis) correlates were evaluated in the PrL pyramidal neurons. Furthermore, transcriptomic cell-specific RNA-analyses (differential gene expression profiling and functional enrichment analyses) were performed in amygdala pyramidal neurons. Our results show that the optogenetic activation of PrL pyramidal neurons in fear-conditioned mice induces fear extinction deficits, reflected in an increase of cellular excitability, excitatory neurotransmission, and spinogenesis of PrL pyramidal neurons, and associated to strong modifications of the transcriptome of amygdala pyramidal neurons. Understanding the electrophysiological, morphological, and transcriptomic architecture of fear extinction may facilitate the comprehension of fear-related disorders.

Functional annotation of human long noncoding RNAs via molecular phenotyping.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

RADICL-seq identifies general and cell type-specific principles of genome-wide RNA-chromatin interactions.

Mammalian genomes encode tens of thousands of noncoding RNAs. Most noncoding transcripts exhibit nuclear localization and several have been shown to play a role in the regulation of gene expression and chromatin remodeling. To investigate the function of such RNAs, methods to massively map the genomic interacting sites of multiple transcripts have been developed; however, these methods have some limitations. Here, we introduce RNA And DNA Interacting Complexes Ligated and sequenced (RADICL-seq), a technology that maps genome-wide RNA-chromatin interactions in intact nuclei. RADICL-seq is a proximity ligation-based methodology that reduces the bias for nascent transcription, while increasing genomic coverage and unique mapping rate efficiency compared with existing methods. RADICL-seq identifies distinct patterns of genome occupancy for different classes of transcripts as well as cell type-specific RNA-chromatin interactions, and highlights the role of transcription in the establishment of chromatin structure.

Microarray-based sketches of the HERV transcriptome landscape.

Human endogenous retroviruses (HERVs) are spread throughout the genome and their long terminal repeats (LTRs) constitute a wide collection of putative regulatory sequences. Phylogenetic similarities and the profusion of integration sites, two inherent characteristics of transposable elements, make it difficult to study individual locus expression in a large-scale approach, and historically apart from some placental and testis-regulated elements, it was generally accepted that HERVs are silent due to epigenetic control. Herein, we have introduced a generic method aiming to optimally characterize individual loci associated with 25-mer probes by minimizing cross-hybridization risks. We therefore set up a microarray dedicated to a collection of 5,573 HERVs that can reasonably be assigned to a unique genomic position. We obtained a first view of the HERV transcriptome by using a composite panel of 40 normal and 39 tumor samples. The experiment showed that almost one third of the HERV repertoire is indeed transcribed. The HERV transcriptome follows tropism rules, is sensitive to the state of differentiation and, unexpectedly, seems not to correlate with the age of the HERV families. The probeset definition within the U3 and U5 regions was used to assign a function to some LTRs (i.e. promoter or polyA) and revealed that (i) autonomous active LTRs are broadly subjected to operational determinism (ii) the cellular gene density is substantially higher in the surrounding environment of active LTRs compared to silent LTRs and (iii) the configuration of neighboring cellular genes differs between active and silent LTRs, showing an approximately 8 kb zone upstream of promoter LTRs characterized by a drastic reduction in sense cellular genes. These gathered observations are discussed in terms of virus/host adaptive strategies, and together with the methods and tools developed for this purpose, this work paves the way for further HERV transcriptome projects.

Custom human endogenous retroviruses dedicated microarray identifies self-induced HERV-W family elements reactivated in testicular cancer upon methylation control.

Endogenous retroviruses (ERVs) are an inherited part of the eukaryotic genomes, and represent approximately 400,000 loci in the human genome. Human endogenous retroviruses (HERVs) can be divided into distinct families, composed of phylogenetically related but structurally heterogeneous elements. The majority of HERVs are silent in most physiological contexts, whereas a significant expression is observed in pathological contexts, such as cancers. Owing to their repetitive nature, few of the active HERV elements have been accurately identified. In addition, there are no criteria defining the active promoters among HERV long-terminal repeats (LTRs). Hence, it is difficult to understand the HERV (de)regulation mechanisms and their implication on the physiopathology of the host. We developed a microarray to specifically detect the LTR-containing transcripts from the HERV-H, HERV-E, HERV-W and HERV-K(HML-2) families. HERV transcriptome was analyzed in the placenta and seven normal/tumoral match-pair samples. We identified six HERV-W loci overexpressed in testicular cancer, including a usually placenta-restricted transcript of ERVWE1. For each locus, specific overexpression was confirmed by quantitative RT-PCR, and comparison of the activity of U3 versus U5 regions suggested a U3-promoted transcription coupled with 5'R initiation. The analysis of DNA from tumoral versus normal tissue revealed that hypomethylation of U3 promoters in tumors is a prerequisite for their activation.

Comparative methylation of ERVWE1/syncytin-1 and other human endogenous retrovirus LTRs in placenta tissues.

Human endogenous retroviruses (HERVs) are globally silent in somatic cells. However, some HERVs display high transcription in physiological conditions. In particular, ERVWE1, ERVFRDE1 and ERV3, three proviruses of distinct families, are highly transcribed in placenta and produce envelope proteins associated with placenta development. As silencing of repeated elements is thought to occur mainly by DNA methylation, we compared the methylation of ERVWE1 and related HERVs to appreciate whether HERV methylation relies upon the family, the integration site, the tissue, the long terminal repeat (LTR) function or the associated gene function. CpG methylation of HERV-W LTRs in placenta-associated tissues was heterogeneous but a joint epigenetic control was found for ERVWE1 5'LTR and its juxtaposed enhancer, a mammalian apparent LTR retrotransposon. Additionally, ERVWE1, ERVFRDE1 and ERV3 5'LTRs were all essentially hypomethylated in cytotrophoblasts during pregnancy, but showed distinct and stage-dependent methylation profiles. In non-cytotrophoblastic cells, they also exhibited different methylation profiles, compatible with their respective transcriptional activities. Comparative analyses of transcriptional activity and LTR methylation in cell lines further sustained a role for methylation in the control of functional LTRs. These results suggest that HERV methylation might not be family related but copy-specific, and related to the LTR function and the tissue. In particular, ERVWE1 and ERV3 could be developmentally epigenetically regulated HERVs.