Explore the vast range of titles available.

Despite the prevalence of sequencing data in biomedical research, the methylome remains underrepresented. Given the importance of DNA methylation in gene regulation and disease, it is crucial to address the need for reliable differential methylation methods. This work presents a novel, transferable approach for extracting information from DNA methylation data. Our agnostic, graph-based pipeline overcomes the limitations of commonly used differential methylation techniques and addresses the “small n, big k” problem. Pheochromocytoma and Paraganglioma (PPGL) tumours with known genetic aetiologies experience extreme hypermethylation genome wide. To highlight the effectiveness of our method in candidate discovery, we present the first phenotypic classifier of PPGLs based on DNA methylation achieving 0.7 ROC-AUC. Each sample is represented by an optimised parenclitic network, a graph representing the deviation of the sample’s DNA methylation from the expected non-aggressive patterns. By extracting meaningful topological features, the dimensionality and, hence, the risk of overfitting is reduced, and the samples can be classified effectively. By using an explainable classification method, in this case logistic regression, the key CG loci influencing the decision can be identified. Our work provides insights into the molecular signature of aggressive PPGLs and we propose candidates for further research. Our optimised parenclitic network implementation improves the potential utility of DNA methylation data and offers an effective and complete pipeline for studying such datasets.

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

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.

Dopa decarboxylase (DDC) synthesizes serotonin in the developing mouse heart where it is encoded by Ddc_exon1a , a tissue-specific paternally expressed imprinted gene. Ddc_exon1a shares an imprinting control region (ICR) with the imprinted, maternally expressed (outside of the central nervous system) Grb10 gene on mouse chromosome 11, but little else is known about the tissue-specific imprinted expression of Ddc_exon1a . Fluorescent immunostaining localizes DDC to the developing myocardium in the pre-natal mouse heart, in a region susceptible to abnormal development and implicated in congenital heart defects in human. Ddc_exon1a and Grb10 are not co-expressed in heart nor in brain where Grb10 is also paternally expressed, despite sharing an ICR, indicating they are mechanistically linked by their shared ICR but not by Grb10 gene expression. Evidence from a Ddc_exon1a gene knockout mouse model suggests that it mediates the growth of the developing myocardium and a thinning of the myocardium is observed in a small number of mutant mice examined, with changes in gene expression detected by microarray analysis. Comparative studies in the human developing heart reveal a paternal expression bias with polymorphic imprinting patterns between individual human hearts at DDC_EXON1a , a finding consistent with other imprinted genes in human.

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.

Nucleotide sequences along a gene provide instructions to transcriptional and co-transcriptional machinery allowing genome expansion into the transcriptome. Interestingly, nucleotide sequence can often be shared between two genes and in some occurrences, a gene is located completely within a different gene, these are known as host/nested genes pairs. In these instances, if both genes are transcribed, overlap can result in a transcriptional crosstalk where genes regulate each other. Despite this, a comprehensive annotation of where such genes are located, and their expression patterns is lacking. To address this, we provide an up-to-date catalogue of host/nested gene pairs in mouse and human, showing that over a tenth of all genes contain a nested gene. We discovered that transcriptional co-occurrence is often tissue-specific. This co-expression was especially prevalent within the transcriptionally permissive tissue, testis. We used this developmental system and scRNA-seq analysis to demonstrate that co-expression of pairs can occur in single cells and transcription in the same place at the same time can enhance transcript diversity of the host gene. In agreement, host genes are more transcript diverse than the rest of the transcriptome and we propose that nested gene expression drives this observed diversity. Given that host/nested gene configurations were common in both human and mouse genomes, the interplay between pairs is therefore likely selected for, highlighting the relevance of transcriptional crosstalk between genes which share nucleic acid sequence. The results and analysis are available on an Rshiny application.

As cells exit the pluripotent state and begin to commit to a specific lineage they must activate genes appropriate for that lineage while silencing genes associated with pluripotency and preventing activation of lineage-inappropriate genes. The Nucleosome Remodelling and Deacetylation (NuRD) complex is essential for pluripotent cells to successfully undergo lineage commitment. NuRD controls nucleosome density at regulatory sequences to facilitate transcriptional responses, and also has been shown to prevent unscheduled transcription (transcriptional noise) in undifferentiated pluripotent cells. How these activities combine to ensure cells engage a gene expression program suitable for successful lineage commitment has not been determined. Here we show that while NuRD is not required to silence all genes, its activity is important to restrict expression of genes primed for activation upon exit from the pluripotent state, and that NuRD activity facilitates their subsequent transcriptional activation. We further show that NuRD coordinates gene expression changes, which acts to maintain a barrier between different stable states. Thus NuRD-mediated chromatin remodelling serves multiple functions, including reducing transcriptional noise, priming genes for activation and coordinating the transcriptional response to facilitate lineage commitment.Competing Interest StatementSara-Jane Dunn was an employee at Microsoft Research during this study and is currently employed at DeepMind. Neither Microsoft Research nor DeepMind have directed any aspect of the study nor exerted any commercial rights over the results. The authors declare no conflicts of interest.

How constitutive allelic methylation at imprinting control regions (ICRs) interacts with other levels of regulation to drive timely parental allele-specific expression along large imprinted domains remains partially understood. To gain insight into the regulation of the Peg13-Kcnk9 domain, an imprinted domain with important brain functions, during neural commitment, we performed an integrative analysis of the epigenetic, transcriptomic and cis-spatial organisation in an allele-specific manner in a mouse stem cell-based model of corticogenesis that recapitulates the control of imprinted gene expression during neurodevelopment. We evidence that despite an allelic higher-order chromatin structure associated with the paternally CTCF-bound Peg13 ICR, the enhancer-Kcnk9 promoter contacts can occur on both alleles, although they are only productive on the maternal allele. This observation challenges the canonical model in which CTCF binding isolates the enhancer and its target gene on either side, and suggests a more nuanced role for allelic CTCF binding at some ICRs.

Enhancers are genomic DNA sequences that bind transcription factors, chromatin regulators and non-coding transcripts to modulate the expression of target genes. They have been found to act from different locations relative to a gene and to modulate the activity of promoters up to ~1 Mb away. Here we report the first 3D genome structures of single mouse ES cells as they are induced to exit pluripotency, transition through a formative stage and undergo neuroectodermal differentiation. In order to directly study how interactions between enhancers and promoters are reconfigured genome wide we have determined 3D structures of haploid cells using single cell Hi-C, where we can unambiguously map the contacts to particular chromosomes. We find that there is a remarkable reorganisation of 3D genome structure where inter-chromosomal intermingling increases dramatically in the formative state. This intermingling is associated with the formation of a large number of multiway hubs that bring together enhancers and promoters with similar chromatin states from typically 5-8 distant chromosomal sites that are often separated by many Mb from each other. Genes important for pluripotency exit establish contacts with emerging enhancers within multiway chromatin hubs as cells enter the formative state, consistent with these structural changes playing an important role in establishing new cell identities. Furthermore, we find that different multiway chromatin hubs, and thus different enhancer-promoter interactions, are formed in different individual cells. Our results suggest that studying genome structure in single cells will be important to identify key changes in enhancer-promoter interactions that occur as cells undergo developmental transitions.Competing Interest StatementThe authors have declared no competing interest.

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 lifespan. 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.