Comparison of Leaf Sheath Transcriptome Profiles with Physiological Traits of Bread Wheat Cultivars under Salinity Stress

Salinity stress has significant negative effects on plant biomass production and crop yield. Salinity tolerance is controlled by complex systems of gene expression and ion transport. The relationship between specific features of mild salinity stress adaptation and gene expression was analyzed using four commercial varieties of bread wheat (Triticum aestivum) that have different levels of salinity tolerance. The high-throughput phenotyping system in The Plant Accelerator at the Australian Plant Phenomics Facility revealed variation in shoot relative growth rate and salinity tolerance among the four cultivars. Comparative analysis of gene expression in the leaf sheaths identified genes whose functions are potentially linked to shoot biomass development and salinity tolerance. Early responses to mild salinity stress through changes in gene expression have an influence on the acquisition of stress tolerance and improvement in biomass accumulation during the early “osmotic” phase of salinity stress. In addition, results revealed transcript profiles for the wheat cultivars that were different from those of usual stress-inducible genes, but were related to those of plant growth. These findings suggest that, in the process of breeding, selection of specific traits with various salinity stress-inducible genes in commercial bread wheat has led to adaptation to mild salinity conditions.


Introduction
Bread wheat (Triticum aestivum) is now grown on a greater area than any other cereal crop, but is only moderately tolerant to salinity stress [1]. With soil salinity affecting one-fifth of irrigated agricultural land worldwide [2], there is an urgent need to improve the salinity tolerance of bread wheat to meet future global food requirements.
Shoot Na + exclusion is only one component of plant salinity tolerance, and some studies have found no clear correlation between plant salinity tolerance and Na + exclusion under conditions of moderate salinity stress [3][4][5][6]. Another important, but not well characterized, tolerance trait is the shoot ion accumulation independent tolerance mechanism (also called "osmotic tolerance"), which is the ability of plants to maintain growth rates after initial exposure to salinity stress that is independent of the extent of shoot ion accumulation [1,7]. This osmotic tolerance mechanism is important for maintaining tissue expansion and tillering during the early stages of salinity stress before Na + and Clions accumulate to toxic concentrations in the shoot. Plants with poor osmotic tolerance will show significant reduction in growth immediately after salt application. This reduction in growth will continue throughout the time of exposure to salinity, resulting in plants with reduced biomass and yield [1,5,7]. Unfortunately, the mechanisms behind osmotic tolerance are poorly understood, and no gene for increased tolerance has yet been linked specifically to this trait. To advance breeding for salinity tolerance in bread wheat, it is imperative to identify the genes important for osmotic stress tolerance.
It is well known that thousands of genes are differentially regulated in response to a range of abiotic stresses including salinity and drought [8,9], due in part to plant stress-tolerance mechanisms in the early phase of stress responses. In wheat, several genes have been reported to be induced rapidly after salinity stress. The salinity stress-inducible transcription factor (TF), TaWRKY10, has been found to be important in regulating drought responses and the accumulation of reactive oxygen species [10]. Wheat calcineurin B-like protein-interacting protein kinase 29 (TaCIPK29) and sucrose nonfermenting 1-related protein kinase 2 (TaSnRK2) are also induced by abiotic stress and are important in plant stress signaling pathways that allow rapid control of multiple tolerance mechanisms [11][12][13][14][15][16][17][18]. However, it is still unclear whether these salinity stress-inducible genes are associated with salt or osmotic stress tolerance.
In the present study, four commercial bread wheat cultivars were exposed to salinity stress conditions to characterize the transcriptional response to early stages of osmotic stress in the key tissue affecting growth over this time period, the leaf sheaths. The high-throughput phenotyping system of The Plant Accelerator at the Australian Plant Phenomics Facility in South Australia allowed the continuous measurement of growth rates to determine the osmotic stress tolerance component of all four cultivars. A comprehensive evaluation based on gene expression profiles of leaf sheaths and phenomics highlighted gene sets associated with individual physiological traits. We discuss the correlation between the salinity stress tolerance of these commercial salinity-tolerant cultivars and their leaf sheath transcriptomes during the early response of shoots to salinity stress. permissions were required for the cultivation in the greenhouse of The Plant Accelerator at the Australian Plant Phenomics Facility described in the text because the non-transformant crops were grown and analyzed in the hermetically-closed greenhouse approved for academic plant science. The field studied did not involve endangered or protected species.
Seeds were soaked in reverse osmosis (RO) water at room temperature for 2 h. The water was drained and the seeds placed in the dark at 4°C for 2 days in sealed 50 mL polypropylene tubes. Evenly sized imbibed seeds were selected for sowing and four seeds per pot were planted 2 cm deep into 145 cm diameter by 190 cm height free-draining pots. The soil substrate was a potting mix consisting of 50% (v/v) University of California mix, 35% (v/v) peat mix and 15% (v/v) clay loam, designed specifically for use on The Plant Accelerator's conveyor system. Plants were watered to a gravimetric water content of 25% (w/w). The seedlings were grown in natural light in a temperature-controlled greenhouse, with a daytime temperature of 22°C and 15°C at night. The pots were randomly placed in the greenhouse. At the time of emergence of leaf 2, seedlings were thinned, leaving one seedling per pot. Two experiments were performed with different concentrations of NaCl to measure the effects of mild salinity stress. The first experiment measured the growth response of wheat seedlings treated with 0 and 75 mM NaCl. The second experiment measured the growth response of seedlings treated with 0 and 100 mM NaCl, with RNA also extracted from leaf sheath tissue for microarray analysis. The experiments were conducted consecutively during the winter months of 2011 (from late May until early August) in Adelaide, South Australia (34°58'17.00"S, 138°38'23.00"E).

