Next Article in Journal
Investigating Human Mitochondrial Genomes in Single Cells
Next Article in Special Issue
A New Intra-Specific and High-Resolution Genetic Map of Eggplant Based on a RIL Population, and Location of QTLs Related to Plant Anthocyanin Pigmentation and Seed Vigour
Previous Article in Journal
Sparsity-Penalized Stacked Denoising Autoencoders for Imputing Single-Cell RNA-seq Data
Previous Article in Special Issue
Molecular Analysis of the Official Algerian Olive Collection Highlighted a Hotspot of Biodiversity in the Central Mediterranean Basin
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 5
10.3390/genes11050533
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Natural Variation Uncovers Candidate Genes for Barley Spikelet Number and Grain Yield under Drought Stress

by
Samar G. Thabet
1,
Yasser S. Moursi
1,
Mohamed A. Karam
1,
Andreas Börner
2 and
Ahmad M. Alqudah
2,*
1
Department of Botany, Faculty of Science, University of Fayoum, Fayoum 63514, Egypt
2
Research Group Resources Genetics and Reproduction, Department Genebank, Leibniz Institute of Plant Genetics and Crop Plant Research, 06466 Seeland OT Gatersleben, Germany
*
Author to whom correspondence should be addressed.
Genes 2020, 11(5), 533; https://0-doi-org.brum.beds.ac.uk/10.3390/genes11050533
Submission received: 22 March 2020 / Revised: 29 April 2020 / Accepted: 5 May 2020 / Published: 11 May 2020
(This article belongs to the Special Issue Genetic Diversity Assessment and Marker-Assisted Selection in Crops)
Browse Figures
Review Reports Versions Notes

Abstract

:
Drought stress can occur at any growth stage and can affect crop productivity, which can result in large yield losses all over the world. In this respect, understanding the genetic architecture of agronomic traits under drought stress is essential for increasing crop yield potential and harvest. Barley is considered the most abiotic stress-tolerant cereal, particularly with respect to drought. In the present study, worldwide spring barley accessions were exposed to drought stress beginning from the early reproductive stage with 35% field capacity under field conditions. Drought stress had significantly reduced the agronomic and yield-related traits such as spike length, awn length, spikelet per spike, grains per spike and thousand kernel weight. To unravel the genetic factors underlying drought tolerance at the early reproductive stage, genome-wide association scan (GWAS) was performed using 121 spring barley accessions and a 9K single nucleotide polymorphisms (SNPs) chip. A total number of 101 significant SNPs, distributed over all seven barley chromosomes, were found to be highly associated with the studied traits, of which five genomic regions were associated with candidate genes at chromosomes 2 and 3. On chromosome 2H, the region between 6469300693-647258342 bp includes two candidate drought-specific genes (HORVU2Hr1G091030 and HORVU2Hr1G091170), which are highly associated with spikelet and final grain number per spike under drought stress conditions. Interestingly, the gene expression profile shows that the candidate genes were highly expressed in spikelet, grain, spike and leaf organs, demonstrating their pivotal role in drought tolerance. To the best of our knowledge, we reported the first detailed study that used GWAS with bioinformatic analyses to define the causative alleles and putative candidate genes underlying grain yield-related traits under field drought conditions in diverse barley germplasm. The identified alleles and candidate genes represent valuable resources for future functional characterization towards the enhancement of barley cultivars for drought tolerance.

1. Introduction

Barley (Hordeum vulgare L.) ranks as one of the most important cereal crops worldwide. Globally, barley is the fourth most important cereal crop in terms of production after maize (Zea mays L.), rice (Oryza sativa L.), and wheat (Triticum spp.) (Faostat 2017, http://www.fao.org/faostat/en/#home). One limitation in achieving the production target is abiotic stress which limits the quality and nutritional value of the grain in cereal crops worldwide [1]. Among all abiotic stresses, drought is the most important environmental stress which limits crop production and yield [2,3,4], and is becoming more common particularly in the arid and semiarid regions [2].
Crops can be exposed to drought during their entire life cycle from vegetative to reproductive stages [5]. Drought stress affects crop growth and yield during all developmental stages [6]. Water shortage at early growth stages can cause severe problems for seedlings, restricting the emergence, growth and development of seedlings, and thus affecting grain yield [7]. Furthermore, the developing plants will have poor tillering capacity, leading to fewer tillers per unit area and thus lower yield potential. Moreover, drought in the period of stem elongation causes a decrease in the number of grains per unit area because it has a negative impact on floret formation and fertility [8].
At post-anthesis, water insufficiency reduces the grain filling rate and duration leading to shriveled grains [2]. Moreover, the effect of drought on yield is highly complex and involves processes as diverse as reproductive organs, gametogenesis, fertilization, embryogenesis and seed development [6,9]. Reproductive and seed development phases are especially sensitive to drought stress [2,10]. In barley, a reduction in the number of grains per spike, grain filling duration and dry matter accumulation have already been reported to decrease grain yield [2,10]. Several studies have reported that early growth stage parameters (e.g., tiller number, biomass formation, etc.) are highly correlated with yield potential and grain quality at harvest under both normal and drought stress conditions in various cereal crops, including barley [11,12]. Accordingly, understanding the genetic basis for drought tolerance in crop plants by identifying the genetic loci and the candidate genes associated with these traits is useful for developing new varieties with more drought-tolerant characters.
The Genome-Wide Association Scan (GWAS) approach is widely employed to reveal associations between genomic loci and advantageous traits in a given population based on linkage disequilibrium [13]. These loci then become targets for improving new genotypes by the breeder. The GWAS is very effective in identifying major candidate genes regulating mono- or oligogenic agricultural traits. Recently, GWAS has been successfully used to identify genes for yield-related traits [14,15]. Barley has a high-level population structure such as two-rowed and six-rowed cultivars, spring and winter barley [16]. This may lead to spurious marker–trait associations in GWAS [17]. Therefore, it is important to use strong statistical methods and strategies to control the population structure [17]. A mixed-linear model (MLM) approach has been developed to control spurious associations through account multiple levels of relatedness leading to better performance [18].
GWAS has successfully yielded genomic locations for quantitative trait loci (QTL) in crop cereals [13]. Identification of genomic locations for QTL using linked segregating markers is considered to be highly useful for marker-assisted breeding. Nevertheless, the ultimate goal of GWAS in crop species is to detect new QTL through genetic dissection of complex traits.
The present study aimed at identifying the genetic basis of drought tolerance at different developmental and growth phases in 121 spring barley accessions under field conditions using GWAS. In total, 101 SNPs showed association with different traits that were distributed across the seven chromosomes of barley. The identified QTL colocalized with several genes that are exclusively distributed on chromosomes 2H and 3H. The annotation and expression of these genes demonstrated their roles in drought tolerance.

2. Materials and Methods

2.1. Experimental Setup and Phenotyping

In total, 121 diverse spring barley accessions from different geographical origins were grown under the field conditions in the 2017/2018 growing season at the Experimental Station of University of Fayoum. The collection included 83 cultivars, 29 landraces, and 9 breeding lines. They originated from Europe (EU, 62), West Asia and North Africa (WANA, 24), East Asia (EA, 22), and the Americas (AM, 13). The row types were two-rowed (72) and six-rowed (49). The population structure and genetic diversity using the genotypic information of the accession in the collection are shown in Figure S1 that demonstrated there is no clear structure on our collection.
Five seeds from each accession were sown in plastic pots (40 cm × 26 cm × 26 cm) filled with field-soil on the 1st of December 2017 under field conditions. The soil texture was classified as clay-loam with pH = 8.0, total N% = 9.4 and available P = 58.0 ppm. Manual irrigation was performed as required, and 5 g (17:11:10/N:P:K) fertilizer was added to each pot. In field-grown plants, each accession was replicated in three pots of each treatment (control and drought). At the beginning of early reproductive (spikelet development phase) phase (~25 days after sowing), the plants were thinned into three plants per pot, standing with a border to eliminate positional and environmental effects on growth and development. Weeds were controlled manually.
Plants were irrigated until the onset of early reproductive, then the plants were then exposed to two watering treatments: (1) well-watered treatment (soil maintained at ~75% of field capacity (FC)); and (2) severe drought stress (at 35% FC). To maintain the targeted (~75% FC) and drought (~35% FC), eight randomly selected accessions were used as a reference. Before irrigation, eight reference pots were weighed and watered to adjust the corresponding field capacity, and the rest of the experiment was watered accordingly. Irrigation and drought treatment continued until maturity after that, irrigation withholds until harvest. Nine morphological, developmental and grain yield-related traits were measured from at least three biological replicates for each accession under each treatment. More information about the phenotypic trait measurements are explained in Table 1. Respective drought tolerance indices were calculated from the recorded data as described in (Table 1).
The average of temperature and humidity during the growing season of this experiment (2017–2018) was 16.62 °C and 55.16%, respectively. The maximum temperature was recorded in November and April (24.33 and 27.66 °C, respectively) with the minimum temperature from January to February (8.29 and 9.12 °C, respectively). The minimum humidity was 30.37% and 30.41% in November and February, respectively, whereas the maximum humidity was recorded in April (78.33%) (Figure S2).

2.2. Data Analysis

Analysis of variance [19] was conducted to compare the controlled and drought stress conditions at p < 0.05 for all measured traits using GENSTAT 18 [20]. Data were analyzed as a randomized complete block design (RCBD) with three replications. The treatments were considered as main plots and the accessions were considered as sub-plots. Broad-sense heritability (H2) for the measured traits under each condition separately was calculated using GENSTAT 18. The phenotypic data were subjected to Residual Maximum Likelihood (REML) to analyze it in mixed linear model (MLM). Mean estimation of each measured trait in each accession under each treatment was calculated as Best Linear Unbiased Estimates (BLUEs) using GENSTAT 18. Correlation matrix analysis among the traits in each treatment was separately calculated by GENSTAT 18. Comparison between treatments at p < 0.05 for each trait including boxplots were calculated using R-studio [21].

2.3. Genome-Wide Marker–Trait Associations

The accessions were genotyped by a 9K IlluminaTM SNPs chip. In the analysis, we only used the markers which passed the quality control as minor allele frequency (MAF) ≥ 0.05 with their physical positions. A mixed linear model (MLM) was performed to determine marker–trait associations between the estimated phenotypic traits (BLUEs) and genotypic data. Different statistical models, e.g., general-linear model (GLM), mixed linear model (MLM), and compressed MLM (CMLM)) were tested in GWAS using GAPIT R package [22]. Finally, we used MLM as a powerful model considering the population structure, including kinship and PCA, to control the population structure influence. False discovery rate (FDR) at 0.001 was calculated for each trait under each treatment separately and association signals passed the threshold of FDR at 0.001 (−log10 p-values ≥ FDR) were used in further analyses. To be sure that our associations were true, we followed the GWAS and post-GWAS protocol published recently [13].

