We propose a statistically valid two-stage procedure for selecting genes with expression levels correlated with that of a 'seed' gene in microarray experiments: Stage I: perform array-normal-scores (ANS) transformation of the raw microarray data and calculate the Pearson correlation coefficients using ANS-transformed data, and Stage II: calculate a resampling-based BH (rsBH) least significant-false discovery rate (LS-FDR), named LS-FDRANSrsBH to select most correlated genes based on a given LS-FDR threshold. By using simulated data sets, we show that the proposed ANS transformation is needed even after the usual data normalization process. The ANS transformation improves the detection of correlated genes, particularly the negatively correlated genes. We demonstrated the proposed procedure by searching for genes involved in insulin resistance and type 2 diabetes mellitus on a dataset the from Human Gene Expression Index (HugeIndex) database.