Selection of reference genes for reverse transcription-qPCR analysis in the biomonitor macrophyte L.

The RT-qPCR has been the method used to analyze gene expression in plants but its benefits have not been completely exploited in the field of plants ecotoxicology when used as molecular biomarkers. The correct use of RT-qPCR demands to establish a certain number of reference genes (RG) which are expected to be invariable in their expression although it does not always happen. The main goals of this work were to: (1) analyze the stability of six potential RG, (2) establish the optimum number of RG, (3) select the most suitable RG to be applied in under different test conditions and tissues and (4) confirm its convenience by normalizing the expression of one gene of interest under three different challenges. When all data were pooled together, the geNorm algorithm pointed out beta-actin and beta-tubulin (TUB) as the optimal RG pair while NormFinder algorithm selected nicotinamide adenine dinucleotide dehydrogenase (NADHD) and histone 3 (H3) as possessing the most invariable levels of expression. On the other hand, when data were grouped by tissues, ANOVA test selected H3 and TUB, while data grouped by conditions indicated that H3 and NADHD were the most stable RG under this analysis. Therefore, for a general-purpose set of RG, the overall analysis showed that a set of three RG would be optimum, and H3, TUB and NADHD were the selected ones. On the other hand, as RG can vary depending on the tissues or conditions, results achieved with ANOVA would be more reliable. Thus, appropriate normalization process would clearly need more than one RG.