Selection of Appropriate Reference Genes for Gene Expression Analysis under Abiotic Stresses in Salix viminalis

Salix viminalis is a fast growing willow species with potential as a plant used for biomass feedstock or for phytoremediation. However, few reference genes (RGs) for quantitative real-time polymerase chain reaction (qPCR) are available in S. viminalis, thereby limiting gene expression studies. Here, we investigated the expression stability of 14 candidate reference genes (RGs) across various organs exposed to five abiotic stresses (cold, heat, drought, salt, and poly-metals). Four RGs ranking algorithms, namely geNorm(PLUS), BestKeeper, NormFinder, and GrayNorm were applied to analyze the qPCR data and the outputs were merged into consensus lists with RankAggreg, a rank aggregation algorithm. In addition, the optimal RG combinations were determined with geNorm(PLUS) and GrayNorm. The genes that were the most stable in the roots were TIP41 and CDC2. In the leaves, TIP41 was the most stable, followed by EF1b and ARI8, depending on the condition tested. Conversely, GAPDH and beta -TUB, two genes commonly used for qPCR data normalization were the least stable across all organs. Nevertheless, both geNorm(PLUS) and GrayNorm recommended the use of a combination of genes rather than a single one. These results are valuable for research of transcriptomic responses in different S. viminalis organs.