used as the gene expression value. Genes with lower levels of expression; i.e. with average counts across all runs <30, were omitted. After conversion to a base-2 logarithm with a pseudo-count of 1, quantile normalization was applied to the data from each experiment, and the average expression levels were subtracted from each gene for each experiment. Using all experiments at once, Pearson's correlation coefficients for each gene pair were calculated, and these values were transferred to the Mutual Rank (MR) value (16), which is the geometric average of asymmetric ranks in coexpressed gene lists.