The miRNA and mRNA array data were analyzed for data summarization, normalization, and quality control using GeneSpring software V13 (Agilent). The default 90th percentile normalization method was performed for data preprocessing. To select the differentially expressed genes, we used fold change threshold values of ≥2 and ≤−2 and a Benjamini–Hochberg corrected P = 0.05. The data were Log2-transformed and median centered by genes using the Adjust Data function in CLUSTER 3.0 software (Stanford University, Palo Alto, CA, USA) and were then further analyzed via hierarchical clustering with average linkage. To better understand the roles of the differentially expressed mRNAs, Gene Ontology (GO) categories derived from the GO database (www.geneontology.org) were determined, and pathway analysis was performed. Three families of GO terms were identified: biological process, cellular component, and molecular function.20 (link) Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg/) analysis allowed us to determine the biological pathways that were significantly enriched with differentially expressed mRNAs.21 (link) Disease enrichment was analyzed by a web server, KOBAS 2.0 (http://kobas.cbi.pku.edu.cn/), which annotates an input set of genes with putative pathways and disease relationships based on mapping to genes with known annotations.22 (link)