GlcNAc modulates mast cell DNA accessibility in human KIT D816 MCs. (A) Mean-difference plot showing the log-fold change and logCPM (log counts par millions) reads of each peak. Peaks with fold changes significantly greater than 2 (false discovery rate <0.05) are highlighted. Significantly open regions (up peaks = 199) and closed regions (down peaks = 356) are highlighted in red and green, respectively. (B) Venn diagram showing the number of common opened regions in KITD816 vs KIT WT and KITD816V vs KITD816V + GlcNAc. (C) Distributions of ATAC-seq peaks in different regions and the proportion of the genomic regions in each category using functional state annotations from Roadmap cell data sets genomic feature in OCR were shown. (D) Bubble plot of Gene Ontology (GO) terms. The z-scores are shown on the x-axis, and the negative log P values are shown on the y-axis. The bubble area is directly proportional to the number of genes associated with a certain GO term. The MAPK signaling pathway had the largest number of genes and the lowest P value. (E) Altered gene expression by GlcNAc. ROSA KIT D816V cells (n = 3) were either left untreated (control) or treated with 20 mM GlcNAc for 24 hours. Selected candidate genes from Figure 6D were validated by reverse transcriptase-quantitative PCR. Bar graph shows relative quantity (RQ) of transcripts in respect to HPRT transcripts. (F) Rankings of the most differentially enriched motifs between +GlcNAc and −GlcNAc samples. The top motifs were the most enriched motifs and the bottom motifs were the less enriched motifs in the +GlcNAc samples. The values in the graphs are presented as the means ± SD (nonsignificant, *P < .05, **P < .01, ***P < .001; 2-tailed, unpaired Student t test).