Assay for growth under salt stress conditions
Following emergence of the third leaf, the free-draining pots were placed in a deep saucer and loaded onto the conveyers of The Plant Accelerator's automated imaging system [19], using a randomized block design (Brien et al, 2013), as shown in S1 and S2 Figs. Daily imaging was performed as described by Brien et al. (2013). Seedlings were watered to 25% (w/w) gravimetric water content. At the emergence of leaf four, the salinity treatment was applied as described by Berger et al. (2012). The salinity treatment solution (in a volume of 210 mL) was added to the saucer for uptake through the base of the pot, a process requiring approximately 2 h. After treatment, the soil water content was measured to be in the order of 35% (w/w).
For the first experiment, after application of the treatments, the pots contained 0 or approximately 50 mM NaCl which, over a period of 6 to 7 days was maintained at 0 or increased to 75 mM by evaporation and transpiration to reach the target gravimetric water content of 25% (w/w). For the second experiment, after treatment the soil contained 0 and approximately 70 mM NaCl and was maintained at 0 or increased to 100 mM NaCl over 6 to 7 days, as described above. Watering recommenced when the soil water content reached 25% (w/w) and NaCl in the soil solution was at the designated concentration. The method was designed to impose gradual salinity stress on the seedlings.
The LemnaTec Scanalyzer 3D (LemnaTec GmbH, Aachen, Germany) at The Plant Accelerator was used for daily image acquisition [20]. Two red-green-blue (RGB) images were captured from the side at a 90°rotation to each other and one RGB image from above. LemnaGrid software (LemnaTec GmbH, Aachen, Germany) was used for image processing to determine the projected shoot area-the sum of pixels identified as being part of the plant in each image [19,21].
To measure shoot growth during the plant osmotic response period, the projected shoot area from days 0 to 8 after salinity application was recorded. At this stage, plants were still young and the use of an exponential growth model is still appropriate. There is also a linear correlation between projected shoot area and shoot biomass, for Triticum monococcum [22] and for bread wheat (B. Berger, unpublished results), so that projected shoot area can be used directly to measure relative growth rates (RGRs). The projected shoot area was plotted over time and exponential growth curves were fitted to determine RGR [19,22]. The salinity tolerance index (STI) was determined: Over the time course of this experiment, this approximates to a measure of the tolerance of plants to the osmotic phase of salinity toxicity, the shoot ion accumulation independent tolerance (Roy et al., 2014).

