Differential contribution to gene expression prediction of histone modifications at enhancers or promoters

Mitotic bookmarking redundancy by nuclear receptors in pluripotent cells

Mitotic bookmarking transcription factors (TFs) are thought to mediate rapid and accurate reactivation after mitotic gene silencing. However, the loss of individual bookmarking TFs often leads to the deregulation of only a small proportion of their mitotic targets, raising doubts on the biological significance and importance of their bookmarking function. Here we used targeted proteomics of the mitotic bookmarking TF ESRRB, an orphan nuclear receptor, to discover a large redundancy in mitotic binding among members of the protein super-family of nuclear receptors. Focusing on the nuclear receptor NR5A2, which together with ESRRB is essential in maintaining pluripotency in mouse embryonic stem cells, we demonstrate conjoint bookmarking activity of both factors on promoters and enhancers of a large fraction of active genes, particularly those most efficiently reactivated in G1. Upon fast and simultaneous degradation of both factors during mitotic exit, hundreds of mitotic targets of ESRRB/NR5A2, including key players of the pluripotency network, display attenuated transcriptional reactivation. We propose that redundancy in mitotic bookmarking TFs, especially nuclear receptors, confers robustness to the reestablishment of gene regulatory networks after mitosis.

Thymine DNA glycosylase regulates cell-cycle-driven p53 transcriptional control in pluripotent cells

Cell cycle progression is linked to transcriptome dynamics and variations in the response of pluripotent cells to differentiation cues, mostly through unknown determinants. Here, we characterized the cell-cycle-associated transcriptome and proteome of mouse embryonic stem cells (mESCs) in naive ground state. We found that the thymine DNA glycosylase (TDG) is a cell-cycle-regulated co-factor of the tumor suppressor p53. Furthermore, TDG and p53 co-bind ESC-specific cis-regulatory elements and thereby control transcription of p53-dependent genes during self-renewal. We determined that the dynamic expression of TDG is required to promote the cell-cycle-associated transcriptional heterogeneity. Moreover, we demonstrated that transient depletion of TDG influences cell fate decisions during the early differentiation of mESCs. Our findings reveal an unanticipated role of TDG in promoting molecular heterogeneity during the cell cycle and highlight the central role of protein dynamics for the temporal control of cell fate during development.

A computational pipeline to learn gene expression predictive models from epigenetic information at enhancers or promoters

Here, we present a computational pipeline to obtain quantitative models that characterize the relationship of gene expression with the epigenetic marking at enhancers or promoters in mouse embryonic stem cells. Our protocol consists of (i) generating predictive models of gene expression from epigenetic information (such as histone modification ChIP-seq) at enhancers or promoters and (ii) assessing the performance of these predictive models. This protocol could be applied to other biological scenarios or other types of epigenetic data. For complete details on the use and execution of this protocol, please refer to Gonzalez-Ramirez et al. (2021).1.

Screening and identification of hub gene and differential gene and mutation sequence analysis of related genes in colorectal cancer based on bioinformatics analysis

Background

At present, the research of genomics is in ascendency, and using bioinformatics analysis methods to systematically explore the pathogenic genes and their regulatory mechanisms will play a great role in promoting the research of cancer. This study was to search The Cancer Genome Atlas (TCGA) database and extract inflammation-related non-coding RNA to construct a prognosis model of colon cancer and search for new immunotherapeutic targets.

Methods

The transcriptome sequencing data and clinical data of 396 colon cancer patients were downloaded from TCGA database, and the inflammation-related non-coding RNA was obtained from the non-coding RNAs in Inflammation (ncRI) database. The prognostic model was constructed by univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression, and the optimal grouping threshold of risk score was determined by X-Tile software. The patients were risk stratified to further explore the differences in immune cell infiltration and biological function between the high- and low-risk groups.

Results

The TCGA dataset of colon cancer was included to screen out 120 differentially expressed genes (DEGs) that overlapped in the 2 datasets, among which 29 genes were up-regulated and 91 genes were down-regulated. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the above 120 DEGs showed that proximal tubule sodium bicarbonate recovery, nitrogen metabolism, pancreatic fluid secretion, and PPAR signaling pathways were closely related to the occurrence of colon cancer. The expression of copper death-related genes was significantly correlated with the correlation coefficient of colon cancer (P<0.01). Gene Ontology analysis showed that the DEGs were mainly enriched in messenger RNA processing, RNA splicing, small G protein-mediated signal transduction, adhesion junction, mitochondrial matrix, mitochondrial protein complex, chromatin binding, small G protein binding, and Ras G protein binding, among others. KEGG analysis showed that the DEGs were enriched in the following pathways: herpes simplex virus type 1 infection, pathways of neurodegenerative diseases, Huntington's disease, prion disease, Parkinson's disease, the Ras signaling pathway, and so on.

Conclusions

The key genes closely related to colon cancer were effectively screened by the bioinformatics method, which provided a theoretical basis for further study of its mechanism.

OCT4 activates a Suv39h1-repressive antisense lncRNA to couple histone H3 Lysine 9 methylation to pluripotency

Histone H3 Lysine 9 (H3K9) methylation, a characteristic mark of heterochromatin, is progressively implemented during development to contribute to cell fate restriction as differentiation proceeds. Accordingly, in undifferentiated and pluripotent mouse Embryonic Stem (ES) cells the global levels of H3K9 methylation are rather low and increase only upon differentiation. How global H3K9 methylation levels are coupled with the loss of pluripotency remains largely unknown. Here, we identify SUV39H1, a major H3K9 di- and tri-methylase, as an indirect target of the pluripotency network of Transcription Factors (TFs). We find that pluripotency TFs, principally OCT4, activate the expression of Suv39h1as, an antisense long non-coding RNA to Suv39h1. In turn, Suv39h1as downregulates Suv39h1 transcription in cis via a mechanism involving the modulation of the chromatin status of the locus. The targeted deletion of the Suv39h1as promoter region triggers increased SUV39H1 expression and H3K9me2 and H3K9me3 levels, affecting all heterochromatic regions, particularly peri-centromeric major satellites and retrotransposons. This increase in heterochromatinization efficiency leads to accelerated and more efficient commitment into differentiation. We report, therefore, a simple genetic circuitry coupling the genetic control of pluripotency with the global efficiency of H3K9 methylation associated with a major cell fate restriction, the irreversible loss of pluripotency.

Argonaute proteins regulate a specific network of genes through KLF4 in mouse embryonic stem cells

The Argonaute proteins (AGOs) are well known for their role in post-transcriptional gene silencing in the microRNA (miRNA) pathway. Here we show that in mouse embryonic stem cells, AGO1&2 serve additional functions that go beyond the miRNA pathway. Through the combined deletion of both Agos, we identified a specific set of genes that are uniquely regulated by AGOs but not by the other miRNA biogenesis factors. Deletion of Ago2&1 caused a global reduction of the repressive histone mark H3K27me3 due to downregulation at protein levels of Polycomb repressive complex 2 components. By integrating chromatin accessibility, prediction of transcription factor binding sites, and chromatin immunoprecipitation sequencing data, we identified the pluripotency factor KLF4 as a key modulator of AGO1&2-regulated genes. Our findings revealed a novel axis of gene regulation that is mediated by noncanonical functions of AGO proteins that affect chromatin states and gene expression using mechanisms outside the miRNA pathway.

•

AGO1&2 regulate a specific set of genes in mESCs, independently of the miRNA pathway

•

PRC2 proteins are downregulated in Ago2&1_KO mESCs, leading to H3K27me3 global loss

•

AGO1&2 regulate gene expression through the pluripotency factor KLF4