Reference gene expression stability obtained with the three algorithms
| geNorm M | Rank | NormFinder | Rank | BestKeeper r+ | Rank | Final rank | |
| GNB2L1 | 0.717 | 1 | 0.275 | 4 | 0.95 | 1 | 1 |
| GUSB | 0.725 | 2 | 0.260 | 3 | 0.91 | 6 | 3 |
| HPRT1 | 0.739 | 3 | 0.246 | 2 | 0.94 | 2 | 2 |
| RPL13A | 0.887 | 4 | 0.279 | 5 | 0.93 | 4 | 4 |
| B2M | 1.01 | 5 | 0.405 | 7 | 0.92 | 5 | 7 |
| ACTB | 1.119 | 6 | 0.309 | 6 | 0.94 | 3 | 5 |
| TBP | 1.175 | 7 | 0.230 | 1 | 0.91 | 7 | 6 |
| GAPDH | 1.227 | 8 | 0.411 | 8 | 0.89 | 8 | 8 |
| RNA18S | 1.336 | 9 | 0.423 | 9 | 0.73 | 9 | 9 |
GNB2L1: guanine nucleotide-binding protein, b-peptide 2-like 1; GUSB: β-glucuronidase; HPRT1: hypoxanthine ribosyltransferase 1; RPL13A: ribosomal protein L13A; B2M: β2 microglobulin; ACTB: β-actin; TBP: TATA-box binding protein; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; RNA 18S: 18 s ribosomal RNA.