Geometric mean was calculated for the pair of best reference genes and multiple reference gene set, and used for normalization. The relative expression levels are depicted as the mean ± SD, which was calculated from three biological replicates.