Differences between treatments were visualized by principal component analysis (PCA) plotting with MeV/TM4. Data were initially filtered out for missing values and then subjected to a CLEAR test that combines differential expression and variability using the GEPAS web server at http://www.gepas.org. In our case, the false discovery rate was set to a stringent level of 5 per cent.