Stress treatment and extraction of RNA for microarrays
Separate replicate pots of each cultivar were grown under conditions identical to those of the second experiment (0 and 100 mM NaCl) for destructive sampling of shoot material for RNA extraction at four separate time points. The leaf sheaths of seedlings of each cultivar were collected immediately before (day 0) or 1, 2, and 3 days after application of 0 or 100 mM NaCl. Samples from three replicate seedlings were taken at each time point under both control and salt stress conditions, and total RNA from the sheaths was isolated using Trizol-modified reagent [23], and then used for the preparation of Cyanin-3 (Cy3)-labelled cRNA probes. Dye incorporation and cRNA yield were checked with the NanoDrop ND-1000 Spectrometer. The Cy3-labelled cRNA was fragmented and hybridized to the Agilent wheat microarray (Design ID: 0222297) according to the manufacturer's guidelines. All microarray data have been deposited in the ArrayExpress database (https://www.ebi.ac.uk/arrayexpress/; accession number E-MTAB-3270).

Data mining of gene expression with Agilent wheat gene expression microarray
Microarray experiments were performed with the four bread wheat cultivars, Berkut, Krichauff, Gladius and Drysdale, using a 44K wheat gene expression array with a one-color gene expression microarray solution (Agilent Technologies). Each analysis was repeated in triplicate as biological replicates-a total of 84 microarray experiments were performed according to the manufacturer's manual. Genes displaying a signal value of >1000 were selected for analysis. Feature extraction and image analysis software (GeneSpring GX version 12.5; Agilent Technologies) was used to locate and delineate each spot in the array and to perform a statistical analysis of the intensity of each spot using the Lowess method (p < 0.01). The statistical significance of gene expression was tested using the ANOVA analysis of variance test combined with a Benjamini and Hochberg false-discovery-rate multiple correction algorithm with a corrected q value of <0.01. The criterion for the fold change threshold selection was a greater than fivefold change of expression. The changes in gene expression were calculated as follows: For expression analysis under control conditions: To identify functional descriptions of genes, EST or probe sequence data were searched against the following nucleotide datasets using the EST data of the NCBI (http://www.ncbi. nlm.nih.gov/), the cDNA data of rice derived from RAP-DB v.2 (http://rapdb.dna.affrc.go.jp/) and the cDNA data present in TAIR release 10 (ftp://ftp.arabidopsis.org/home/tair/Sequences/ blast_dataset/). The top-scoring hit for each query was selected for gene annotation with a threshold E-value of <1E-5, and the Load ID was used for the functional classification using the gene ontology (GO) of agriGO (http://bioinfo.cau.edu.cn/agriGO/), amiGO 2 (http:// amigo.geneontology.org/amigo) and the Gene Ontology (http://www.geneontology.org/). The percentage of identified genes was calculated as follows: percentage ð%Þ ¼ ðnumber of the genes classified by the GO termÞ ðtotal number of genes used for the classificationÞ Â 100 Gene expression profile charts under stress conditions were searched against the following datasets using the eFP Browser (http://bbc.botany.utoronto.ca/efp/cgi-bin/efpWeb.cgi) for Arabidopsis and the RiceXPro (http://ricexpro.dna.affrc.go.jp/) for rice.

