Benchmarking the impact different stress assays have on Arabidopsis gene expression.

A: 22-day-old Arabidopsis growth on plates under either 1.67ξ hard agar (HA), 20 % PEG, 100 mM mannitol, or 75 mM NaCl treatments. B: Water potential measurements of treatment media (n = 3 - 4). C: Dry weight of 22-day-old Arabidopsis seedlings under different doses of each stress treatment (n = 11 - 12). D: Number and intersect of differentially expressed genes (DEGs) that are dose-responsive to each stress treatment within root and shoot tissue. E - I: Heatmaps displaying the top 50 most significant upregulated or downregulated genes in response to (E) HA, (F) PEG, (G) mannitol, (H) NaCl and (I) vermiculite drying in the Arabidopsis root (n = 2 - 3 biological replicates). Membership of Gene Ontology (GO) Terms for ‘response to stress’, ‘response to chemical stimulus’ or ‘response to ABA stimulus’ are indicated.

– Comparing HA, PEG, mannitol, and NaCl gene expression responses to vermiculite drying.

A: Heatmap displaying genes differentially expressed in response to vermiculite drying in shoot or root tissue compared to their dose-responsive expression within each plate-based assay. Level of ‘directional agreement’ (i.e., differentially expressed in the same direction) found within each assay reported. B: Overlap analysis of genes found differentially expressed due to vermiculite drying, compared to those found differentially expressed within each dose of PEG, mannitol, NaCl or HA agar assays in either shoot and root. C - D: Expression patterns of HOMEOBOX12 (HB12) and DROUGHT HYPERSENSITIVE 2 (DRY2) across each assay in root tissue. E: Shoot area of seedlings grown under increasing doses of HA, agar or nutrient concentrations (n = 19). F: Number and percent overlap of genes found differentially expressed in response to increasing doses of HA, agar or nutrient concentrations with those differentially expressed in response to vermiculite drying. G: Total shoot area of Arabidopsis accessions grown under either 1ξ or 2ξ HA agar treatment (n = 5 - 12). H: Images of Arabidopsis Trs-0 or UKSE06-325 accessions grown on either 1ξ or 2ξ HA treatment.

Plant growth responses to stress assays.

A – C: Images of 22-day-old Arabidopsis seedlings grown under different doses of each agar stress assay.

Shoot gene expression responses to each stress assay are dose responsive.

Heatmap displaying the top 50 most significant upregulated or downregulated genes in shoots in response to (A) HA, (B) PEG, (C) mannitol, (D) NaCl and (E) vermiculite drying (n = 2 - 3 biological replicates). Membership of Gene Ontology (GO) Terms for ‘response to stress’, ‘response to chemical stimulus’ or ‘response to ABA stimulus’ indicated.

Overlapping differentially expressed genes (DEGs) responsive to different assay types.

Overlap of dose-responsive differentially expressed genes in shoot (A) and root (B) in response to either HA, PEG, mannitol or NaCl (replicated from Figure 1). Number of upregulated or downregulated dose-responsive genes in response to each treatment type in shoot (C) and root (D). Overlapping gene sets in (E) shoot or (F) root tissue (permutation test, p < 0.001).

Treating vermiculite-grown Arabidopsis plants to mild drought stress.

A: Field capacity measurements of vermiculite as water evaporated over a 5-day period (n = 6 - 12). B: Shoot dry weight of Arabidopsis rosettes as they grew either under well-watered conditions or drought conditions over a 5-day period (t-test, n = 12). C: Images of plants after 5 days of water stress. D: Seed yield resulting from Arabidopsis plants after drought recovery (t-test, n = 50).

Comparing gene expression response to vermiculite drying and PEG treatment with previous studies.

A: Intersect analysis of root or shoot genes found differentially expressed in response to vermiculite drying within this study, and genes found differentially expressed in response to soil drying by Lozano-Elena et al. (16) or by Wilkins et al. (36) (permutation test, p < 0.001). B - C: Intersect analysis of genes found differentially expressed in response to PEG treatment in shoot (B) or root (C) in this study, with those found differentially expressed in response to PEG treatment by Wong et al. (37) and Wang et al. (38) (permutation test, p < 0.001). D - E: Heatmap displaying direction of shoot (D) or root (E) differentially expressed in response to vermiculite drying (this study) and Wilkins et al. (36) or Lozano-Elena et al. (16) respectively. Directional agreement with this study’s vermiculite drying response indicated. F - G : Heatmap displaying direction of genes differentially expressed in response to PEG treatment across each study. We note that both Wong et al. and Wang et al. assess transcriptomic responses of whole seedlings (both root and shoot), and thus we compare our shoot (F) and root (G) data separately. E: Examining the 119 genes that were differentially expressed in response to drought (this study), PEG treatment (this study) and PEG treatment reported in Wong et al.

Comparing ABA-induced differential expression to vermiculite drying and HA induced gene expression patterns.

A: Overlap analysis of genes found differentially expressed in response to vermiculite drying, compared to those within each dose of either transient ABA treatment, PEG, mannitol, NaCl, or HA agar assays in both root and shoot (Fisher exact test, adj. p < 0.05). B: Heatmap displaying genes differentially expressed under vermiculite drying in root tissue compared to their dose-responsive expression within each stress assay. Direction of gene expression agreement with vermiculite-drying responsive gene expression (i.e., ‘directional agreement’) indicated.

Gene expression profiles of individual genes.

A - E: Expression patterns of individual genes under doses of each assay in shoot tissue: (A) RD29B, (B) RD20, (C) HB7, (D) NCED3 and (E) P5CS1. F - H: Expression patterns of individual genes under doses of each assay in root tissue: (F) GCL1, (G) RD21 and (H) RHS18.

Physiological measurements of Arabidopsis seedlings in response to HA treatment.

A: Measurement of primary root growth rate across 8 days of growth under 1ξHA (no treatment) and 2.5ξ HA conditions (n= 16, t-test p). B: Shoot water potential measurements of seedlings grown under different HA media doses (n = 3, Pearson p = 0.009). C: Measurement of maximum quantum yield of PSII (Fv / Fm) under different HA media doses (n=4, t-test p).

Comparing the separate effects of nutrient concentration and agar concentration on seedling growth.

Image of Arabidopsis seedlings grown on either (A) 1ξ HA agar (i.e. 1X LS, 2% agar) or (B) 2.5ξ HA, which increased both nutrient and agar concentrations to 2.5X and 5%, respectively. C: Image of seedlings grown on an increased 2.5X nutrient concentration (without a change in agar concentration). D: Image of seedlings grown on an increased 5% agar concentration (without a change in nutrient concentration). E: Water potential measurements of media presented in A – D (n = 3). F: Intersection of differentially expressed genes responsive to either agar concentration, nutrient concentration, HA treatment or drought stress (permutation test, p < 0.001).

Associating hard agar’s (HA) impact on shoot size with plant fitness.

Comparing the impact HA treatment has on shoot size of 20 different Arabidopsis accessions to the change in their fitness found under drought conditions in the field, as reported in (27).

The volume of HA agar has minimal impact on gene expression.

Image of Arabidopsis seedlings grown on either (A) 75 mL of 1ξ HA agar, (B) 30 mL of 2.5ξ HA agar, or (C) 75 mL of 2.5ξ HA agar. Number of genes found differentially expressed is reported for comparisons (A) and (B), as well as (B) and (C) (DESeq, adj. p < 0.01). Of the 29 genes found differentially expressed between (B) and (C), 13 are found in comparison (A) and (B). (E) Heatmap of genes found differentially expressed in either comparison.