Fig. 3

Dynamic changes of H3K18la reflect transcriptional adaptations. A Box plots showing H3K18la log₂FC changes from MT versus MB over different genomic features. B Scatterplot showing the correlation between significant (FDR < 0.05) H3K18la log₂FC (>0.5) in promoters and their corresponding gene expression log₂FC (>0.5) based on the overlapping genes from MT versus MB differential analysis. Pearson correlation coefficient R is indicated. C IGV genome browser [50] snapshot of H3K18la profile at the Myh1 promoter region from various mouse samples and the corresponding SEACR-called peak regions. H3K18la level is depicted (rpm). Genomic regions are indicated on the top, as well as RefSeq gene names. D Fold enrichment of significant differentially H3K18la-marked peaks (FDR < 0.05, |log₂FC| > 1.5) in ENCODE cCREs. E Scatterplot showing the correlation between significant (FDR < 0.05) H3K18la log₂FC (>0.5) in dELS and their closest gene expression log₂FC (>0.5) based on the overlapping genes from MT versus MB differential analysis. Pearson correlation coefficient R is indicated. F Top 10 GO terms (category “Biological Process”) based on the GO analysis of the overlapping upregulated genes in MT from E (first quadrant red dots). G Box plots showing H3K18la log₂FC of peaks overlapping with MB- or MT-specific enhancers and of peaks not overlapping with these enhancers. Wilcoxon test p-values are indicated for each pair of groups. H Gene expression changes (log₂FC) after treatment of MB with 10 mM sodium-L-lactate of genes containing a H3K18la peak in MB, in MT, in both or in none of both. Wilcoxon test p-values are indicated for each pair of groups. *p value < 0.05, **p value < 0.01, ***p value < 0.001, ****p value < 0.0001