Relative expression levels of the late embryogenesis abundant (<i>LEA</i>) and heat shock protein 70 (<i>HSP70</i>) genes under different stress conditions when the selected reference genes were used for normalization.

<p>(A) <i>LEA</i> and <i>HSP70</i> expression in the cold experimental set of <i>S</i>. <i>uncinata</i> was normalized using a combination of the two most stable genes <i>POB1</i> and <i>UFD2</i>, and the least stable genes <i>HIS</i> and <i>tIF</i>. Red bar indicates fold-change values of gene expression, which were calculated based on the normalized FPKM from the RNA-Seq analysis of the cold experimental set (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0199356#pone.0199356.s006" target="_blank">S2 Table</a>). (B) The expression of <i>LEA</i> and <i>HSP70</i> in the drought experimental set of <i>S</i>. <i>uncinata</i> was normalized using a combination of the two most stable genes <i>POB1</i> and <i>AKB</i>, and the least stable genes <i>UFD2</i> and <i>TUB250</i>. (C) The expression of <i>LEA</i> and <i>HSP70</i> in the WGF experimental set of <i>S</i>. <i>uncinata</i> was normalized using a combination of the two most stable genes <i>UFD2</i> and <i>AKB</i>, and the least stable genes <i>SuTub</i> and <i>ACT5</i>. Error bars represent the standard deviation of the mean of three biological replicates.