2.4. SNP-Gene-Based Association and Haplotype Analysis

At each chromosome, linkage disequilibrium (LD, r2) among the significant SNPs within highly associated genomic region was calculated and presented as heatmap plot. This allowed us to define the most important physical position that had been used for candidate gene identification. The physical positions of SNPs exceeding FDR within the linkage disequilibrium interval were used in annotation for high-confidence (HC) candidate gene with other respective information using the barley genome explorer web-based with recent barley genome dataset (BARLEX; http://apex.ipk-gatersleben.de).
SNPs within the candidate gene physical position were used for further validation of SNP-Gene-based haplotype analyses and expression analyses. T-test at p < 0.05 was used to calculate the significant differences between alleles on the associated trait(s) [13]. RNA-seq datasets were derived from 16 different tissues of barley cv. ‘Morex’ cultivar, each with three biological replicates. In total 48 samples were used for generation of RNA-seq data.
From seven vegetative, six inflorescence, two developing grain and one germinating grain tissues, more details about RNA-seq experiments was published by Mascher et al. [23]. We used BARLEX; an expression database for barley that presented as FPPM (fragments per kilobase per million mapped reads).

3. Results

3.1. Phenotypic Characteristics and Natural Variation

In total, nine traits were recorded under control and drought treatments. Additionally, respective drought tolerance indices were calculated and used as a derived trait for GWAS. A wide range of phenotypic variation with normal distribution was detected for all traits (Figure 1 and Figure S3). The means of most of the studied traits showed a significant reduction under drought treatment compared to control conditions (Table 2 and Figure S3). There were no significant differences between the treatments in phase transition i.e., AT, SH and AE developmental stages (Table 2 and Figure S3). Drought treatment influenced significantly other developmental and yield traits such as plant height and spikelet number per spike (Table 2 and Figure S3). Notably, the genotypes showed wide range of variation in all drought tolerance indices (Table 2, Figure 2 and Figure S4).
Furthermore, heritability values were relatively high under drought, ranging from 0.64 for TGW to 0.84 for AT, whereas they ranged from 0.72 to 0.80 for AT and SH, respectively under control conditions. Additionally, the heritability values varied for tolerance indices from 0.68 for AE_DTI and NGS_DTI to 0.78 for AT_DTI (Table 2).

3.2. Correlations Analysis

Significant positive correlations were observed among various traits within both treatments. For example, AE showed a significantly high positive correlation with AT and SH (r = 0.98 *** and r = 1.00 ***, respectively) under control conditions and (r = 0.93 *** and r = 1.00 ***, respectively) under drought treatment (Figure 3A,B and Figure S5). On the contrary, some traits showed high significant negative correlations under both conditions. Interestingly, TGW showed negative correlation with NGS under control and drought conditions (r = −0.21 and −0.44 ***, respectively; Figure 3A,B).
For DTI, AT_DTI showed high significant positive correlation with SH_DTI and AE_DTI (r = 0.92 *** and r = 0.91 ***, respectively). Additionally, SH_DTI exhibited high positive correlation with AE_DTI (r = 0.99 ***). The SL_DTI with AL_DTI showed a significant negative correlation (r = −0.27 ***; Figure S5).

3.3. Natural Genetic Variation and Candidate Genes Potentially Underlying Drought Tolerance

GWAS analysis of 121 selected accessions was performed to find out the natural genetic variation of the studied traits. We detected a total number of 101 significant marker–trait associations (with −log10 p-value ≥ 3) distributed over the seven barley chromosomes (Figure 4 and Table S1). There was plenty of natural genetic variation of all studied traits in this collection. Through GWAS analysis, we found twelve interesting genomic regions for genetic variation of all studied traits distributed on only two chromosomes (2H and 3H).
Highly significant LD was found among significant SNPs within these genomic regions (Figure 4b), indicating that these significant SNPs are potentially harboring important candidate genes in addition to their useful for marker-assisted selection. On chromosome 2H, two adaptive genes were identified, whereas on chromosome 3H, ten genes that represent a combination of both adaptive and constitutive genes were identified. Adaptive genes might be control-specific, i.e., genes that regulate trait variation under control only, or drought-specific (genes that regulate trait variation under drought only). Constitutive genes regulate trait variation under both control and drought conditions (Table 3).
In total, eight SNPs were associated (−log10 p-value ≥ 3) with AT parameters. Out of these, five SNPs were adaptive: three were control-specific and two were drought-specific; the remaining three SNPs were constitutive. The constitutive SNPs were identified on chromosomes 1H, 3H and 4H. The most significant one (with −log10 p-value = 8.7) was observed on chromosome 3H at 126.69 cM. Only one constitutive gene was identified HORVU3Hr1G098200 (Chr. 3; 126.69 cM) (Table 3).
Sixteen SNPs were associated (−log10 p-value ≥ 3) with SH parameters. Out of these, seven SNPs were adaptive: four SNPs were control-specific and three drought-specific; the remaining nine SNPs were constitutive. The constitutive SNPs were identified on chromosomes 1H, 2H, 3H, 4H, 5H and 6H where the most significant one (with −log10 p-value = 4.7) was observed on chromosome 6H at 53.54 cM. Only seven SNPs showed association with candidate genes. Four genes were control-specific; HORVU3Hr1G088300, HORVU3Hr1G089160, HORVU3Hr1G089080 and HORVU3Hr1G098200. Three constitutive genes were identified: HORVU3Hr1G098200, HORVU3Hr1G116790 and HORVU3Hr1G115810 (Table 3).
In total, 18 SNPs showed association (with −log10 p-value ≥ 3) with AE parameters. Ten SNPs were adaptive: six were control-specific and four were drought-specific; the remaining SNPs were constitutive. These constitutive SNPs were mapped on chromosomes 2H, 3H, 4H, 5H and 6H. The most significant one (with −log10 p-value = 4.6) was detected on chromosomes 6H at 53.54 cM. Out of these, twelve SNPs showed association with candidate genes. Eight genes were control-specific: HORVU3Hr1G018650, HORVU3Hr1G020430, HORVU3Hr1G020660, HORVU3Hr1G019590, HORVU3Hr1G088300, HORVU3Hr1G089160, HORVU3Hr1G089080 and HORVU3Hr1G098200 and four genes are constitutive: HORVU3Hr1G116790, HORVU3Hr1G098200 and HORVU3Hr1G115810 (Table 3).
For the trait PH, six SNPs were associated. Out of these, five SNPs were adaptive: four SNPs were control-specific and one SNP was drought-specific. Only one constitutive SNP was found on chromosome 7H at 131.59 cM. There were no SNPs associated with candidate genes for PH (Table 3).
GWAS analysis showed that six SNPs were associated with SL parameters. Only one SNP was associated with a candidate gene that was control-specific, HORVU3Hr1G098200 (Table 3).
In total, eight SNPs showed association (p-value ≤ 0.001) with AL. Out of these, six SNPs were adaptive: two were control-specific and four were drought-specific; the remaining two SNPs were constitutive and mapped on chromosomes 2H and 3H. Six SNPs were associated with candidate genes. Only one gene was control-specific, HORVU3Hr1G098200; and one was constitutive, HORVU3Hr1G098200 (Table 3).
For NSS, 15 SNPs showed significant association with −log10 p-value ≥ 3. Six of them were associated with adaptive candidate genes. Three SNPs were control-specific: HORVU3Hr1G018650, HORVU3Hr1G020430 and HORVU3Hr1G020660 and another three were drought-specific: HORVU2Hr1G091030, HORVU2Hr1G091170 and HORVU3Hr1G019590. No constitutive genes were identified (Table 3).
In total, 15 SNPs were associated with NGS parameters. Seven SNPs showed association with candidate genes. Five genes were control-specific: HORVU2Hr1G091030, HORVU3Hr1G018650, HORVU3Hr1G020430, HORVU3Hr1G020660 and HORVU3Hr1G019590 and two drought-specific genes: HORVU2Hr1G091030 and HORVU2Hr1G091170 (Table 3).
For TGW parameters, seven SNPs showed significant association with −log10 p-value ≥ 3. Only two SNPs showed association with candidate genes. One was control-specific, HORVU3Hr1G098200, and the other drought-specific one was HORVU3Hr1G098200 (Table 3).
Overall, all the identified genes revealed a pleiotropic effect, i.e., each gene controlled more than one trait. The gene HORVU3Hr1G098200, for instance, regulates the variation of ten traits (Table 3). On the contrary, they differ in their mode of action, and some of them are adaptive genes such as HORVU2Hr1G091170 and HORVU3Hr1G018650. The first gene modulates the variation of NGS and NSS under drought; whereas the second controls them under control. Other genes are constitutive, such as HORVU3Hr1G116790 and HORVU3Hr1G115810 (Table 3). Surprisingly, the gene HORVU3Hr1G098200 showed a constitutive/adaptive mode of action. It controls the variation of (TGW) constitutively, while regulating the variation of (AE, AL, SH and SH) in an adaptive manner, i.e., under control only.

3.4. SNP-Gene-Based Analysis

Twelve SNPs at chromosomes 2H and 3H were physically co-located inside the candidate genes (Table 3). Two SNPs at 2H (SCRI_RS_166540 and SCRI_RS_157347) were detected within the genes HORVU2Hr1G091030 and HORVU2Hr1G091170, respectively. Meanwhile, ten SNPs at 3H were co-located within the physical position of candidate genes. For example, SNP numbers BOPA1_2391-566 and SCRI_RS_177313 were located within HORVU3Hr1G019590 and HORVU3Hr1G098200 genes, respectively. Interestingly, these SNPs were mostly associated with NGS, NSS, and TKW under drought stress. The rest of the genes at 3H were associated with traits under drought stress or with DTI for SH and AE.
The allelic analysis of the SNPs that are associated with traits under drought showed that alleles A, G and A from markers SCRI_RS_166540, SCRI_RS_157347 and BOPA1_2391-566, respectively, have a highly significant impact on NSS (Figure 5a and Figure 6). The genes HORVU2Hr1G091030 and HORVU2Hr1G091170 at 2H controlled NGS under drought via markers SCRI_RS_166540 and SCRI_RS_157347, where alleles G and A, respectively, increased NGS significantly (Figure 5b and Figure 6). Only one marker (SCRI_RS_177313 from HORVU3Hr1G098200 gene) showed a significant effect on TKW by allele A that increased the value under drought (Figure 5b).

3.5. Expression Analysis of Candidate Genes

The expression analysis of candidate genes in different organs showed a wide range of expression for the genes (Figure 6). Notably, the associated genes with spikelet and grain number per spike under drought stress (Figure 5) showed a high expression for most of the organs (Figure 6). Gene HORVU2Hr1G091170 at 2H in particular that had high impact on spikelet and grain numbers under drought showed highest expression in the respective grain organs, e.g., lodicule and rachis in addition to spike at 1–1.5 cm length (Figure 6). The second highest expressed gene from the highly associated ones was HORVU2Hr1G091030. This gene was highly expressed in developing grain and lodicule (Figure 6). The expression of these genes demonstrated their biological roles in the spikelet and grain development under drought conditions. Other genes, e.g., HORVU3Hr1G020660 and HORVU3Hr1G018650, showed high expression particularly in senescing leaves, suggesting their roles in leaf development (Figure 6).

4. Discussion

Studying drought stress tolerance under field conditions in cereals such as barley is very limited, as it requires complex and laborious experiments for population characterization, in addition to being influenced by environmental factors [3,4]. Nevertheless, the present study focused on investigating the natural variation in diverse spring barley collections and on identifying candidate genes associated with the traits of interest under field conditions.
In the present study, there was a considerable reduction in most traits under drought stress compared to control conditions. These results indicated that drought stress reduced grain yield by decreasing NSS, NGS, and TGW. Our results are supported by the findings of [2,24], who examined the response of spring barley to pre- and post-anthesis drought and reported a yield reduction due to pre-anthesis water deficit on several fertile spikes and grains per plant. Our results are in agreement with those of Samarah, Alqudah, Amayreh and McAndrews [2], who found that drought stress reduced grain yield by reducing the number of tillers, spikes and grains per plant and individual grain weight in barley (Hordeum vulgare L.) as a result of early maturity and shortened grain filling duration at 25% field capacity compared to control. In conclusion, drought stress negatively influenced barley yield through impairing grain development, size and grain filling duration.
Agronomic traits such as grain yield and its components (NSS, NGS, and TGW) are the major selection criteria for drought tolerance in barley breeding [25]. Therefore, understanding the interplay among these parameters is of high importance. The correlation between TGW and NGS was always negative and significant under control and drought, respectively. The positive and significant correlation between TGW and SL indicates that yield mainly depends on spike length. Zhou et al. [26] suggested that the grain yield traits interacting with each other, increase in one of them (e.g., TGW) can be correlated with a reduction in another (e.g., NGS), which is in agreement with our results. In support of our findings, in a European spring barley collection [27], it was found that NGS showed negative correlations with all other yield parameters except grain length, concluding that the yield was mainly dependent on grain size and SL rather than NGS. Several reports showed that drought-tolerant genotypes implement high productivity under both stressed and unstressed conditions [2,28]. Therefore, the comparative analysis of the yield components under drought and well-watered conditions can be used for predicting stress tolerance of genotypes, and then in selection of more tolerant barley lines for future breeding purposes [24,29].

The Role of Putative Candidate Genes in Drought Tolerance

Promising genomic regions are located at position 75.56 cM on chromosome 2H, harboring two important candidate genes. The first one is HORVU2Hr1G091030, which encodes RNA polymerase II C-terminal domain phosphatase-like 1 (CPL1), for NGS_C, NGS_D, and NSS_D. CPL1 is a negative regulator of stress-responsive gene transcription, ABA, and stress responses [30]. In Arabidopsis, CPL1 regulates gene expression under various osmotic stresses through ABA signaling [30]. Loss of CPL1 function in the mutants enhances tolerance to oxidative stress including drought and salt stresses [31]. In rice, OsCPL1 is expressed in the young spikelets [32]. Most likely, this gene is expressed during the development of barley spikelets under drought. Additionally, it controls the NGS under drought and control conditions in a constitutive manner. The second candidate, HORVU2Hr1G091170, encodes expansin B3 for NSS_D and NGS_D. Expansins have been implicated in the responses of various plant species to water stress. In maize, increased expansin activity was found to be involved in maintaining the growth of primary roots at a low water potential [33]. The expression of a root β-expansin gene, GmEXPB2, is remarkably participated in root system architecture responses to several abiotic stresses, such as Fe, P, and water deficiency [34]. RhEXPA4 is a rose expansin gene that is up-regulated in rose petals after dehydration [35], and it confers salt and drought tolerance to transgenic Arabidopsis [36]. In potato, most of expansin-like B genes have a potential role in multistress tolerance and upregulated under stresses including drought [37]. These two genes—HORVU2Hr1G091030 and HORVU2Hr1G091170—showed different expression profiles (Figure 6). The noisy expression profile of HORVU2Hr1G091170, suggesting a stress-responsive gene, exclusively regulates the variation of NGS_D and NSS_D under drought (Table 3). While the slightly constant expression profile of HORVU2Hr1G091030 indicates a constitutive gene that regulates the variation of NGS under control and drought. These findings are in accordance with the results of several authors who found that the stress-responsive genes exhibiting noisy expression profiles, while the constitutive ones showing constant expression patterns reviewed in Lopez-Maury et al. [38].
Notably, the allelic diversity analysis shows that allele A and G from HORVU2Hr1G091030 and HORVU2Hr1G091170 genes, respectively, have a significant positive impact on NSS and then NGS under drought stress. The expression of these genes during spike, spikelet and grain developmental stages demonstrates their influence on agronomic traits under drought stress conditions. In terms of molecular breeding, the above-mentioned alleles were the highest alleles explained the natural variation under drought stress suggesting their usefulness in breeding programs. Taken together, this provides evidence that these genes are drought-specific and involved in the drought stress-tolerance pathway. These findings indicate that both genes are of high importance for enhancing barley grain yield under drought stress.
Interestingly, ten candidate genes were detected on chromosome 3H. Out of these, three genes at positions 44.26 and 45.55 cM were identified as candidates for various traits such as AE, NGS, and NSS under control conditions. The first one is HORVU3Hr1G018650 encoding pyruvate decarboxylase-2 (PDC2), which belongs to pyruvate decarboxylase (PDC) gene family. Pyruvate has been involved in the ethanolic fermentation pathway that is associated with flooding tolerance when plant cells switch from respiration to anaerobic fermentation [39]. Additionally, fermentation has important functions in the presence of oxygen, mainly in germinating pollen and during abiotic stress. This indicates the interdependency between pyruvate decarboxylase (PDC) and AE, NGS, and NSS under control conditions. Furthermore, PDC, which catalyzes the first step in this pathway, is thought to be a main regulatory enzyme [40]. In Arabidopsis, the expression of PDC genes during abiotic stresses has been reported [41]. In maize and Arabidopsis, strong induction of fermentation genes takes place in anaerobic conditions [42]. Thus, it is conceivable that ethanolic fermentation is part of a general response to environmental stress, e.g., drought stress.
For traits such as AE_C, NGS_C, and NSS_D at position 46.25 cM, HORVU3Hr1G019590 encodes myb domain protein 37. MYB (myeloblastosis) has a regulatory role in ABA signaling by activating some stress-inducible genes [43]. In Arabidopsis, MYB, namely AtMYB60, AtMYB44, and AtMYB15, have been implicated in the regulation of stomatal closure and ABA-mediated response to drought and salt stresses [44]. Agarwal et al. [45] detected that AtMYB15 was expressed in both vegetative and reproductive organs and up-regulated by cold and salt stresses. The differential expression of numerous MYB TFs in the Triticeae was shown to be involved in the response of abiotic stress conditions such as drought and salt stresses [46,47]. These findings suggest that the genes in question are constitutive genes and may have a role in drought tolerance in barley during heading and maturation.
For SH_C and AE_C, HORVU3Hr1G089160 encodes AP2-like ethylene-responsive transcription factor at position 104.32 cM. The AP2/ERF superfamily regulates diverse developmental responses such as flower pedicel abscission [48], leaf senescence [49], cell proliferation and shoot branching [50]. Houston et al. [51] observed that mutants of HvAP2 internode lead to the reduction of elongation in both the culm and the spike in barley, suggesting that the HvAP2 alleles increase grain yield by controlling spike density. In the current study, the allelic variation was only observed under control condition indicating that this gene might not be involved directly in drought tolerance.
Additionally, HORVU3Hr1G116790 encodes Aquaporin-like superfamily protein which is a candidate for SH_DTI and AE_DTI. Aquaporins (AQPs) are a class of water channel proteins that belong to the major intrinsic protein (MIP) superfamily of membrane proteins [52]. These proteins regulate the movement of water and other small molecules across plant vacuolar and plasma membranes [53]. Aquaporins have also been suggested to have an essential role in plant tolerance of biotic and abiotic stresses [54], and extension growth [52]. Furthermore, it was reported in [55] that various aquaporin homologs are involved in plant stress responses against a variety of environmental stresses that disturb plant cell osmotic balance and nutrient homeostasis.
The most effective gene HORVU3Hr1G098200 was mapped at 126.69 cM. This gene orchestrates the variation of ten traits both in constitutive and adaptive manner (Table 3). This explains the significant correlations between these traits, either under control or drought (Table 3 and Figure 3A,B). Shi et al. [56] reported that Chromosome 3B harbors genes that may be significant in controlling agronomical important traits, such as yield and resistance to biotic and abiotic stress in wheat. The QTL for these traits maps quite close to semidwarf1 (sdw1) on chromosome 3H at 126.69 cM. The semi dwarf 1 (sdw1) gene has previously been found to control the most desired agronomic traits barley reviewed by Hedden [57]. The pleiotropic effect of sdw1 gene was evidenced in wheat [58,59].
Additionally, the gene HORVU3Hr1G098200 showed an interplay between the constitutive and adaptive control pattern. For example, it constitutively controls the variation of thousand grain weight (TGW) under both control and drought stress. At the same time, it controls the spike heading only under control in an adaptive manner (Table 3). This pattern indicates that HORVU3Hr1G098200 is partially constitutive and partially stress-responsive gene. Therefore, we considered it an adaptive gene when it controlled a trait(s) under control or drought. On the other hand, we considered it constitutive when it controlled trait(s) under both conditions. This gene, nevertheless, exhibited constant expression profile in different plant tissues, spanning the reproductive period from anther extrusion until the seed set, suggesting a key role in controlling different traits and more likely to be constitutive (Figure 6). This finding is consistent with the results of [60], who found that cells express some genes constantly to maintain the concentration of some proteins tuned with the cell physiological needs. The constant expression pattern characterizes the constitutive genes rather than the stress responsive (adaptive) ones (reviewed in Lopez-Maury, Marguerat and Bahler [38]. According to Blum [61], drought stress when expressed as a final yield is affected by constitutive and adaptive plant traits (genes). These constitutive QTL/genes represent an instrumental tool for selection as they show stability across different environments compared to the adaptive ones. Moreover, the selection for drought tolerance based on these QTL/genes does not require drought stress.
The last candidate gene is HORVU3Hr1G115810, which encodes Kinetochore protein spc25 and affects AE_DTI and SH_DTI. Kinetochore proteins may have a pivotal role for centromere and kinetochore functioning [62,63,64], and chromosome segregation mediating [65]. Specifically, kinetochore protein MIS12 is required for the co-orientation of sister kinetochores during meiosis I in maize [66]. NDC80 kinetochore protein serves as a contact point for chromosome–spindle interaction [67]. Interestingly, QTL at 3H 126 and 154 cM have previously been reported to be associated with grain number and yield in barley under drought stress conditions [4,68,69].
On chromosome 3H, three genes (HORVU3Hr1G020660, HORVU3Hr1G088300, and HORVU3Hr1G098200) are counterparts of the ortholog in Arabidopsis, AT2G36270. The corresponding genes are expressed during the reproductive stage in the different flower and seed organs (Figure 6), indicating their significance in flowering and seed set. Our results are in accordance with the findings of Klepikova et al. [70], who compared the spatiotemporal expression and stability of a lot of genes in 79 organs and developmental stages. Additionally, they found that AT2G36270 had been significantly expressed in senescent organs (leaves and silique); this is similar to expression profile of the gene HORVU3Hr1G020660 that highly expressed in senescent leave (Figure 6). Taking these findings together, the similarity of spatiotemporal expression patterns, as well as functions of these genes in both barley and Arabidopsis, suggests that they might be are homologous for AT2G36270.
The significant role of AT2G36270 (ABI5, AtABI5 and BZIP39) during different growth stages, as well as in drought tolerance, has been evidenced in several studies. Finkelstein et al. [71] reported that it encodes a member of the basic leucine zipper transcription factor family, involved in ABA-regulated gene expression during germination, seed development and subsequent reproductive stage. In particular, ABI5 regulates a set of the late embryogenesis-abundant genes during both seed and ABA-inducible vegetative gene expression in wild-type and abi5-1 plants [71]. Mittal et al. [72] reported that overexpression of AtABI5 in transgenic cotton (Gossypium hirsutum) showed resistance to the imposed drought stress through ROS scavenging and osmotic adjustment, enhancing photosynthesis, as well as traits of drought avoidance (bigger root and leaf systems) and tolerance (longer internode length and higher stem weight) leading to better establishment under water shortage. In rice (Oryza sativa), overexpression of OsbZIP46CA1 significantly increased tolerance to drought and osmotic stresses at flowering stage, and suggested that OsbZIP46 is a positive regulator of ABA signaling and drought stress tolerance by modulating many stress-related genes [73].
In summary, the present study showed the value of using field experiments to investigate natural phenotypic and genetic variation in a worldwide panel of barley accessions underlying agronomic traits such as grain yield and its components which can be exploited for crop improvement. Drought stress negatively influences most of the studied traits. In addition, we observed significant positive correlations among various traits within both control and drought treatments. Candidate genes associated with drought response were detected on two chromosomes, notably 2H and 3H. Interestingly, most the candidate genes are described to be involved in responses to abiotic stresses such as drought and salt. Interestingly, the genomic regions at 75 cM on 2H and 126.69 cM on 3H harbor three candidate gene HORVU2Hr1G091030 HORVU2Hr1G091170 and HORVU3Hr1G098200, which are highly associated with spikelet and final grain number per spike, suggesting the crucial role in controlling grain yield under drought conditions. The discovered SNPs and candidate genes for drought response will be helpful for breeding drought tolerant barley cultivars.

5. Conclusions

Conclusively, the present study showed that drought negatively affects the yield-related traits. Despite of the reduction in most traits under drought, the heritability estimates for all respective traits were high, indicating the potential of this collection to conduct a GWAS analysis looking for drought-controlling alleles/genes. Additionally, the present study revealed that combining GWAS and bioinformatics is a very instrumental approach to identify candidate genes even for polygenic traits such as the yield-related components. Our results confirmed that the yield-related components are under polygenetic control; under contrasting growth conditions (control and drought). The candidate genes exhibited different patterns of traits control; some genes were adaptive (HORVU2Hr1G091170), while other genes were constitutive (HORVU2Hr1G091030 and HORVU3Hr1G098200). The constitutive pleiotropic genes are of high importance to improve drought tolerance because they can be employed to improve several traits at a time with no need to test under drought. The causative genes showed different expression patterns; the constitutive genes showed constant expression profiles, while the adaptive genes showed a noisy expression profile. Most of the causative genes were expressed in spikelet organs (palea, lema and lodicules), as well as in grain, spike and leaf, indicating their potential role in drought tolerance, in particular, during the reproductive stage. To get more comprehensive and clear answers, a new detailed experiment is underway to study the gene expression of the candidate genes identified in this study especially, the gene HORVU3Hr1G098200 because it regulates the variation of ten traits, and because of its constitutive/adaptive mode of action.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/5/533/s1, Figure S1: Population structure based on their row type (two-rowed vs. six-rowed), geographical of origin; Figure S2: Weather data at Fayoum station that indicates Max-Min Temperature and Humidity; Figure S3: Boxplot analysis of variation of the traits; Figure S4: Boxplot analysis of drought tolerance index variation of the traits in barley genotypes under control and drought stress; Figure S5: Correlations matrix among drought tolerant index; Table S1: Summary of the most significant SNP associated with studied traits.

Author Contributions

Conceptualization, S.G.T., Y.S.M., M.A.K. and A.M.A.; methodology, S.G.T., Y.S.M., M.A.K. and A.M.A.; software, A.M.A.; validation, S.G.T., Y.S.M., M.A.K. and A.M.A.; formal analysis, A.M.A.; investigation, S.G.T., Y.S.M. and A.M.A.; resources, A.B.; data curation, S.G.T., Y.S.M. and A.M.A.; writing—original draft preparation, S.G.T.; writing—review and editing, S.G.T., Y.S.M. and A.M.A.; visualization, S.G.T., Y.S.M. and A.M.A.; supervision, Y.S.M., M.A.K. and A.M.A.; software, A.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank all who contributed to this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lobell, D.B.; Schlenker, W.; Costa-Roberts, J. Climate trends and global crop production since 1980. Science 2011, 333, 616–620. [Google Scholar] [CrossRef] [Green Version]
  2. Samarah, N.H.; Alqudah, A.M.; Amayreh, J.A.; McAndrews, G.M. The effect of late-terminal drought stress on yield components of four barley cultivars. J. Agron. Crop Sci. 2009, 195, 427–441. [Google Scholar] [CrossRef]
  3. Thabet, S.G.; Alqudah, A.M. Crops and drought. In eLS; John Wiley & Sons, Ltd., Ed.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2019; pp. 1–8. [Google Scholar] [CrossRef]
  4. Sallam, A.; Alqudah, A.M.; Dawood, M.F.A.; Baenziger, P.S.; Borner, A. Drought stress tolerance in wheat and barley: Advances in physiology, breeding and genetics research. Int. J. Mol. Sci. 2019, 20, 3137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Shavrukov, Y.; Kurishbayev, A.; Jatayev, S.; Shvidchenko, V.; Zotova, L.; Koekemoer, F.; de Groot, S.; Soole, K.; Langridge, P. Early flowering as a drought escape mechanism in plants: How can it aid wheat production? Front. Plant Sci. 2017, 8, 1950. [Google Scholar] [CrossRef] [PubMed]
  6. Alqudah, A.M.; Samarah, N.H.; Mullen, R.E. Drought stress effect on crop pollination, seed set, yield and quality. In Alternative Farming Systems, Biotechnology, Drought Stress and Ecological Fertilization; Springer: Dordrecht, The Netherlands, 2011; pp. 193–213. [Google Scholar] [CrossRef]
  7. Thabet, S.G.; Moursi, Y.S.; Karam, M.A.; Graner, A.; Alqudah, A.M. Genetic basis of drought tolerance during seed germination in barley. PLoS ONE 2018, 13, e0206682. [Google Scholar] [CrossRef]
  8. Fischer, R.; Turner, N.C. Plant productivity in the arid and semiarid zones. Annu. Rev. Plant Physiol. 1978, 29, 277–317. [Google Scholar] [CrossRef]
  9. Barnabas, B.; Jager, K.; Feher, A. The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 2008, 31, 11–38. [Google Scholar] [CrossRef]
  10. Samarah, N.; Alqudah, A. Effects of late-terminal drought stress on seed germination and vigor of barley (Hordeum vulgareL.). Arch. Agron. Soil Sci. 2011, 57, 27–32. [Google Scholar] [CrossRef]
  11. Botwright, T.; Condon, A.; Rebetzke, G.; Richards, R. Field evaluation of early vigour for genetic improvement of grain yield in wheat. Aust. J. Agric. Res. 2002, 53, 1137–1145. [Google Scholar] [CrossRef]
  12. Cowley, R.; Luckett, D.; Moroni, J.; Diffey, S. Use of remote sensing to determine the relationship of early vigour to grain yield in canola (Brassica napus L.) germplasm. Crop Pasture Sci. 2014, 65, 1288–1299. [Google Scholar] [CrossRef]
  13. Alqudah, A.M.; Sallam, A.; Stephen Baenziger, P.; Borner, A. GWAS: Fast-forwarding gene identification and characterization in temperate Cereals: Lessons from Barley—A review. J. Adv. Res. 2020, 22, 119–135. [Google Scholar] [CrossRef] [PubMed]
  14. Reig-Valiente, J.L.; Marques, L.; Talon, M.; Domingo, C. Genome-wide association study of agronomic traits in rice cultivated in temperate regions. BMC Genom. 2018, 19, 706. [Google Scholar] [CrossRef] [PubMed]
  15. Mazaheri, M.; Heckwolf, M.; Vaillancourt, B.; Gage, J.L.; Burdo, B.; Heckwolf, S.; Barry, K.; Lipzen, A.; Ribeiro, C.B.; Kono, T.J.Y.; et al. Genome-wide association analysis of stalk biomass and anatomical traits in maize. BMC Plant Biol. 2019, 19, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Teulat, B.; Merah, O.; Sirault, X.; Borries, C.; Waugh, R.; This, D. QTLs for grain carbon isotope discrimination in field-grown barley. Theor. Appl. Genet. 2002, 106, 118–126. [Google Scholar] [CrossRef] [PubMed]
  17. Hou, S.; Zhu, G.; Li, Y.; Li, W.; Fu, J.; Niu, E.; Li, L.; Zhang, D.; Guo, W. Genome-wide association studies reveal genetic variation and candidate genes of drought stress related traits in cotton (Gossypium hirsutum L.). Front. Plant Sci. 2018, 9, 1276. [Google Scholar] [CrossRef] [PubMed]
  18. De Ronde, J.A.; Cress, W.A.; Kruger, G.H.; Strasser, R.J.; Van Staden, J. Photosynthetic response of transgenic soybean plants, containing an Arabidopsis P5CR gene, during heat and drought stress. J. Plant Physiol. 2004, 161, 1211–1224. [Google Scholar] [CrossRef]
  19. Haberle, J.; Svoboda, P.; Raimanová, I. The effect of post-anthesis water supply on grain nitrogen concentration and grain nitrogen šeld of winter wheat. Plant Soil Environ. 2018, 54, 304–312. [Google Scholar] [CrossRef] [Green Version]
  20. VSN-International. VSN GenStat for Windows, 18th ed.; VSN International: Hemel Hempstead, UK, 2016; Available online: https://genstat.kb.vsni.co.uk/knowledge-base/hcitegen/ (accessed on 15 August 2019).
  21. RStudio-Team. RStudio: Integrated Development for R; RStudio, Inc.: Boston, MA, USA, 2015; Available online: http://www.rstudio.com/ (accessed on 15 August 2019).
  22. Lipka, A.E.; Tian, F.; Wang, Q.; Peiffer, J.; Li, M.; Bradbury, P.J.; Gore, M.A.; Buckler, E.S.; Zhang, Z. GAPIT: Genome association and prediction integrated tool. Bioinformatics 2012, 28, 2397–2399. [Google Scholar] [CrossRef] [Green Version]
  23. Mascher, M.; Gundlach, H.; Himmelbach, A.; Beier, S.; Twardziok, S.O.; Wicker, T.; Radchuk, V.; Dockter, C.; Hedley, P.E.; Russell, J.; et al. A chromosome conformation capture ordered sequence of the barley genome. Nature 2017, 544, 427–433. [Google Scholar] [CrossRef] [Green Version]
  24. Al-Ajlouni, Z.; Al-Abdallat, A.; Al-Ghzawi, A.; Ayad, J.; Abu Elenein, J.; Al-Quraan, N.; Baenziger, P. Impact of pre-anthesis water deficit on yield and yield components in barley (Hordeum vulgare L.) plants grown under controlled conditions. Agronomy 2016, 6, 33. [Google Scholar] [CrossRef] [Green Version]
  25. Hossain, A.; Teixeira da Silva, J.A.; Lozovskaya, M.V.; Zvolinsky, V.P. High temperature combined with drought affect rainfed spring wheat and barley in South-Eastern Russia: I. Phenology and growth. Saudi J. Biol. Sci. 2012, 19, 473–487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zhou, H.; Liu, S.; Liu, Y.; Liu, Y.; You, J.; Deng, M.; Ma, J.; Chen, G.; Wei, Y.; Liu, C.; et al. Mapping and validation of major quantitative trait loci for kernel length in wild barley (Hordeum vulgare ssp. spontaneum). BMC Genet. 2016, 17, 130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Xu, X.; Sharma, R.; Tondelli, A.; Russell, J.; Comadran, J.; Schnaithmann, F.; Pillen, K.; Kilian, B.; Cattivelli, L.; Thomas, W.T.B.; et al. Genome-wide association analysis of grain yield-associated traits in a pan-european barley cultivar collection. Plant Genome 2018, 11, 0073. [Google Scholar] [CrossRef]
  28. Sharafi, S.; Ghassemi-Golezani, K.; Mohammadi, S.; Lak, S.; Sorkhy, B. Evaluation of drought tolerance and yield potential in winter barley (Hordeum vulgare) genotypes. J. Food Agric. Environ. 2011, 9, 419–422. [Google Scholar]
  29. Perlikowski, D.; Kosmala, A.; Rapacz, M.; Koscielniak, J.; Pawlowicz, I.; Zwierzykowski, Z. Influence of short-term drought conditions and subsequent re-watering on the physiology and proteome of Lolium multiflorum/Festuca arundinacea introgression forms, with contrasting levels of tolerance to long-term drought. Plant Biol. Stuttg 2014, 16, 385–394. [Google Scholar] [CrossRef] [PubMed]
  30. Koiwa, H.; Barb, A.W.; Xiong, L.; Li, F.; McCully, M.G.; Lee, B.H.; Sokolchik, I.; Zhu, J.; Gong, Z.; Reddy, M.; et al. C-terminal domain phosphatase-like family members (AtCPLs) differentially regulate Arabidopsis thaliana abiotic stress signaling, growth, and development. Proc. Natl. Acad. Sci. USA 2002, 99, 10893–10898. [Google Scholar] [CrossRef] [Green Version]
  31. Thatcher, L.F.; Foley, R.; Casarotto, H.J.; Gao, L.L.; Kamphuis, L.G.; Melser, S.; Singh, K.B. The Arabidopsis RNA Polymerase II Carboxyl Terminal Domain (CTD) Phosphatase-Like1 (CPL1) is a biotic stress susceptibility gene. Sci. Rep. 2018, 8, 13454. [Google Scholar] [CrossRef]
  32. Ji, H.; Kim, S.R.; Kim, Y.H.; Kim, H.; Eun, M.Y.; Jin, I.D.; Cha, Y.S.; Yun, D.W.; Ahn, B.O.; Lee, M.C.; et al. Inactivation of the CTD phosphatase-like gene OsCPL1 enhances the development of the abscission layer and seed shattering in rice. Plant J. 2010, 61, 96–106. [Google Scholar] [CrossRef]
  33. Wu, Y.; Thorne, E.T.; Sharp, R.E.; Cosgrove, D.J. Modification of expansin transcript levels in the maize primary root at low water potentials. Plant Physiol. 2001, 126, 1471–1479. [Google Scholar] [CrossRef] [Green Version]
  34. Guo, W.; Zhao, J.; Li, X.; Qin, L.; Yan, X.; Liao, H. A soybean beta-expansin gene GmEXPB2 intrinsically involved in root system architecture responses to abiotic stresses. Plant J. 2011, 66, 541–552. [Google Scholar] [CrossRef]
  35. Dai, F.; Nevo, E.; Wu, D.; Comadran, J.; Zhou, M.; Qiu, L.; Chen, Z.; Beiles, A.; Chen, G.; Zhang, G. Tibet is one of the centers of domestication of cultivated barley. Proc. Natl. Acad. Sci. USA 2012, 109, 16969–16973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lu, P.; Kang, M.; Jiang, X.; Dai, F.; Gao, J.; Zhang, C. RhEXPA4, a rose expansin gene, modulates leaf growth and confers drought and salt tolerance to Arabidopsis. Planta 2013, 237, 1547–1559. [Google Scholar] [CrossRef] [PubMed]
  37. Chen, Y.; Li, C.; Yi, J.; Yang, Y.; Lei, C.; Gong, M. Transcriptome response to drought, rehydration and re-dehydration in potato. Int. J. Mol. Sci. 2019, 21, 159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Lopez-Maury, L.; Marguerat, S.; Bahler, J. Tuning gene expression to changing environments: From rapid responses to evolutionary adaptation. Nat. Rev. Genet. 2008, 9, 583–593. [Google Scholar] [CrossRef]
  39. Kürsteiner, O.; Dupuis, I.; Kuhlemeier, C. The pyruvate decarboxylase1 gene of Arabidopsis is required during anoxia but not other environmental stresses. Plant Physiol. 2003, 132, 968–978. [Google Scholar] [CrossRef] [Green Version]
  40. Uga, Y.; Kitomi, Y.; Ishikawa, S.; Yano, M. Genetic improvement for root growth angle to enhance crop production. Breed Sci. 2015, 65, 111–119. [Google Scholar] [CrossRef] [Green Version]
  41. Dolferus, R.; Ellis, M.; De Bruxelles, G.; Trevaskis, B.; Hoeren, F.; Dennis, E.; Peacock, W. Strategies of gene action in Arabidopsis during hypoxia. Ann. Bot. 1997, 79, 21–31. [Google Scholar] [CrossRef] [Green Version]
  42. Wahid, A.; Gelani, S.; Ashraf, M.; Foolad, M.R. Heat tolerance in plants: An overview. Environ. Exp. Bot. 2007, 61, 199–223. [Google Scholar] [CrossRef]
  43. Abe, H.; Urao, T.; Ito, T.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 2003, 15, 63–78. [Google Scholar] [CrossRef] [Green Version]
  44. Jaradat, M.R.; Feurtado, J.A.; Huang, D.; Lu, Y.; Cutler, A.J. Multiple roles of the transcription factor AtMYBR1/AtMYB44 in ABA signaling, stress responses, and leaf senescence. BMC Plant Biol. 2013, 13, 192. [Google Scholar] [CrossRef] [Green Version]
  45. Agarwal, M.; Hao, Y.; Kapoor, A.; Dong, C.-H.; Fujii, H.; Zheng, X.; Zhu, J.-K. A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance. J. Biol. Chem. 2006, 281, 37636–37645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Rahaie, M.; Xue, G.P.; Naghavi, M.R.; Alizadeh, H.; Schenk, P.M. A MYB gene from wheat (Triticum aestivum L.) is up-regulated during salt and drought stresses and differentially regulated between salt-tolerant and sensitive genotypes. Plant Cell Rep. 2010, 29, 835–844. [Google Scholar] [CrossRef] [PubMed]
  47. Zhang, L.; Zhao, G.; Jia, J.; Liu, X.; Kong, X. Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress. J. Exp. Bot. 2012, 63, 203–214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Nakano, T.; Fujisawa, M.; Shima, Y.; Ito, Y. The AP2/ERF transcription factor SlERF52 functions in flower pedicel abscission in tomato. J. Exp. Bot. 2014, 65, 3111–3119. [Google Scholar] [CrossRef] [Green Version]
  49. Koyama, T.; Nii, H.; Mitsuda, N.; Ohta, M.; Kitajima, S.; Ohme-Takagi, M.; Sato, F. A regulatory cascade involving class II ETHYLENE RESPONSE FACTOR transcriptional repressors operates in the progression of leaf senescence. Plant Physiol. 2013, 162, 991–1005. [Google Scholar] [CrossRef] [Green Version]
  50. Mehrnia, M.; Balazadeh, S.; Zanor, M.I.; Mueller-Roeber, B. EBE, an AP2/ERF transcription factor highly expressed in proliferating cells, affects shoot architecture in Arabidopsis. Plant Physiol. 2013, 162, 842–857. [Google Scholar] [CrossRef] [Green Version]
  51. Houston, K.; McKim, S.M.; Comadran, J.; Bonar, N.; Druka, I.; Uzrek, N.; Cirillo, E.; Guzy-Wrobelska, J.; Collins, N.C.; Halpin, C.; et al. Variation in the interaction between alleles of HvAPETALA2 and microRNA172 determines the density of grains on the barley inflorescence. Proc. Natl. Acad. Sci. USA 2013, 110, 16675–16680. [Google Scholar] [CrossRef] [Green Version]
  52. Johansson, I.; Karlsson, M.; Johanson, U.; Larsson, C.; Kjellbom, P. The role of aquaporins in cellular and whole plant water balance. Biochim. Biophys. Acta 2000, 1465, 324–342. [Google Scholar] [CrossRef] [Green Version]
  53. Maurel, C.; Verdoucq, L.; Luu, D.T.; Santoni, V. Plant aquaporins: Membrane channels with multiple integrated functions. Annu. Rev. Plant Biol. 2008, 59, 595–624. [Google Scholar] [CrossRef] [Green Version]
  54. Li, J.; Ban, L.; Wen, H.; Wang, Z.; Dzyubenko, N.; Chapurin, V.; Gao, H.; Wang, X. An aquaporin protein is associated with drought stress tolerance. Biochem. Biophys. Res. Commun. 2015, 459, 208–213. [Google Scholar] [CrossRef] [Green Version]
  55. Afzal, Z.; Howton, T.C.; Sun, Y.; Mukhtar, M.S. The roles of aquaporins in plant stress responses. J. Dev. Biol. 2016, 4, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Shi, X.; Ling, H.-Q. Current advances in genome sequencing of common wheat and its ancestral species. Crop J. 2018, 6, 15–21. [Google Scholar] [CrossRef]
  57. Hedden, P. The genes of the green revolution. Trends Genet. 2003, 19, 5–9. [Google Scholar] [CrossRef]
  58. Borner, A.; Worland, A.; Plaschke, J.; Schumann, E.; Law, C. Pleiotropic effects of genes for reduced height (Rht) and day-length insensitivity (Ppd) on yield and its components for wheat grown in middle Europe. Plant Breed. 1993, 111, 204–216. [Google Scholar] [CrossRef]
  59. Youssefian, S.; Kirby, E.; Gale, M. Pleiotropic effects of the GA-insensitive Rht dwarfing genes in wheat. 2. Effects on leaf, stem, ear and floret growth. Field Crops Res. 1992, 28, 191–210. [Google Scholar] [CrossRef]
  60. Geisel, N. Constitutive versus responsive gene expression strategies for growth in changing environments. PLoS ONE 2011, 6, e27033. [Google Scholar] [CrossRef] [Green Version]
  61. Blum, A. Drought resistance and its improvement. In Plant Breeding for Water-Limited Environments; Springer: Berlin/Heidelberg, Germany, 2011; pp. 53–152. [Google Scholar]
  62. Dawe, R.K.; Reed, L.M.; Yu, H.-G.; Muszynski, M.G.; Hiatt, E.N. A maize homolog of mammalian CENPC is a constitutive component of the inner kinetochore. Plant Cell 1999, 11, 1227–1238. [Google Scholar] [CrossRef]
  63. Zhong, C.X.; Marshall, J.B.; Topp, C.; Mroczek, R.; Kato, A.; Nagaki, K.; Birchler, J.A.; Jiang, J.; Dawe, R.K. Centromeric retroelements and satellites interact with maize kinetochore protein CENH3. Plant Cell 2002, 14, 2825–2836. [Google Scholar] [CrossRef]
  64. Cheeseman, I.M.; Desai, A. Molecular architecture of the kinetochore-microtubule interface. Nat. Rev. Mol. Cell. Biol. 2008, 9, 33–46. [Google Scholar] [CrossRef]
  65. Zheng, H.; Wu, H.; Pan, X.; Jin, W.; Li, X. Aberrant meiotic modulation partially contributes to the lower germination rate of pollen grains in maize (Zea mays L.) under low nitrogen supply. Plant Cell Physiol. 2017, 58, 342–353. [Google Scholar] [CrossRef]
  66. Li, X.; Dawe, R.K. Fused sister kinetochores initiate the reductional division in meiosis I. Nat. Cell Biol. 2009, 11, 1103–1108. [Google Scholar] [CrossRef] [PubMed]
  67. Powers, A.F.; Franck, A.D.; Gestaut, D.R.; Cooper, J.; Gracyzk, B.; Wei, R.R.; Wordeman, L.; Davis, T.N.; Asbury, C.L. The Ndc80 kinetochore complex forms load-bearing attachments to dynamic microtubule tips via biased diffusion. Cell 2009, 136, 865–875. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Jabbari, M.; Fakheri, B.A.; Aghnoum, R.; Mahdi Nezhad, N.; Ataei, R. GWAS analysis in spring barley (Hordeum vulgare L.) for morphological traits exposed to drought. PLoS ONE 2018, 13, e0204952. [Google Scholar] [CrossRef]
  69. Varshney, R.K.; Paulo, M.J.; Grando, S.; van Eeuwijk, F.A.; Keizer, L.C.P.; Guo, P.; Ceccarelli, S.; Kilian, A.; Baum, M.; Graner, A. Genome wide association analyses for drought tolerance related traits in barley (Hordeum vulgare L.). Field Crops Res. 2012, 126, 171–180. [Google Scholar] [CrossRef] [Green Version]
  70. Klepikova, A.V.; Kasianov, A.S.; Gerasimov, E.S.; Logacheva, M.D.; Penin, A.A. A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling. Plant J. 2016, 88, 1058–1070. [Google Scholar] [CrossRef] [PubMed]
  71. Finkelstein, R.; Gampala, S.S.L.; Lynch, T.J.; Thomas, T.L.; Rock, C.D. Redundant and distinct functions of the ABA response loci ABA-INSENSITIVE(ABI)5 and ABRE-BINDING FACTOR (ABF)3. Plant Mol. Biol. 2005, 59, 253–267. [Google Scholar] [CrossRef]
  72. Mittal, A.; Gampala, S.S.L.; Ritchie, G.L.; Payton, P.; Burke, J.J.; Rock, C.D. Related to ABA-Insensitive3(ABI3)/Viviparous1 and AtABI5 transcription factor coexpression in cotton enhances drought stress adaptation. Plant Biotechnol. J. 2014, 12, 578–589. [Google Scholar] [CrossRef] [Green Version]
  73. Tang, N.; Zhang, H.; Li, X.; Xiao, J.; Xiong, L. Constitutive activation of transcription factor OsbZIP46 improves drought tolerance in rice. Plant Physiol. 2012, 158, 1755–1768. [Google Scholar] [CrossRef] [Green Version]
Genes 11 00533 g001 550
Figure 1. Histogram of phenotypic values distribution analysis of the studied traits; (a) Awn Tipping, (b) Spike Heading, (c) Anther Extrusion, (d) Plant Height, (e) Spike length, (f) Awn Length, (g) Spikelet per Spike, (h) Grain per Spike and (i) Thousand Grain Weight (TGW) in 121 spring barley accessions under control and drought stress.
Figure 1. Histogram of phenotypic values distribution analysis of the studied traits; (a) Awn Tipping, (b) Spike Heading, (c) Anther Extrusion, (d) Plant Height, (e) Spike length, (f) Awn Length, (g) Spikelet per Spike, (h) Grain per Spike and (i) Thousand Grain Weight (TGW) in 121 spring barley accessions under control and drought stress.
Genes 11 00533 g001
Genes 11 00533 g002 550
Figure 2. Histogram of phenotypic values distribution analysis of drought tolerance index variation of the studied traits in 121 spring barley accessions under control and drought stress. (a) Awn Tipping, (b) Spike Heading, (c) Anther Extrusion, (d) Plant Height, (e) Spike length, (f) Awn Length, (g) Spikelet per Spike, (h) Grain per Spike and (i) Thousand Grain Weight (TGW).
Figure 2. Histogram of phenotypic values distribution analysis of drought tolerance index variation of the studied traits in 121 spring barley accessions under control and drought stress. (a) Awn Tipping, (b) Spike Heading, (c) Anther Extrusion, (d) Plant Height, (e) Spike length, (f) Awn Length, (g) Spikelet per Spike, (h) Grain per Spike and (i) Thousand Grain Weight (TGW).
Genes 11 00533 g002
Genes 11 00533 g003 550
Figure 3. Correlations matrix among the studied traits in barley genotypes (A) under control, and (B) under drought stress. The degree of significance is indicated as * p, 0.05; ** p, 0.01; *** p, 0.001; ns: not significant.
Figure 3. Correlations matrix among the studied traits in barley genotypes (A) under control, and (B) under drought stress. The degree of significance is indicated as * p, 0.05; ** p, 0.01; *** p, 0.001; ns: not significant.
Genes 11 00533 g003
Genes 11 00533 g004a 550Genes 11 00533 g004b 550
Figure 4. (a) The significant SNPs (101 SNPs, −log10 ≥ 3) associated with all traits under control and drought stress conditions. The x-axis shows the chromosomes and the SNP positions. The y-axis shows the −log10 (P-value) for each SNP marker. (b) Heatmap linkage disequilibrium to detect the region of candidate genes, which show a consistent effect on traits and associated with SNPs passing ≥ FDR using their physical position.
Figure 4. (a) The significant SNPs (101 SNPs, −log10 ≥ 3) associated with all traits under control and drought stress conditions. The x-axis shows the chromosomes and the SNP positions. The y-axis shows the −log10 (P-value) for each SNP marker. (b) Heatmap linkage disequilibrium to detect the region of candidate genes, which show a consistent effect on traits and associated with SNPs passing ≥ FDR using their physical position.
Genes 11 00533 g004aGenes 11 00533 g004b
Genes 11 00533 g005 550
Figure 5. SNP-gene-based analysis for (a) Spikelet number per spike, (b) Grain number per spike in barley genotype, and (c) TKW (gram). The degree of significance indicated as * p, 0.05; ** p, 0.01; *** p, 0.001; **** p, 0.0001.
Figure 5. SNP-gene-based analysis for (a) Spikelet number per spike, (b) Grain number per spike in barley genotype, and (c) TKW (gram). The degree of significance indicated as * p, 0.05; ** p, 0.01; *** p, 0.001; **** p, 0.0001.
Genes 11 00533 g005
Genes 11 00533 g006 550
Figure 6. Expression analysis of the candidate genes in different plant organs at different developmental stages.
Figure 6. Expression analysis of the candidate genes in different plant organs at different developmental stages.
Genes 11 00533 g006
Table
Table 1. The name and abbreviation of measured traits and respective description of measurements.
Table 1. The name and abbreviation of measured traits and respective description of measurements.
TraitAbbreviationDescription
Control Drought
Awn TippingAT_CAT_DThe number of days from planting up to awn tipping.
Spike HeadingSH_CSH_DThe number of days from planting up to spike heading
Anther ExtrusionAE_CAE_DThe number of days from planting up to anther extrusion.
Plant HeightPH_CPH_DThe distance between the ground level to the tip of the highest spikelet (excluding awns) in cm.
Spike LengthSL_CSL_DDistance from the base of the spike to the tip of the highest spikelet (excluding awns) in cm.
Awn LengthAL_CAL_DDistance from the tip of the spike to the end of the awn in cm.
No of Spikelets per SpikeNSS_CNSS_DThe actual count of the number of spikelets.
No of Grains per SpikeNGS_CNGS_DThe actual count of the number of the grains.
Thousand Grain WeightTGW_CTGW_DThe weight of 1000 grains randomly taken from each plot in gram (g).
Drought tolerance index
(Awn Tipping)		DTI_ATDTI ( AT ) = AT   under   drought AT   under   control × 100
Drought tolerance index
(Spike Heading)		DTI_SHDTI ( SH ) = SH   under   drought SH   under   control × 100
Drought tolerance index
(Anther Extrusion)		DTI_AEDTI ( AE ) = AE   under   drought AE   under   control × 100
Drought tolerance index
(Plant Height)		DTI_PHDTI ( PH ) = PH   under   drought PH   under   control × 100
Drought tolerance index
(Spike Length)		DTI_SLDTI ( SL ) = SL   under   drought SL   under   control × 100
Drought tolerance index
(Awn Length)		DTI_ALDTI ( AL ) = AL   under   drought AL   under   control × 100
Drought tolerance index
(No of Spikelet per Spike)		DTI_NSSDTI ( NSS ) = NSS   under   drought NSS   under   control × 100
Drought tolerance index
(No of Grain per Spike)		DTI_NGSDTI ( NGS ) = NGS   under   drought NGS   under   control × 100
Drought tolerance index
(Thousand Grain Weight)		DTI_TGWDTI ( TGW ) = TGW   under   drought TGW   under   control × 100
Table
Table 2. Analysis of variance and heritability for the measured traits under control and drought treatments.
Table 2. Analysis of variance and heritability for the measured traits under control and drought treatments.
TraitControlH2DroughtH2
TGT × GTGT × G
Awn Tippingns*****0.72ns*****0.84
Spike Headingns*****0.80ns*****0.78
Anther Extrusionns****0.79ns****0.81
Plant Height*********0.75*********0.71
Spike Length*********0.74*********0.69
Awn Length*********0.79*********0.76
No of Spikelet per Spike*********0.76*********0.65
No of Grain per Spike*********0.77*********0.68
Thousand Grain Weight*********0.72*********0.64
Awn Tipping_DTI̶***̶0.78
Spike Heading_DTI̶***̶0.71
Anther Extrusion_DTI̶***̶0.68
Plant Height_DTI̶***̶0.73
Spike Length_DTI̶***̶0.74
Awn Length_DTI̶***̶0.71
No of Spikelet per Spike _DTI̶***̶0.74
No of Grain per Spike _DTI̶***̶0.68
Thousand Grain Weight_DTI̶***̶0.77
H2—Heritability; T—Treatment; G—Genotype; T × G—Treatment by Genotype interaction; DTI—Drought Tolerant Index; The degree of significance is indicated as * p, 0.05; ** p, 0.01; *** p, 0.001; ns: not significant.
Table
Table 3. The functional annotation of the putative candidate genes associated with the studied traits under drought and control growth conditions.
Table 3. The functional annotation of the putative candidate genes associated with the studied traits under drought and control growth conditions.
TraitMarkerChrPos.SNP Pos.GeneGene Pos.GO IDAnnotation in BarleyOrthologs
NGS_C
NGS_D
NSS_D		SCRI_RS_1665402H75.56646934425HORVU2Hr1G091030646930069-646939693NoneRNA polymerase II C-terminal domain phosphatase-like 1AT4G21670
NGS_D
NSS_D		SCRI_RS_1573472H75.56647255135HORVU2Hr1G091170647252087-647258342Noneexpansin B3 AT4G28250
AE_C
NGS_C NSS_C		SCRI_RS_2296933H44.2648643412HORVU3Hr1G01865048643412-48647662GO:0000287
GO:0003824
GO:0030976		pyruvate decarboxylase-2AT5G54960
AE_C
NGS_C NSS_C		BOPA2_12_307373H45.5563626447HORVU3Hr1G02043063623393-63623651NoneProtein HASTY 1AT3G05040
AE_C
NGS_C
NSS_C		BOPA1_7045-9503H45.5564618663HORVU3Hr1G02066064617802-64621872GO:0005515Chromosome 3B, genomic scaffold, cultivar Chinese SpringAT2G36270
AE_C
NGS_C
NSS_D		BOPA1_2391-5663H46.2555590274HORVU3Hr1G01959055915087-55917096GO:0003677myb domain protein 37AT5G23000
AE_C
SH_C		SCRI_RS_2300753H103.27624309908HORVU3Hr1G088300624309585-624311169NoneChromosome 3B, genomic scaffold, cultivar Chinese SpringAT2G36270
AE_C
SH_C		BOPA2_12_302233H104.32627767267HORVU3Hr1G089160627749990-627754917GO:0003677
GO:0003700
GO:0006355		AP2-like ethylene-responsive transcription factor AT2G28550
AE_C
SH_C		SCRI_RS_1923603H104.26627260224HORVU3Hr1G089080627258327-627265762Noneundescribed protein
AE_DTI
AL_DTI
AE_C
AL_C
SH_C
SL_C
TGW_C
TGW_D
SH_DTI
AT_DTI		SCRI_RS_1773133H126.69657713242HORVU3Hr1G098200657712024-657713975GO:0005515Chromosome 3B, genomic scaffold, cultivar Chinese SpringAT2G36270
AE_DTI
SH_DTI		BOPA2_12_109813H154.15696452271HORVU3Hr1G116790696450874-696453390GO:0005215
GO:0006810
GO:0016020		Aquaporin-like superfamily proteinAT2G36830
AE_DTI
SH_DTI		SCRI_RS_2377383H154.6694105354HORVU3Hr1G115810694103839-694106776NoneKinetochore protein spc25

Share and Cite

MDPI and ACS Style

Thabet, S.G.; Moursi, Y.S.; Karam, M.A.; Börner, A.; Alqudah, A.M. Natural Variation Uncovers Candidate Genes for Barley Spikelet Number and Grain Yield under Drought Stress. Genes 2020, 11, 533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11050533

AMA Style

Thabet SG, Moursi YS, Karam MA, Börner A, Alqudah AM. Natural Variation Uncovers Candidate Genes for Barley Spikelet Number and Grain Yield under Drought Stress. Genes. 2020; 11(5):533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11050533

Chicago/Turabian Style

Thabet, Samar G., Yasser S. Moursi, Mohamed A. Karam, Andreas Börner, and Ahmad M. Alqudah. 2020. "Natural Variation Uncovers Candidate Genes for Barley Spikelet Number and Grain Yield under Drought Stress" Genes 11, no. 5: 533. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11050533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop