Identification of genes under dynamic post-transcriptional regulation from time-series epigenomic data

2-Hydroxyglutarate modulates histone methylation at specific loci and alters gene expression via Rph1 inhibition

2-Hydroxyglutarate (2-HG) is an oncometabolite that accumulates in certain cancers. Gain-of-function mutations in isocitrate dehydrogenase lead to 2-HG accumulation at the expense of alpha-ketoglutarate. Elevated 2-HG levels inhibit histone and DNA demethylases, causing chromatin structure and gene regulation changes with tumorigenic consequences. We investigated the effects of elevated 2-HG levels in Saccharomyces cerevisiae, a yeast devoid of DNA methylation and heterochromatin-associated histone methylation. Our results demonstrate genetic background-dependent gene expression changes and altered H3K4 and H3K36 methylation at specific loci. Analysis of histone demethylase deletion strains indicated that 2-HG inhibits Rph1 sufficiently to induce extensive gene expression changes. Rph1 is the yeast homolog of human KDM4 demethylases and, among the yeast histone demethylases, was the most sensitive to the inhibitory effect of 2-HG in vitro. Interestingly, Rph1 deficiency favors gene repression and leads to further down-regulation of already silenced genes marked by low H3K4 and H3K36 trimethylation, but abundant in H3K36 dimethylation. Our results provide novel insights into the genome-wide effects of 2-HG and highlight Rph1 as its preferential demethylase target.

Bioinformatics prediction and experimental verification of key biomarkers for diabetic kidney disease based on transcriptome sequencing in mice

Background

Diabetic kidney disease (DKD) is the leading cause of death in people with type 2 diabetes mellitus (T2DM). The main objective of this study is to find the potential biomarkers for DKD.

Materials and Methods

Two datasets (GSE86300 and GSE184836) retrieved from Gene Expression Omnibus (GEO) database were used, combined with our RNA sequencing (RNA-seq) results of DKD mice (C57 BLKS-32w db/db) and non-diabetic (db/m) mice for further analysis. After processing the expression matrix of the three sets of data using R software "Limma", differential expression analysis was performed. The significantly differentially expressed genes (DEGs) (-logFC- > 1, p-value < 0.05) were visualized by heatmaps and volcano plots respectively. Next, the co-expression genes expressed in the three groups of DEGs were obtained by constructing a Venn diagram. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were further analyzed the related functions and enrichment pathways of these co-expression genes. Then, qRT-PCR was used to verify the expression levels of co-expression genes in the kidney of DKD and control mice. Finally, protein-protein interaction network (PPI), GO, KEGG analysis and Pearson correlation test were performed on the experimentally validated genes, in order to clarify the possible mechanism of them in DKD.

Results

Our RNA-seq results identified a total of 125 DEGs, including 59 up-regulated and 66 down-regulated DEGs. At the same time, 183 up-regulated and 153 down-regulated DEGs were obtained in GEO database GSE86300, and 76 up-regulated and 117 down-regulated DEGs were obtained in GSE184836. Venn diagram showed that 13 co-expression DEGs among the three groups of DEGs. GO analysis showed that biological processes (BP) were mainly enriched inresponse to stilbenoid, response to fatty acid, response to nutrient, positive regulation of macrophage derived foam cell differentiation, triglyceride metabolic process. KEGG pathway analysis showed that the three major enriched pathways were cholesterol metabolism, drug metabolism-cytochrome P450, PPAR signaling pathway. After qRT-PCR validation, we obtained 11 genes that were significant differentially expressed in the kidney tissues of DKD mice compared with control mice. (The mRNA expression levels of Aacs, Cpe, Cd36, Slc22a7, Slc1a4, Lpl, Cyp7b1, Akr1c14 and Apoh were declined, whereas Abcc4 and Gsta2 were elevated).

Conclusion

Our study, based on RNA-seq results, GEO databases and qRT-PCR, identified 11 significant dysregulated DEGs, which play an important role in lipid metabolism and the PPAR signaling pathway, which provide novel targets for diagnosis and treatment of DKD.