Quantitative RT-PCR analysis
Five microgram of total RNA was used for cDNA synthesis with random hexamer primers and SuperScript III reverse transcriptase (Invitrogen). SYBR Premis Ex Taq (TaKaRa) and genespecific primers provided for the Primer Express Software v3.0.1 (Life technologies) were used for the reactions. Quantitative RT-PCR (qRT-PCR) analyses were performed using 7500 Fast Real-Time PCR systems (Applied Biosystems). The relative expression levels in each transcript were obtained by normalization to the actin gene (AB181991). The gene-specific primer sets were listed in S1 Table. Results Continuous monitoring of shoot biomass production reveals differences in growth responses to salinity between the bread wheat cultivars Four commercial bread wheat cultivars were investigated for variation in salinity tolerance: Berkut, Krichauff, Gladius and Drysdale. Previous studies indicate that the doubled-haploid (DH) mapping population from a cross between Berkut and Krichauff were analyzed as salinity tolerance (ST) quantitative trait loci (QTL) to improve the grain yield-independent ST in field [24,25]. A high-throughput phenotyping system was used to measure the growth of these bread wheat cultivars, and to quantify the effects of salinity on their growth. To evaluate the shoot ion independent component of salinity stress, the relative growth rate (RGR) of plants was estimated using the images of the projected shoot area for 8 days after treatment, and exponential growth curves fitted through the projected shoot area data [21,26,27]. For the relatively short time period of this experiment, the assumption of exponential growth was found to be a satisfactory approximation (Fig 1 and S3 Fig).
Over the time course of the experiments, Berkut showed the highest mean RGR under control conditions (Figs 1 and 2). Treatment with 75 mM and 100 mM NaCl caused little reduction in RGR for both Drysdale and Krichauff, and a greater reduction in RGR under salinity was observed for Gladius and Berkut (Fig 2 and S3E Fig). These differences can be best summarized by calculating the salinity tolerance index (STI: Table 1), which indicated that Krichauff had the highest salinity tolerance under those salinity stress conditions (strictly, the greatest ability to maintain RGR over this time period in saline conditions relative to nonsaline conditions).
Cluster analysis of variation in gene expression in leaf sheaths with respect to shoot growth and salinity stress tolerance Microarray profiling of gene transcripts from leaf sheaths was performed to identify genes responding to the early stages of salinity stress. We selected the treatment of wheat plants with 100 mM NaCl for microarray analysis because the treatment with 100 mM NaCl affected the further widening of the difference of STI in comparison with the treatment of wheat with 75 mM NaCl among four cultivars (Table 1) Although the patterns of leaf sheath transcriptome regulation varied sharply among the cultivars, an attempt was made to identify common responsive genes governing shoot biomass development and salinity tolerance. We focused first on RGRs under control conditions. Berkut showed the highest mean RGR under control conditions in both experiments (  . Cultivars were treated with RO water (open symbols) or 100 mM NaCl (filled symbols) at the emergence of the fourth leaf (vertical line indicates day 0 or day of treatment). Projected shoot area (pixels) was determined by image analysis from digital images taken with RGB cameras. Exponential curves were fitted to data from day 0 to day 7 after treatment (rectangle within graphs A to D) to calculate relative growth rates for (E) Berkut, (F) Krichauff, (G) Gladius and (H) Drysdale seedlings under control and saline conditions. doi:10.1371/journal.pone.0133322.g001  (Fig 3A).
The timing of expression of the 39 genes exhibited three tendencies. The first, shown as cluster I, comprises the early response that occurred within 24 h after control conditions ( Fig  3B). Cluster II comprises expressed genes in the second day under control conditions. Cluster III includes the late response that occurred on the 3 rd day. All 39 genes were searched against the annotated genome datasets of both Arabidopsis and rice (TAIR 10 and RAP-DB) to identify their function. A comparative gene annotation analysis of wheat with rice and Arabidopsis indicated that most homologous genes have known functions (Table 2). Functional category classification using the gene ontology (GO) term revealed that most of the 39 genes were assigned to cellular and metabolic processes, developmental processes, and cell organization and biogenesis, suggesting that Berkut has a specific ability to maintain a high growth rate in the early stages of growth by transcription of particular genes involved in plant growth ( Fig 3C).
Drysdale and Krichauff showed higher STI with the 100 mM NaCl treatment in comparison with Gladius and Berkut (Table 1), and they showed specific temporal variations in gene expression in response to salinity stress (Fig 4). Most upregulated and downregulated genes were expressed within 24 h of salinity treatment in both cultivars ( Fig 4A). Of the 47 genes in Drysdale and 96 genes in Krichauff that were upregulated under salinity stress compared with control (RO water) conditions, 42 (89.4%) and 91 (94.8%), respectively, were observed to be upregulated one day after treatment. However, only five genes were similarly regulated in both cultivars ( Fig 5A and Table 3). In contrast, the more salt sensitive Gladius and Berkut varieties did not exhibit dramatic changes in the numbers of upregulated genes during the three days after stress treatment, with only about 10 genes upregulated by mild salinity stress in each of the varieties. The same tendency was observed in the gene sets which were downregulated under mild salinity stress conditions (Figs 4 and 5A and S2 Table).
Krichauff consistently showed the highest STI compared with the other three varieties under moderate salinity stress conditions (Table 1). Transcriptome analysis showed that 87 genes were upregulated only in the leaf sheaths of Krichauff under salinity conditions (Fig 5A). The annotated nucleotide database of RAP-DB contained rice homologs of 48 of the 87 genes (Table 4). Using gene ontology analysis, most of the 48 genes were assigned as encoding proteins involved in cellular and metabolic processes (Fig 5B). The Krichauff-specific genes were accordingly classified as genes associated with the persistence of vigorous growth under mild salinity stress conditions. Drysdale showed second highest STI compared with the other three varieties (Table 1). GO analysis showed that rice homologs of 21 of the 39 genes which were upregulated only in Drysdale were also classified as genes associated with the cellular, metabolic and biological processes ( Fig 5C and Table 5). To identify differences between gene expression profiles associated with shoot biomass development under control and saline conditions, gene sets showing expression only in Berkut under control conditions (Table 2 and S3 Table), only in Krichauff under saline conditions (Table 4 and S4 Table) and only in Drysdale under saline conditions (Table 5 and S5 Table) were compared. However, the Venn diagram shows no overlap among the 39 genes upregulated in Berkut, the 87 genes upregulated in Krichauff and the 39 genes upregulated in Drysdale (Fig 6A). The 58 genes downregulated in Berkut, the 124 genes downregulated in Krichauff and the 4 genes downregulated in Drysdale also showed no overlap (Fig 6B).
Quantitative RT-PCR (qRT-PCR) was performed to validate the gene expression profiles obtained by microarray analysis. We randomly selected several genes which were expressed in  either control or 100 mM NaCl treated plants (Figs 3A and 5A). The scatter plots showed that the results of the qRT-PCR analysis correlate highly with those of the microarray analysis (S5 Fig). The determination coefficient (R 2 = 0.8418) indicates that the microarray analysis provides the reliable data to identify genes responding to the early stages of salinity stress.

Comparative analysis of gene expression profiles of bread wheat with those of Arabidopsis and rice revealed physiological differences in salinity stress responses
To determine whether the upregulated genes identified only in Krichauff and only in Drysdale under salinity treatment are also expressed in other plants, the expression profiles of the homologous genes in Arabidopsis and rice were examined using public gene expression databases (Tables 3-5). Using the Arabidopsis eFP Browser and Rice Expression Profile Database (RiceXPro), it was possible to determine the expression patterns of the 48 homologous of the upregulated Krichauff genes in Arabidopsis or rice after they had been exposed to either abiotic stress or hormone treatment. Only 12 rice and/or Arabidopsis homologs showed a tendency to respond to salinity stress and/or abscisic acid-inducible expression. In the 21 homologous genes upregulated only in Drysdale, the 7 rice and/or Arabidopsis homologs showed a tendency to response to stress-inducible expression. These results suggest that the gene expression  responses in commercial bread wheat under mild salinity stress may be different from those in Arabidopsis and rice under severe salinity stress, and that bread wheat varieties with high RGR and high salinity tolerance have a unique transcriptome pattern.

Discussion
Several genes resolved in our microarray results have been reported as regulators of stress tolerance in other plants. Dehydrin/responsive to ABA (RAB), the homolog of CK213509 and CV761504, is a well-known family of hydrophilic proteins, and strongly responsive to water stresses and ABA treatment [28]. The cell wall-associated receptor kinase (WAK) and WAKlike kinase (WAKL) gene family are involved in biotic and abiotic stress responses and are required for cell elongation and development [29,30]. The expression of WAKL4, the homolog of CJ851704, is induced by several mineral responses including to salt. This finding indicates that these membrane-to-cytoplasm communication signals are present in wheat under salinity stress as well as in other plants [31]. Lipid transfer protein (LTP) is involved in stress responses, and overexpression of LTP mediates salinity tolerance [32]. The LTP homolog AK332187 is expressed under salinity stress in stress-resistant wheat cultivars, suggesting a process of regulation similar to that for equilibrating metabolite accumulation. In addition, the results from the temporal microarray analysis performed in our study indicate that the experimental system properly captured the early response of gene expressions during the osmotic phase. Our analysis sheds light on the unique molecular mechanisms of gene expression in commercial bread wheat under mild salinity conditions. The results from the 3D PCA revealed a wide diversity in leaf sheath transcriptome profiles for the four commercial cultivars studied (S4 Fig). Our results indicate that the salt stress-inducible gene sets generated during the osmotic phase are distinct between the four cultivars (Fig 5). The comparative analysis of variation in gene expression in relation to physiological function revealed the causal components associated with development of biomass and tolerance to mild salinity stress (Fig 6). In addition, most genes expressed in the leaf sheath of the four bread wheat cultivars tested do not  correspond with salinity stress-inducible genes described in other well-characterized plants, such as rice and Arabidopsis (Tables 3-5). These results suggest that genetic variation resulting from selection during breeding exerts a strong influence on the development of genome architecture and adaptation to salinity conditions. An integrated analysis of the transcriptomes and epigenomes in two rice subspecies revealed a high correlation of allelic bias of epigenetic modification with gene expression [33]. Bioinformatics approaches have shown that cis-regulatory mutations are the major source of evolutionary innovations and that these alterations lead to changes in gene regulatory information and phenotypic defects [34]. A recent genome-wide association study identified sequence variants between cultivars and wild rice in promoter regions [35]. Comparative transcriptome analysis between domesticated and wild tomatoes revealed that thousands of gene expression patterns were different and that most of them were associated with environmental responses and stress tolerance [36]. These studies are consistent with the present study in suggesting that genomic variation affects salinity tolerance through changes in gene expression.
In particular, gene ontology analysis revealed that the genes involved in a number of cellular and metabolic processes were more upregulated for shoot biomass development than stress response genes (Fig 5), indicating that salinity-tolerant cultivars have the ability of improving biomass development with the alterations in metabolism-mediated gene expressions compared with stress tolerance under mild salinity conditions. For Australian bread wheat, improvements in salinity tolerance and the ability to continue to increase biomass under mild salinity stress conditions are particularly important. The breeding of salinity tolerant cultivars can exploit this study for knowledge of an understanding of sequence variants associated with the domestication of salinity tolerant Australian bread wheat.
In this study, we evaluated, under control and mild saline conditions in a greenhouse, the physiological differences between four commercial bread wheat cultivars. The high-throughput phenotyping system at The Plant Accelerator of the Australian Plant Phenomics Facility allowed the measurement through time of shoot biomass, thus enabling calculations of RGR  (Fig 1, S3 Fig and Table 1). This phenotyping analysis in The Plant Accelerator was more quantitative than is practically achievable in most other experiments in both controlled conditions and in the field. Therefore, the phenotyping system facilitates the analysis of useful germplasm using detailed molecular analyses, including microarrays. This study established that a highthroughput greenhouse phenotyping system, such as The Plant Accelerator, allows a powerful approach, combining physiological trait evaluation with molecular analysis of gene expression in bread wheat cultivars under mild saline conditions.
Our microarray analysis indicated that changes in the expression of salt stress-inducible genes reproducibly decreases within 2 days (Fig 4), suggesting that rapidly-induced gene expression, occurring shortly after the perception of stress, has a significant influence on adaptation to mild salinity conditions in the early osmotic phase. When combined with previous studies on Na + and Clions accumulation which gradually occurs under stress conditions, these results indicate that two phases, osmotic and ionic, are involved in responses and adaptation to mild salinity stress conditions in the bread wheat cultivars tested, and that the physiological characteristics of each phase are different from the perspective of acquiring salinity tolerance. Combined analysis of relative shoot growth, stress tolerance and Na + exclusion constitutes a detailed analysis of mechanisms underlying salinity tolerance in both the osmotic and ionic phases during mild salinity treatment.
In conclusion, our integrated systematic analyses of phenomics and transcriptomics of various cultivars can provide useful information for a better understanding of the physiological and molecular processes involved in mild salinity tolerance in Australian bread wheat cultivars. The high-throughput phenotyping system at The Plant Accelerator of the Australian Plant Phenomics Facility has permitted a systematic analysis of the transcriptome profiles of bread wheat plants under well-controlled conditions and shed light on mild salinity-tolerance mechanisms represented by gene expression in leaf sheath during the early "osmotic" phase of salinity response and adaptation.
Supporting Information S1 Table. Primer sets for quantitative RT-PCR analysis. (XLSX) S2 Table. Differentially downregulated genes during salinity stress in Drysdale and Krichauff when compared to Berkut and Gladius. Listed are genes whose expression was reduced more than five-fold in Drysdale and Krichauff compared with the other cultivars under salinity stress conditions. p values of <0.01 were included. Asterisk indicates stress-and/ or abscisic acid treatment-inducible gene in rice and Arabidopsis. a Probe name as given by Agilent 44K wheat gene expression array (Design ID: 022297). b Locus indicates GenBank accession. c Load ID as given by RAP-DB. d AGI, Arabidopsis Genome Initiative. e Description as given by The Institute for Genomic Research database. (XLSX) S3 Table. Differentially downregulated genes during biomass accumulation under control conditions in Berkut. Listed are genes whose expression was reduced more than five-fold in Berkut compared with the other cultivars under control conditions. p values of <0.01 were included. Asterisk indicates stress-and/or abscisic acid treatment-inducible gene in rice and Arabidopsis. a Probe name as given by Agilent 44K wheat gene expression array (Design ID: 022297). b Locus indicates GenBank accession. c Load ID as given by RAP-DB. d AGI, Arabidopsis Genome Initiative. e Description as given by The Institute for Genomic Research database. (XLSX) S4 Table. Differentially downregulated genes during salinity stress in Krichauff. Listed are genes whose expression was reduced more than five-fold in Krichauff compared with the other cultivars under salinity stress conditions. p values of <0.01 were included. Asterisk indicates stress-and/or abscisic acid treatment-inducible gene in rice and Arabidopsis. a Probe name as given by Agilent 44K wheat gene expression array (Design ID: 022297). b Locus indicates Gen-Bank accession. c Load ID as given by RAP-DB. d AGI, Arabidopsis Genome Initiative. e Description as given by The Institute for Genomic Research database. (XLSX) S5 Table. Differentially downregulated genes during salinity stress in Drysdale. Listed are genes whose expression was reduced more than five-fold in Drysdale compared with the other cultivars under salinity stress conditions. p values of <0.01 were included. Asterisk indicates stress-and/or abscisic acid treatment-inducible gene in rice and Arabidopsis. a Probe name as given by Agilent 44K wheat gene expression array (Design ID: 022297). b Locus indicates Gen-Bank accession. c Load ID as given by RAP-DB. d