L1 chimeric transcripts are expressed in healthy brain and their deregulation in glioma follows that of their host locus.

Besides the consequences of retrotransposition, long interspersed element 1 (L1) retrotransposons can affect the host genome through their antisense promoter. In addition to the sense promoter, the evolutionarily recent L1 retrotransposons, which are present in several thousand copies, also possess an anti-sense promoter that can produce L1 chimeric transcripts (LCT) composed of the L1 5' UTR followed by the adjacent genomic sequence. The full extent to which LCT expression occurs in a given tissue and whether disruption of the defense mechanisms that normally control L1 retrotransposons affects their expression and function in cancer cells, remain to be established. By using CLIFinder, a dedicated bioinformatics tool, we found that LCT expression was widespread in normal brain and aggressive glioma samples, and that approximately 17% of recent L1 retrotransposons, from the L1PA1 to L1PA7 subfamilies, were involved in their production. Importantly, the transcriptional activities of the L1 antisense promoters and of their host loci were coupled. Accordingly, we detected LCT-producing L1 retrotransposons mainly in transcriptionally active genes and genomic loci. Moreover, changes in the host genomic locus expression level in glioma were associated with a similar change in LCT expression level, regardless of the L1 promoter methylation status. Our findings support a model in which the host genomic locus transcriptional activity is the main driving force of LCT expression. We hypothesize that this model is more applicable when host gene and LCT are transcribed from the same strand.

Transcriptional alterations in glioma result primarily from DNA methylation-independent mechanisms.

In cancer cells, aberrant DNA methylation is commonly associated with transcriptional alterations, including silencing of tumor suppressor genes. However, multiple epigenetic mechanisms, including polycomb repressive marks, contribute to gene deregulation in cancer. To dissect the relative contribution of DNA methylation-dependent and -independent mechanisms to transcriptional alterations at CpG island/promoter-associated genes in cancer, we studied 70 samples of adult glioma, a widespread type of brain tumor, classified according to their isocitrate dehydrogenase (IDH1) mutation status. We found that most transcriptional alterations in tumor samples were DNA methylation-independent. Instead, altered histone H3 trimethylation at lysine 27 (H3K27me3) was the predominant molecular defect at deregulated genes. Our results also suggest that the presence of a bivalent chromatin signature at CpG island promoters in stem cells predisposes not only to hypermethylation, as widely documented, but more generally to all types of transcriptional alterations in transformed cells. In addition, the gene expression strength in healthy brain cells influences the choice between DNA methylation- and H3K27me3-associated silencing in glioma. Highly expressed genes were more likely to be repressed by H3K27me3 than by DNA methylation. Our findings support a model in which altered H3K27me3 dynamics, more specifically defects in the interplay between polycomb protein complexes and the brain-specific transcriptional machinery, is the main cause of transcriptional alteration in glioma cells. Our study provides the first comprehensive description of epigenetic changes in glioma and their relative contribution to transcriptional changes. It may be useful for the design of drugs targeting cancer-related epigenetic defects.

Chromatin immunoprecipitation indirect peaks highlight long-range interactions of insulator proteins and Pol II pausing.

Eukaryotic chromosomes are partitioned into topologically associating domains (TADs) that are demarcated by distinct insulator-binding proteins (IBPs) in Drosophila. Whether IBPs regulate specific long-range contacts and how this may impact gene expression remains unclear. Here we identify "indirect peaks" of multiple IBPs that represent their distant sites of interactions through long-range contacts. Indirect peaks depend on protein-protein interactions among multiple IBPs and their common cofactors, including CP190, as confirmed by high-resolution analyses of long-range contacts. Mutant IBPs unable to interact with CP190 impair long-range contacts as well as the expression of hundreds of distant genes that are specifically flanked by indirect peaks. Regulation of distant genes strongly correlates with RNAPII pausing, highlighting how this key transcriptional stage may trap insulator-based long-range interactions. Our data illustrate how indirect peaks may decipher gene regulatory networks through specific long-range interactions.

TET3 controls the expression of the H3K27me3 demethylase Kdm6b during neural commitment.

