Discovering local structure in gene expression data: the order-preserving submatrix problem