Fig. 1

Skeletal muscle group distribution and Six1 dependence of gene expression. DNA microarray gene expression profiling from Sakakibara et al. was analyzed to identify differentially expressed genes in the gastrocnemius muscles of wild-type versus Six1-cKO mice. The cut-offs used were abs (log₂FC) > 0.58 and Benjamini-Hochberg adjusted p value < 0.05. Genes were annotated based on whether their transcription start sites are within 50 kb of a Six1-binding site in primary myotubes (red bars). Pearson correlation hierarchical clustering was applied to rows (genes) and columns (samples). To further annotate these Six1-dependent genes, their expression from RNA-seq performed in several muscle tissue groups from the muscleDB database (Terry et al.) was also plotted. The rows retained the clustering solution order from the Sakakibara study, and columns were independently clustered by Pearson correlation (clustering performed separately from the Sakakibara data). The log₂ fold changes of each gene in the two experiments are shown in the vertical histograms, blue and yellow bars indicating negative and positive values, respectively. A red arrowhead indicates where the Slc16a10 gene appears in these graphs. The gene symbols in this figure are also listed in supplementary Table S7