The acquisition of cell identity is associated with developmentally regulated changes in the cellular histone methylation signatures. For instance, commitment to neural differentiation relies on the tightly controlled gain or loss of H3K27me3, a hallmark of polycomb-mediated transcriptional gene silencing, at specific gene sets. The KDM6B demethylase, which removes H3K27me3 marks at defined promoters and enhancers, is a key factor in neurogenesis. Therefore, to better understand the epigenetic regulation of neural fate acquisition, it is important to determine how Kdm6b expression is regulated. Here, we investigated the molecular mechanisms involved in the induction of Kdm6b expression upon neural commitment of mouse embryonic stem cells. We found that the increase in Kdm6b expression is linked to a rearrangement between two 3D configurations defined by the promoter contact with two different regions in the Kdm6b locus. This is associated with changes in 5-hydroxymethylcytosine (5hmC) levels at these two regions, and requires a functional ten-eleven-translocation (TET) 3 protein. Altogether, our data support a model whereby Kdm6b induction upon neural commitment relies on an intronic enhancer the activity of which is defined by its TET3-mediated 5-hmC level. This original observation reveals an unexpected interplay between the 5-hmC and H3K27me3 pathways during neural lineage commitment in mammals. It also questions to which extent KDM6B-mediated changes in H3K27me3 level account for the TET-mediated effects on gene expression.

The Long Non-Coding RNA HOXA-AS2 Promotes Proliferation of Glioma Stem Cells and Modulates Their Inflammation Pathway Mainly through Post-Transcriptional Regulation.

Glioblastomas represent approximatively half of all gliomas and are the most deadly and aggressive form. Their therapeutic resistance and tumor relapse rely on a subpopulation of cells that are called Glioma Stem Cells (GSCs). Here, we investigated the role of the long non-coding RNA HOXA-AS2 in GSC biology using descriptive and functional analyses of glioma samples classified according to their isocitrate dehydrogenase (IDH) gene mutation status, and of GSC lines. We found that HOXA-AS2 is overexpressed only in aggressive (IDHwt) glioma and GSC lines. ShRNA-based depletion of HOXA-AS2 in GSCs decreased cell proliferation and altered the expression of several hundreds of genes. Integrative analysis revealed that these expression changes were not associated with changes in DNA methylation or chromatin signatures at the promoter of the majority of genes deregulated following HOXA-AS2 silencing in GSCs, suggesting a post-transcriptional regulation. In addition, transcription factor binding motif enrichment and correlation analyses indicated that HOXA-AS2 affects, directly or indirectly, the expression of key transcription factors implicated in GCS biology, including E2F8, E2F1, STAT1, and ATF3, thus contributing to GCS aggressiveness by promoting their proliferation and modulating the inflammation pathway.

High-salt-recovered sequences are associated with the active chromosomal compartment and with large ribonucleoprotein complexes including nuclear bodies.

The mammalian cell nucleus contains numerous discrete suborganelles named nuclear bodies. While recruitment of specific genomic regions into these large ribonucleoprotein (RNP) complexes critically contributes to higher-order functional chromatin organization, such regions remain ill-defined. We have developed the high-salt-recovered sequences-sequencing (HRS-seq) method, a straightforward genome-wide approach whereby we isolated and sequenced genomic regions associated with large high-salt insoluble RNP complexes. By using mouse embryonic stem cells (ESCs), we showed that these regions essentially correspond to the most highly expressed genes, and to cis-regulatory sequences like super-enhancers, that belong to the active A chromosomal compartment. They include both cell-type-specific genes, such as pluripotency genes in ESCs, and housekeeping genes associated with nuclear bodies, such as histone and snRNA genes that are central components of Histone Locus Bodies and Cajal bodies. We conclude that HRSs are associated with the active chromosomal compartment and with large RNP complexes including nuclear bodies. Association of such chromosomal regions with nuclear bodies is in agreement with the recently proposed phase separation model for transcription control and might thus play a central role in organizing the active chromosomal compartment in mammals.

CLIFinder: identification of LINE-1 chimeric transcripts in RNA-seq data.

Summary: L1 Chimeric Transcripts (LCTs) are initiated by repeated LINE-1 element antisense promoters and include the L1 5'UTR sequence in antisense orientation followed by the adjacent genomic region. LCTs have been characterized mainly using bioinformatics approaches to query dbEST. To take advantage of NGS data to unravel the transcriptome composition, we developed Chimeric LIne Finder (CLIFinder), a new bioinformatics tool. Using stranded paired-end RNA-seq data, we demonstrated that CLIFinder can identify genome-wide transcribed chimera sequences corresponding to potential LCTs. Moreover, CLIFinder can be adapted to study transcription from other repeat types. Availability and implementation: The code is available at: https://github.com/GReD-Clermont/CLIFinder; and for Galaxy users, it is directly accessible in the tool shed at: https://toolshed.g2.bx.psu.edu/view/clifinder/clifinder/. Contact: catherine.barriere@uca.fr. Supplementary information: Supplementary data are available at Bioinformatics online.

Human Oocyte-Derived Methylation Differences Persist in the Placenta Revealing Widespread Transient Imprinting.

Thousands of regions in gametes have opposing methylation profiles that are largely resolved during the post-fertilization epigenetic reprogramming. However some specific sequences associated with imprinted loci survive this demethylation process. Here we present the data describing the fate of germline-derived methylation in humans. With the exception of a few known paternally methylated germline differentially methylated regions (DMRs) associated with known imprinted domains, we demonstrate that sperm-derived methylation is reprogrammed by the blastocyst stage of development. In contrast a large number of oocyte-derived methylation differences survive to the blastocyst stage and uniquely persist as transiently methylated DMRs only in the placenta. Furthermore, we demonstrate that this phenomenon is exclusive to primates, since no placenta-specific maternal methylation was observed in mouse. Utilizing single cell RNA-seq datasets from human preimplantation embryos we show that following embryonic genome activation the maternally methylated transient DMRs can orchestrate imprinted expression. However despite showing widespread imprinted expression of genes in placenta, allele-specific transcriptional profiling revealed that not all placenta-specific DMRs coordinate imprinted expression and that this maternal methylation may be absent in a minority of samples, suggestive of polymorphic imprinted methylation.

Imprinting control regions (ICRs) are marked by mono-allelic bivalent chromatin when transcriptionally inactive.

Parental allele-specific expression of imprinted genes is mediated by imprinting control regions (ICRs) that are constitutively marked by DNA methylation imprints on the maternal or paternal allele. Mono-allelic DNA methylation is strictly required for the process of imprinting and has to be faithfully maintained during the entire life-span. While the regulation of DNA methylation itself is well understood, the mechanisms whereby the opposite allele remains unmethylated are unclear. Here, we show that in the mouse, at maternally methylated ICRs, the paternal allele, which is constitutively associated with H3K4me2/3, is marked by default by H3K27me3 when these ICRs are transcriptionally inactive, leading to the formation of a bivalent chromatin signature. Our data suggest that at ICRs, chromatin bivalency has a protective role by ensuring that DNA on the paternal allele remains unmethylated and protected against spurious and unscheduled gene expression. Moreover, they provide the proof of concept that, beside pluripotent cells, chromatin bivalency is the default state of transcriptionally inactive CpG island promoters, regardless of the developmental stage, thereby contributing to protect cell identity.

Maternal mutations of FOXF1 cause alveolar capillary dysplasia despite not being imprinted.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare cause of pulmonary hypertension in newborns. Maternally inherited point mutations in Forkhead Box F1 gene (FOXF1), deletions of the gene, or its long-range enhancers on the maternal allele are responsible for this neonatal lethal disorder. Here, we describe monozygotic twins and one full-term newborn with ACD and gastrointestinal malformations caused by de novo mutations of FOXF1 on the maternal-inherited alleles. Since this parental transmission is consistent with genomic imprinting, the parent-of-origin specific monoallelic expression of genes, we have undertaken a detailed analysis of both allelic expression and DNA methylation. FOXF1 and its neighboring gene FENDRR were both biallelically expressed in a wide range of fetal tissues, including lung and intestine. Furthermore, detailed methylation screening within the 16q24.1 regions failed to identify regions of allelic methylation, suggesting that disrupted imprinting is not responsible for ACDMPV.